OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
c3evec3.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!|| c3evec3 ../engine/source/elements/sh3n/coque3n/c3evec3.F
25!||--- called by ------------------------------------------------------
26!|| c3forc3 ../engine/source/elements/sh3n/coque3n/c3forc3.F
27!|| c3forc3_crk ../engine/source/elements/xfem/c3forc3_crk.F
28!||--- calls -----------------------------------------------------
29!|| clskew3 ../engine/source/elements/sh3n/coquedk/cdkcoor3.F
30!||--- uses -----------------------------------------------------
31!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
32!||====================================================================
33 SUBROUTINE c3evec3(ELBUF_STR,DIR_A ,DIR_B ,JFT ,JLT ,
34 . IREP ,E1X0 ,E1Y0 ,E1Z0 ,E2X0 ,
35 . E2Y0 ,E2Z0 ,E3X0 ,E3Y0 ,E3Z0 ,
36 . E1X ,E1Y ,E1Z ,E2X ,
37 . E2Y ,E2Z ,E3X ,E3Y ,E3Z ,
38 . NLAY ,OFFG ,ECOS ,ESIN ,ISH3NFRAM,
39 . NEL ,AREA ,X21 ,Y21 ,Z21 ,
40 . X31 ,Y31 ,Z31 ,
41 . X1 ,X2 ,X3 ,Y1 ,Y2 ,
42 . Y3 ,Z1 ,Z2 ,Z3 )
43C-----------------------------------------------
44C M o d u l e s
45C-----------------------------------------------
46 USE elbufdef_mod
47C-----------------------------------------------
48C I m p l i c i t T y p e s
49C-----------------------------------------------
50#include "implicit_f.inc"
51C-----------------------------------------------
52C G l o b a l P a r a m e t e r s
53C-----------------------------------------------
54#include "mvsiz_p.inc"
55C-----------------------------------------------
56C D u m m y A r g u m e n t s
57C-----------------------------------------------
58 INTEGER JFT,JLT,IREP,NLAY,ISH3NFRAM,NEL
59 my_real
60 . DIR_A(NEL,*), DIR_B(NEL,*),E1X0(*), E1Y0(*), E1Z0(*),
61 . E2X0(*), E2Y0(*), E2Z0(*),E3X0(*), E3Y0(*), E3Z0(*),OFFG(*),
62 . E1X(MVSIZ), E1Y(MVSIZ), E1Z(MVSIZ),
63 . E2X(MVSIZ), E2Y(MVSIZ), E2Z(MVSIZ),
64 . E3X(MVSIZ), E3Y(MVSIZ), E3Z(MVSIZ),
65 . X21(MVSIZ), Y21(MVSIZ), Z21(MVSIZ),
66 . X31(MVSIZ), Y31(MVSIZ), Z31(MVSIZ),
67 . ecos(*) ,esin(*), area(mvsiz)
68! SP issue :
69 REAL(kind=8), dimension(mvsiz), INTENT(in) ::x1,x2,x3
70 REAL(kind=8), dimension(mvsiz), INTENT(in) ::y1,y2,y3
71 REAL(kind=8), dimension(mvsiz), INTENT(in) ::z1,z2,z3
72
73
74 TYPE (ELBUF_STRUCT_) :: ELBUF_STR
75C-----------------------------------------------
76C L o c a l V a r i a b l e s
77C-----------------------------------------------
78 INTEGER I,II,J,N,ILAW, IDRAPE,IGTYP,NPTT,IT,IPT_ALL,IPT
79C REAL
80 my_real x32(mvsiz), y32(mvsiz), z32(mvsiz),sum(mvsiz),
81 . e11(mvsiz),e12(mvsiz),e13(mvsiz),
82 . e21(mvsiz),e22(mvsiz),e23(mvsiz)
83 my_real v1,v2,v3,vr,vs,aa,bb,suma
84 my_real, DIMENSION(:) ,POINTER :: dir1, dir2
85C=======================================================================
86 idrape = elbuf_str%IDRAPE
87 igtyp = elbuf_str%IGTYP
88 DO i=jft,jlt
89 x21(i)=x2(i)-x1(i)
90 y21(i)=y2(i)-y1(i)
91 z21(i)=z2(i)-z1(i)
92 x31(i)=x3(i)-x1(i)
93 y31(i)=y3(i)-y1(i)
94 z31(i)=z3(i)-z1(i)
95 x32(i)=x3(i)-x2(i)
96 y32(i)=y3(i)-y2(i)
97 z32(i)=z3(i)-z2(i)
98 ENDDO
99C
100 IF (irep > 0) THEN
101 DO i=jft,jlt
102 e11(i) = x21(i)
103 e12(i) = y21(i)
104 e13(i) = z21(i)
105 e21(i) = x31(i)
106 e22(i) = y31(i)
107 e23(i) = z31(i)
108 ENDDO
109 ENDIF
110C
111 DO i=jft,jlt
112 e1x(i)= x21(i)
113 e1y(i)= y21(i)
114 e1z(i)= z21(i)
115 sum(i) = sqrt(e1x(i)*e1x(i)+e1y(i)*e1y(i)+e1z(i)*e1z(i))
116 e1x(i)=e1x(i)/sum(i)
117 e1y(i)=e1y(i)/sum(i)
118 e1z(i)=e1z(i)/sum(i)
119 ENDDO
120C
121 DO i=jft,jlt
122 e3x(i)=y31(i)*z32(i)-z31(i)*y32(i)
123 e3y(i)=z31(i)*x32(i)-x31(i)*z32(i)
124 e3z(i)=x31(i)*y32(i)-y31(i)*x32(i)
125 sum(i) = sqrt(e3x(i)*e3x(i)+e3y(i)*e3y(i)+e3z(i)*e3z(i))
126 e3x(i)=e3x(i)/sum(i)
127 e3y(i)=e3y(i)/sum(i)
128 e3z(i)=e3z(i)/sum(i)
129 area(i) = half * sum(i)
130 ENDDO
131C
132 DO i=jft,jlt
133 e2x(i)=e3y(i)*e1z(i)-e3z(i)*e1y(i)
134 e2y(i)=e3z(i)*e1x(i)-e3x(i)*e1z(i)
135 e2z(i)=e3x(i)*e1y(i)-e3y(i)*e1x(i)
136 sum(i) = sqrt(e2x(i)*e2x(i)+e2y(i)*e2y(i)+e2z(i)*e2z(i))
137 e2x(i)=e2x(i)/sum(i)
138 e2y(i)=e2y(i)/sum(i)
139 e2z(i)=e2z(i)/sum(i)
140 ENDDO
141C---
142 DO i=jft,jlt
143 e1x0(i) = e1x(i)
144 e1y0(i) = e1y(i)
145 e1z0(i) = e1z(i)
146 e2x0(i) = e2x(i)
147 e2y0(i) = e2y(i)
148 e2z0(i) = e2z(i)
149 e3x0(i) = e3x(i)
150 e3y0(i) = e3y(i)
151 e3z0(i) = e3z(i)
152 ENDDO
153 IF(ish3nfram ==0 ) THEN
154 ii = 0
155 CALL clskew3(jft,jlt,ii,
156 . x21, y21, z21,
157 . x31, y31, z31,
158 . e1x,e2x,e3x,e1y,e2y,e3y,e1z,e2z,e3z,sum,offg)
159 DO i=jft,jlt
160 ecos(i) = e1x(i)*e1x0(i)+e1y(i)*e1y0(i)+e1z(i)*e1z0(i)
161 aa = max(zero,one-ecos(i)*ecos(i))
162 esin(i) = sqrt(aa)
163 bb = e2x(i)*e1x0(i)+e2y(i)*e1y0(i)+e2z(i)*e1z0(i)
164 IF (bb >zero) esin(i) = -esin(i)
165 ENDDO
166 END IF
167C-----------------------------------------
168 IF(idrape > 0 .AND. (igtyp == 51 .OR. igtyp == 52)) THEN
169 ipt_all = 0
170 IF (irep == 1) THEN
171 DO n=1,nlay
172 nptt = elbuf_str%BUFLY(n)%NPTT
173! J = (N-1)*NEL
174 DO it = 1,nptt
175 dir1 => elbuf_str%BUFLY(n)%LBUF_DIR(it)%DIRA
176 ipt = ipt_all + it
177 j = 2*(ipt - 1)
178 DO i=jft,jlt
179 aa = dir1(i)
180 bb = dir1(i + nel)
181 v1 = aa*e11(i) + bb*e21(i)
182 v2 = aa*e12(i) + bb*e22(i)
183 v3 = aa*e13(i) + bb*e23(i)
184 vr = v1*e1x(i) + v2*e1y(i) + v3*e1z(i)
185 vs = v1*e2x(i) + v2*e2y(i) + v3*e2z(i)
186 suma=sqrt(vr*vr + vs*vs)
187 dir_a(i,j+1) = vr/suma
188 dir_a(i,j+2) = vs/suma
189 ENDDO
190 ENDDO
191 ipt_all = ipt_all + nptt
192 ENDDO
193 ELSEIF (irep == 2) THEN
194 DO n=1,nlay
195 nptt = elbuf_str%BUFLY(n)%NPTT
196! J = (N-1)*NEL
197 DO it = 1,nptt
198 dir1 => elbuf_str%BUFLY(n)%LBUF_DIR(it)%DIRA
199 dir2 => elbuf_str%BUFLY(n)%LBUF_DIR(it)%DIRB
200 ipt = ipt_all + it
201 j = 2*(ipt - 1)
202 DO i=jft,jlt
203C--- Axe I
204 aa = dir1(i)
205 bb = dir1(i + nel)
206 v1 = aa*e11(i) + bb*e21(i)
207 v2 = aa*e12(i) + bb*e22(i)
208 v3 = aa*e13(i) + bb*e23(i)
209 vr = v1*e1x(i) + v2*e1y(i) + v3*e1z(i)
210 vs = v1*e2x(i) + v2*e2y(i) + v3*e2z(i)
211 suma = one / max( sqrt(vr*vr + vs*vs), em20)
212 dir_a(i,j+1) = vr*suma
213 dir_a(i,j+2) = vs*suma
214C--- Axe II
215 aa = dir2(i)
216 bb = dir2(i + nel)
217 v1 = aa*e11(i) + bb*e21(i)
218 v2 = aa*e12(i) + bb*e22(i)
219 v3 = aa*e13(i) + bb*e23(i)
220 vr = v1*e1x(i) + v2*e1y(i) + v3*e1z(i)
221 vs = v1*e2x(i) + v2*e2y(i) + v3*e2z(i)
222 suma = one / max( sqrt(vr*vr + vs*vs), em20)
223 dir_b(i,j+1) = vr*suma
224 dir_b(i,j+2) = vs*suma
225 ENDDO
226 ENDDO
227 ipt_all = ipt_all + nptt
228 ENDDO
229 ELSEIF (irep == 3) THEN
230C mi xing law58 with other user laws with IREP = 0 within PID51
231 DO n=1,nlay
232 ilaw = elbuf_str%BUFLY(n)%ILAW
233 nptt = elbuf_str%BUFLY(n)%NPTT
234 IF (ilaw == 58) THEN
235! J = (N-1)*NEL
236 DO it = 1,nptt
237 dir1 => elbuf_str%BUFLY(n)%LBUF_DIR(it)%DIRA
238 dir2 => elbuf_str%BUFLY(n)%LBUF_DIR(it)%DIRB
239 ipt = ipt_all + it
240 j = 2*(ipt - 1)
241 DO i=jft,jlt
242C--- Axe I
243 aa = dir1(i)
244 bb = dir1(i + nel)
245 v1 = aa*e11(i) + bb*e21(i)
246 v2 = aa*e12(i) + bb*e22(i)
247 v3 = aa*e13(i) + bb*e23(i)
248 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
249 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
250 suma = one / max( sqrt(vr*vr + vs*vs), em20)
251 dir_a(i,j+1) = vr*suma
252 dir_a(i,j+2) = vs*suma
253C--- Axe II
254 aa = dir2(i)
255 bb = dir2(i + nel)
256 v1 = aa*e11(i) + bb*e21(i)
257 v2 = aa*e12(i) + bb*e22(i)
258 v3 = aa*e13(i) + bb*e23(i)
259 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
260 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
261 suma = one / max( sqrt(vr*vr + vs*vs), em20)
262 dir_b(i,j+1) = vr*suma
263 dir_b(i,j+2) = vs*suma
264 ENDDO
265 ENDDO
266 ipt_all = ipt_all + nptt
267 ELSE ! IREP = 0 within PID51
268 DO it = 1,nptt
269 dir1 => elbuf_str%BUFLY(n)%LBUF_DIR(it)%DIRA
270 ipt = ipt_all + it
271 j = 2*(ipt - 1)
272 DO i=jft,jlt
273 dir_a(i,j+1) = dir1(i)
274 dir_a(i,j+2) = dir1(i+nel)
275 ENDDO
276 ENDDO
277 ipt_all = ipt_all + nptt
278 ENDIF
279 ENDDO
280 ELSEIF (irep == 4) THEN
281C mi xing law58 with other user laws with IREP = 1 within PID51
282 DO n=1,nlay
283 ilaw = elbuf_str%BUFLY(n)%ILAW
284 nptt = elbuf_str%BUFLY(n)%NPTT
285! J = (N-1)*NEL
286 IF (ilaw == 58) THEN
287 DO it = 1,nptt
288 dir1 => elbuf_str%BUFLY(n)%LBUF_DIR(it)%DIRA
289 dir2 => elbuf_str%BUFLY(n)%LBUF_DIR(it)%DIRB
290 ipt = ipt_all + it
291 j = 2*(ipt - 1)
292 DO i=jft,jlt
293C--- Axe I
294 aa = dir1(i)
295 bb = dir1(i + nel)
296 v1 = aa*e11(i) + bb*e21(i)
297 v2 = aa*e12(i) + bb*e22(i)
298 v3 = aa*e13(i) + bb*e23(i)
299 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
300 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
301 suma = one / max( sqrt(vr*vr + vs*vs), em20)
302 dir_a(i,j+1) = vr*suma
303 dir_a(i,j+2) = vs*suma
304C--- Axe II
305 aa = dir2(i)
306 bb = dir2(i + nel)
307 v1 = aa*e11(i) + bb*e21(i)
308 v2 = aa*e12(i) + bb*e22(i)
309 v3 = aa*e13(i) + bb*e23(i)
310 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
311 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
312 suma = one / max( sqrt(vr*vr + vs*vs), em20)
313 dir_b(i,j+1) = vr*suma
314 dir_b(i,j+2) = vs*suma
315 ENDDO
316 ENDDO
317 ipt_all = ipt_all + nptt
318 ELSE ! IREP = 1 within PID51
319
320 DO it = 1,nptt
321 dir1 => elbuf_str%BUFLY(n)%LBUF_DIR(it)%DIRA
322 ipt = ipt_all + it
323 j = 2*(ipt - 1)
324 DO i=jft,jlt
325 aa = dir1(i)
326 bb = dir1(i + nel)
327 v1 = aa*e11(i) + bb*e21(i)
328 v2 = aa*e12(i) + bb*e22(i)
329 v3 = aa*e13(i) + bb*e23(i)
330 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
331 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
332 suma=sqrt(vr*vr + vs*vs)
333 dir_a(i,j+1) = vr/suma
334 dir_a(i,j+2) = vs/suma
335 ENDDO
336 ENDDO
337 ipt_all = ipt_all + nptt
338 ENDIF
339 ENDDO
340 ENDIF
341 ELSE ! idrape > 0
342 IF (irep == 1) THEN
343 DO n=1,nlay
344 dir1 => elbuf_str%BUFLY(n)%DIRA
345 j = 2*(n-1)
346 DO i=jft,jlt
347 aa = dir1(i)
348 bb = dir1(i + nel)
349 v1 = aa*e11(i) + bb*e21(i)
350 v2 = aa*e12(i) + bb*e22(i)
351 v3 = aa*e13(i) + bb*e23(i)
352 vr = v1*e1x(i) + v2*e1y(i) + v3*e1z(i)
353 vs = v1*e2x(i) + v2*e2y(i) + v3*e2z(i)
354 suma=sqrt(vr*vr + vs*vs)
355 dir_a(i,j+1) = vr/suma
356 dir_a(i,j+2) = vs/suma
357 ENDDO
358 ENDDO
359 ELSEIF (irep == 2) THEN
360 DO n=1,nlay
361 dir1 => elbuf_str%BUFLY(n)%DIRA
362 dir2 => elbuf_str%BUFLY(n)%DIRB
363 j = (n-1)*nel
364 j = 2*(n-1)
365 DO i=jft,jlt
366! Axe I
367 aa = dir1(i)
368 bb = dir1(i + nel)
369 v1 = aa*e11(i) + bb*e21(i)
370 v2 = aa*e12(i) + bb*e22(i)
371 v3 = aa*e13(i) + bb*e23(i)
372 vr = v1*e1x(i) + v2*e1y(i) + v3*e1z(i)
373 vs = v1*e2x(i) + v2*e2y(i) + v3*e2z(i)
374 suma = one / max( sqrt(vr*vr + vs*vs), em20)
375 dir_a(i,j+1) = vr*suma
376 dir_a(i,j+2) = vs*suma
377! Axe II
378 aa = dir2(i)
379 bb = dir2(i + nel)
380 v1 = aa*e11(i) + bb*e21(i)
381 v2 = aa*e12(i) + bb*e22(i)
382 v3 = aa*e13(i) + bb*e23(i)
383 vr = v1*e1x(i) + v2*e1y(i) + v3*e1z(i)
384 vs = v1*e2x(i) + v2*e2y(i) + v3*e2z(i)
385 suma = one / max( sqrt(vr*vr + vs*vs), em20)
386 dir_b(i,j+1) = vr*suma
387 dir_b(i,j+2) = vs*suma
388 ENDDO
389 ENDDO
390 ELSEIF (irep == 3) THEN
391C mi xing law58 with other user laws with IREP = 0 within PID51
392 DO n=1,nlay
393 ilaw = elbuf_str%BUFLY(n)%ILAW
394 j = (n-1)*nel
395 j = 2*(n-1)
396 IF (ilaw == 58) THEN
397 dir1 => elbuf_str%BUFLY(n)%DIRA
398 dir2 => elbuf_str%BUFLY(n)%DIRB
399 DO i=jft,jlt
400! Axe I
401 aa = dir1(i)
402 bb = dir1(i + nel)
403 v1 = aa*e11(i) + bb*e21(i)
404 v2 = aa*e12(i) + bb*e22(i)
405 v3 = aa*e13(i) + bb*e23(i)
406 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
407 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
408 suma = one / max( sqrt(vr*vr + vs*vs), em20)
409 dir_a(i,j+1) = vr*suma
410 dir_a(i,j+2) = vs*suma
411! Axe II
412 aa = dir2(i)
413 bb = dir2(i + nel)
414 v1 = aa*e11(i) + bb*e21(i)
415 v2 = aa*e12(i) + bb*e22(i)
416 v3 = aa*e13(i) + bb*e23(i)
417 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
418 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
419 suma = one / max( sqrt(vr*vr + vs*vs), em20)
420 dir_b(i,j+1) = vr*suma
421 dir_b(i,j+2) = vs*suma
422 ENDDO
423 ELSE ! IREP = 0 within PID51
424 dir1 => elbuf_str%BUFLY(n)%DIRA
425 DO i=jft,jlt
426 dir_a(i,j+1) = dir1(i)
427 dir_a(i,j+2) = dir1(i+nel)
428 ENDDO
429 ENDIF
430 ENDDO
431 ELSEIF (irep == 4) THEN
432! mi xing law58 with other user laws with IREP = 1 within PID51
433 DO n=1,nlay
434 ilaw = elbuf_str%BUFLY(n)%ILAW
435 j = (n-1)*nel
436 j = 2*(n-1)
437 IF (ilaw == 58) THEN
438 dir1 => elbuf_str%BUFLY(n)%DIRA
439 dir2 => elbuf_str%BUFLY(n)%DIRB
440 DO i=jft,jlt
441! Axe I
442 aa = dir1(i)
443 bb = dir1(i + nel)
444 v1 = aa*e11(i) + bb*e21(i)
445 v2 = aa*e12(i) + bb*e22(i)
446 v3 = aa*e13(i) + bb*e23(i)
447 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
448 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
449 suma = one / max( sqrt(vr*vr + vs*vs), em20)
450 dir_a(i,j+1) = vr*suma
451 dir_a(i,j+2) = vs*suma
452! Axe II
453 aa = dir2(i)
454 bb = dir2(i + nel)
455 v1 = aa*e11(i) + bb*e21(i)
456 v2 = aa*e12(i) + bb*e22(i)
457 v3 = aa*e13(i) + bb*e23(i)
458 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
459 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
460 suma = one / max( sqrt(vr*vr + vs*vs), em20)
461 dir_b(i,j+1) = vr*suma
462 dir_b(i,j+2) = vs*suma
463 ENDDO
464 ELSE ! IREP = 1 within PID51
465 dir1 => elbuf_str%BUFLY(n)%DIRA
466 DO i=jft,jlt
467 aa = dir1(i)
468 bb = dir1(i + nel)
469 v1 = aa*e11(i) + bb*e21(i)
470 v2 = aa*e12(i) + bb*e22(i)
471 v3 = aa*e13(i) + bb*e23(i)
472 vr = v1*e1x(i)+ v2*e1y(i) + v3*e1z(i)
473 vs = v1*e2x(i)+ v2*e2y(i) + v3*e2z(i)
474 suma=sqrt(vr*vr + vs*vs)
475 dir_a(i,j+1) = vr/suma
476 dir_a(i,j+2) = vs/suma
477 ENDDO
478 ENDIF
479 ENDDO
480 ENDIF
481 ENDIF
482C------ restore old local sys compute cos, sin relative to old
483 IF(ish3nfram ==0 ) THEN
484 DO i=jft,jlt
485 e11(i) = e1x(i)
486 e12(i) = e1y(i)
487 e13(i) = e1z(i)
488 e21(i) = e2x(i)
489 e22(i) = e2y(i)
490 e23(i) = e2z(i)
491 ENDDO
492 DO i=jft,jlt
493 e1x(i) = e1x0(i)
494 e1y(i) = e1y0(i)
495 e1z(i) = e1z0(i)
496 e2x(i) = e2x0(i)
497 e2y(i) = e2y0(i)
498 e2z(i) = e2z0(i)
499 e3x(i) = e3x0(i)
500 e3y(i) = e3y0(i)
501 e3z(i) = e3z0(i)
502 ENDDO
503C-------------Ej0 is used in mulawc for user'laws (99)
504 DO i=jft,jlt
505 e1x0(i) = e11(i)
506 e1y0(i) = e12(i)
507 e1z0(i) = e13(i)
508 e2x0(i) = e21(i)
509 e2y0(i) = e22(i)
510 e2z0(i) = e23(i)
511 ENDDO
512 END IF
513c-----------
514 RETURN
515 END
516!||====================================================================
517!|| c3newve3 ../engine/source/elements/sh3n/coque3n/c3evec3.F
518!||====================================================================
519 SUBROUTINE c3newve3(JFT ,JLT ,TG3 ,ECOS ,ESIN )
520C-----------------------------------------------
521C I m p l i c i t T y p e s
522C-----------------------------------------------
523#include "implicit_f.inc"
524C-----------------------------------------------
525C G l o b a l P a r a m e t e r s
526C-----------------------------------------------
527#include "mvsiz_p.inc"
528C-----------------------------------------------
529C D u m m y A r g u m e n t s
530C-----------------------------------------------
531 INTEGER JFT, JLT
532 my_real
533 . ECOS(MVSIZ),ESIN(MVSIZ),TG3(MVSIZ)
534C-----------------------------------------------
535C L o c a l V a r i a b l e s
536C-----------------------------------------------
537 INTEGER I, II
538C REAL
539 my_real
540 . theta
541C=======================================================================
542 DO i=jft,jlt
543 theta = half*(atan(tg3(i)))
544 ecos(i)= cos(theta)
545 esin(i)= sin(theta)
546 ENDDO
547c-----------
548 RETURN
549 END
550!||====================================================================
551!|| shroto3 ../engine/source/elements/sh3n/coque3n/c3evec3.F
552!||--- called by ------------------------------------------------------
553!|| c3forc3 ../engine/source/elements/sh3n/coque3n/c3forc3.F
554!||====================================================================
555 SUBROUTINE shroto3(JFT,JLT,ECOS,ESIN,EXX,
556 . EYY,EXY,EXZ,EYZ,KXX,
557 . KYY,KXY)
558C-----------------------------------------------
559C I m p l i c i t T y p e s
560C-----------------------------------------------
561#include "implicit_f.inc"
562C-----------------------------------------------
563C G l o b a l P a r a m e t e r s
564C-----------------------------------------------
565#include "mvsiz_p.inc"
566C-----------------------------------------------
567C D u m m y A r g u m e n t s
568C-----------------------------------------------
569 INTEGER JFT, JLT
570C REAL
571 my_real
572 . EXX(MVSIZ), EYY(MVSIZ), EXY(MVSIZ), EXZ(MVSIZ), EYZ(MVSIZ),
573 . KXX(MVSIZ), KYY(MVSIZ), KXY(MVSIZ),ECOS(MVSIZ),ESIN(MVSIZ)
574C-----------------------------------------------
575C L o c a l V a r i a b l e s
576C-----------------------------------------------
577 INTEGER I,J
578C REAL
579 my_real
580 . m2(mvsiz),n2(mvsiz),mn(mvsiz),cm(4),mn2(mvsiz),t1,txy,tyx
581C-----------------------------------------------
582 DO i=jft,jlt
583 m2(i)= ecos(i)*ecos(i)
584 n2(i)= esin(i)*esin(i)
585 mn(i)= ecos(i)*esin(i)
586 mn2(i)= two*mn(i)
587 ENDDO
588C------Exy :2*exy
589 DO i=jft,jlt
590 t1 = exy(i)*mn(i)
591 cm(1)=m2(i)*exx(i)+n2(i)*eyy(i)-t1
592 cm(2)=n2(i)*exx(i)+m2(i)*eyy(i)+t1
593 cm(3)=(exx(i)-eyy(i))*mn2(i)+exy(i)*(m2(i)-n2(i))
594 exx(i)=cm(1)
595 eyy(i)=cm(2)
596 exy(i)=cm(3)
597 ENDDO
598 DO i=jft,jlt
599 t1 = kxy(i)*mn(i)
600 cm(1)=m2(i)*kxx(i)+n2(i)*kyy(i)-t1
601 cm(2)=n2(i)*kxx(i)+m2(i)*kyy(i)+t1
602 cm(3)=(kxx(i)-kyy(i))*mn2(i)+kxy(i)*(m2(i)-n2(i))
603 kxx(i)=cm(1)
604 kyy(i)=cm(2)
605 kxy(i)=cm(3)
606 ENDDO
607C
608 DO i=jft,jlt
609 cm(1)= exz(i)*ecos(i) - eyz(i)*esin(i)
610 cm(2)= exz(i)*esin(i) + eyz(i)*ecos(i)
611 exz(i) = cm(1)
612 eyz(i) = cm(2)
613 ENDDO
614C
615 RETURN
616 END
617!||====================================================================
618!|| shtroto3 ../engine/source/elements/sh3n/coque3n/c3evec3.F
619!||--- called by ------------------------------------------------------
620!|| c3forc3 ../engine/source/elements/sh3n/coque3n/c3forc3.F
621!||====================================================================
622 SUBROUTINE shtroto3(JFT,JLT,ECOS,ESIN,GSTR,
623 . F_DEF,ISMSTR,NEL)
624C-----------------------------------------------
625C I m p l i c i t T y p e s
626C-----------------------------------------------
627#include "implicit_f.inc"
628C-----------------------------------------------
629C G l o b a l P a r a m e t e r s
630C-----------------------------------------------
631#include "mvsiz_p.inc"
632C-----------------------------------------------
633C D u m m y A r g u m e n t s
634C-----------------------------------------------
635 INTEGER JFT, JLT,ISMSTR,NEL
636C REAL
637 my_real
638 . GSTR(NEL,8), F_DEF(MVSIZ,8),ECOS(MVSIZ),ESIN(MVSIZ)
639C-----------------------------------------------
640C L o c a l V a r i a b l e s
641C-----------------------------------------------
642 INTEGER I,J
643C REAL
644 my_real
645 . M2(MVSIZ),N2(MVSIZ),MN(MVSIZ),CM(4),MN2(MVSIZ),T1,T2,T3
646C-----------------------------------------------
647 DO i=jft,jlt
648 m2(i)= ecos(i)*ecos(i)
649 n2(i)= esin(i)*esin(i)
650 mn(i)= ecos(i)*esin(i)
651 mn2(i)= two*mn(i)
652 ENDDO
653C------Exy :2*exy
654 IF (ismstr==10) THEN
655 DO i=jft,jlt
656 t1 = (f_def(i,3)+f_def(i,4))*mn(i)
657 cm(1)=m2(i)*f_def(i,1)+n2(i)*f_def(i,2)-t1
658 cm(2)=n2(i)*f_def(i,1)+m2(i)*f_def(i,2)+t1
659 t2 = (f_def(i,1)-f_def(i,2))*mn2(i)
660 cm(3)=t2+f_def(i,3)*m2(i)-f_def(i,4)*n2(i)
661 cm(4)=t2+f_def(i,4)*m2(i)-f_def(i,3)*n2(i)
662 f_def(i,1)=cm(1)
663 f_def(i,2)=cm(2)
664 f_def(i,3)=cm(3)
665 f_def(i,4)=cm(4)
666 ENDDO
667 ELSE
668 DO i=jft,jlt
669 t1 = gstr(i,3)*mn(i)
670 cm(1)=m2(i)*gstr(i,1)+n2(i)*gstr(i,2)-t1
671 cm(2)=n2(i)*gstr(i,1)+m2(i)*gstr(i,2)+t1
672 cm(3)=(gstr(i,1)-gstr(i,2))*mn2(i)+gstr(i,3)*(m2(i)-n2(i))
673 gstr(i,1)=cm(1)
674 gstr(i,2)=cm(2)
675 gstr(i,3)=cm(3)
676 ENDDO
677 END IF
678C
679 RETURN
680 END
681!||====================================================================
682!|| c3sroto3 ../engine/source/elements/sh3n/coque3n/c3evec3.f
683!||--- called by ------------------------------------------------------
684!|| c3forc3 ../engine/source/elements/sh3n/coque3n/c3forc3.F
685!||====================================================================
686 SUBROUTINE c3sroto3(JFT ,JLT ,ECOS ,ESIN ,FOR,
687 . MOM ,NFOR ,NMOM ,IFRAM_OLD ,NEL)
688C-----------------------------------------------
689C I m p l i c i t T y p e s
690C-----------------------------------------------
691#include "implicit_f.inc"
692C-----------------------------------------------
693C G l o b a l P a r a m e t e r s
694C-----------------------------------------------
695#include "mvsiz_p.inc"
696C-----------------------------------------------
697C D u m m y A r g u m e n t s
698C-----------------------------------------------
699 INTEGER JFT, JLT,IFRAM_OLD,NEL
700C REAL
701 my_real
702 . FOR(NEL,5), MOM(NEL,3),NFOR(NEL,5), NMOM(NEL,3),
703 . ECOS(MVSIZ),ESIN(MVSIZ)
704C-----------------------------------------------
705C L o c a l V a r i a b l e s
706C-----------------------------------------------
707 INTEGER I,J
708C REAL
709 my_real
710 . M2(MVSIZ),N2(MVSIZ),MN(MVSIZ),CM(5),MN2(MVSIZ),T1,T2,T3
711C-----------------------------------------------
712 IF(ifram_old ==0 ) THEN
713 DO i=jft,jlt
714 m2(i)= ecos(i)*ecos(i)
715 n2(i)= esin(i)*esin(i)
716 mn(i)= ecos(i)*esin(i)
717 mn2(i)= two*mn(i)
718 ENDDO
719C------
720 DO i=jft,jlt
721 t1 = for(i,3)*mn2(i)
722 nfor(i,1)=m2(i)*for(i,1)+n2(i)*for(i,2)+t1
723 nfor(i,2)=n2(i)*for(i,1)+m2(i)*for(i,2)-t1
724 nfor(i,3)=(-for(i,1)+for(i,2))*mn(i)+for(i,3)*(m2(i)-n2(i))
725 nfor(i,4)=ecos(i)*for(i,4)-esin(i)*for(i,5)
726 nfor(i,5)=ecos(i)*for(i,5)+esin(i)*for(i,4)
727 ENDDO
728 DO i=jft,jlt
729 t1 = mom(i,3)*mn2(i)
730 nmom(i,1)=m2(i)*mom(i,1)+n2(i)*mom(i,2)+t1
731 nmom(i,2)=n2(i)*mom(i,1)+m2(i)*mom(i,2)-t1
732 nmom(i,3)=(-mom(i,1)+mom(i,2))*mn(i)+mom(i,3)*(m2(i)-n2(i))
733 ENDDO
734 ELSE
735 DO i=jft,jlt
736 nfor(i,1) = for(i,1)
737 nfor(i,2) = for(i,2)
738 nfor(i,3) = for(i,3)
739 nfor(i,4) = for(i,4)
740 nfor(i,5) = for(i,5)
741 nmom(i,1) = mom(i,1)
742 nmom(i,2) = mom(i,2)
743 nmom(i,3) = mom(i,3)
744 ENDDO
745 END IF
746C
747 RETURN
748 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
shtroto3
subroutine shtroto3(jft, jlt, ecos, esin, gstr, f_def, ismstr, nel)
Definition c3evec3.F:624
c3newve3
subroutine c3newve3(jft, jlt, tg3, ecos, esin)
Definition c3evec3.F:520
c3evec3
subroutine c3evec3(elbuf_str, dir_a, dir_b, jft, jlt, irep, e1x0, e1y0, e1z0, e2x0, e2y0, e2z0, e3x0, e3y0, e3z0, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, nlay, offg, ecos, esin, ish3nfram, nel, area, x21, y21, z21, x31, y31, z31, x1, x2, x3, y1, y2, y3, z1, z2, z3)
Definition c3evec3.F:43
c3sroto3
subroutine c3sroto3(jft, jlt, ecos, esin, for, mom, nfor, nmom, ifram_old, nel)
Definition c3evec3.F:688
shroto3
subroutine shroto3(jft, jlt, ecos, esin, exx, eyy, exy, exz, eyz, kxx, kyy, kxy)
Definition c3evec3.F:558
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
max
#define max(a, b)
Definition macros.h:21