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

Go to the source code of this file.

Functions/Subroutines

subroutine s4orth3 (jhbe, gama, nel, r11, r12, r13, r21, r22, r23, r31, r32, r33, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4)
subroutine s4ortht3 (jhbe, gama, nel, r11, r12, r13, r21, r22, r23, r31, r32, r33, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4)

Function/Subroutine Documentation

◆ s4orth3()

subroutine s4orth3 ( integer jhbe,
gama,
integer nel,
r11,
r12,
r13,
r21,
r22,
r23,
r31,
r32,
r33,
double precision, dimension(*) x1,
double precision, dimension(*) x2,
double precision, dimension(*) x3,
double precision, dimension(*) x4,
double precision, dimension(*) y1,
double precision, dimension(*) y2,
double precision, dimension(*) y3,
double precision, dimension(*) y4,
double precision, dimension(*) z1,
double precision, dimension(*) z2,
double precision, dimension(*) z3,
double precision, dimension(*) z4 )

Definition at line 28 of file s4orth3.F.

33C-----------------------------------------------
34C I m p l i c i t T y p e s
35C-----------------------------------------------
36#include "implicit_f.inc"
37C-----------------------------------------------
38C G l o b a l P a r a m e t e r s
39C-----------------------------------------------
40#include "mvsiz_p.inc"
41C-----------------------------------------------
42C D u m m y A r g u m e n t s
43C-----------------------------------------------
44 INTEGER JHBE,NEL
45C REAL
46 double precision
47 . x1(*), x2(*), x3(*), x4(*),
48 . y1(*), y2(*), y3(*), y4(*),
49 . z1(*), z2(*), z3(*), z4(*)
50 my_real
51 . gama(nel,6),
52 . r11(*), r12(*), r13(*),
53 . r21(*), r22(*), r23(*),
54 . r31(*), r32(*), r33(*)
55C-----------------------------------------------
56C C o m m o n B l o c k s
57C-----------------------------------------------
58#include "vect01_c.inc"
59C-----------------------------------------------
60C L o c a l V a r i a b l e s
61C-----------------------------------------------
62 INTEGER I
63C REAL
64 my_real
65 . g11(mvsiz),g12(mvsiz),g13(mvsiz),
66 . g21(mvsiz),g22(mvsiz),g23(mvsiz),
67 . g31(mvsiz),g32(mvsiz),g33(mvsiz)
68 my_real
69 . xl(mvsiz),yl(mvsiz),zl(mvsiz)
70C-----------------------------------------------
71C PASSAGE AU REPERE ORTHOTROPE.
72 IF (jhbe/=24 .AND. jhbe/=222 .AND. jhbe/=14)THEN
73 DO i=lft,llt
74C Extration de G tq Xortho=Transpose(G) Xcvt
75C =Transpose(G) Transpose(R) Xglobal.
76 g11(i)=gama(i,1)
77 g21(i)=gama(i,2)
78 g31(i)=gama(i,3)
79 g12(i)=gama(i,4)
80 g22(i)=gama(i,5)
81 g32(i)=gama(i,6)
82 g13(i)=g21(i)*g32(i)-g31(i)*g22(i)
83 g23(i)=g31(i)*g12(i)-g11(i)*g32(i)
84 g33(i)=g11(i)*g22(i)-g21(i)*g12(i)
85C
86 xl(i)=g11(i)*x1(i)+g21(i)*y1(i)+g31(i)*z1(i)
87 yl(i)=g12(i)*x1(i)+g22(i)*y1(i)+g32(i)*z1(i)
88 zl(i)=g13(i)*x1(i)+g23(i)*y1(i)+g33(i)*z1(i)
89 x1(i)=xl(i)
90 y1(i)=yl(i)
91 z1(i)=zl(i)
92 xl(i)=g11(i)*x2(i)+g21(i)*y2(i)+g31(i)*z2(i)
93 yl(i)=g12(i)*x2(i)+g22(i)*y2(i)+g32(i)*z2(i)
94 zl(i)=g13(i)*x2(i)+g23(i)*y2(i)+g33(i)*z2(i)
95 x2(i)=xl(i)
96 y2(i)=yl(i)
97 z2(i)=zl(i)
98 xl(i)=g11(i)*x3(i)+g21(i)*y3(i)+g31(i)*z3(i)
99 yl(i)=g12(i)*x3(i)+g22(i)*y3(i)+g32(i)*z3(i)
100 zl(i)=g13(i)*x3(i)+g23(i)*y3(i)+g33(i)*z3(i)
101 x3(i)=xl(i)
102 y3(i)=yl(i)
103 z3(i)=zl(i)
104 xl(i)=g11(i)*x4(i)+g21(i)*y4(i)+g31(i)*z4(i)
105 yl(i)=g12(i)*x4(i)+g22(i)*y4(i)+g32(i)*z4(i)
106 zl(i)=g13(i)*x4(i)+g23(i)*y4(i)+g33(i)*z4(i)
107 x4(i)=xl(i)
108 y4(i)=yl(i)
109 z4(i)=zl(i)
110C
111C MATRICE DE PASSAGE GLOBAL -> ORTHOTROPE.
112 r11(i)=r11(i)*g11(i)+r12(i)*g21(i)+r13(i)*g31(i)
113 r12(i)=r11(i)*g12(i)+r12(i)*g22(i)+r13(i)*g32(i)
114 r13(i)=r11(i)*g13(i)+r12(i)*g23(i)+r13(i)*g33(i)
115 r21(i)=r21(i)*g11(i)+r22(i)*g21(i)+r23(i)*g31(i)
116 r22(i)=r21(i)*g12(i)+r22(i)*g22(i)+r23(i)*g32(i)
117 r23(i)=r21(i)*g13(i)+r22(i)*g23(i)+r23(i)*g33(i)
118 r31(i)=r31(i)*g11(i)+r32(i)*g21(i)+r33(i)*g31(i)
119 r32(i)=r31(i)*g12(i)+r32(i)*g22(i)+r33(i)*g32(i)
120 r33(i)=r31(i)*g13(i)+r32(i)*g23(i)+r33(i)*g33(i)
121 END DO
122 END IF
123c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
124 RETURN
my_real
#define my_real
Definition cppsort.cpp:32

