First Problem

This is the shortest end-to-end linear elasticity workflow.

using LowLevelFEM

structured_rect_mesh() # uses Gmsh to create mesh

mat = Material("body", E=2e5, ν=0.3)
prob = Problem([mat], type=:PlaneStress)

bc = displacementConstraint("supp", ux=0, uy=0)
ld = load("load", fy=-1)

u = solveDisplacement(prob, load=[ld], support=[bc])
S = solveStress(u)

showDoFResults(u, :ux)
showStressResults(S, :s)

openPostProcessor()

Next Steps

• Continue with Mesh and Physical Groups
• See full worked models in Tutorials