Glassvision/JuliaFEM.jl
1
0
Fork 0
mirror of https://github.com/JuliaFEM/JuliaFEM.jl.git synced 2026-09-09 20:43:33 +00:00
Code Issues Packages Projects Releases Wiki Activity

weighted gap calculation

Browse Source
This commit is contained in:
Jukka Aho
2015-12-31 07:41:09 +02:00
parent 2667c641ab
commit c0f3fec0dd
1 changed files with 0 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/mortar.jl
-11
View File
@@ -708,16 +708,6 @@ function assemble!{E<:MortarElements2D}(assembly::BoundaryAssembly, problem::Bou
N2 = master_element(xi_projected, time)
S = w*kron(N1', N1)
M = w*kron(N1', N2)
X1 = slave_element("geometry", xi_gauss, time)
u1 = slave_element("displacement", xi_gauss, time)
X2 = master_element("geometry", xi_projected, time)
u2 = master_element("displacement", xi_projected, time)
x1 = X1+u1
x2 = X2+u2
g = x1-x2 # gap in strong form
# info("gap in strong form", g)
gh = w*g*N1 # weighted gap (gap in weak form)
# info("weighted gap", gh)
for i=1:field_dim
sd = slave_dofs[i:field_dim:end]
md = master_dofs[i:field_dim:end]
@@ -725,7 +715,6 @@ function assemble!{E<:MortarElements2D}(assembly::BoundaryAssembly, problem::Bou
add!(assembly.C1, sd, md, -M)
add!(assembly.C2, sd, sd, S)
add!(assembly.C2, sd, md, -M)
add!(assembly.g, sd, gh[i,:])
end
end