OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
shell_rota.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!|| shell_rota ../engine/source/output/dynain/shell_rota.F
25!||--- called by ------------------------------------------------------
26!|| dynain_c_strsg ../engine/source/output/dynain/dynain_c_strsg.F
27!||--- uses -----------------------------------------------------
28!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
29!||====================================================================
30 SUBROUTINE shell_rota(I ,ILAY ,NEL ,IORTH ,ITY ,
31 2 IGTYP ,MLW ,JDIR ,SIG ,ELBUF_STR,
32 3 RX ,RY ,RZ ,SX ,SY ,
33 4 SZ ,E1X ,E2X , E3X ,E1Y ,
34 5 E2Y ,E3Y ,E1Z , E2Z ,E3Z ,
35 6 DIR_A ,DIR_B )
36C-----------------------------------------------
37C M o d u l e s
38C-----------------------------------------------
39 USE elbufdef_mod
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C G l o b a l P a r a m e t e r s
46C-----------------------------------------------
47#include "mvsiz_p.inc"
48C-----------------------------------------------
49C D u m m y A r g u m e n t s
50C-----------------------------------------------
51 INTEGER NEL, ILAY,I , NS1,ITY, IORTH, IGTYP,MLW ,JDIR
52 my_real SIG(6)
53 my_real
54 . RX(MVSIZ), RY(MVSIZ), RZ(MVSIZ),
55 . SX(MVSIZ), SY(MVSIZ), SZ(MVSIZ),
56 . e1x(mvsiz), e1y(mvsiz), e1z(mvsiz),
57 . e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
58 . e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),
59 . dir_a(*) ,dir_b(*)
60
61 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
62C-----------------------------------------------
63C L o c a l V a r i a b l e s
64C-----------------------------------------------
65 INTEGER ILAW , II
66 my_real
67 . dir(2),
68 . v1,v2,v3,vr,vs,aa,bb,suma,
69 . r1,s1,s2,r2,rs1,rs2,r12a,r22a,s12b,s22b,
70 . rs3,r3r3,s3s3,t1,t2,t3,s3,s4,s5,
71 . l11,l12,l13,l22,l23,l33,
72 . r11,r12,r13,r21,r22,r23,r31,r32,r33,
73 . s11,s12,s21,s13,s31,s22,s23,s32,s33
74 my_real,
75 . DIMENSION(:) , POINTER :: dir1
76C=======================================================================
77C*******************************************
78C Rotate to global system
79C*******************************************
80
81 IF(ilay > 0 ) THEN
82 ilaw = elbuf_str%BUFLY(ilay)%ILAW
83 ELSE
84 ilaw = mlw
85 ENDIF
86C-----------------------------------
87C Rotation orthotrop > Element
88C------------------------------------
89 IF (iorth > 0) THEN
90
91 IF (igtyp == 16) THEN
92 ii = jdir + i-1
93 r1 = dir_a(ii)
94 s1 = dir_a(ii+nel)
95 r2 = dir_b(ii)
96 s2 = dir_b(ii+nel)
97
98 rs1= r1*s1
99 rs2= r2*s2
100 r12a = r1*r1
101 r22a = r2*r2
102 s12b = s1*s1
103 s22b = s2*s2
104
105 rs3 = s1*s2-r1*r2
106 r3r3= one+s1*r2+r1*s2
107 r3r3= half*r3r3
108 s3s3= one-s1*r2-r1*s2
109 s3s3= half*s3s3
110 t1 = sig(1)
111 t2 = sig(2)
112 t3 = sig(3)
113 sig(1) = r12a*t1 + r22a*t2 - rs3*t3
114 sig(2) = s12b*t1 + s22b*t2 + rs3*t3
115 sig(3) = rs1*t1 + rs2*t2 + (r3r3 - s3s3)*t3
116c
117 ELSEIF ((igtyp == 51 .OR. igtyp == 52) .AND. ilaw == 58) THEN
118 ii = jdir + i-1
119 r1 = dir_a(ii)
120 s1 = dir_a(ii+nel)
121 r2 = dir_b(ii)
122 s2 = dir_b(ii+nel)
123c
124 rs1= r1*s1
125 rs2= r2*s2
126 r12a = r1*r1
127 r22a = r2*r2
128 s12b = s1*s1
129 s22b = s2*s2
130 rs3 = s1*s2-r1*r2
131 r3r3= one+s1*r2+r1*s2
132 r3r3= half*r3r3
133 s3s3= one-s1*r2-r1*s2
134 s3s3= half*s3s3
135 t1 = sig(1)
136 t2 = sig(2)
137 t3 = sig(3)
138c
139 sig(1) = r12a*t1 + r22a*t2 - rs3*t3
140 sig(2) = s12b*t1 + s22b*t2 + rs3*t3
141 sig(3) = rs1*t1 + rs2*t2 + (r3r3 - s3s3)*t3
142 ELSE
143 IF (ilaw /= 1 .and. ilaw /= 2 .and. ilaw /= 19 .and. ilaw /= 27 .and. ilaw /= 32) THEN
144 ii = jdir + i-1
145 dir(1) = dir_a(ii)
146 dir(2) = dir_a(ii+nel)
147 s1 = dir(1)*dir(1)*sig(1)
148 . + dir(2)*dir(2)*sig(2)
149 . -two*dir(1)*dir(2)*sig(3)
150 s2 = dir(2)*dir(2)*sig(1)
151 . + dir(1)*dir(1)*sig(2)
152 . +two*dir(2)*dir(1)*sig(3)
153 s3 = dir(1)*dir(2)*sig(1)
154 . - dir(2)*dir(1)*sig(2)
155 . +(dir(1)*dir(1)-dir(2)*dir(2))*sig(3)
156 s4 = dir(2)*sig(5)+dir(1)*sig(4)
157 s5 = dir(1)*sig(5)-dir(2)*sig(4)
158
159 sig(1)=s1
160 sig(2)=s2
161 sig(3)=s3
162 sig(4)=s4
163 sig(5)=s5
164 ENDIF
165 ENDIF
166
167 ENDIF
168C-----------------------------------
169C Rotation Element > Global
170C TENSOR ROTATION.
171C R passage Global -> Local.
172C R S Transpose(R)
173C------------------------------------
174
175 l11 =sig(1)
176 l22 =sig(2)
177 l33 =zero
178 l12 =sig(3)
179 l23 =sig(4)
180 l13 =sig(5)
181
182 r11 = e1x(i)
183 r12 = e1y(i)
184 r13 = e1z(i)
185
186 r21 = e2x(i)
187 r22 = e2y(i)
188 r23 = e2z(i)
189
190 r31 = e3x(i)
191 r32 = e3y(i)
192 r33 = e3z(i)
193
194 s11 =l11*r11+l12*r12+l13*r13
195 s12 =l11*r21+l12*r22+l13*r23
196 s13 =l11*r31+l12*r32+l13*r33
197 s21 =l12*r11+l22*r12+l23*r13
198 s22 =l12*r21+l22*r22+l23*r23
199 s23 =l12*r31+l22*r32+l23*r33
200 s31 =l13*r11+l23*r12
201 s32 =l13*r21+l23*r22
202 s33 =l13*r31+l23*r32
203
204 sig(1) =r11*s11+r12*s21+r13*s31
205 sig(2) =r21*s12+r22*s22+r23*s32
206 sig(3) =r31*s13+r32*s23+r33*s33
207 sig(4) =r11*s12+r12*s22+r13*s32
208 sig(5) =r21*s13+r22*s23+r23*s33
209 sig(6) =r11*s13+r12*s23+r13*s33
210C
211C-----------------------------------------------
212 RETURN
213 END
214!||====================================================================
215!|| shell_local_frame ../engine/source/output/dynain/shell_rota.F
216!||--- called by ------------------------------------------------------
217!|| dynain_c_strsg ../engine/source/output/dynain/dynain_c_strsg.F
218!||--- calls -----------------------------------------------------
219!|| clskew3 ../engine/source/elements/sh3n/coquedk/cdkcoor3.F
220!||--- uses -----------------------------------------------------
221!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
222!||====================================================================
223 SUBROUTINE shell_local_frame(JFT , JLT ,ITY ,IHBE ,IGTYP ,
224 2 IXC ,IXTG ,NFT ,X ,OFFG ,
225 3 RX ,RY ,RZ ,SX ,SY ,
226 4 SZ ,E1X ,E2X ,E3X ,E1Y ,
227 5 E2Y ,E3Y ,E1Z ,E2Z ,E3Z )
228C-----------------------------------------------
229C M o d u l e s
230C-----------------------------------------------
231 USE elbufdef_mod
232C-----------------------------------------------
233C I m p l i c i t T y p e s
234C-----------------------------------------------
235#include "implicit_f.inc"
236C-----------------------------------------------
237C G l o b a l P a r a m e t e r s
238C-----------------------------------------------
239#include "mvsiz_p.inc"
240#include "com01_c.inc"
241C-----------------------------------------------
242C D u m m y A r g u m e n t s
243C-----------------------------------------------
244 INTEGER NEL, ITY, NFT ,JFT, JLT , IHBE, IGTYP ,
245 . IXC(NIXC,*), IXTG(NIXTG,*)
246 my_real X(3,*), OFFG(*)
247 my_real
248 . E1X(MVSIZ), E1Y(MVSIZ), E1Z(MVSIZ),
249 . E2X(MVSIZ), E2Y(MVSIZ), E2Z(MVSIZ),
250 . E3X(MVSIZ), E3Y(MVSIZ), E3Z(MVSIZ),
251 . RX(MVSIZ), RY(MVSIZ), RZ(MVSIZ),
252 . SX(MVSIZ), SY(MVSIZ), SZ(MVSIZ)
253C-----------------------------------------------
254C L o c a l V a r i a b l e s
255C-----------------------------------------------
256 INTEGER I , K, I1,I2,I3,I4,IREL
257 my_real
258 . SUMA(MVSIZ)
259
260C--------------------------------------
261C LOCAL ELEMENT FRAME COMPUTING
262C--------------------------------------
263
264 IF (ITY == 3) then
265C------------------------
266C shells 4 nodes
267C-------------------------
268
269 IF (ihbe>10.OR.igtyp==16.OR.ishfram ==0) THEN
270 irel=0
271 ELSEIF (ishfram ==1) THEN
272 irel=2
273 ELSE
274 irel=1
275 END IF
276
277 DO i=jft,jlt
278 i1=ixc(2,i+nft)
279 i2=ixc(3,i+nft)
280 i3=ixc(4,i+nft)
281 i4=ixc(5,i+nft)
282
283 rx(i)=x(1,i2)+x(1,i3)-x(1,i1)-x(1,i4)
284 sx(i)=x(1,i3)+x(1,i4)-x(1,i1)-x(1,i2)
285 ry(i)=x(2,i2)+x(2,i3)-x(2,i1)-x(2,i4)
286 sy(i)=x(2,i3)+x(2,i4)-x(2,i1)-x(2,i2)
287 rz(i)=x(3,i2)+x(3,i3)-x(3,i1)-x(3,i4)
288 sz(i)=x(3,i3)+x(3,i4)-x(3,i1)-x(3,i2)
289 ENDDO
290C----------------------------
291C LOCAL SYSTEM
292C----------------------------
293 CALL clskew3(jft,jlt,irel,
294 . rx, ry, rz,
295 . sx, sy, sz,
296 . e1x,e1y,e1z,e2x,e2y,e2z,e3x,e3y,e3z,suma,offg )
297
298 ELSE
299
300C---------------------
301C shells 3 nodes
302C---------------------
303
304 IF (ihbe>=30) THEN
305 irel=0
306 ELSE
307 irel=2
308 END IF
309
310
311 DO i=jft,jlt
312 i1=ixtg(2,i+nft)
313 i2=ixtg(3,i+nft)
314 i3=ixtg(4,i+nft)
315
316 rx(i)=x(1,i2)-x(1,i1)
317 ry(i)=x(2,i2)-x(2,i1)
318 rz(i)=x(3,i2)-x(3,i1)
319 sx(i)=x(1,i3)-x(1,i1)
320 sy(i)=x(2,i3)-x(2,i1)
321 sz(i)=x(3,i3)-x(3,i1)
322 ENDDO
323C----------------------------
324C LOCAL SYSTEM
325C----------------------------
326 CALL clskew3(jft,jlt,irel,
327 . rx, ry, rz,
328 . sx, sy, sz,
329 . e1x,e1y,e1z,e2x,e2y,e2z,e3x,e3y,e3z,suma,offg )
330
331 ENDIF
332C
333C
334C-----------------------------------------------
335 RETURN
336 END
clskew3
subroutine clskew3(jft, jlt, irep, rx, ry, rz, sx, sy, sz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, det)
Definition clskew.F:34
shell_local_frame
subroutine shell_local_frame(jft, jlt, ity, ihbe, igtyp, ixc, ixtg, nft, x, offg, rx, ry, rz, sx, sy, sz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition shell_rota.F:228
shell_rota
subroutine shell_rota(i, ilay, nel, iorth, ity, igtyp, mlw, jdir, sig, elbuf_str, rx, ry, rz, sx, sy, sz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, dir_a, dir_b)
Definition shell_rota.F:36