MeshFile non_regression.msh

CutMeshPlane 0 1 0 0
CutMeshCylinder 0 0.25 0 0 0 1 0.1
CutMeshSphere 0 -0.25 0 0.1

#_Sphere,_cylindre,_bas,_haut
CutElasticDomain 53 1e9 0.01
#CutElasticDomain 55 1e9 0.01
CutElasticDomain 56 1e9 0.49
CutElasticDomain 57 1e10 0.25

#UnCutElasticDomain 51 1e9 0.3

FaceDisplacement 35 0 0
FaceDisplacement 35 1 0
FaceDisplacement 35 2 0


#_Displacement_function
#----------------------------------------
#
#_Nouvelle_numerotation_du_noeud_1->3
#NodeDisplacementFunction 3 0 x x y z
#NodeDisplacementFunction 3 1 y x y z
#NodeDisplacementFunction 3 2 z x y z
#
#_Arrete_supperieur_laterale.
EdgeDisplacementFunction 11 0 0.1*z x y z
#EdgeDisplacementFunction 11 1 y x y z
#EdgeDisplacementFunction 11 2 z x y z
#
#_Face_originale
#FaceDisplacementFunction 35 0 x x y z
#FaceDisplacementFunction 35 1 y x y z
#FaceDisplacementFunction 35 2 z x y z


#_Force_function
#----------------------------------------
NodeForceFunction 3 -x*10^6 -y*10^6 -z*10^6 x y z
#
EdgeForceFunction 11 x*10^6 y*10^6 z*10^6 x y z
#
FaceForceFunction 33 x*10^9 y*10^9 z*10^9 x y z
#
VolumeForceFunction 53 x*10^11 -y*10^11 z*10^11 x y z


# Lagrange multipliers function
#----------------------------------------
StableLagMultField 40 1 0.1
StableLagMultFieldFunction 40 2 -y*0.2 x y z
#LagrangeMultipliersFunction 1 0.01 1 0 -1 0 1*0.01 x y z


Solver 3

#--------------------------------
# Valeur reference
# elastic energy=368673525.198862
#