Hallo zusammen,
ich habe mittlerweile eine Subroutine erstellt (Interaktion von UMAT und MPC), womit die besagte Rissmodellierung möglich ist. Momentan läuft das ganze noch etwas instabil, aber das bekomme ich sicherlich auch noch in den Griff.
Die MPC-Subroutine für eine direkte Interaktion (TIE-Constraint) zwischen zwei Knoten sieht bei mir jedenfalls wie folgt aus:
IF (Bedingung für Constraint) THEN
IF (MDOF.EQ.2) THEN
A(1,1,1) = 1.D0
A(2,2,1) = 1.D0
A(1,1,2) = -1.D0
A(2,2,2) = -1.D0
JDOF(1,1) = 1
JDOF(2,1) = 2
JDOF(1,2) = 1
JDOF(2,2) = 2
ENDIF
IF (MDOF.EQ.3) THEN
A(1,1,1) = 1.D0
A(2,2,1) = 1.D0
A(3,3,1) = 1.D0
A(1,1,2) = -1.D0
A(2,2,2) = -1.D0
A(3,3,2) = -1.D0
JDOF(1,1) = 1
JDOF(2,1) = 2
JDOF(3,1) = 6
JDOF(1,2) = 1
JDOF(2,2) = 2
JDOF(3,2) = 6
ENDIF
ELSE
LMPC = 0
ENDIF
Nur für den Fall, das jemand da nochmal Bedarf für haben sollte. Nachdem man sich das entsprechende ABAQUS-Beispiel durchgelesen hat kommt man aber auch recht schnell selber drauf.
Grüße,
C.G.
[Diese Nachricht wurde von DrReinerKlimpke am 28. Feb. 2017 editiert.]
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP