add plastic strain for type3 beam

OpenRadioss · Jan 11, 2023 · 96f63c1 · 96f63c1
1 parent 4108e0c
commit 96f63c1
Showing 1 changed file with 7 additions and 4 deletions.
diff --git a/engine/source/output/th/thpout.F b/engine/source/output/th/thpout.F
@@ -196,7 +196,6 @@ SUBROUTINE THPOUT(IPARG , NTHGRP2  , ITHGRP ,GEO,  IXP,
                                              LBUF  => BUFLY%LBUF(1,1,IPT)
                                              WA(IJK) = WA(IJK)+ LBUF%SIG(KK(3)+I) * AREAPT/AREA
                                             ENDDO 
-
                                            END IF !(BUFLY%L_SIG > 0) 
                                         END IF! (IGTYP) 
                                     ELSEIF (ITHBUF(L) > 11 .AND.ITHBUF(L) <= 254 ) THEN
@@ -209,8 +208,12 @@ SUBROUTINE THPOUT(IPARG , NTHGRP2  , ITHGRP ,GEO,  IXP,
                                             WA(IJK) = LBUF%SIG(KK(IK)+I)
                                         ENDIF
                                     ELSEIF (ITHBUF(L) == 255) THEN
-                                        WA(IJK)=ZERO 
-                                        IF(IGTYP == 18 ) THEN
+                                        WA(IJK)=ZERO
+                                        IF(IGTYP == 3 ) THEN 
+                                          IF(GBUF%G_PLA>0)THEN
+                                            WA(IJK)=GBUF%PLA(I) 
+                                          ENDIF
+                                        ELSEIF(IGTYP == 18 ) THEN
                                            IF (BUFLY%L_PLA > 0) THEN                                            
                                             DO IPT = 1,NPT
                                              AREAPT  = GEO(IPA+IPT,PID) 
@@ -224,7 +227,7 @@ SUBROUTINE THPOUT(IPARG , NTHGRP2  , ITHGRP ,GEO,  IXP,
                                            IF (BUFLY%L_PLA > 0) THEN                                                                                     
                                             IPT = ITHBUF(L) - 255 
                                             LBUF  => BUFLY%LBUF(1,1,IPT) 
-                                            WA(IJK) = LBUF%SIG(KK(IK)+I)
+                                            WA(IJK) = LBUF%PLA(KK(IK)+I)
                                            ENDIF
                                         ENDIF