hallo,
mit dem angehängten code sollen zwei übereinander stehende, durch eine schräge fuge getrennte, körper simuliert werden. leider bewegt sich in der fuge nichts. nicht mit dem auskommentierten reibungsmodell noch mit dem kohäsivzonenmodell. es scheint als wäre die fuge eine starre verbindung. hat jemand eine idee was ich falsch gemacht habe? ich sehe es leider nicht.
gruß
!------------------------------------------------------------------------------------------------------
!-------------STARTEN DES PREPROCESSOR-----------------------------------------------------------------
!------------------------------------------------------------------------------------------------------
/PREP7
!----Geometrie-----------------------------------------------------------------------------------------
*AFUN,DEG
b=150 ! [mm]
h=600 ! [mm]
Druckstrebenwinkel=50 ! [°]
Elementlaenge=15 !Definition der Elementlänge
Last=10 ![mm]
!----Materialparameter
MU=0.66 !Reibkennwert
!----Elementwahl---------------------------------------------------------------------------------------
!----unterer Betonkörper----------
ET,1,PLANE182
KEYOPT,1,1,2 !Enhanced Strain
kEYOPT,1,3,0 !Ebener Spannungszustand
!----oberer Betonkörper----------
ET,2,PLANE182
KEYOPT,2,1,2 !Enhanced Strain
kEYOPT,2,3,0 !Ebener Spannungszustand
!----Grenzfläche----
ET,3,CONTA171
KEYOPT,3,12,5 !Bonded always in contact
ET,4,TARGE169
!----Eingabe der Materialdaten-------------------------------------------------------------------------
!----Betondaten-----
!MP,EX,1,47538.70
!MP,PRXY,1,0.2
!TB,FRIC,3,,,,ISO
!Materialmodell, Coulomb-Reibung, Mat#,Art d. Reibung
!TBDATA,1,MU
!TBDATA,Startpunkt, Reibung
!----Kohäsivzone----
TB, CZM , 3 , 1 , 1 , CBDD
!Materialmodell, Cohesive zone, Mat#3, für 1 Temp, Anahl d. Punkte pro Temp, TBOption
TBDATA , 1 , 0 , 0 ,50.8, 0.22, 0.001 , 0
!TBDATA, Startpunkt, \sigma_n, u_n, \tau_max, u^c_t , dämpfung, u_t bei \sigma_n
!----Real Constants------------------------------------------------------------------------------------
R,4,,,
RMODIF,4,12,-317.5 !springen an Eingabezeile #12 FKT (Tangent penalty stiffness factor)
!----Koordinateneingabe der Keypoints------------------------------------------------------------------
!K,#,x,y
K,1,0,0
K,2,b,0
K,3,0,h
K,4,b,h
K,5,0,h/2+(TAN(Druckstrebenwinkel)*(b/2))
K,6,b,h/2-(TAN(Druckstrebenwinkel)*(b/2))
!----Eingabe von Linien--------------------------------------------------------------------------------
!L,vonK #,nach K#
L,1,2 !L# 1
L,3,4 !L# 2
L,1,5 !L# 3
L,5,3 !L# 4
L,2,6 !L# 5
L,6,4 !L# 6
L,5,6 !L# 7
!----Eingabe von Flächen-------------------------------------------------------------------------------
A,1,2,6,5 !unterer Betonkörper
A,3,4,6,5 !oberer Betonkörper
!----Vernetzen-----------------------------------------------------------------------------------------
!----unterer Betonkörper----------
MAT,1
REAL,1
TYPE,1
LESIZE,1,Elementlaenge
LESIZE,3,Elementlaenge
LESIZE,5,Elementlaenge
LESIZE,7,Elementlaenge
AMESH,1
!----oberer Betonkörper----------
MAT,1
REAL,2
TYPE,2
LESIZE,2,Elementlaenge
LESIZE,4,Elementlaenge
LESIZE,6,Elementlaenge
AMESH,2
!----Contact--------
TYPE,3
REAL,3
MAT,3
ESEL,S,TYPE,,1 !Selektiere alle Elemente des unteren Betonkörper für Kontaktelemente
NSLE,S !Selektiere alle Knoten der selektierten Elemente
LSEL,S,LINE,,7,7,1 !Selektiere Linennummer 7
NSLL,R,1 !Reselektiere alle die Knoten mit der selektierten Liniennummer
ESURF !Erzeuge Kontaktelemente
ALLSEL
!----Target---------
TYPE,4
REAL,4
MAT,3
ESEL,S,TYPE,,2 !Selektiere alle Elemente des oberen Betonkörper für Kontaktelemente
NSLE,S !Selektiere alle Knoten der selektierten Elemente
LSEL,S,LINE,,7,7,1 !Selektiere Linennummer 7
NSLL,R,1 !Reselektiere alle die Knoten mit der selektierten Liniennummer
ESURF
ALLSEL
!------------------------------------------------------------------------------------------------------
!-------------STARTEN DER LOESUNGSPHASE----------------------------------------------------------------
!------------------------------------------------------------------------------------------------------
/SOLU
!----Randbedingungen-----------------------------------------------------------------------------------
NSEL,S,LOC,Y,0
D,ALL,UY,0 !Vertikales Auflager unten
ALLSEL
!----Lasten--------------------------------------------------------------------------------------------
NSEL,S,LOC,Y,h
D,ALL,UY,-Last
ALLSEL
!----Lösungsbedingungen--------------------------------------------------------------------------------
NLGEOM,ON !Nonlinear geometry on
SOLVE
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP
Foren auf CAD.de (alle Foren)
| 
