OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sta_c_get_t3lsys.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!|| get_t3lsys ../engine/source/output/sta/sta_c_get_t3lsys.F
25!||--- called by ------------------------------------------------------
26!|| stat_c_strsfg ../engine/source/output/sta/stat_c_strsfg.f
27!||--- calls -----------------------------------------------------
28!|| clskew3 ../engine/source/elements/sh3n/coquedk/cdkcoor3.F
29!|| cortdir3 ../engine/source/elements/shell/coque/cortdir3.F
30!||--- uses -----------------------------------------------------
31!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
32!|| element_mod ../common_source/modules/elements/element_mod.F90
33!||====================================================================
34 SUBROUTINE get_t3lsys(JFT ,JLT ,IXTG ,X ,OFFG ,
35 . IREL ,VQ ,NLAY ,IREP ,NEL ,
36 . DIR_A ,DIR_B,ELBUF_STR)
37C-----------------------------------------------
38C M o d u l e s
39C-----------------------------------------------
40 USE elbufdef_mod
41 use element_mod , only : nixtg
42C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
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 IXTG(NIXTG,*),JFT,JLT,IREL,NLAY,IREP,NEL
52 my_real
53 . x(3,*), offg(*),vq(3,3,mvsiz),dir_a(*),dir_b(*)
54 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
55C-----------------------------------------------
56C L O C A L V A R I A B L E S
57C-----------------------------------------------
58 INTEGER I
59 INTEGER I2,I3,I1
60 my_real
61 . rx(mvsiz),ry(mvsiz),rz(mvsiz),sx(mvsiz),sy(mvsiz),
62 . r11(mvsiz),r12(mvsiz),r13(mvsiz),r21(mvsiz),r22(mvsiz),
63 . r23(mvsiz),r31(mvsiz),r32(mvsiz),r33(mvsiz),
64 . sz(mvsiz),deta1(mvsiz)
65C----------------------------------------------
66 DO i=jft,jlt
67 i1=ixtg(2,i)
68 i2=ixtg(3,i)
69 i3=ixtg(4,i)
70
71 rx(i)=x(1,i2)-x(1,i1)
72 ry(i)=x(2,i2)-x(2,i1)
73 rz(i)=x(3,i2)-x(3,i1)
74 sx(i)=x(1,i3)-x(1,i1)
75 sy(i)=x(2,i3)-x(2,i1)
76 sz(i)=x(3,i3)-x(3,i1)
77 ENDDO
78C----------------------------
79C LOCAL SYSTEM
80C----------------------------
81 CALL clskew3(jft,jlt,irel,
82 . rx, ry, rz,
83 . sx, sy, sz,
84 . r11,r12,r13,r21,r22,r23,r31,r32,r33,deta1,offg )
85 DO i=jft,jlt
86 vq(1,1,i)=r11(i)
87 vq(2,1,i)=r21(i)
88 vq(3,1,i)=r31(i)
89 vq(1,2,i)=r12(i)
90 vq(2,2,i)=r22(i)
91 vq(3,2,i)=r32(i)
92 vq(1,3,i)=r13(i)
93 vq(2,3,i)=r23(i)
94 vq(3,3,i)=r33(i)
95 ENDDO
96 CALL cortdir3(elbuf_str,dir_a ,dir_b ,jft ,jlt ,
97 . nlay ,irep ,rx ,ry ,rz ,
98 . sx ,sy ,sz ,r11 ,r21 ,
99 . r31 ,r12 ,r22 ,r32 ,nel )
100C
101 RETURN
102 END
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
subroutine cortdir3(elbuf_str, dir_a, dir_b, jft, jlt, nlay, irep, rx, ry, rz, sx, sy, sz, e1x, e1y, e1z, e2x, e2y, e2z, nel)
Definition cortdir3.F:45
subroutine get_t3lsys(jft, jlt, ixtg, x, offg, irel, vq, nlay, irep, nel, dir_a, dir_b, elbuf_str)
subroutine stat_c_strsfg(elbuf_tab, x, iparg, ipm, igeo, ixc, ixtg, wa, wap0, ipartc, iparttg, ipart_state, stat_indxc, stat_indxtg, thke, sizp0, geo, stack, drape_sh4n, drape_sh3n, drapeg)