OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
iqela1.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!|| iqela1 ../engine/source/ale/inter/iqela1.F
25!||--- called by ------------------------------------------------------
26!|| intal2 ../engine/source/ale/inter/intal2.F
27!||--- calls -----------------------------------------------------
28!|| shapeh ../engine/source/ale/inter/shapeh.F
29!||--- uses -----------------------------------------------------
30!|| anim_mod ../common_source/modules/output/anim_mod.F
31!|| h3d_mod ../engine/share/modules/h3d_mod.F
32!||====================================================================
33 SUBROUTINE iqela1(A,E,MSM,IRECT,CRST,MSR,
34 1 NSV,ILOC,IRTL,MS,NOR,LCODE,ISKEW,FSAV,
35 2 FCONT,FNCONT,H3D_DATA,NSN,NMN)
36C-----------------------------------------------
37C M o d u l e s
38C-----------------------------------------------
39 USE h3d_mod
40 USE anim_mod
41C-----------------------------------------------
42C D e s c r i p t i o n
43C-----------------------------------------------
44C This subroutine is related to option /INTER/TYPE1
45C and computes forces. See E(1:3)=FSN.N(1:3)
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C C o m m o n B l o c k s
52C-----------------------------------------------
53#include "com01_c.inc"
54#include "com04_c.inc"
55#include "com06_c.inc"
56#include "com08_c.inc"
57#include "scr07_c.inc"
58#include "scr08_a_c.inc"
59#include "scr14_c.inc"
60#include "scr16_c.inc"
61#include "param_c.inc"
62#include "comlock.inc"
63#include "tabsiz_c.inc"
64C-----------------------------------------------
65C D u m m y A r g u m e n t s
66C-----------------------------------------------
67 INTEGER,INTENT(INOUT) :: IRECT(4,*), MSR(*), NSV(*), ILOC(*), IRTL(*),LCODE(*), ISKEW(*)
68 my_real,INTENT(INOUT) :: A(SA), E(*), MSM(*), CRST(2,*), MS(*),NOR(3,*),FSAV(*)
69 my_real,INTENT(INOUT) :: fcont(3,*),fncont(3,*)
70 INTEGER, INTENT(IN) :: NSN,NMN
71 TYPE(h3d_database) :: H3D_DATA
72C-----------------------------------------------
73C L o c a l V a r i a b l e s
74C-----------------------------------------------
75 INTEGER NIR, I, J, I3, J3, I2, J2, I1, J1, II, L, JJ, NN, JJ3, JJ2, JJ1, ISK, LCOD
76 my_real H(4), N1, N2, N3, AA(3), SSS, TTT, XMSS, FXI, FYI, FZI, FSN
77 my_real :: fsn_sav, fxi_sav, fyi_sav, fzi_sav,impx,impy,impz
78 LOGICAL ICONT, IPCONT, IANIM
79C-----------------------------------------------
80C S o u r c e L i n e s
81C-----------------------------------------------
82 icont = .false.
83 ipcont = .false.
84 ianim = .false.
85 icont = (anim_v(4)+outp_v(4) > 0+h3d_data%N_VECT_CONT)
86 IF(anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT > 0)THEN
87 IF( (tt>=tanim .AND. tt<=tanim_stop) .OR.tt >= toutp.OR.tt >= h3d_data%TH3D.OR.
88 . (manim >= 4.AND.manim <= 15).OR. h3d_data%MH3D /= 0)THEN
89 ipcont = .true.
90 ENDIF
91 ENDIF
92 IF(icont .OR. ipcont)ianim=.true.
93
94 fsn_sav = zero
95 fxi_sav = zero
96 fyi_sav = zero
97 fzi_sav = zero
98
99 nir=2
100 IF(n2d == 0)nir=4
101 !main nodes
102 DO i=1,nmn
103 j=msr(i)
104 i3=3*i
105 i2=i3-1
106 i1=i2-1
107 msm(i)=ms(j)
108 e(i1)=zero
109 e(i2)=zero
110 e(i3)=zero
111 ENDDO
112
113 !secnd nodes
114 DO ii=1,nsn
115 i=nsv(ii)
116 j=iloc(ii)
117 IF(j >= 1) THEN
118 l=irtl(ii)
119 DO jj=1,nir
120 nn=irect(jj,l)
121 iy(jj)=nn
122 ENDDO
123 !parametric coordinates on main face
124 sss=crst(1,ii)
125 ttt=crst(2,ii)
126 !normal
127 n1=nor(1,ii)
128 n2=nor(2,ii)
129 n3=nor(3,ii)
130 !A(1:3,I) <-> A(3*I - 1:3)
131 i3=3*i
132 i2=i3-1
133 i1=i2-1
134 CALL shapeh(h,sss,ttt)
135 DO jj=1,nir
136 j3=3*iy(jj)
137 j2=j3-1
138 j1=j2-1
139 jj3=3*msr(iy(jj))
140 jj2=jj3-1
141 jj1=jj2-1
142 aa(1)=a(i1) !Accel for itab(I), where I=NSV(II)
143 aa(2)=a(i2)
144 aa(3)=a(i3)
145 isk=iskew(iy(jj))
146 lcod=lcode(iy(jj))
147 xmss=ms(i)*h(jj)
148 fxi=aa(1)-a(jj1)
149 fyi=aa(2)-a(jj2)
150 fzi=aa(3)-a(jj3)
151 fsn=(fxi*n1+fyi*n2+fzi*n3)*xmss
152 fsn_sav = fsn_sav + fsn
153 fxi_sav = fxi_sav + fxi*xmss
154 fyi_sav = fyi_sav + fyi*xmss
155 fzi_sav = fzi_sav + fzi*xmss
156 e(j1)=e(j1)+fsn*n1
157 e(j2)=e(j2)+fsn*n2
158 e(j3)=e(j3)+fsn*n3
159 msm(iy(jj))=msm(iy(jj))+xmss
160 enddo!next JJ
161 ENDIF
162 enddo!next II
163
164 !SUM(E(1:3,*)) is NF(1:3) on main segment
165
166 !---------------------------------
167 ! /TH/INTER
168 ! NORMAL IMPULSE BACKUP
169 !---------------------------------
170 impx = fxi_sav*dt12
171 impy = fyi_sav*dt12
172 impz = fzi_sav*dt12
173 fsn_sav = fsn_sav*dt12
174#include "lockon.inc"
175 fsav(1)=fsav(1) + impx
176 fsav(2)=fsav(2) + impy
177 fsav(3)=fsav(3) + impz
178 fsav(8)=fsav(8) + abs(impx)
179 fsav(9)=fsav(9) + abs(impy)
180 fsav(10)=fsav(10)+ abs(impz)
181 fsav(11)=fsav(11)+ fsn_sav
182#include "lockoff.inc"
183 !---------------------------------
184
185
186 IF(ianim)THEN
187 !---------------------------------
188 ! /ANIM/VECT/CONT
189 !---------------------------------
190 IF(icont)THEN
191#include "lockon.inc"
192 DO i=1,nmn
193 j = msr(i)
194 i3 = 3*i
195 i2 = i3-1
196 i1 = i2-1
197 fcont(1,j) = fcont(1,j)+e(i1)
198 fcont(2,j) = fcont(2,j)+e(i2)
199 fcont(3,j) = fcont(3,j)+e(i3)
200 ENDDO
201#include "lockoff.inc"
202 ENDIF
203 !---------------------------------
204 ! /ANIM/VECT/PCONT
205 !---------------------------------
206 IF(ipcont)THEN
207#include "lockon.inc"
208 DO i=1,nmn
209 j = msr(i)
210 i3 = 3*i
211 i2 = i3-1
212 i1 = i2-1
213 fncont(1,j) = fncont(1,j)+e(i1)
214 fncont(2,j) = fncont(2,j)+e(i2)
215 fncont(3,j) = fncont(3,j)+e(i3)
216 ENDDO
217#include "lockoff.inc"
218 ENDIF
219 ENDIF
220
221
222
223C-----------------------------------------------------
224
225 RETURN
226 END
iqela1
subroutine iqela1(a, e, msm, irect, crst, msr, nsv, iloc, irtl, ms, nor, lcode, iskew, fsav, fcont, fncont, h3d_data, nsn, nmn)
Definition iqela1.F:36
anim_mod
Definition anim_mod.F:96
h3d_mod
Definition h3d_mod.F:308
shapeh
subroutine shapeh(h, s, t)
Definition shapeh.F:34
h3d_mod::h3d_database
Definition h3d_mod.F:435