OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sms_rgwall.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!|| sms_rgwall_impact ../engine/source/ams/sms_rgwall.F
25!||--- called by ------------------------------------------------------
26!|| sms_rgwal_0 ../engine/source/ams/sms_rgwal0.F
27!||====================================================================
28 SUBROUTINE sms_rgwall_impact
29 1 (x ,a ,v ,rwl ,nsw ,
30 2 nsn ,itied ,msr ,ms ,weight ,
31 3 rwsav ,nimpact,impact ,nsms ,nrwl_sms)
32C-----------------------------------------------
33C I m p l i c i t T y p e s
34C-----------------------------------------------
35#include "implicit_f.inc"
36#include "comlock.inc"
37C-----------------------------------------------
38C C o m m o n B l o c k s
39C-----------------------------------------------
40#include "com08_c.inc"
41#include "sms_c.inc"
42C-----------------------------------------------
43C D u m m y A r g u m e n t s
44C-----------------------------------------------
45 INTEGER NSN, ITIED, MSR, NIMPACT, NSMS
46 INTEGER NSW(*), WEIGHT(*), IMPACT(*), NRWL_SMS(*)
47C REAL
48 my_real
49 . x(*), a(*), v(*), rwl(*), ms(*), rwsav(*)
50C-----------------------------------------------
51C L o c a l V a r i a b l e s
52C-----------------------------------------------
53 INTEGER I, N, N3, N2, N1,IFQ, J, K, M1, M2, M3
54
55C REAL
56 my_real
57 . xwl, ywl, zwl, vxw, vyw, vzw,
58 . vx, vy, vz, ux, uy, uz, xc, yc, zc, dp
59C-----------------------------------------------
60C
61C fixed wall only.
62C
63 IF(msr==0)THEN
64 xwl=rwl(4)
65 ywl=rwl(5)
66 zwl=rwl(6)
67 vxw=zero
68 vyw=zero
69 vzw=zero
70 ELSE
71 m3=3*msr
72 m2=m3-1
73 m1=m2-1
74 vxw=v(m1)+a(m1)*dt12
75 vyw=v(m2)+a(m2)*dt12
76 vzw=v(m3)+a(m3)*dt12
77 xwl=x(m1)+vxw*dt2
78 ywl=x(m2)+vyw*dt2
79 zwl=x(m3)+vzw*dt2
80 ENDIF
81C
82 nimpact=0
83C
84 ifq = nint(rwl(15))
85 IF (ifq>0) THEN
86 DO 20 j=1,nsms
87 i=nrwl_sms(j)
88 n=nsw(i)
89 n3=3*n
90 n2=n3-1
91 n1=n2-1
92 vx=v(n1)+a(n1)*dt12
93 vy=v(n2)+a(n2)*dt12
94 vz=v(n3)+a(n3)*dt12
95 ux=x(n1)+vx*dt2
96 uy=x(n2)+vy*dt2
97 uz=x(n3)+vz*dt2
98 xc=ux-xwl
99 yc=uy-ywl
100 zc=uz-zwl
101 dp=xc*rwl(1)+yc*rwl(2)+zc*rwl(3)
102 IF(dp>zero)THEN
103 k = 3*(i-1)
104 rwsav(k+1) = zero
105 rwsav(k+2) = zero
106 rwsav(k+3) = zero
107 GOTO 20
108 ENDIF
109C--- test pour noeuds penetres
110 IF((vx-vxw)*rwl(1)+(vy-vyw)*rwl(2)+(vz-vzw)*rwl(3)>zero)
111 . GOTO 20
112 nimpact = nimpact+1
113 impact(nimpact) = i
114 20 CONTINUE
115 ELSE
116 DO 21 j=1,nsms
117 i=nrwl_sms(j)
118 n=nsw(i)
119 n3=3*n
120 n2=n3-1
121 n1=n2-1
122 vx=v(n1)+a(n1)*dt12
123 vy=v(n2)+a(n2)*dt12
124 vz=v(n3)+a(n3)*dt12
125 ux=x(n1)+vx*dt2
126 uy=x(n2)+vy*dt2
127 uz=x(n3)+vz*dt2
128 xc=ux-xwl
129 yc=uy-ywl
130 zc=uz-zwl
131 dp=xc*rwl(1)+yc*rwl(2)+zc*rwl(3)
132 IF(dp>zero) GOTO 21
133C--- test pour noeuds penetres
134 IF((vx-vxw)*rwl(1)+(vy-vyw)*rwl(2)+(vz-vzw)*rwl(3)>zero)
135 . GOTO 21
136 nimpact = nimpact+1
137 impact(nimpact) = i
138 21 CONTINUE
139 ENDIF
140C
141 IF(nimpact/=0.AND.itied==2)ifricw=1
142C
143 RETURN
144 END
145!||====================================================================
146!|| sms_rgwall_fric ../engine/source/ams/sms_rgwall.F
147!||--- called by ------------------------------------------------------
148!|| sms_rgwal_0 ../engine/source/ams/sms_rgwal0.F
149!||====================================================================
150 SUBROUTINE sms_rgwall_fric
151 1 (x ,a ,v ,rwl ,nsw ,
152 2 nsn ,itied ,msr ,ms ,weight ,
153 3 rwsav ,nimpact,impact ,nsms ,nrwl_sms,
154 4 fsav ,fopt ,res ,r ,frea )
155C-----------------------------------------------
156C I m p l i c i t T y p e s
157C-----------------------------------------------
158#include "implicit_f.inc"
159#include "comlock.inc"
160C-----------------------------------------------
161C C o m m o n B l o c k s
162C-----------------------------------------------
163#include "com08_c.inc"
164C-----------------------------------------------
165C D u m m y A r g u m e n t s
166C-----------------------------------------------
167 INTEGER NSN, ITIED, MSR, NIMPACT, NSMS
168 INTEGER NSW(*), WEIGHT(*), IMPACT(*), NRWL_SMS(*)
169C REAL
170 my_real
171 . x(*), a(*), v(*), rwl(*), ms(*), rwsav(*), fsav(*),
172 . fopt(*), res(*), frea(*), r(*)
173C-----------------------------------------------
174C L o c a l V a r i a b l e s
175C-----------------------------------------------
176 INTEGER I, N, N3, N2, N1,IFQ, J, K, M1, M2, M3
177
178C REAL
179 my_real
180 . xwl, ywl, zwl, vxw, vyw, vzw, vx, vy, vz,
181 . dv, da, dvt, fn, msw,
182 . fnxn, fnyn, fnzn, fnxt, fnyt, fnzt, fndfn, ftdft, fric, fric2,
183 . fcoe, fac, alpha, alphi, fxt, fyt, fzt, fxn, fyn, fzn
184C-----------------------------------------------
185C
186 IF(msr==0)THEN
187 xwl=rwl(4)
188 ywl=rwl(5)
189 zwl=rwl(6)
190 vxw=zero
191 vyw=zero
192 vzw=zero
193 ELSE
194 m3=3*msr
195 m2=m3-1
196 m1=m2-1
197 vxw=v(m1)+a(m1)*dt12
198 vyw=v(m2)+a(m2)*dt12
199 vzw=v(m3)+a(m3)*dt12
200 xwl=x(m1)+vxw*dt2
201 ywl=x(m2)+vyw*dt2
202 zwl=x(m3)+vzw*dt2
203 ENDIF
204C
205 ifq = nint(rwl(15))
206C
207 IF(itied == 2)THEN
208C
209 IF (ifq>0) THEN
210C--- friction filtering
211 fric = rwl(13)
212 alpha= rwl(14)
213 IF (ifq==3) alpha = alpha * dt12
214 alphi = one - alpha
215 fric2 = fric**2
216 fac=one/dt12
217 DO j = 1,nimpact
218 i = abs(impact(j))
219 n=nsw(i)
220 n3=3*n
221 n2=n3-1
222 n1=n2-1
223C---
224 fn=res(n1)*rwl(1)+res(n2)*rwl(2)+res(n3)*rwl(3)
225 fn=fn*dt12
226 fnxn=fn*rwl(1)
227 fnyn=fn*rwl(2)
228 fnzn=fn*rwl(3)
229 fnxt=res(n1)*dt12-fnxn
230 fnyt=res(n2)*dt12-fnyn
231 fnzt=res(n3)*dt12-fnzn
232C--- filter
233 k = 3*(i-1)+1
234 fnxt = fnxt * alpha + rwsav(k) * alphi
235 fnyt = fnyt * alpha + rwsav(k+1) * alphi
236 fnzt = fnzt * alpha + rwsav(k+2) * alphi
237C---
238 fndfn=fnxn**2+fnyn**2+fnzn**2
239 ftdft=fnxt**2+fnyt**2+fnzt**2
240 IF (fndfn==0) THEN
241 rwsav(k) = zero
242 rwsav(k+1) = zero
243 rwsav(k+2) = zero
244 ELSE
245 rwsav(k) = fnxt
246 rwsav(k+1) = fnyt
247 rwsav(k+2) = fnzt
248 ENDIF
249C---
250 IF(ftdft <= fric2*fndfn)THEN
251 ELSE
252C glissement
253 fcoe=fric*sqrt(fndfn/ftdft)
254 fnxt=fcoe*fnxt
255 fnyt=fcoe*fnyt
256 fnzt=fcoe*fnzt
257C
258C apply Ft
259 fxt=fnxt*fac
260 fyt=fnyt*fac
261 fzt=fnzt*fac
262 r(n1)=r(n1)-fxt
263 r(n2)=r(n2)-fyt
264 r(n3)=r(n3)-fzt
265C
266 frea(n1) = fxt
267 frea(n2) = fyt
268 frea(n3) = fzt
269C
270 impact(j)=-impact(j)
271 END IF
272 ENDDO
273 ELSE
274C--- no friction filtering
275 fric=rwl(13)
276 fric2=fric**2
277 fac=one/dt12
278 DO j = 1,nimpact
279 i = impact(j)
280 n=nsw(i)
281 n3=3*n
282 n2=n3-1
283 n1=n2-1
284C---
285 fn=res(n1)*rwl(1)+res(n2)*rwl(2)+res(n3)*rwl(3)
286 fn=fn*dt12
287 fnxn=fn*rwl(1)
288 fnyn=fn*rwl(2)
289 fnzn=fn*rwl(3)
290 fnxt=res(n1)*dt12-fnxn
291 fnyt=res(n2)*dt12-fnyn
292 fnzt=res(n3)*dt12-fnzn
293C
294C---
295 fndfn=fnxn**2+fnyn**2+fnzn**2
296 ftdft=fnxt**2+fnyt**2+fnzt**2
297 IF(ftdft <= fric2*fndfn)THEN
298C adherence
299 ELSE
300C glissement
301 fcoe=fric*sqrt(fndfn/ftdft)
302 fnxt=fcoe*fnxt
303 fnyt=fcoe*fnyt
304 fnzt=fcoe*fnzt
305C
306C apply Ft
307 fxt=fnxt*fac
308 fyt=fnyt*fac
309 fzt=fnzt*fac
310 r(n1)=r(n1)-fxt
311 r(n2)=r(n2)-fyt
312 r(n3)=r(n3)-fzt
313C
314 fxn=fnxn*fac
315 fyn=fnyn*fac
316 fzn=fnzn*fac
317 frea(n1) = fxt
318 frea(n2) = fyt
319 frea(n3) = fzt
320C
321 impact(j)=-impact(j)
322 END IF
323 ENDDO
324 ENDIF
325 ENDIF
326C
327 RETURN
328 END
329!||====================================================================
330!|| sms_rgwall_bcs_0 ../engine/source/ams/sms_rgwall.F
331!||--- called by ------------------------------------------------------
332!|| sms_rgwal_0 ../engine/source/ams/sms_rgwal0.F
333!||====================================================================
334 SUBROUTINE sms_rgwall_bcs_0
335 1 (x ,a ,v ,rwl ,nsw ,
336 2 nsn ,itied ,msr ,diag_sms,weight ,
337 3 rwsav ,nimpact,impact ,nsms ,nrwl_sms)
338C-----------------------------------------------
339C I m p l i c i t T y p e s
340C-----------------------------------------------
341#include "implicit_f.inc"
342#include "comlock.inc"
343C-----------------------------------------------
344C C o m m o n B l o c k s
345C-----------------------------------------------
346#include "com08_c.inc"
347C-----------------------------------------------
348C D u m m y A r g u m e n t s
349C-----------------------------------------------
350 INTEGER NSN, ITIED, MSR, NIMPACT, NSMS
351 INTEGER NSW(*), WEIGHT(*), IMPACT(*), NRWL_SMS(*)
352C REAL
353 my_real
354 . x(*), a(*), v(*), rwl(*), diag_sms(*), rwsav(*)
355C-----------------------------------------------
356C L o c a l V a r i a b l e s
357C-----------------------------------------------
358 INTEGER I, N, N3, N2, N1,IFQ, J, K, M1, M2, M3, M
359
360C REAL
361 my_real
362 . xwl, ywl, zwl, vxw, vyw, vzw,
363 . dv, da, dvt, msw, fnxn, fnyn, fnzn
364C-----------------------------------------------
365C
366C fixed wall only.
367 IF(msr==0)THEN
368 xwl=rwl(4)
369 ywl=rwl(5)
370 zwl=rwl(6)
371 vxw=zero
372 vyw=zero
373 vzw=zero
374 ELSE
375 m3=3*msr
376 m2=m3-1
377 m1=m2-1
378 vxw=v(m1)+a(m1)*dt12
379 vyw=v(m2)+a(m2)*dt12
380 vzw=v(m3)+a(m3)*dt12
381 xwl=x(m1)+vxw*dt2
382 ywl=x(m2)+vyw*dt2
383 zwl=x(m3)+vzw*dt2
384 ENDIF
385C
386 ifq = nint(rwl(15))
387C
388 IF(itied==0)THEN
389C
390 DO 40 j = 1,nimpact
391 i = impact(j)
392 n=nsw(i)
393 n3=3*n
394 n2=n3-1
395 n1=n2-1
396 dv =(v(n1)-vxw)*rwl(1)+(v(n2)-vyw)*rwl(2)+(v(n3)-vzw)*rwl(3)
397 da =a(n1)*rwl(1)+a(n2)*rwl(2)+a(n3)*rwl(3)
398 da =dv/dt12+da
399C
400 a(n1)=a(n1)-da*rwl(1)
401 a(n2)=a(n2)-da*rwl(2)
402 a(n3)=a(n3)-da*rwl(3)
403C
404 40 CONTINUE
405C
406 ELSEIF(itied==1)THEN
407C
408 DO 60 j = 1,nimpact
409 i = impact(j)
410 n=nsw(i)
411 n3=3*n
412 n2=n3-1
413 n1=n2-1
414C
415 a(n1)=-(v(n1)-vxw)/dt12
416 a(n2)=-(v(n2)-vyw)/dt12
417 a(n3)=-(v(n3)-vzw)/dt12
418 60 CONTINUE
419C
420 ELSE
421C
422C--- friction
423 DO j = 1,nimpact
424 i = abs(impact(j))
425 n=nsw(i)
426 n3=3*n
427 n2=n3-1
428 n1=n2-1
429C---
430C adherence
431 a(n1)=-(v(n1)-vxw)/dt12
432 a(n2)=-(v(n2)-vyw)/dt12
433 a(n3)=-(v(n3)-vzw)/dt12
434 ENDDO
435 ENDIF
436C
437 RETURN
438 END
439!||====================================================================
440!|| sms_rgwall_bcs_1 ../engine/source/ams/sms_rgwall.F
441!||--- called by ------------------------------------------------------
442!|| sms_rgwal_0 ../engine/source/ams/sms_rgwal0.F
443!||====================================================================
444 SUBROUTINE sms_rgwall_bcs_1
445 1 (x ,a ,v ,rwl ,nsw ,
446 2 nsn ,itied ,msr ,ms ,weight ,
447 3 rwsav ,nimpact,impact ,nsms ,nrwl_sms)
448C-----------------------------------------------
449C I m p l i c i t T y p e s
450C-----------------------------------------------
451#include "implicit_f.inc"
452#include "comlock.inc"
453C-----------------------------------------------
454C D u m m y A r g u m e n t s
455C-----------------------------------------------
456 INTEGER NSN, ITIED, MSR, NIMPACT, NSMS
457 INTEGER NSW(*), WEIGHT(*), IMPACT(*), NRWL_SMS(*)
458C REAL
459 my_real
460 . x(*), a(*), v(*), rwl(*), ms(*), rwsav(*)
461C-----------------------------------------------
462C L o c a l V a r i a b l e s
463C-----------------------------------------------
464 INTEGER I, N, N3, N2, N1,IFQ, J
465
466C REAL
467 my_real
468 . da
469C-----------------------------------------------
470C
471 ifq = nint(rwl(15))
472C
473 IF(itied==0)THEN
474C
475 DO 40 j = 1,nimpact
476 i = impact(j)
477 n=nsw(i)
478 n3=3*n
479 n2=n3-1
480 n1=n2-1
481 da =a(n1)*rwl(1)+a(n2)*rwl(2)+a(n3)*rwl(3)
482C
483 a(n1)=a(n1)-da*rwl(1)
484 a(n2)=a(n2)-da*rwl(2)
485 a(n3)=a(n3)-da*rwl(3)
486 40 CONTINUE
487C
488 ELSEIF(itied==1)THEN
489C
490 DO 60 j = 1,nimpact
491 i = impact(j)
492 n=nsw(i)
493 n3=3*n
494 n2=n3-1
495 n1=n2-1
496 a(n1)=zero
497 a(n2)=zero
498 a(n3)=zero
499 60 CONTINUE
500C
501 ELSE
502C
503C--- friction
504 DO j = 1,nimpact
505 i = abs(impact(j))
506 n=nsw(i)
507 n3=3*n
508 n2=n3-1
509 n1=n2-1
510 da=a(n1)*rwl(1)+a(n2)*rwl(2)+a(n3)*rwl(3)
511C---
512 IF(impact(j) > 0)THEN
513C adherence
514 a(n1)=zero
515 a(n2)=zero
516 a(n3)=zero
517 ELSE
518C glissement
519 a(n1)=a(n1)-da*rwl(1)
520 a(n2)=a(n2)-da*rwl(2)
521 a(n3)=a(n3)-da*rwl(3)
522 END IF
523 ENDDO
524 ENDIF
525C
526 RETURN
527 END
528!||====================================================================
529!|| sms_rgwall_bilan ../engine/source/ams/sms_rgwall.F
530!||--- called by ------------------------------------------------------
531!|| sms_rgwal_0 ../engine/source/ams/sms_rgwal0.F
532!||--- calls -----------------------------------------------------
533!|| sum_6_float ../engine/source/system/parit.F
534!||====================================================================
535 SUBROUTINE sms_rgwall_bilan
536 1 (x ,frea ,v ,rwl ,nsw ,
537 2 nsn ,itied ,msr ,ms ,weight ,
538 3 rwsav ,nimpact,impact ,nsms ,nrwl_sms,
539 4 fsav ,fopt ,frwl6 ,a ,wfext)
540C-----------------------------------------------
541C I m p l i c i t T y p e s
542C-----------------------------------------------
543#include "implicit_f.inc"
544#include "comlock.inc"
545C-----------------------------------------------
546C C o m m o n B l o c k s
547C-----------------------------------------------
548#include "com06_c.inc"
549#include "com08_c.inc"
550C-----------------------------------------------
551C D u m m y A r g u m e n t s
552C-----------------------------------------------
553 INTEGER NSN, ITIED, MSR, NIMPACT, NSMS
554 INTEGER NSW(*), WEIGHT(*), IMPACT(*), NRWL_SMS(*)
555C REAL
556 my_real
557 . x(*), v(*), rwl(*), ms(*), rwsav(*), fsav(*), frea(3,*),
558 . fopt(*), a(*)
559 DOUBLE PRECISION FRWL6(7,6)
560 DOUBLE PRECISION,INTENT(INOUT) :: WFEXT
561C-----------------------------------------------
562C L o c a l V a r i a b l e s
563C-----------------------------------------------
564 INTEGER I, N, N3, N2, N1,IFQ, J, K, M1, M2, M3
565
566 my_real
567 . vxw, vyw, vzw, vx, vy, vz,
568 . wfextt,
569 . fnxn, fnyn, fnzn, fnxt, fnyt, fnzt, fn, msw,
570 . fxn, fyn, fzn, fxt, fyt, fzt,
571 . f1(nsn), f2(nsn), f3(nsn), f4(nsn), f5(nsn), f6(nsn), f7(nsn)
572C-----------------------------------------------
573C
574C fixed wall only.
575 IF(msr==0)THEN
576 vxw=zero
577 vyw=zero
578 vzw=zero
579 ELSE
580 m3=3*msr
581 m2=m3-1
582 m1=m2-1
583C WFEXT only <=> dt12/2.
584 vxw=v(m1)+half*a(m1)*dt12
585 vyw=v(m2)+half*a(m2)*dt12
586 vzw=v(m3)+half*a(m3)*dt12
587 ENDIF
588C
589 wfextt= zero
590C
591 ifq = nint(rwl(15))
592C
593 IF(itied==0)THEN
594C
595 DO 40 j = 1,nimpact
596 i = impact(j)
597 n=nsw(i)
598C
599 fn=frea(1,n)*rwl(1)+frea(2,n)*rwl(2)+frea(3,n)*rwl(3)
600 fn=weight(n)*fn
601 fxn=fn*rwl(1)
602 fyn=fn*rwl(2)
603 fzn=fn*rwl(3)
604C
605 f1(j) = fxn
606 f2(j) = fyn
607 f3(j) = fzn
608 f4(j) = ms(n)
609 f5(j) = zero
610 f6(j) = zero
611 f7(j) = zero
612C
613C 1er impact (WFEXT avec decalage 1/2 cycle)
614c N3=3*N
615c N2=N3-1
616c N1=N2-1
617c VX=V(N1)+HALF*A(N1)*DT12
618c VY=V(N2)+HALF*A(N2)*DT12
619c VZ=V(N3)+HALF*A(N3)*DT12
620c WFEXTT = WFEXTT -DT12*((VX-VXW)*FXN+(VY-VYW)*FYN+(VZ-VZW)*FZN)
621 40 CONTINUE
622C
623 ELSE
624C
625 DO 60 j = 1,nimpact
626 i = abs(impact(j))
627 n=nsw(i)
628C
629 fn=frea(1,n)*rwl(1)+frea(2,n)*rwl(2)+frea(3,n)*rwl(3)
630 fn=weight(n)*fn
631 fxn=fn*rwl(1)
632 fyn=fn*rwl(2)
633 fzn=fn*rwl(3)
634C
635 f1(j) = fxn
636 f2(j) = fyn
637 f3(j) = fzn
638 f4(j) = ms(n)
639C
640 n3=3*n
641 n2=n3-1
642 n1=n2-1
643 vx=v(n1)+half*a(n1)*dt12
644 vy=v(n2)+half*a(n2)*dt12
645 vz=v(n3)+half*a(n3)*dt12
646C 1er impact (WFEXT avec decalage 1/2 cycle)
647c WFEXTT = WFEXTT -DT12*((VX-VXW)*FXN+(VY-VYW)*FYN+(VZ-VZW)*FZN)
648C
649 fxt=weight(n)*frea(1,n)-fxn
650 fyt=weight(n)*frea(2,n)-fyn
651 fzt=weight(n)*frea(3,n)-fzn
652 f5(j) = fxt
653 f6(j) = fyt
654 f7(j) = fzt
655 wfextt = wfextt -dt12*((vx-vxw)*fxt+(vy-vyw)*fyt+(vz-vzw)*fzt)
656C
657 60 CONTINUE
658 ENDIF
659C
660#include "lockon.inc"
661 wfext= wfext+wfextt
662#include "lockoff.inc"
663C
664 CALL sum_6_float(1, nimpact, f1, frwl6(1,1), 7)
665 CALL sum_6_float(1, nimpact, f2, frwl6(2,1), 7)
666 CALL sum_6_float(1, nimpact, f3, frwl6(3,1), 7)
667 CALL sum_6_float(1, nimpact, f4, frwl6(4,1), 7)
668 CALL sum_6_float(1, nimpact, f5, frwl6(5,1), 7)
669 CALL sum_6_float(1, nimpact, f6, frwl6(6,1), 7)
670 CALL sum_6_float(1, nimpact, f7, frwl6(7,1), 7)
671C
672 RETURN
673 END
my_real
#define my_real
Definition cppsort.cpp:32
alpha
#define alpha
Definition eval.h:35
sum_6_float
subroutine sum_6_float(jft, jlt, f, f6, n)
Definition parit.F:64
sms_rgwall_bcs_0
subroutine sms_rgwall_bcs_0(x, a, v, rwl, nsw, nsn, itied, msr, diag_sms, weight, rwsav, nimpact, impact, nsms, nrwl_sms)
Definition sms_rgwall.F:338
sms_rgwall_bcs_1
subroutine sms_rgwall_bcs_1(x, a, v, rwl, nsw, nsn, itied, msr, ms, weight, rwsav, nimpact, impact, nsms, nrwl_sms)
Definition sms_rgwall.F:448
sms_rgwall_fric
subroutine sms_rgwall_fric(x, a, v, rwl, nsw, nsn, itied, msr, ms, weight, rwsav, nimpact, impact, nsms, nrwl_sms, fsav, fopt, res, r, frea)
Definition sms_rgwall.F:155
sms_rgwall_bilan
subroutine sms_rgwall_bilan(x, frea, v, rwl, nsw, nsn, itied, msr, ms, weight, rwsav, nimpact, impact, nsms, nrwl_sms, fsav, fopt, frwl6, a, wfext)
Definition sms_rgwall.F:540
sms_rgwall_impact
subroutine sms_rgwall_impact(x, a, v, rwl, nsw, nsn, itied, msr, ms, weight, rwsav, nimpact, impact, nsms, nrwl_sms)
Definition sms_rgwall.F:32