OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
iqela2.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!|| iqela2 ../engine/source/ale/inter/iqela2.F
25!||--- called by ------------------------------------------------------
26!|| intal2 ../engine/source/ale/inter/intal2.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| bcs2 ../engine/source/constraints/general/bcs/bcs2.F
30!|| shapeh ../engine/source/ale/inter/shapeh.F
31!||--- uses -----------------------------------------------------
32!|| message_mod ../engine/share/message_module/message_mod.F
33!||====================================================================
34 SUBROUTINE iqela2(SKEW,A,E,MSM,IRECT,CRST,MSR,NSV,IRTL,LCODE,V,ISKEW,NOR,ITAB,NSN,NMN)
35C-----------------------------------------------
36C M o d u l e s
37C-----------------------------------------------
38 USE message_mod
39C-----------------------------------------------
40C I m p l i c i t T y p e s
41C-----------------------------------------------
42#include "implicit_f.inc"
43C-----------------------------------------------
44C C o m m o n B l o c k s
45C-----------------------------------------------
46#include "com01_c.inc"
47#include "com04_c.inc"
48#include "scr08_a_c.inc"
49#include "param_c.inc"
50#include "warn_c.inc"
51#include "tabsiz_c.inc"
52C-----------------------------------------------
53C D u m m y A r g u m e n t s
54C-----------------------------------------------
55 INTEGER IRECT(4,*), MSR(*), NSV(*), IRTL(*), LCODE(*), ISKEW(*), ITAB(*)
56 my_real skew(lskew,*), a(sa), e(*), msm(*), crst(2,*), v(sv),nor(3,*)
57 INTEGER, INTENT(IN) :: NSN,NMN
58C-----------------------------------------------
59C L o c a l V a r i a b l e s
60C-----------------------------------------------
61 INTEGER JBC(3), NIR, I, J, I3, J3, I2, J2, I1, J1, ISK, LCOD, II, L, JJ, NN, LK, IBC
62 my_real xn(3), yn(3), zn(3), h(4), n1, n2, n3, ss, tt, ax, ay, az,
63 . vx, vy, vz, amn, vmn, amod, vmod, bvz, baz, bvx, bvy, bax, bay,
64 . a11, a12, a13, a21, a22, a23, a31, a32, a33, det
65C-----------------------------------------------
66C S o u r c e L i n e s
67C-----------------------------------------------
68 nir=2
69 IF(n2d == 0)nir=4
70C-------------------------------------
71C ACCELERATION OF MAIN NODES
72C-------------------------------------
73 DO i=1,nmn
74 j=msr(i)
75 i3=3*i
76 j3=3*j
77 i2=i3-1
78 j2=j3-1
79 i1=i2-1
80 j1=j2-1
81 IF(msm(i) > zero)THEN
82 a(j1)=a(j1)+e(i1)/msm(i)
83 a(j2)=a(j2)+e(i2)/msm(i)
84 a(j3)=a(j3)+e(i3)/msm(i)
85 ENDIF
86 isk = iskew(j)
87 lcod = lcode(j)
88 IF(lcod>0)CALL bcs2(a(j1),skew(1,isk),isk,lcod)
89 ENDDO
90C------------------------------------
91C ACCELERATION ON SECND NODES
92C------------------------------------
93 DO ii=1,nsn
94 i=nsv(ii)
95 l=irtl(ii)
96 DO jj=1,nir
97 nn=irect(jj,l)
98 ix(jj)=msr(nn)
99 ENDDO !NEXT JJ
100 ss=crst(1,ii)
101 tt=crst(2,ii)
102 n1=nor(1,ii)
103 n2=nor(2,ii)
104 n3=nor(3,ii)
105 CALL shapeh(h,ss,tt)
106 i3=3*i
107 i2=i3-1
108 i1=i2-1
109 ax=zero
110 ay=zero
111 az=zero
112 vx=zero
113 vy=zero
114 vz=zero
115 DO jj=1,nir
116 j3=3*ix(jj)
117 j2=j3-1
118 j1=j2-1
119 ax=ax+a(j1)*h(jj)
120 ay=ay+a(j2)*h(jj)
121 az=az+a(j3)*h(jj)
122 vx=vx+v(j1)*h(jj)
123 vy=vy+v(j2)*h(jj)
124 vz=vz+v(j3)*h(jj)
125 ENDDO !NEXT JJ
126
127 amn = n1*ax+n2*ay+n3*az
128 vmn = n1*vx+n2*vy+n3*vz
129 amod= amn-n1*a(i1)-n2*a(i2)-n3*a(i3)
130 vmod= vmn-n1*v(i1)-n2*v(i2)-n3*v(i3)
131
132 lcod=lcode(i) ! LCOD \in [1,7] (3-bit-integer representing XYZ tags)
133
134 IF(lcod /= 0)THEN
135 !--------------------------------
136 ! IN CASE ON BOUNDARY CONDITION DEFINED ON THIS NODE
137 !--------------------------------
138 xn(1)=n1
139 yn(1)=n2
140 zn(1)=n3
141
142 jbc(1:3) = 0
143 IF(lcode(i) > 0) THEN ! LCODE(I) is a 3-bit-integer representing XYZ tags
144 jbc(3) = iand(lcode(i), 1) !1st bit
145 jbc(2) = iand(lcode(i), 2) !2nd bit
146 jbc(1) = iand(lcode(i), 4) !3rd bit
147 ENDIF
148
149 lk=iskew(i)
150 ibc=2
151
152 IF(jbc(1) /= 0)THEN
153 xn(ibc)=skew(1,lk)
154 yn(ibc)=skew(2,lk)
155 zn(ibc)=skew(3,lk)
156 ibc=ibc+1
157 ENDIF
158 IF(jbc(2) /= 0)THEN
159 xn(ibc)=skew(4,lk)
160 yn(ibc)=skew(5,lk)
161 zn(ibc)=skew(6,lk)
162 ibc=ibc+1
163 ENDIF
164 IF(jbc(3) /= 0)THEN
165 IF(ibc == 4)THEN
166 !BCS check also in Starter, but BCS may be updated in Engine file (/BCS, /BCSR)
167 CALL ancmsg(msgid=11,anmode=aninfo,i1=itab(i))
168 ierr=1
169 EXIT
170 ELSE
171 xn(ibc)=skew(7,lk)
172 yn(ibc)=skew(8,lk)
173 zn(ibc)=skew(9,lk)
174 ibc=ibc+1
175 ENDIF
176 ENDIF
177 IF(ibc == 3)THEN
178 !C---------------------------
179 !C SINGLE CONDITION ON NODE
180 !C---------------------------
181 xn(3)=yn(1)*zn(2)-zn(1)*yn(2)
182 yn(3)=zn(1)*xn(2)-xn(1)*zn(2)
183 zn(3)=xn(1)*yn(2)-yn(1)*xn(2)
184 bvz=v(i1)*xn(3)+v(i2)*yn(3)+v(i3)*zn(3)
185 baz=a(i1)*xn(3)+a(i2)*yn(3)+a(i3)*zn(3)
186 ELSE
187 !-----------------------------
188 ! DOUBLE CONDITION ON NODE
189 !-----------------------------
190 bvz=zero
191 baz=zero
192 ENDIF
193
194 bvx=vmn
195 bvy=zero
196 bax=amn
197 bay=zero
198
199 a11=yn(2)*zn(3)-zn(2)*yn(3)
200 a12=zn(2)*xn(3)-xn(2)*zn(3)
201 a13=xn(2)*yn(3)-yn(2)*xn(3)
202 a21=yn(3)*zn(1)-zn(3)*yn(1)
203 a22=zn(3)*xn(1)-xn(3)*zn(1)
204 a23=xn(3)*yn(1)-yn(3)*xn(1)
205 a31=yn(1)*zn(2)-zn(1)*yn(2)
206 a32=zn(1)*xn(2)-xn(1)*zn(2)
207 a33=xn(1)*yn(2)-yn(1)*xn(2)
208 det=xn(1)*a11+yn(1)*a12+zn(1)*a13
209
210 !-----------------------------------
211 ! VELOCITY AND ACCELERATION
212 !-----------------------------------
213 IF(det /= zero)THEN
214 v(i1)=(a11*bvx+a21*bvy+a31*bvz)/det
215 v(i2)=(a12*bvx+a22*bvy+a32*bvz)/det
216 v(i3)=(a13*bvx+a23*bvy+a33*bvz)/det
217 a(i1)=(a11*bax+a21*bay+a31*baz)/det
218 a(i2)=(a12*bax+a22*bay+a32*baz)/det
219 a(i3)=(a13*bax+a23*bay+a33*baz)/det
220 ENDIF
221
222 ELSEIF(lcod == 0)THEN
223 !--------------------------------
224 ! NO BOUNDARY CONDITION ON NODE
225 !--------------------------------
226 a(i1)=a(i1)+amod*n1
227 a(i2)=a(i2)+amod*n2
228 a(i3)=a(i3)+amod*n3
229 v(i1)=v(i1)+vmod*n1
230 v(i2)=v(i2)+vmod*n2
231 v(i3)=v(i3)+vmod*n3
232 ENDIF
233
234 ENDDO !next II
235C--------------------------------
236 RETURN
237 END
bcs2
subroutine bcs2(a, b, j, k)
Definition bcs2.F:32
my_real
#define my_real
Definition cppsort.cpp:32
iqela2
subroutine iqela2(skew, a, e, msm, irect, crst, msr, nsv, irtl, lcode, v, iskew, nor, itab, nsn, nmn)
Definition iqela2.F:35
message_mod
Definition message_mod.F:1249
shapeh
subroutine shapeh(h, s, t)
Definition shapeh.F:34
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889