FEM: Fix the shell to shell contact example (#23533)

* FEM: Update constraint_contact_shell_shell.py

* FEM: Update constraint_contact_shell_shell.inp
This commit is contained in:
FEA-eng
2025-09-08 03:31:38 +02:00
committed by GitHub
parent a1a1849cce
commit 3e5c5ceccd
2 changed files with 18 additions and 18 deletions

View File

@@ -40,7 +40,7 @@ def get_information():
"name": "Constraint Contact Shell Shell",
"meshtype": "face",
"meshelement": "Tria3",
"constraints": ["fixed", "force", "contact"],
"constraints": ["displacement", "force", "contact"],
"solvers": ["ccxtools"],
"material": "solid",
"equations": ["mechanical"],
@@ -171,13 +171,16 @@ def setup(doc=None, solvertype="ccxtools"):
material_obj.Material = mat
analysis.addObject(material_obj)
# constraint fixed
con_fixed = ObjectsFem.makeConstraintFixed(doc, "ConstraintFixed")
con_fixed.References = [
# constraint displacement
con_displacement = ObjectsFem.makeConstraintDisplacement(doc, "ConstraintDisplacement")
con_displacement.References = [
(lower_tube, "Edge2"),
(upper_tube, "Edge3"),
]
analysis.addObject(con_fixed)
con_displacement.xFree = False
con_displacement.yFree = False
con_displacement.zFree = False
analysis.addObject(con_displacement)
# constraint force
con_force = ObjectsFem.makeConstraintForce(doc, "ConstraintForce")

View File

@@ -22998,9 +22998,9 @@ Efaces
Efaces
***********************************************************
** constraints fixed node sets
** ConstraintFixed
*NSET,NSET=ConstraintFixed
** constraints displacement node sets
** ConstraintDisplacement
*NSET,NSET=ConstraintDisplacement
2,
3,
136,
@@ -38383,15 +38383,12 @@ DEPConstraintContact, INDConstraintContact
***********************************************************
** Fixed Constraints
** ConstraintFixed
** Displacement constraint applied
** ConstraintDisplacement
*BOUNDARY
ConstraintFixed,1
ConstraintFixed,2
ConstraintFixed,3
ConstraintFixed,4
ConstraintFixed,5
ConstraintFixed,6
ConstraintDisplacement,1,1,0.0
ConstraintDisplacement,2,2,0.0
ConstraintDisplacement,3,3,0.0
***********************************************************
@@ -38410,8 +38407,8 @@ U
*EL FILE
S, E
** outputs --> dat file
** reaction forces for Constraint fixed
*NODE PRINT, NSET=ConstraintFixed, TOTALS=ONLY
** reaction forces for Constraint displacement constraining translation
*NODE PRINT, NSET=ConstraintDisplacement, TOTALS=ONLY
RF
** contact forces for Constraint contact
*CONTACT PRINT, MASTER=INDConstraintContact, SLAVE=DEPConstraintContact