Note
Go to the end to download the full example code.
HyperElastic2#
A hyper elastic cube in compression.
12 from EasyFEA import Display, ElemType, Models, Simulations, PyVista
13 from EasyFEA.Geoms import Domain
14
15 if __name__ == "__main__":
16 Display.Clear()
17
18 # ----------------------------------------------
19 # Configuration
20 # ----------------------------------------------
21
22 # geom
23 L = 1
24 h = 1
25
26 # ----------------------------------------------
27 # Mesh
28 # ----------------------------------------------
29 meshSize = h / 10
30
31 contour = Domain((0, 0), (L, h), h / 10)
32
33 mesh = contour.Mesh_Extrude(
34 [], [0, 0, h], [h / meshSize], ElemType.HEXA8, isOrganised=True
35 )
36
37 nodesX0 = mesh.Nodes_Conditions(lambda x, y, z: x == 0)
38 nodesXL = mesh.Nodes_Conditions(lambda x, y, z: x == L)
39
40 # ----------------------------------------------
41 # Simulation
42 # ----------------------------------------------
43
44 isot = Models.ElasIsot(3, E=1, v=0.3)
45 lmbda = isot.get_lambda()
46 mu = isot.get_mu()
47 mat = Models.SaintVenantKirchhoff(3, lmbda, mu)
48
49 simu = Simulations.HyperElasticSimu(mesh, mat)
50
51 uc = -0.3
52 simu.add_dirichlet(nodesX0, [0, 0, 0], simu.Get_unknowns())
53 simu.add_dirichlet(nodesXL, [uc, 0, 0], simu.Get_unknowns())
54
55 simu.Solve()
56
57 # ----------------------------------------------
58 # Results
59 # ----------------------------------------------
60
61 PyVista.Plot_BoundaryConditions(simu).show()
62 PyVista.Plot(simu, "ux", 1, plotMesh=True).show()
Total running time of the script: (0 minutes 0.997 seconds)