Hi,
I’m trying to extract how the stiffness across a contact is distributed depending on the mesh density.
I have request *CONTACT FILE CELS but for some reason the result is only shown at one node of the slave surface, and the value is not as expected (looks like the value of the first integration point).
*CONTACT PRINT CELS works fine and the addition of all contact springs energy (slave element+face,energy) for all contact elements shows the expected result.
Any idea what am I doing wrong?
*NODE
1,1.2,0.5,1
2,1.3,-0.6,0
3,1.3,0.6,0
4,0.1,-0.6,0
5,1.3,0.6,-0.2
6,1.3,-0.6,-0.2
7,0,0,-0.15
8,0.1,0.6,0
9,0.1,0.6,-0.2
10,0.1,-0.6,-0.2
11,1.2,-0.5,0
12,1.2,0.5,0
13,0.2,0.5,1
14,0.2,-0.5,1
15,0.2,0.5,0
16,1.2,-0.5,1
17,0.2,-0.5,0
*ELEMENT,TYPE=C3D8I
1,17,14,16,11,15,13,1,12
2,3,2,4,8,5,6,10,9
*NSET,NSET=NSET_CONTACT
11
12
15
17
*NSET,NSET=CELS_NODES
11
12
15
17
*ELSET,ELSET=HEX8
1
2
*ELSET,ELSET=Bases
2
*ELSET,ELSET=Linear
1
*SURFACE,NAME=Slave1
1,S6
*SURFACE,NAME=Master1
2,S1
*MATERIAL,NAME=Steel
*ELASTIC,TYPE=ISOTROPIC
100000000000,0.3
*DENSITY
7850
*SOLID SECTION,ELSET=Bases,MATERIAL=Steel
*SOLID SECTION,ELSET=Linear,MATERIAL=Steel
*BOUNDARY
1,3,,0
2,2,,0
4,1,,0
4,2,,0
13,3,,0
14,1,,0
14,2,,0
14,3,,0
16,2,,0
16,3,,0
*AMPLITUDE,NAME=A_1
0,0
1,5E-05
*SURFACE INTERACTION,NAME=SI_6
*SURFACE BEHAVIOR,PRESSURE-OVERCLOSURE=LINEAR
5E12
*CONTACT PAIR,INTERACTION=SI_6,TYPE=SURFACE TO SURFACE,ADJUST=1
Slave1,Master1
*STEP,NLGEOM=YES,INC=110,AMPLITUDE=STEP
*STATIC,SOLVER=PARDISO
0.01,1,0,0.01
*BOUNDARY,AMPLITUDE=A_1
2,3,,1
*BOUNDARY,AMPLITUDE=A_1
3,3,,1
*BOUNDARY,AMPLITUDE=A_1
4,3,,1
*BOUNDARY,AMPLITUDE=A_1
8,3,,1
*NODE FILE,GLOBAL=YES
U
*EL FILE
S,NOE
*CONTACT FILE
CDIS,CSTR
*CONTACT FILE
CDIS,CSTR,CELS
*CONTACT PRINT
CELS
*END STEP