OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i3for2.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| i3for2 ../engine/source/interfaces/inter2d/i3for2.F
25!||--- called by ------------------------------------------------------
26!|| intvo2 ../engine/source/interfaces/inter2d/intvo2.F
27!||--- calls -----------------------------------------------------
28!|| ibcoff ../engine/source/interfaces/interf/ibcoff.F
29!||--- uses -----------------------------------------------------
30!|| h3d_mod ../engine/share/modules/h3d_mod.F
31!|| output_mod ../common_source/modules/output/output_mod.F90
32!||====================================================================
33 SUBROUTINE i3for2(OUTPUT, X ,E ,IRECT,LMSR ,MSR ,
34 1 NSV ,ILOC ,IRTL ,STF ,NSN ,
35 2 NMN ,NTY ,CST ,IRTLO,FRIC0,
36 3 FRIGAP,STFN ,IBC ,ICODT,IMAST,
37 4 FSAV ,FSKYI,ISKY ,PTMAX,AREAS,
38 5 FCONT ,FNCONT,FTCONT,H3D_DATA)
39C-----------------------------------------------
40C M o d u l e s
41C-----------------------------------------------
42 USE h3d_mod
43 USE output_mod
44C-----------------------------------------------
45C I m p l i c i t T y p e s
46C-----------------------------------------------
47#include "implicit_f.inc"
48#include "comlock.inc"
49C-----------------------------------------------
50C C o m m o n B l o c k s
51C-----------------------------------------------
52#include "parit_c.inc"
53#include "com01_c.inc"
54#include "com08_c.inc"
55#include "scr18_c.inc"
56#include "scr16_c.inc"
57#include "com06_c.inc"
58#include "scr07_c.inc"
59#include "scr14_c.inc"
60C-----------------------------------------------
61C D u m m y A r g u m e n t s
62C-----------------------------------------------
63 TYPE(output_), INTENT(inout) :: OUTPUT
64 INTEGER NSN, NMN, NTY, IBC, IMAST
65 INTEGER IRECT(4,*), LMSR(*), MSR(*), NSV(*), ILOC(*), IRTL(*),
66 . IRTLO(*), ICODT(*), ISKY(*)
67C REAL
68 my_real
69 . x(3,*), e(*), stf(*), cst(2,*), fric0(3,*), frigap(*),
70 . stfn(*), fsav(*),fskyi(lskyi,nfskyi),ptmax, areas(*),
71 . fcont(3,*),fncont(3,*), ftcont(3,*)
72 TYPE(h3d_database) :: H3D_DATA
73C-----------------------------------------------
74C L o c a l V a r i a b l e s
75C-----------------------------------------------
76 INTEGER IX(2), II, I, J, K, L, M, IMP, I3, I2, JJ, J3, J2, LOLD,
77 . niskyl, ipas
78C REAL
79 my_real
80 . h(2), n2, n3, fric, gap, ym1, zm1, ym2, zm2, ys, zs, t2, t3,
81 . xl, ans, ss, stif, fni, fyi, fzi, ss0, fti, ds, anst, fmax,
82 . stfri, ax, fs, ft
83C-----------------------------------------------
84 fric=frigap(1)
85 gap =frigap(2)
86C
87 DO 500 ii=1,nsn
88 ipas = 0
89 i=nsv(ii)
90 j=iloc(ii)
91 k=msr(j)
92 l=irtl(ii)
93 m=msr(irect(1,l))
94 ix(1)=m
95 ym1=x(2,m)
96 zm1=x(3,m)
97 m=msr(irect(2,l))
98 ix(2)=m
99 ym2=x(2,m)
100 zm2=x(3,m)
101 ys =x(2,i)
102 zs =x(3,i)
103 IF(n2d==1)THEN
104 ax=ys
105 ELSE
106 ax=one
107 ENDIF
108 t2=ym2-ym1
109 t3=zm2-zm1
110 xl=sqrt(t2**2+t3**2)
111 t2=t2/xl
112 t3=t3/xl
113 n2= t3
114 n3=-t2
115C
116 imp=0
117 i3=3*i
118 i2=i3-1
119C
120 ans =n2*(ys-ym1)+n3*(zs-zm1)
121 ans =ans-gap
122 IF(ans>zero)GOTO 120
123 h(2)=t2*(ys-ym1)+t3*(zs-zm1)
124 h(2)=h(2)/xl
125 h(1)=one-h(2)
126 ss=h(2)-h(1)
127 IF(ss> onep05)GO TO 120
128 IF(ss<-onep05)GO TO 120
129 ss= max(-one,ss)
130 ss= min( one,ss)
131C
132 IF(nty==5)THEN
133C Add test on Stiffness Second side
134 IF (stfn(ii)<zero) THEN
135 stif = zero
136 ELSE
137 stif=stf(l)
138 ENDIF
139 ELSE
140 stif=stf(l)*stfn(ii)/ max(em20,(stf(l)+stfn(ii)))
141 ENDIF
142 fni=ans*stif
143 fyi=n2*fni
144 fzi=n3*fni
145 imp=1
146C-------------------------------------
147C SAUVEGARDE DE L'IMPULSION TOTALE
148C-------------------------------------
149 fsav(2)=fsav(2)+fyi*imast*dt12*ax
150 fsav(3)=fsav(3)+fzi*imast*dt12*ax
151C
152 IF(iparit==0)THEN
153 DO 100 jj=1,2
154 j3=3*ix(jj)
155 j2=j3-1
156 e(j2)=e(j2)+fyi*h(jj)
157 e(j3)=e(j3)+fzi*h(jj)
158 100 CONTINUE
159 e(i2)=e(i2)-fyi
160 e(i3)=e(i3)-fzi
161 ELSE
162#include "lockon.inc"
163 niskyl = nisky
164 nisky = nisky + 3
165#include "lockoff.inc"
166 ipas = 1
167C
168 IF(kdtint==0)THEN
169 fskyi(niskyl+1,1)= zero
170 fskyi(niskyl+1,2)= fyi*h(1)
171 fskyi(niskyl+1,3)= fzi*h(1)
172 fskyi(niskyl+1,4)= zero
173 isky(niskyl+1) = ix(1)
174C
175 fskyi(niskyl+2,1)= zero
176 fskyi(niskyl+2,2)= fyi*h(2)
177 fskyi(niskyl+2,3)= fzi*h(2)
178 fskyi(niskyl+2,4)= zero
179 isky(niskyl+2) = ix(2)
180C
181 fskyi(niskyl+3,1)= zero
182 fskyi(niskyl+3,2)= -fyi
183 fskyi(niskyl+3,3)= -fzi
184 fskyi(niskyl+3,4)= zero
185 isky(niskyl+3) = i
186 ELSE
187 fskyi(niskyl+1,1)= zero
188 fskyi(niskyl+1,2)= fyi*h(1)
189 fskyi(niskyl+1,3)= fzi*h(1)
190 fskyi(niskyl+1,4)= zero
191 fskyi(niskyl+1,5)= zero
192 isky(niskyl+1) = ix(1)
193C
194 fskyi(niskyl+2,1)= zero
195 fskyi(niskyl+2,2)= fyi*h(2)
196 fskyi(niskyl+2,3)= fzi*h(2)
197 fskyi(niskyl+1,4)= zero
198 fskyi(niskyl+1,5)= zero
199 isky(niskyl+2) = ix(2)
200C
201 fskyi(niskyl+3,1)= zero
202 fskyi(niskyl+3,2)= -fyi
203 fskyi(niskyl+3,3)= -fzi
204 fskyi(niskyl+1,4)= zero
205 fskyi(niskyl+1,5)= zero
206 isky(niskyl+3) = i
207 ENDIF
208 ENDIF
209C
210 IF(anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT>0.AND.
211 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
212 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))THEN
213 fcont(2,ix(1)) =fcont(2,ix(1)) + fyi*h(1)
214 fcont(3,ix(1)) =fcont(3,ix(1)) + fzi*h(1)
215 fcont(2,ix(2)) =fcont(2,ix(2)) + fyi*h(2)
216 fcont(3,ix(2)) =fcont(3,ix(2)) + fzi*h(2)
217c
218 fcont(2,i)=fcont(2,i)- fyi
219 fcont(3,i)=fcont(3,i)- fzi
220 ENDIF
221
222 IF(anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT>0.AND.
223 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP) .OR.
224 . (manim>=4.AND.manim<=15).OR.h3d_data%MH3D/=0))THEN
225
226 fncont(2,ix(1)) =fncont(2,ix(1)) + fyi*h(1)
227 fncont(3,ix(1)) =fncont(3,ix(1)) + fzi*h(1)
228 fncont(2,ix(2)) =fncont(2,ix(2)) + fyi*h(2)
229 fncont(3,ix(2)) =fncont(3,ix(2)) + fzi*h(2)
230C
231 fncont(2,i)=fncont(2,i)- fyi
232 fncont(3,i)=fncont(3,i)- fzi
233 ENDIF
234
235 IF(ibc/=0)CALL ibcoff(ibc,icodt(i))
236C
237 120 CONTINUE
238 IF(fric==zero)GO TO 500
239 IF(imp==0) THEN
240 irtlo(ii)=0
241 fric0(2,ii)=zero
242 fric0(3,ii)=zero
243 GO TO 500
244 ENDIF
245C
246 lold=irtlo(ii)
247 IF(lold==0)THEN
248 irtlo(ii)=l
249 cst(1,ii)=ss
250 GO TO 500
251 ENDIF
252C
253 ss0=cst(1,ii)
254 fti=fric0(1,ii)
255 ds=ss-ss0
256 anst=half*ds*xl
257 fmax=-min(fric*fni,zero)
258 stfri=em01*stif
259 fti=fti + anst*stfri
260C
261 IF(fti>fmax)THEN
262 fti=fmax
263 ELSE
264 IF(fti<-fmax)THEN
265 fti=-fmax
266 ELSE
267 fric0(1,ii)=fti
268 irtlo(ii)=l
269 cst(1,ii)=ss
270 ENDIF
271 ENDIF
272C
273c-------limit tangential force is ON : FT<= YIELD/(S*sqrt(3))
274 fs = ptmax*areas(ii)/sqrt(three)
275 ft =fti
276 IF(fs>zero) THEN
277 IF(fti>fs)THEN
278 ft=fs
279 ELSEIF(fti<-fs)THEN
280 ft=-fs
281 ENDIF
282 ENDIF
283C---------------------------------------------------------------
284 fyi=t2*ft
285 fzi=t3*ft
286C-------------------------------------
287C SAUVEGARDE DE L'IMPULSION TOTALE
288C-------------------------------------
289 fsav(5)=fsav(5)+fyi*imast*dt12*ax
290 fsav(6)=fsav(6)+fzi*imast*dt12*ax
291C
292 IF(iparit==0)THEN
293 DO 400 jj=1,2
294 j3=3*ix(jj)
295 j2=j3-1
296 e(j2)=e(j2)+fyi*h(jj)
297 400 e(j3)=e(j3)+fzi*h(jj)
298 e(i2)=e(i2)-fyi
299 e(i3)=e(i3)-fzi
300 ELSE
301C
302 IF(ipas==0) THEN
303#include "lockon.inc"
304 niskyl = nisky
305 nisky = nisky + 3
306#include "lockoff.inc"
307 IF(kdtint==0)THEN
308 fskyi(niskyl,1)= zero
309 fskyi(niskyl+1,2)= fyi*h(1)
310 fskyi(niskyl+1,3)= fzi*h(1)
311 fskyi(niskyl+1,4)= zero
312 isky(niskyl+1) = ix(1)
313C
314 fskyi(niskyl+2,1)= zero
315 fskyi(niskyl+2,2)= fyi*h(2)
316 fskyi(niskyl+2,3)= fzi*h(2)
317 fskyi(niskyl+2,4)= zero
318 isky(niskyl+2) = ix(2)
319C
320 fskyi(niskyl+3,1)= zero
321 fskyi(niskyl+3,2)= -fyi
322 fskyi(niskyl+3,3)= -fzi
323 fskyi(niskyl+3,4)= zero
324 isky(niskyl+3) = i
325 ELSE
326 fskyi(niskyl,1)= zero
327 fskyi(niskyl+1,2)= fyi*h(1)
328 fskyi(niskyl+1,3)= fzi*h(1)
329 fskyi(niskyl+1,4)= zero
330 fskyi(niskyl+1,5)= zero
331 isky(niskyl+1) = ix(1)
332C
333 fskyi(niskyl+2,1)= zero
334 fskyi(niskyl+2,2)= fyi*h(2)
335 fskyi(niskyl+2,3)= fzi*h(2)
336 fskyi(niskyl+2,4)= zero
337 fskyi(niskyl+1,5)= zero
338 isky(niskyl+2) = ix(2)
339C
340 fskyi(niskyl+3,1)= zero
341 fskyi(niskyl+3,2)= -fyi
342 fskyi(niskyl+3,3)= -fzi
343 fskyi(niskyl+3,4)= zero
344 fskyi(niskyl+1,5)= zero
345 isky(niskyl+3) = i
346 ENDIF
347 ELSE
348C already treated higher
349 fskyi(niskyl+1,2)= fskyi(niskyl+1,2)+fyi*h(1)
350 fskyi(niskyl+1,3)= fskyi(niskyl+1,3)+fzi*h(1)
351C
352 fskyi(niskyl+2,2)= fskyi(niskyl+2,2)+fyi*h(2)
353 fskyi(niskyl+2,3)= fskyi(niskyl+2,3)+fzi*h(2)
354C
355 fskyi(niskyl+3,2)= fskyi(niskyl+3,2)-fyi
356 fskyi(niskyl+3,3)= fskyi(niskyl+3,3)-fzi
357 ENDIF
358 ENDIF
359C
360 IF(anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT>0.AND.
361 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
362 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))THEN
363 fcont(2,ix(1)) =fcont(2,ix(1)) + fyi*h(1)
364 fcont(3,ix(1)) =fcont(3,ix(1)) + fzi*h(1)
365 fcont(2,ix(2)) =fcont(2,ix(2)) + fyi*h(2)
366 fcont(3,ix(2)) =fcont(3,ix(2)) + fzi*h(2)
367c
368 fcont(2,i)=fcont(2,i)- fyi
369 fcont(3,i)=fcont(3,i)- fzi
370 ENDIF
371 IF(anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT>0.AND.
372 . ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
373 . (manim>=4.AND.manim<=15).OR.h3d_data%MH3D/=0))THEN
374 ftcont(2,ix(1)) =ftcont(2,ix(1)) + fyi*h(1)
375 ftcont(3,ix(1)) =ftcont(3,ix(1)) + fzi*h(1)
376 ftcont(2,ix(2)) =ftcont(2,ix(2)) + fyi*h(2)
377 ftcont(3,ix(2)) =ftcont(3,ix(2)) + fzi*h(2)
378C
379 ftcont(2,i)=ftcont(2,i)- fyi
380 ftcont(3,i)=ftcont(3,i)- fzi
381 ENDIF
382
383 500 CONTINUE
384 RETURN
385 END
my_real
#define my_real
Definition cppsort.cpp:32
i3for2
subroutine i3for2(output, x, e, irect, lmsr, msr, nsv, iloc, irtl, stf, nsn, nmn, nty, cst, irtlo, fric0, frigap, stfn, ibc, icodt, imast, fsav, fskyi, isky, ptmax, areas, fcont, fncont, ftcont, h3d_data)
Definition i3for2.F:39
ibcoff
subroutine ibcoff(ibc, icodt)
Definition ibcoff.F:44
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
h3d_mod
Definition h3d_mod.F:311
h3d_mod::h3d_database
Definition h3d_mod.F:438