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

Go to the source code of this file.

Functions/Subroutines

subroutine szhour_ctl (pm, rho, off, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, px1h1, px1h2, px1h3, px2h1, px2h2, px2h3, px3h1, px3h2, px3h3, px4h1, px4h2, px4h3, vol, fhour, mtn, dt1, mat, cxx, eint, npropm, nummat, strhg, jr, js, jt, vol0, dn, sti, nel)

Function/Subroutine Documentation

◆ szhour_ctl()

subroutine szhour_ctl ( intent(in) pm,
intent(in) rho,
intent(in) off,
intent(in) vx1,
intent(in) vx2,
intent(in) vx3,
intent(in) vx4,
intent(in) vx5,
intent(in) vx6,
intent(in) vx7,
intent(in) vx8,
intent(in) vy1,
intent(in) vy2,
intent(in) vy3,
intent(in) vy4,
intent(in) vy5,
intent(in) vy6,
intent(in) vy7,
intent(in) vy8,
intent(in) vz1,
intent(in) vz2,
intent(in) vz3,
intent(in) vz4,
intent(in) vz5,
intent(in) vz6,
intent(in) vz7,
intent(in) vz8,
intent(inout) f11,
intent(inout) f21,
intent(inout) f31,
intent(inout) f12,
intent(inout) f22,
intent(inout) f32,
intent(inout) f13,
intent(inout) f23,
intent(inout) f33,
intent(inout) f14,
intent(inout) f24,
intent(inout) f34,
intent(inout) f15,
intent(inout) f25,
intent(inout) f35,
intent(inout) f16,
intent(inout) f26,
intent(inout) f36,
intent(inout) f17,
intent(inout) f27,
intent(inout) f37,
intent(inout) f18,
intent(inout) f28,
intent(inout) f38,
intent(in) px1h1,
intent(in) px1h2,
intent(in) px1h3,
intent(in) px2h1,
intent(in) px2h2,
intent(in) px2h3,
intent(in) px3h1,
intent(in) px3h2,
intent(in) px3h3,
intent(in) px4h1,
intent(in) px4h2,
intent(in) px4h3,
intent(in) vol,
intent(inout) fhour,
integer, intent(in) mtn,
intent(in) dt1,
integer, dimension(mvsiz) mat,
intent(in) cxx,
intent(inout) eint,
integer, intent(in) npropm,
integer, intent(in) nummat,
intent(inout) strhg,
intent(in) jr,
intent(in) js,
intent(in) jt,
intent(in) vol0,
intent(in) dn,
intent(inout) sti,
integer, intent(in) nel )

Definition at line 33 of file szhour_ctl.F.

