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

some surface debugging ...

Browse Source
This commit is contained in:
Jukka Aho
2016-11-22 06:01:46 +02:00
parent 3641324b10
commit 48e72ea645
1 changed files with 13 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files
geometry/3d_beam_nx/beam_3d.jl
+13 -1
View File
@@ -78,12 +78,24 @@ bc1 = create_bc(mesh, "fixed")
# load
load = Problem(Elasticity, "pressure load", 3)
load.elements = create_surface_elements(mesh, :load)
nload = length(load.elements)
info("$nload elements in load surface")
area = 0.0
time = 0.0
for element in load.elements
for ip in get_integration_points(element)
detJ = element(ip, time, Val{:detJ})
area += ip.weight*detJ
end
end
info("load surface area: $area")
update!(load, "surface pressure", 20.0)
# solution
isfile("results.h5") && rm("results.h5")
isfile("results.xmf") && rm("results.xmf")
solver = Solver(Linear, beam)
solver = Solver(Linear, beam, load)
solver.xdmf = Xdmf("results")
solver()
close(get(solver.xdmf).hdf)