OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fail_snconnect.F File Reference
#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine fail_snconnect (nel, nuparam, nuvar, nfunc, ifunc, npf, tf, time, timestep, uparam, uvar, ngl, ipg, npg, ndamf, epsd, pla, offg, offl, isolid, signzz, signyz, signzx, sym, area, dmg, damt, dfmax, tdele)

Function/Subroutine Documentation

◆ fail_snconnect()

subroutine fail_snconnect ( integer nel,
integer nuparam,
integer nuvar,
integer nfunc,
integer, dimension(*) ifunc,
integer, dimension(*) npf,
tf,
time,
timestep,
uparam,
uvar,
integer, dimension(nel) ngl,
integer ipg,
integer npg,
integer ndamf,
dimension(nel) epsd,
dimension(nel) pla,
dimension(nel) offg,
dimension(nel) offl,
integer isolid,
dimension(nel) signzz,
dimension(nel) signyz,
dimension(nel) signzx,
dimension(nel) sym,
dimension(nel) area,
dimension(nel) dmg,
damt,
dimension(nel) dfmax,
dimension(nel) tdele )

Definition at line 30 of file fail_snconnect.F.

37C-----------------------------------------------
38c SNCONNECT failure model for solid spotwelds
39C-----------------------------------------------
40#include "implicit_f.inc"
41C---------+---------+---+---+--------------------------------------------
42C VAR | SIZE |TYP| RW| DEFINITION
43C---------+---------+---+---+--------------------------------------------
44C NEL | 1 | I | R | SIZE OF THE ELEMENT GROUP NEL
45C NUPARAM | 1 | I | R | SIZE OF THE USER PARAMETER ARRAY
46C NUVAR | 1 | I | R | NUMBER OF FAILURE ELEMENT VARIABLES
47C---------+---------+---+---+--------------------------------------------
48C NFUNC | 1 | I | R | NUMBER FUNCTION USED FOR THIS USER LAW not used
49C IFUNC | NFUNC | I | R | FUNCTION INDEX not used
50C NPF | * | I | R | FUNCTION ARRAY
51C TF | * | F | R | FUNCTION ARRAY
52C---------+---------+---+---+--------------------------------------------
53C TIME | 1 | F | R | CURRENT TIME
54C TIMESTEP| 1 | F | R | CURRENT TIME STEP
55C UPARAM | NUPARAM | F | R | USER FAILURE PARAMETER ARRAY
56C EPSP | NEL | F | W | STRAIN RATE
57C---------+---------+---+---+--------------------------------------------
58C UVAR |NEL*NUVAR| F |R/W| USER ELEMENT VARIABLE ARRAY
59C OFFG | NEL | F |R/W| DELETED ELEMENT FLAG (=1. ON, =0. OFF)
60C-----------------------------------------------
61C Y = FINTER(IFUNC(J),X,NPF,TF,DYDX)
62C Y : y = f(x)
63C X : x
64C DYDX : f'(x) = dy/dx
65C IFUNC(J): FUNCTION INDEX
66C J : FIRST(J=1), SECOND(J=2) .. FUNCTION USED FOR THIS LAW
67C NPF,TF : FUNCTION PARAMETER
68CC-----------------------------------------------
69C C o m m o n B l o c k s
70C-----------------------------------------------
71#include "units_c.inc"
72#include "comlock.inc"
73C-----------------------------------------------
74C D u m m y A r g u m e n t s
75C-----------------------------------------------
76 INTEGER NEL,NUPARAM,NUVAR,NFUNC,IPG,NPG,NDAMF,ISOLID
77 INTEGER NGL(NEL),NPF(*),IFUNC(*)
78 my_real time,timestep
79 my_real uparam(nuparam),uvar(nel,nuvar),damt(nel,ndamf),tf(*)
80 my_real , DIMENSION(NEL) :: offg,offl,epsd,pla,dein,deit,
81 . signzz,signyz,signzx,sym,dmg,dfmax,tdele,area
82C-----------------------------------------------
83C L o c a l V a r i a b l e s
84C-----------------------------------------------
85 INTEGER I,J,IDEL,IDEV,NINDX,NINDXD,NINDXA,ISYM,FUNN,FUNT,IFUN2N,
86 . IFUN2T,IFUN3N,IFUN3T
87 INTEGER INDX(NEL),INDXD(NEL),INDXA(NEL)
88 my_real a2,b2,a3,b3,xscale2,xscale3,svmn,svmt,t1,t2,ttn,tts,fct,
89 . dydx,dama,ssym,cphi,sphi,areascale,defo
90 my_real , DIMENSION(NEL) :: fun2n,fun2t,fun3n,fun3t,dydx2n,
91 . dydx2t,dydx3n,dydx3t,phi,pla1,pla2
92C-----------------------------------------------
93C E x t e r n a l F u n c t i o n s
94C-----------------------------------------------
95 my_real finter
96 EXTERNAL finter
97C=======================================================================
98c Internal VARIABLES :
99c-----------------
100c UVAR1 initial plastic elongation
101c UVAR2 final plastic elongation (at rupture)
102c-----------------
103c Damage output VARIABLES :
104c-----------------
105c DAMT1 damage factor d, where (sig = sig *(1-d))
106c DAMT2 damage function at damage start
107c DAMT3 damage function at rupture
108C=======================================================================
109 a2 = uparam(1)
110 b2 = uparam(2)
111 a3 = uparam(3)
112 b3 = uparam(4)
113 isolid = nint(uparam(5))
114 xscale2 = uparam(6)
115 xscale3 = uparam(7)
116 isym = nint(uparam(8))
117 areascale = uparam(9)
118
119 ifun2n = ifunc(1)
120 ifun2t = ifunc(2)
121 ifun3n = ifunc(3)
122 ifun3t = ifunc(4)
123C-----------------------------------------------
124 IF (uvar(1,3) == zero) THEN
125 DO i=1,nel
126 uvar(i,3)= area(i) ! initial area at time = 0
127 ENDDO
128 ENDIF
129c
130 DO i=1,nel
131 pla1(i) = uvar(i,1)
132 pla2(i) = uvar(i,2)
133 END DO
134C-----------------------------------------------
135C RATE FUNCTIONS INERPOLATION
136C-----------------------------------------------
137 IF (ifun2n == 0) THEN
138 fun2n(1:nel) = one
139 ELSE
140 DO i=1,nel
141 fun2n(i) = finter(ifun2n,epsd(i)*xscale2,npf,tf,dydx2n)
142 ENDDO
143 ENDIF
144 IF (ifun2t == 0) THEN
145 fun2t(1:nel) = one
146 ELSE
147 DO i=1,nel
148 fun2t(i) = finter(ifun2t,epsd(i)*xscale2,npf,tf,dydx2t)
149 ENDDO
150 ENDIF
151 IF (ifun3n == 0) THEN
152 DO i=1,nel
153 fun3n(i) = one
154 ENDDO
155 ELSE
156 DO i=1,nel
157 fun3n(i) = finter(ifun3n,epsd(i)*xscale3,npf,tf,dydx3n)
158 ENDDO
159 ENDIF
160 IF (ifun3t == 0) THEN
161 DO i=1,nel
162 fun3t(i) = one
163 ENDDO
164 ELSE
165 DO i=1,nel
166 fun3t(i) = finter(ifun3t,epsd(i)*xscale3,npf,tf,dydx3t)
167 ENDDO
168 ENDIF
169C-----------------------------------------------
170 nindx = 0
171 nindxd = 0
172 nindxa = 0
173c
174 DO i=1,nel
175 idel = 0
176 dama = zero
177 ssym = sin(sym(i))
178 svmn = abs(signzz(i))
179 svmt = sqrt(signyz(i)**2 + signzx(i)**2)
180 phi(i)= atan(svmn/max(em20,svmt))
181 sphi = sin(phi(i))
182 cphi = cos(phi(i))
183c
184 IF (pla1(i) == zero) THEN ! No damage yet
185 IF (isym == 1 .AND. signzz(i)<= zero) THEN
186 t1 = zero
187 ELSE
188 t1 = sphi/(one-a2*ssym)/fun2n(i)
189 ENDIF
190 t2 = cphi/fun2t(i)
191 ttn = t1*pla(i)
192 tts = t2*pla(i)
193 fct = (ttn**b2 + tts**b2)**(one/b2)
194 damt(i,2) = min(fct, one)
195 IF (fct > one ) THEN
196 pla1(i) = (t1**b2 + t2**b2)**(-one/b2)
197 IF (isym == 1 .AND. signzz(i) <= zero)THEN
198 t1 = zero
199 ELSE
200 t1 = sphi/(one-a3*ssym)/fun3n(i)
201 ENDIF
202 t2 = cphi/fun3t(i)
203 ttn = t1*pla(i)
204 tts = t2*pla(i)
205 fct = (ttn**b3 + tts**b3) **(one / b3)
206 pla2(i) = (t1**b3 + t2**b3)**(-one/b3)
207 damt(i,1) = (pla(i)-pla1(i))/max(em20,(pla2(i) - pla1(i)))
208 damt(i,1) = min(damt(i,1), one)
209 damt(i,3) = min(fct, one)
210 dama = damt(i,1)
211 nindxd = nindxd+1
212 indxd(nindxd) = i
213 ENDIF
214c
215 ELSE ! (PLA1(I) > ZERO)
216 IF (isym == 1 .AND. signzz(i) <= zero) THEN
217 t1 = zero
218 ELSE
219 t1 = sphi/(one-a3*ssym)/fun3n(i)
220 ENDIF
221 t2 = cphi/fun3t(i)
222 ttn = t1*pla(i)
223 tts = t2*pla(i)
224 fct = (ttn**b3 + tts**b3) **(one / b3)
225 pla2(i) = (t1**b3 + t2**b3)**(-one/b3)
226 damt(i,1) = (pla(i)-pla1(i))/max(em20,(pla2(i) - pla1(i)))
227 damt(i,1) = min(damt(i,1), one)
228 damt(i,3) = min(fct, one)
229 dama = damt(i,1)
230c----- check if rupture ...
231 IF (fct > one) THEN
232 idel =1
233 dama = one
234 ENDIF
235 ENDIF
236 dmg(i) = max(dmg(i),dama)
237 dmg(i) = min(dmg(i),one)
238c-----------------------------
239 IF (idel == 1 .AND. offl(i) == one) THEN
240 offl(i) = zero ! local integ point rupture
241 nindx = nindx+1
242 indx(nindx) = i
243 tdele(i) = time
244 ENDIF
245C
246 uvar(i,1) = pla1(i)
247 uvar(i,2) = pla2(i)
248C------------- Maximum Damage storing for output : 0 < DFMAX < 1
249 dfmax(i) = min(one,max(dfmax(i),damt(i,1)))
250c
251c--- deformed elements check
252 IF (areascale > zero ) THEN
253 defo = uvar(i,3) * areascale
254 IF (area(i) > defo .AND. offg(i) == one) THEN
255 offl(i) = zero
256 nindxa = nindxa+1
257 indxa(nindxa) = i
258 tdele(i) = time
259 isolid = 1
260 ENDIF
261 ENDIF
262c-------------------------------
263 ENDDO ! I=1,NEL
264C-----------------------------------------------
265 IF (nindxd > 0) THEN
266 DO j=1,nindxd
267 i = indxd(j)
268#include "lockon.inc"
269 WRITE(iout ,1000) ngl(i),ipg,pla1(i)
270 WRITE(istdo,1100) ngl(i),ipg,pla1(i),time
271#include "lockoff.inc"
272 END DO
273 ELSEIF (nindx > 0) THEN
274 DO j=1,nindx
275 i = indx(j)
276#include "lockon.inc"
277 WRITE(iout ,1200) ngl(i),ipg,pla2(i)
278 WRITE(istdo,1300) ngl(i),ipg,pla2(i),time
279#include "lockoff.inc"
280 END DO
281 ELSEIF (nindxa > 0) THEN
282 DO j=1,nindxa
283 i = indxa(j)
284#include "lockon.inc"
285 WRITE(iout ,1400) ngl(i),area(i)
286 WRITE(istdo,1500) ngl(i),area(i),time
287#include "lockoff.inc"
288 END DO
289 ENDIF
290C-----------------------------------------------
291 1000 FORMAT(5x,'START DAMAGE CONNECTION ELEMENT ',i10,
292 . ' INTEGRATION POINT',i2,', PLASTIC ELONGATION =',1pe16.9)
293 1100 FORMAT(5x,'START DAMAGE CONNECTION ELEMENT ',i10,
294 . ' INTEGRATION POINT',i2,', PLASTIC ELONGATION =',1pe16.9,
295 . ' AT TIME ',1pe16.9)
296 1200 FORMAT(5x,'FAILURE CONNECTION SOLID ELEMENT ',i10,
297 . ' INTEGRATION POINT',i2,', PLASTIC ELONGATION=',1pe16.9)
298 1300 FORMAT(5x,'FAILURE CONNECTION SOLID ELEMENT ',i10,
299 . ' INTEGRATION POINT',i2,', PLASTIC ELONGATION :',1pe16.9,
300 . ' AT TIME ',1pe16.9)
301 1400 FORMAT(5x,'FAILURE CONNECTION SOLID ELEMENT ',i10,
302 . ', AREA (LIMIT REACHED) :',1pe16.9)
303 1500 FORMAT(5x,'FAILURE CONNECTION SOLID ELEMENT ',i10,
304 . ', AREA (LIMIT REACHED) :',1pe16.9,' AT TIME ',1pe16.9)
305C-----------------------------------------------
306 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21