◆ s4ortht3()

subroutine s4ortht3 ( integer jhbe,
gama,
integer nel,
r11,
r12,
r13,
r21,
r22,
r23,
r31,
r32,
r33,
x1,
x2,
x3,
x4,
y1,
y2,
y3,
y4,
z1,
z2,
z3,
z4 )

Definition at line 131 of file s4orth3.F.

136C-----------------------------------------------
137C I m p l i c i t T y p e s
138C-----------------------------------------------
139#include "implicit_f.inc"
140C-----------------------------------------------
141C G l o b a l P a r a m e t e r s
142C-----------------------------------------------
143#include "mvsiz_p.inc"
144C-----------------------------------------------
145C D u m m y A r g u m e n t s
146C-----------------------------------------------
147 INTEGER JHBE,NEL
148C REAL
149 my_real
150 . x1(*), x2(*), x3(*), x4(*),
151 . y1(*), y2(*), y3(*), y4(*),
152 . z1(*), z2(*), z3(*), z4(*),
153 . gama(nel,6),
154 . r11(*), r12(*), r13(*),
155 . r21(*), r22(*), r23(*),
156 . r31(*), r32(*), r33(*)
157C-----------------------------------------------
158C C o m m o n B l o c k s
159C-----------------------------------------------
160#include "vect01_c.inc"
161C-----------------------------------------------
162C L o c a l V a r i a b l e s
163C-----------------------------------------------
164 INTEGER I
165C REAL
166 my_real
167 . g11(mvsiz),g12(mvsiz),g13(mvsiz),
168 . g21(mvsiz),g22(mvsiz),g23(mvsiz),
169 . g31(mvsiz),g32(mvsiz),g33(mvsiz)
170 my_real
171 . xl(mvsiz),yl(mvsiz),zl(mvsiz)
172C-----------------------------------------------
173C PASSAGE AU REPERE ORTHOTROPE.
174 IF (jhbe/=24 .AND. jhbe/=222 .AND. jhbe/=14)THEN
175 DO i=lft,llt
176C Extration de G tq Xortho=Transpose(G) Xcvt
177C =Transpose(G) Transpose(R) Xglobal.
178 g11(i)=gama(i,1)
179 g21(i)=gama(i,2)
180 g31(i)=gama(i,3)
181 g12(i)=gama(i,4)
182 g22(i)=gama(i,5)
183 g32(i)=gama(i,6)
184 g13(i)=g21(i)*g32(i)-g31(i)*g22(i)
185 g23(i)=g31(i)*g12(i)-g11(i)*g32(i)
186 g33(i)=g11(i)*g22(i)-g21(i)*g12(i)
187C
188 xl(i)=g11(i)*x1(i)+g21(i)*y1(i)+g31(i)*z1(i)
189 yl(i)=g12(i)*x1(i)+g22(i)*y1(i)+g32(i)*z1(i)
190 zl(i)=g13(i)*x1(i)+g23(i)*y1(i)+g33(i)*z1(i)
191 x1(i)=xl(i)
192 y1(i)=yl(i)
193 z1(i)=zl(i)
194 xl(i)=g11(i)*x2(i)+g21(i)*y2(i)+g31(i)*z2(i)
195 yl(i)=g12(i)*x2(i)+g22(i)*y2(i)+g32(i)*z2(i)
196 zl(i)=g13(i)*x2(i)+g23(i)*y2(i)+g33(i)*z2(i)
197 x2(i)=xl(i)
198 y2(i)=yl(i)
199 z2(i)=zl(i)
200 xl(i)=g11(i)*x3(i)+g21(i)*y3(i)+g31(i)*z3(i)
201 yl(i)=g12(i)*x3(i)+g22(i)*y3(i)+g32(i)*z3(i)
202 zl(i)=g13(i)*x3(i)+g23(i)*y3(i)+g33(i)*z3(i)
203 x3(i)=xl(i)
204 y3(i)=yl(i)
205 z3(i)=zl(i)
206 xl(i)=g11(i)*x4(i)+g21(i)*y4(i)+g31(i)*z4(i)
207 yl(i)=g12(i)*x4(i)+g22(i)*y4(i)+g32(i)*z4(i)
208 zl(i)=g13(i)*x4(i)+g23(i)*y4(i)+g33(i)*z4(i)
209 x4(i)=xl(i)
210 y4(i)=yl(i)
211 z4(i)=zl(i)
212C
213C MATRICE DE PASSAGE GLOBAL -> ORTHOTROPE.
214 r11(i)=r11(i)*g11(i)+r12(i)*g21(i)+r13(i)*g31(i)
215 r12(i)=r11(i)*g12(i)+r12(i)*g22(i)+r13(i)*g32(i)
216 r13(i)=r11(i)*g13(i)+r12(i)*g23(i)+r13(i)*g33(i)
217 r21(i)=r21(i)*g11(i)+r22(i)*g21(i)+r23(i)*g31(i)
218 r22(i)=r21(i)*g12(i)+r22(i)*g22(i)+r23(i)*g32(i)
219 r23(i)=r21(i)*g13(i)+r22(i)*g23(i)+r23(i)*g33(i)
220 r31(i)=r31(i)*g11(i)+r32(i)*g21(i)+r33(i)*g31(i)
221 r32(i)=r31(i)*g12(i)+r32(i)*g22(i)+r33(i)*g32(i)
222 r33(i)=r31(i)*g13(i)+r32(i)*g23(i)+r33(i)*g33(i)
223 END DO
224 END IF
225c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
226 RETURN