OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fail_tsaiwu_c.F File Reference
#include "implicit_f.inc"
#include "units_c.inc"
#include "comlock.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine fail_tsaiwu_c (nel, nuparam, nuvar, uparam, uvar, time, ngl, ipg, ilay, ipt, signxx, signyy, signxy, signyz, signzx, off, foff, dmg_flag, dmg_scale, dfmax, lf_dammx, tdel, timestep)

Function/Subroutine Documentation

◆ fail_tsaiwu_c()

subroutine fail_tsaiwu_c ( integer, intent(in) nel,
integer, intent(in) nuparam,
integer, intent(in) nuvar,
intent(in) uparam,
intent(inout) uvar,
intent(in) time,
integer, dimension(nel), intent(in) ngl,
integer, intent(in) ipg,
integer, intent(in) ilay,
integer, intent(in) ipt,
intent(in) signxx,
intent(in) signyy,
intent(in) signxy,
intent(in) signyz,
intent(in) signzx,
intent(in) off,
integer, dimension(nel), intent(inout) foff,
integer, intent(out) dmg_flag,
intent(out) dmg_scale,
intent(inout) dfmax,
integer, intent(in) lf_dammx,
intent(out) tdel,
intent(in) timestep )

Definition at line 29 of file fail_tsaiwu_c.F.

35C-----------------------------------------------
36C modified Puck model ------
37C-----------------------------------------------
38C I m p l i c i t T y p e s
39C-----------------------------------------------
40#include "implicit_f.inc"
41C-----------------------------------------------
42C G l o b a l P a r a m e t e r s
43C-----------------------------------------------
44#include "units_c.inc"
45#include "comlock.inc"
46C-----------------------------------------------
47C I N P U T A r g u m e n t s
48C-----------------------------------------------
49 INTEGER ,INTENT(IN) :: NEL,NUPARAM,NUVAR,IPG,ILAY,IPT,LF_DAMMX
50 INTEGER ,DIMENSION(NEL) ,INTENT(IN) :: NGL
51 my_real ,INTENT(IN) :: time,timestep
52 my_real ,DIMENSION(NEL) ,INTENT(IN) :: off,
53 . signxx,signyy,signxy,signyz,signzx
54 my_real,DIMENSION(NUPARAM) ,INTENT(IN) :: uparam
55C-----------------------------------------------
56C I N P U T O U T P U T A r g u m e n t s
57C-----------------------------------------------
58 INTEGER ,INTENT(OUT) ::DMG_FLAG
59 INTEGER ,DIMENSION(NEL) ,INTENT(INOUT) :: FOFF
60 my_real ,DIMENSION(NEL,LF_DAMMX),INTENT(INOUT) :: dfmax
61 my_real ,DIMENSION(NEL) ,INTENT(OUT) :: tdel,dmg_scale
62 my_real ,DIMENSION(NEL,NUVAR) ,INTENT(INOUT) :: uvar
63C-----------------------------------------------
64C L o c a l V a r i a b l e s
65C-----------------------------------------------
66 INTEGER :: I,J,NINDX,IFAIL_SH
67 INTEGER ,DIMENSION(NEL) :: INDX
68 my_real
69 . tmax,fcut
70 my_real
71 . f1,f2,f11,f22,f66,f12,asrate,a,b,
72 . sxx(nel),syy(nel),sxy(nel),findex,rfactr
73C!----------------------------------------------
74 !=======================================================================
75 ! - INITIALISATION OF COMPUTATION ON TIME STEP
76 !=======================================================================
77 ! Recovering model parameters
78 f1 = uparam(1)
79 f2 = uparam(2)
80 f11 = uparam(3)
81 f22 = uparam(4)
82 f66 = uparam(5)
83 f12 = uparam(6)
84 tmax = uparam(7)
85 fcut = uparam(8)
86 ifail_sh = int(uparam(9))
87c
88 ! Stress softening activation
89 dmg_flag = 1
90c
91 ! Stress tensor filtering
92 IF (fcut > zero) THEN
93 asrate = two*pi*fcut*timestep
94 asrate = asrate/(one+asrate)
95 DO i = 1,nel
96 sxx(i) = asrate*signxx(i) + (one - asrate)*uvar(i,2)
97 syy(i) = asrate*signyy(i) + (one - asrate)*uvar(i,3)
98 sxy(i) = asrate*signxy(i) + (one - asrate)*uvar(i,4)
99 uvar(i,2) = sxx(i)
100 uvar(i,3) = syy(i)
101 uvar(i,4) = sxy(i)
102 ENDDO
103 ELSE
104 DO i = 1,nel
105 sxx(i) = signxx(i)
106 syy(i) = signyy(i)
107 sxy(i) = signxy(i)
108 ENDDO
109 ENDIF
110c
111 !====================================================================
112 ! - COMPUTATION OF THE DAMAGE VARIABLE EVOLUTION
113 !====================================================================
114 ! Initialization of element failure index
115 nindx = 0
116 indx(1:nel) = 0
117c
118 ! Loop over the elements
119 DO i=1,nel
120c
121 ! If damage has not been reached yet
122 IF (dfmax(i,1)<one) THEN
123c
124 ! Compute parameters A and B
125 a = f11*(sxx(i)**2) + f22*(syy(i)**2) +
126 . f66*(sxy(i)**2) + two*f12*sxx(i)*syy(i)
127 b = f1*sxx(i) + f2*syy(i)
128c
129 ! Compute failure index
130 findex = a + b
131 findex = max(zero,findex)
132c
133 ! Compute reserve factor
134 rfactr = (-b + sqrt((b**2)+four*a))/max((two*a),em20)
135 dfmax(i,2) = rfactr
136c
137 ! Damage variable update
138 dfmax(i,1) = min(one ,max(findex,dfmax(i,1)))
139 IF (dfmax(i,1) >= one) THEN
140 nindx = nindx+1
141 indx(nindx) = i
142 IF (ifail_sh > 0) THEN
143 uvar(i,1) = time
144 ENDIF
145 ENDIF
146 ENDIF
147c
148 ! Stress relaxation in case of damage reached
149 IF ((uvar(i,1) > zero).AND.(foff(i) /= 0).AND.(ifail_sh > 0)) THEN
150 dmg_scale(i) = exp(-(time - uvar(i,1))/tmax)
151 IF (dmg_scale(i) < em02) THEN
152 foff(i) = 0
153 tdel(i) = time
154 dmg_scale(i) = zero
155 ENDIF
156 ENDIF
157 ENDDO
158c
159 !====================================================================
160 ! - PRINTOUT DATA ABOUT FAILED ELEMENTS
161 !====================================================================
162 IF (nindx > 0) THEN
163 DO j=1,nindx
164 i = indx(j)
165#include "lockon.inc"
166 WRITE(iout, 2000) ngl(i),ipg,ilay,ipt
167 WRITE(istdo,2100) ngl(i),ipg,ilay,ipt,time
168#include "lockoff.inc"
169 END DO
170 END IF
171c------------------------
172 2000 FORMAT(1x,'FAILURE (TSAIWU) OF SHELL ELEMENT ',i10,1x,',GAUSS PT',
173 . i2,1x,',LAYER',i3,1x,',INTEGRATION PT',i3)
174 2100 FORMAT(1x,'FAILURE (TSAIWU) OF SHELL ELEMENT ',i10,1x,',GAUSS PT',
175 . i2,1x,',LAYER',i3,1x,',INTEGRATION PT',i3,1x,'AT TIME :',1pe12.4)
176c------------------------
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21