54C-----------------------------------------------
55C M o d u l e s
56C-----------------------------------------------
57 USE shour_ctl_mod, ONLY : shour_ctl
58C-----------------------------------------------
59C I m p l i c i t T y p e s
60C-----------------------------------------------
61#include "implicit_f.inc"
62C-----------------------------------------------
63C G l o b a l P a r a m e t e r s
64C-----------------------------------------------
65#include "mvsiz_p.inc"
66C-----------------------------------------------
67C C o m m o n B l o c k s
68C-----------------------------------------------
69#include "scr14_c.inc"
70C-----------------------------------------------
71C D u m m y A r g u m e n t s
72C-----------------------------------------------
73 INTEGER, INTENT(IN) :: NEL
74 INTEGER, INTENT(IN) :: MTN,NPROPM,NUMMAT
75 INTEGER, DIMENSION(MVSIZ) :: MAT
76 my_real, DIMENSION(NPROPM,NUMMAT) ,INTENT(IN):: pm
77 my_real, DIMENSION(NEL,3,4) ,INTENT(INOUT):: fhour
78 my_real, DIMENSION(NEL,18) ,INTENT(INOUT):: strhg
79 my_real, DIMENSION(NEL) ,INTENT(IN):: rho,vol0
80 my_real, DIMENSION(MVSIZ) ,INTENT(IN):: off,vol,
81 . vx1,vx2,vx3,vx4,vx5,vx6,vx7,vx8,
82 . vy1,vy2,vy3,vy4,vy5,vy6,vy7,vy8,
83 . vz1,vz2,vz3,vz4,vz5,vz6,vz7,vz8,
84 . px1h1, px1h2, px1h3,
85 . px2h1, px2h2, px2h3,
86 . px3h1, px3h2, px3h3,
87 . px4h1, px4h2, px4h3,cxx,jr,js,jt
88 my_real, DIMENSION(MVSIZ) ,INTENT(INOUT)::
89 . f11,f21,f31,f12,f22,f32,
90 . f13,f23,f33,f14,f24,f34,
91 . f15,f25,f35,f16,f26,f36,
92 . f17,f27,f37,f18,f28,f38,eint,sti
93 my_real, INTENT(IN):: dt1,dn
94C-----------------------------------------------
95C L o c a l V a r i a b l e s
96C-----------------------------------------------
97 INTEGER I, MX, J, II, IC
98C REAL
99 my_real
100 . caq(mvsiz), fcl(mvsiz), edt(mvsiz),
101 . hx1(mvsiz), hx2(mvsiz), hx3(mvsiz), hx4(mvsiz),
102 . hy1(mvsiz), hy2(mvsiz), hy3(mvsiz), hy4(mvsiz),
103 . hz1(mvsiz), hz2(mvsiz), hz3(mvsiz), hz4(mvsiz),
104 . hgx1(mvsiz), hgx2(mvsiz), hgx3(mvsiz), hgx4(mvsiz),
105 . hgy1(mvsiz), hgy2(mvsiz), hgy3(mvsiz), hgy4(mvsiz),
106 . hgz1(mvsiz), hgz2(mvsiz), hgz3(mvsiz), hgz4(mvsiz),
107 . g11(mvsiz),g21(mvsiz),g31(mvsiz),g41(mvsiz),
108 . g51(mvsiz),g61(mvsiz),g71(mvsiz),g81(mvsiz),
109 . g12(mvsiz),g22(mvsiz),g32(mvsiz),g42(mvsiz),
110 . g52(mvsiz),g62(mvsiz),g72(mvsiz),g82(mvsiz),
111 . g13(mvsiz),g23(mvsiz),g33(mvsiz),g43(mvsiz),
112 . g53(mvsiz),g63(mvsiz),g73(mvsiz),g83(mvsiz),nu2(mvsiz),
113 . jr_1(mvsiz),js_1(mvsiz),jt_1(mvsiz),nu4(mvsiz),
114 . e0,g0,c1,nu,ehourt ,qh,lamg ,stif ,ll ,fvl,nu1
115C-----------------------------------------------
116 DO i=1,nel
117 f11(i)=zero
118 f12(i)=zero
119 f13(i)=zero
120 f14(i)=zero
121 f15(i)=zero
122 f16(i)=zero
123 f17(i)=zero
124 f18(i)=zero
125C
126 f21(i)=zero
127 f22(i)=zero
128 f23(i)=zero
129 f24(i)=zero
130 f25(i)=zero
131 f26(i)=zero
132 f27(i)=zero
133 f28(i)=zero
134C
135 f31(i)=zero
136 f32(i)=zero
137 f33(i)=zero
138 f34(i)=zero
139 f35(i)=zero
140 f36(i)=zero
141 f37(i)=zero
142 f38(i)=zero
143 ENDDO
144 CALL shour_ctl(
145 . pm, rho, off, vx1,
146 . vx2, vx3, vx4, vx5,
147 . vx6, vx7, vx8, vy1,
148 . vy2, vy3, vy4, vy5,
149 . vy6, vy7, vy8, vz1,
150 . vz2, vz3, vz4, vz5,
151 . vz6, vz7, vz8, f11,
152 . f21, f31, f12, f22,
153 . f32, f13, f23, f33,
154 . f14, f24, f34, f15,
155 . f25, f35, f16, f26,
156 . f36, f17, f27, f37,
157 . f18, f28, f38, px1h1,
158 . px1h2, px1h3, px2h1, px2h2,
159 . px2h3, px3h1, px3h2, px3h3,
160 . px4h1, px4h2, px4h3, vol,
161 . fhour, mtn, dt1 , mat,
162 . cxx, eint, npropm, nummat,
163 . vol0, dn , sti , nel )
164!
165 IF(((anim_n(iad_gps+400+1) == 1) .OR. (anim_n(iad_gps+400+2) == 1) .OR.
166 . (anim_n(iad_gps+400+3) == 1) .OR. (anim_n(iad_gps+400+4) == 1) .OR.
167 . (anim_n(iad_gps+400+5) == 1) .OR. (anim_n(iad_gps+400+6) == 1)) )THEN
168 mx = mat(1)
169 nu=pm(21,mx)
170 nu1 =two/(one-nu)
171 DO i=1,nel
172 jr_1(i) = one/max(em20,jr(i))
173 js_1(i) = one/max(em20,js(i))
174 jt_1(i) = one/max(em20,jt(i))
175 nu2(i) =half*nu*nu1
176 nu4(i) =nu
177 ENDDO
178 CALL szstrainhg(
179 1 jr_1, js_1, jt_1, strhg,
180 2 nel, hgx1, hgx2, hgx3,
181 3 hgx4, hgy1, hgy2, hgy3,
182 4 hgy4, hgz1, hgz2, hgz3,
183 5 hgz4, nu4, nu2)
184 ENDIF
185C
186 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
max
#define max(a, b)
Definition macros.h:21
szstrainhg
subroutine szstrainhg(jr_1, js_1, jt_1, strhg, nel, hgx1, hgx2, hgx3, hgx4, hgy1, hgy2, hgy3, hgy4, hgz1, hgz2, hgz3, hgz4, nu, nu1)
Definition szstrainhg.F:37