OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i15ass.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!|| i15ass ../engine/source/interfaces/int15/i15ass.F
25!||--- called by ------------------------------------------------------
26!|| i15cmp ../engine/source/interfaces/int15/i15cmp.F
27!||--- uses -----------------------------------------------------
28!|| anim_mod ../common_source/modules/output/anim_mod.F
29!|| groupdef_mod ../common_source/modules/groupdef_mod.F
30!|| h3d_mod ../engine/share/modules/h3d_mod.F
31!|| output_mod ../common_source/modules/output/output_mod.F90
32!||====================================================================
33 SUBROUTINE i15ass(OUTPUT,AF ,X ,V ,KSURF ,IGRSURF ,
34 2 BUFSF ,STIFN ,FS ,FCONT ,FSKYI ,
35 3 ISKY ,DE ,WNF ,WTF ,WNS ,
36 4 FNORMX ,FNORMY ,FNORMZ ,FTANGX ,FTANGY ,
37 5 FTANGZ ,NNC ,KNC ,H3D_DATA)
38C-----------------------------------------------
39C M o d u l e s
40C-----------------------------------------------
41 USE h3d_mod
42 USE groupdef_mod
43 USE anim_mod
44 USE output_mod, ONLY : output_
45C-----------------------------------------------
46C I m p l i c i t T y p e s
47C-----------------------------------------------
48#include "implicit_f.inc"
49#include "comlock.inc"
50C-----------------------------------------------
51C C o m m o n B l o c k s
52C-----------------------------------------------
53#include "com04_c.inc"
54#include "com06_c.inc"
55#include "com08_c.inc"
56#include "scr07_c.inc"
57#include "scr14_c.inc"
58#include "scr16_c.inc"
59#include "parit_c.inc"
60#include "param_c.inc"
61#include "scr18_c.inc"
62C-----------------------------------------------
63C D u m m y A r g u m e n t s
64C-----------------------------------------------
65 TYPE(output_), INTENT(INOUT) :: OUTPUT
66 INTEGER KSURF,ISKY(*),NNC,KNC(*)
67 my_real
68 . AF(*) , X(3,*), V(3,*),BUFSF(*),
69 . stifn(*), fs(nthvki),
70 . fcont(3,*),fskyi(lskyi,nfskyi), de,
71 . wnf(3,*) ,wtf(3,*) ,wns(*) ,
72 . fnormx,fnormy,fnormz,ftangx,ftangy,ftangz
73 TYPE(h3d_database) :: H3D_DATA
74 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
75C-----------------------------------------------
76C L o c a l V a r i a b l e s
77C-----------------------------------------------
78 INTEGER ADRBUF, I, IN, I3, I2, I1
79 INTEGER NFORC , NISKYL
80 INTEGER NDEB, NREST
81 my_real
82 . rot(9), xk , yk, zk, fx, fy, fz, am1, am2, am3,
83 . fn1, fn2, fn3, ft1, ft2, ft3,
84 . stf, dd
85C-----------------------------------------------
86 adrbuf=igrsurf(ksurf)%IAD_BUFR
87 DO i=1,9
88 rot(i)=bufsf(adrbuf+7+i-1)
89 END DO
90C---------------------------------
91C SORTIES dans le rep. global.
92C---------------------------------
93 fn1=rot(1)*fnormx+rot(4)*fnormy+rot(7)*fnormz
94 fn2=rot(2)*fnormx+rot(5)*fnormy+rot(8)*fnormz
95 fn3=rot(3)*fnormx+rot(6)*fnormy+rot(9)*fnormz
96 fs(1)=fs(1)-fn1*dt1
97 fs(2)=fs(2)-fn2*dt1
98 fs(3)=fs(3)-fn3*dt1
99 ft1=rot(1)*ftangx+rot(4)*ftangy+rot(7)*ftangz
100 ft2=rot(2)*ftangx+rot(5)*ftangy+rot(8)*ftangz
101 ft3=rot(3)*ftangx+rot(6)*ftangy+rot(9)*ftangz
102 fs(4)=fs(4)-ft1*dt1
103 fs(5)=fs(5)-ft2*dt1
104 fs(6)=fs(6)-ft3*dt1
105C------------------------------------------------------------
106C RETOUR EN GLOBAL.
107C------------------------------------------------------------
108#include "vectorize.inc"
109 DO i=1,nnc
110 in=knc(i)
111 fx=wnf(1,in)+wtf(1,in)
112 fy=wnf(2,in)+wtf(2,in)
113 fz=wnf(3,in)+wtf(3,in)
114 wnf(1,in)=rot(1)*fx+rot(4)*fy+rot(7)*fz
115 wnf(2,in)=rot(2)*fx+rot(5)*fy+rot(8)*fz
116 wnf(3,in)=rot(3)*fx+rot(6)*fy+rot(9)*fz
117 ENDDO
118C------------------------------------------------------------
119C ASSEMBLAGE AU NOEUD MAIN DE LA SURFACE.
120C------------------------------------------------------------
121 DO i=1,nnc
122 in=knc(i)
123 xk=x(1,in)-bufsf(adrbuf+16)
124 yk=x(2,in)-bufsf(adrbuf+17)
125 zk=x(3,in)-bufsf(adrbuf+18)
126 fx =wnf(1,in)
127 fy =wnf(2,in)
128 fz =wnf(3,in)
129 am1=yk*fz-zk*fy
130 am2=zk*fx-xk*fz
131 am3=xk*fy-yk*fx
132C-----
133 bufsf(adrbuf+25)=bufsf(adrbuf+25)-fx
134 bufsf(adrbuf+26)=bufsf(adrbuf+26)-fy
135 bufsf(adrbuf+27)=bufsf(adrbuf+27)-fz
136 bufsf(adrbuf+28)=bufsf(adrbuf+28)-am1
137 bufsf(adrbuf+29)=bufsf(adrbuf+29)-am2
138 bufsf(adrbuf+30)=bufsf(adrbuf+30)-am3
139C-----
140 stf=wns(in)
141 bufsf(adrbuf+31)=bufsf(adrbuf+31)+stf
142 dd = xk**2+yk**2+zk**2
143 bufsf(adrbuf+32)=bufsf(adrbuf+32)+dd*stf
144 ENDDO
145C---------------------------------
146C ASSEMBLAGE des FORCES aux NOEUDS SECONDS.
147C---------------------------------
148 IF (iparit/=0) THEN
149#include "lockon.inc"
150 niskyl = nisky
151 nisky = nisky+nnc
152#include "lockoff.inc"
153 END IF
154 IF (iparit==0) THEN
155#include "vectorize.inc"
156 DO 300 i=1,nnc
157 in=knc(i)
158 fx=wnf(1,in)
159 fy=wnf(2,in)
160 fz=wnf(3,in)
161 i3=3*in
162 i2=i3-1
163 i1=i2-1
164 af(i1)=af(i1)+fx
165 af(i2)=af(i2)+fy
166 af(i3)=af(i3)+fz
167 stifn(in)=stifn(in)+wns(in)
168 300 CONTINUE
169 ELSE
170 IF(kdtint==0)THEN
171 DO 350 i=1,nnc
172 in=knc(i)
173 fx=wnf(1,in)
174 fy=wnf(2,in)
175 fz=wnf(3,in)
176 niskyl = niskyl + 1
177 fskyi(niskyl,1)=fx
178 fskyi(niskyl,2)=fy
179 fskyi(niskyl,3)=fz
180 fskyi(niskyl,4)=wns(in)
181 isky(niskyl) =in
182 350 CONTINUE
183 ELSE
184 DO i=1,nnc
185 in=knc(i)
186 fx=wnf(1,in)
187 fy=wnf(2,in)
188 fz=wnf(3,in)
189 niskyl = niskyl + 1
190 fskyi(niskyl,1)=fx
191 fskyi(niskyl,2)=fy
192 fskyi(niskyl,3)=fz
193 fskyi(niskyl,4)=wns(in)
194 fskyi(niskyl,5)=zero
195 isky(niskyl) =in
196 ENDDO
197 ENDIF
198 ENDIF
199C------------------------------------------------------------
200C ANIM (FORCES DE CONTACT).
201C------------------------------------------------------------
202 IF(anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT>0.AND.
203 . ((tt>=tanim .AND. tt<=tanim_stop).OR.tt>=toutp.OR.(tt>=h3d_data%TH3D.AND.tt<=h3d_data%TH3D_STOP).OR.
204 . (manim>=4.AND.manim<=15).OR. h3d_data%MH3D /= 0))THEN
205#include "lockon.inc"
206#include "vectorize.inc"
207 DO 400 i=1,nnc
208 in=knc(i)
209 fcont(1,in) =fcont(1,in) + wnf(1,in)
210 fcont(2,in) =fcont(2,in) + wnf(2,in)
211 fcont(3,in) =fcont(3,in) + wnf(3,in)
212 400 CONTINUE
213#include "lockoff.inc"
214 ENDIF
215C---------------------------------
216C Pour Travail des forces sur noeuds seconds
217C 1ere partie : ici
218C 2eme partie : apres calcul de DT2.
219C---------------------------------
220 DO 450 i=1,nnc
221 in=knc(i)
222 fx=wnf(1,in)
223 fy=wnf(2,in)
224 fz=wnf(3,in)
225 de=de+fx*v(1,in)+fy*v(2,in)+fz*v(3,in)
226 450 CONTINUE
227C---------------------------------
228C Working force at interface (Madymo)
229C---------------------------------
230 fs(7)=fs(7)+de*dt1*half
231 IF (igrsurf(ksurf)%TYPE==100) THEN
232C Madymo Ellipsoids
233!$OMP ATOMIC
234 output%TH%WFEXT=output%TH%WFEXT+de*dt1*half
235 ENDIF
236C----------------------------------
237 RETURN
238 END
i15ass
subroutine i15ass(output, af, x, v, ksurf, igrsurf, bufsf, stifn, fs, fcont, fskyi, isky, de, wnf, wtf, wns, fnormx, fnormy, fnormz, ftangx, ftangy, ftangz, nnc, knc, h3d_data)
Definition i15ass.F:38
anim_mod
Definition anim_mod.F:96
groupdef_mod
Definition groupdef_mod.F:662
h3d_mod
Definition h3d_mod.F:308
h3d_mod::h3d_database
Definition h3d_mod.F:435