OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cgshell.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine cgshell4 (elbuf_str, jft, jlt, pm, ixc, x, mas, xc, yc, zc)
subroutine cgshell3 (elbuf_str, jft, jlt, pm, ixtg, x, mas, xc, yc, zc)

Function/Subroutine Documentation

◆ cgshell3()

subroutine cgshell3 ( type (elbuf_struct_), target elbuf_str,
integer jft,
integer jlt,
pm,
integer, dimension(nixtg,*) ixtg,
x,
mas,
xc,
yc,
zc )

Definition at line 104 of file cgshell.F.

106C-----------------------------------------------
107 USE elbufdef_mod
108 use element_mod , only : nixtg
109C-----------------------------------------------
110C I m p l i c i t T y p e s
111C-----------------------------------------------
112#include "implicit_f.inc"
113C-----------------------------------------------
114C C o m m o n B l o c k s
115C-----------------------------------------------
116#include "mvsiz_p.inc"
117#include "param_c.inc"
118C-----------------------------------------------
119C D u m m y A r g u m e n t s
120C-----------------------------------------------
121 INTEGER IXTG(NIXTG,*), JFT, JLT
122C REAL
123 my_real
124 . pm(npropm,*),x(3,*),mas,xc,yc,zc
125 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
126C-----------------------------------------------
127C L o c a l V a r i a b l e s
128C-----------------------------------------------
129 INTEGER I,MAT
130C REAL
131 my_real
132 . xx,yy,zz,xmas25,xmas(mvsiz),area(mvsiz),
133 . rx, ry, rz,sx,sy,sz,e3x,e3y,e3z
134 TYPE(G_BUFEL_) ,POINTER :: GBUF
135C-----------------------------------------------
136 gbuf => elbuf_str%GBUF
137 DO i=jft,jlt
138 rx=x(1,ixtg(3,i))-x(1,ixtg(2,i))
139 sx=x(1,ixtg(4,i))-x(1,ixtg(2,i))
140 ry=x(2,ixtg(3,i))-x(2,ixtg(2,i))
141 sy=x(2,ixtg(4,i))-x(2,ixtg(2,i))
142 rz=x(3,ixtg(3,i))-x(3,ixtg(2,i))
143 sz=x(3,ixtg(4,i))-x(3,ixtg(2,i))
144 e3x = ry * sz - rz * sy
145 e3y = rz * sx - rx * sz
146 e3z = rx * sy - ry * sx
147 area(i) =half*sqrt(e3x*e3x + e3y*e3y + e3z*e3z)
148 END DO
149 DO i=jft,jlt
150 mat = ixtg(1,i)
151 xmas(i)=pm(1,mat)*area(i)*gbuf%THK(i)
152 END DO
153C
154C
155 DO i=jft,jlt
156 IF (gbuf%OFF(i) ==zero) cycle
157 xx= x(1,ixtg(2,i))+x(1,ixtg(3,i))+x(1,ixtg(4,i))
158 yy= x(2,ixtg(2,i))+x(2,ixtg(3,i))+x(2,ixtg(4,i))
159 zz= x(3,ixtg(2,i))+x(3,ixtg(3,i))+x(3,ixtg(4,i))
160 xmas25 = third*xmas(i)
161 mas = mas+xmas(i)
162 xc = xc + xmas25*xx
163 yc = yc + xmas25*yy
164 zc = zc + xmas25*zz
165 ENDDO
166C
167 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567

◆ cgshell4()

subroutine cgshell4 ( type (elbuf_struct_), target elbuf_str,
integer jft,
integer jlt,
pm,
integer, dimension(nixc,*) ixc,
x,
mas,
xc,
yc,
zc )

Definition at line 31 of file cgshell.F.

33C-----------------------------------------------
34 USE elbufdef_mod
35 use element_mod , only : nixc
36C-----------------------------------------------
37C I m p l i c i t T y p e s
38C-----------------------------------------------
39#include "implicit_f.inc"
40C-----------------------------------------------
41C C o m m o n B l o c k s
42C-----------------------------------------------
43#include "mvsiz_p.inc"
44#include "param_c.inc"
45C-----------------------------------------------
46C D u m m y A r g u m e n t s
47C-----------------------------------------------
48 INTEGER IXC(NIXC,*), JFT, JLT
49C REAL
50 my_real
51 . pm(npropm,*),x(3,*),mas,xc,yc,zc
52 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
53C-----------------------------------------------
54C L o c a l V a r i a b l e s
55C-----------------------------------------------
56 INTEGER I,MAT
57C REAL
58 my_real
59 . xx,yy,zz,xmas25,xmas(mvsiz),area(mvsiz),
60 . rx, ry, rz,sx,sy,sz,e3x,e3y,e3z
61 TYPE(G_BUFEL_) ,POINTER :: GBUF
62C-----------------------------------------------
63 gbuf => elbuf_str%GBUF
64 DO i=jft,jlt
65 rx=x(1,ixc(3,i))+x(1,ixc(4,i))-x(1,ixc(2,i))-x(1,ixc(5,i))
66 sx=x(1,ixc(4,i))+x(1,ixc(5,i))-x(1,ixc(2,i))-x(1,ixc(3,i))
67 ry=x(2,ixc(3,i))+x(2,ixc(4,i))-x(2,ixc(2,i))-x(2,ixc(5,i))
68 sy=x(2,ixc(4,i))+x(2,ixc(5,i))-x(2,ixc(2,i))-x(2,ixc(3,i))
69 rz=x(3,ixc(3,i))+x(3,ixc(4,i))-x(3,ixc(2,i))-x(3,ixc(5,i))
70 sz=x(3,ixc(4,i))+x(3,ixc(5,i))-x(3,ixc(2,i))-x(3,ixc(3,i))
71 e3x = ry * sz - rz * sy
72 e3y = rz * sx - rx * sz
73 e3z = rx * sy - ry * sx
74 area(i) =fourth*sqrt(e3x*e3x + e3y*e3y + e3z*e3z)
75 END DO
76 DO i=jft,jlt
77 mat = ixc(1,i)
78 xmas(i)=pm(1,mat)*area(i)*gbuf%THK(i)
79 END DO
80C
81C
82 DO i=jft,jlt
83 IF (gbuf%OFF(i) ==zero) cycle
84 xx= x(1,ixc(2,i))+x(1,ixc(3,i))+x(1,ixc(4,i))+x(1,ixc(5,i))
85 yy= x(2,ixc(2,i))+x(2,ixc(3,i))+x(2,ixc(4,i))+x(2,ixc(5,i))
86 zz= x(3,ixc(2,i))+x(3,ixc(3,i))+x(3,ixc(4,i))+x(3,ixc(5,i))
87 xmas25 = fourth*xmas(i)
88 mas = mas+xmas(i)
89 xc = xc + xmas25*xx
90 yc = yc + xmas25*yy
91 zc = zc + xmas25*zz
92 ENDDO
93C
94 RETURN