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

Go to the source code of this file.

Functions/Subroutines

subroutine sigeps01g (jft, jlt, g, thk, off, gs, a1, a2, nu, thk0, nel, f_def, ismstr, depsxx, depsyy, depsxy, depsyz, depszx, depbxx, depbyy, depbxy, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, momoxx, momoyy, momoxy, signxx, signyy, signxy, signyz, signzx, momnxx, momnyy, momnxy, zshift)

Function/Subroutine Documentation

◆ sigeps01g()

subroutine sigeps01g ( integer jft,
integer jlt,
g,
thk,
off,
gs,
a1,
a2,
nu,
thk0,
integer nel,
f_def,
integer ismstr,
depsxx,
depsyy,
depsxy,
depsyz,
depszx,
depbxx,
depbyy,
depbxy,
sigoxx,
sigoyy,
sigoxy,
sigoyz,
sigozx,
momoxx,
momoyy,
momoxy,
signxx,
signyy,
signxy,
signyz,
signzx,
momnxx,
momnyy,
momnxy,
intent(in) zshift )

Definition at line 28 of file sigeps01g.F.

36C-----------------------------------------------
37C I m p l i c i t T y p e s
38C-----------------------------------------------
39#include "implicit_f.inc"
40C-----------------------------------------------
41C G l o b a l P a r a m e t e r s
42C-----------------------------------------------
43#include "mvsiz_p.inc"
44C-----------------------------------------------
45C D u m m y A r g u m e n t s
46C-----------------------------------------------
47 INTEGER JFT,JLT,NEL,ISMSTR
48C REAL
49 my_real
50 . off(*),gs(*)
51 my_real
52 . nu(mvsiz),g(mvsiz),a1(mvsiz),a2(mvsiz),f_def(mvsiz,8),
53 . depsxx(nel),depsyy(nel),depsxy(nel),
54 . depbxx(nel),depbyy(nel),depbxy(nel),
55 . depsyz(nel),depszx(nel),
56 . sigoxx(nel),sigoyy(nel),sigoxy(nel),
57 . sigoyz(nel),sigozx(nel),
58 . momoxx(nel),momoyy(nel),momoxy(nel)
59 my_real, intent(in) :: zshift
60C-----------------------------------------------
61C O U T P U T A R G U M E N T S
62C-----------------------------------------------
63 my_real
64 . thk(*),
65 . signxx(nel),signyy(nel),signxy(nel),
66 . momnxx(nel),momnyy(nel),momnxy(nel),
67 . signyz(nel),signzx(nel)
68C-----------------------------------------------
69C L o c a l V a r i a b l e s
70C-----------------------------------------------
71 INTEGER I
72C REAL
73 my_real
74 . b1(mvsiz),b2(mvsiz),thk0(mvsiz),
75 . b3(mvsiz),
76 . thk08,ezz(mvsiz)
77 my_real
78 . eigv(mvsiz,3,2),trav(mvsiz),rootv(mvsiz),
79 . eh1(mvsiz),eh2(mvsiz),eb1(mvsiz),eb2(mvsiz),
80 . epsxx(mvsiz),epsyy(mvsiz),epsxy(mvsiz),
81 . for1(mvsiz),for2(mvsiz),mom1(mvsiz),mom2(mvsiz)
82c . TAN_A,TAN_A2,SIN2_A(MVSIZ),COS2_A(MVSIZ),SINCOS_A(MVSIZ),ANG
83C-----------------------------------------------
84!---
85 IF (ismstr == 10) THEN
86C---------[F]=[F_DEF]+[1]; [B]=[F][F]^t strain----
87 DO i=jft,jlt
88 epsxx(i)=f_def(i,1)*(two+f_def(i,1))+
89 . f_def(i,3)*f_def(i,3) + one
90 epsyy(i)=f_def(i,2)*(two+f_def(i,2))+
91 . f_def(i,4)*f_def(i,4) + one
92c EPSZZ(I)=MFZZ(I)*(TWO+MFZZ(I))+
93c . MFZX(I)*MFZX(I)+MFZY(I)*MFZY(I)
94 epsxy(i)=f_def(i,3)+f_def(i,4)+f_def(i,1)*f_def(i,4)+
95 . f_def(i,3)*f_def(i,2)
96 ENDDO
97C principal stretch (def gradient eigenvalues)
98 DO i=jft,jlt
99 trav(i) = epsxx(i)+epsyy(i)
100 rootv(i) = sqrt((epsxx(i)-epsyy(i))*(epsxx(i)-epsyy(i))
101 . + epsxy(i)*epsxy(i)*four)
102 eb1(i) = half*(trav(i)+rootv(i))
103 eb2(i) = half*(trav(i)-rootv(i))
104 eh1(i) = half*log(eb1(i))
105 eh2(i) = half*log(eb2(i))
106 for1(i)=a1(i)*eh1(i)+a2(i)*eh2(i)
107 for2(i)=a1(i)*eh2(i)+a2(i)*eh1(i)
108 ENDDO
109C rot matrix (eigenvectors), ROOTV= 2R
110 DO i=jft,jlt
111 IF (abs(eb2(i)-eb1(i)) < em15) THEN
112 eigv(i,1,1)=one
113 eigv(i,2,1)=zero
114 eigv(i,3,1)=zero
115 eigv(i,1,2)=zero
116 eigv(i,2,2)=one
117 eigv(i,3,2)=zero
118 ELSE
119 eigv(i,1,1) = (epsxx(i)-eb2(i)) / rootv(i)
120 eigv(i,2,1) = (epsyy(i)-eb2(i)) / rootv(i)
121 eigv(i,3,1) = (epsxy(i)) / rootv(i)
122 eigv(i,1,2) = (eb1(i)-epsxx(i)) / rootv(i)
123 eigv(i,2,2) = (eb1(i)-epsyy(i)) / rootv(i)
124 eigv(i,3,2) =-(epsxy(i)) / rootv(i)
125 ENDIF
126 ENDDO
127 ENDIF ! IF (ISMSTR == 10) THEN
128!---
129C
130 DO i=jft,jlt
131 thk08 = thk0(i)*(one_over_12+zshift*zshift)
132 b1(i) = a1(i)*thk08
133 b2(i) = a2(i)*thk08
134 b3(i) = g(i) *thk08
135 ENDDO
136C
137 IF (ismstr /= 10) THEN
138 DO i=jft,jlt
139 signxx(i)=sigoxx(i)+a1(i)*depsxx(i)+a2(i)*depsyy(i)
140 signyy(i)=sigoyy(i)+a2(i)*depsxx(i)+a1(i)*depsyy(i)
141 signxy(i)=sigoxy(i)+g(i) *depsxy(i)
142 signyz(i)=sigoyz(i)+gs(i)*depsyz(i)
143 signzx(i)=sigozx(i)+gs(i)*depszx(i)
144C
145 momnxx(i)=momoxx(i)+b1(i)*depbxx(i)+b2(i)*depbyy(i)
146 momnyy(i)=momoyy(i)+b1(i)*depbyy(i)+b2(i)*depbxx(i)
147 momnxy(i)=momoxy(i)+b3(i)*depbxy(i)
148 ENDDO
149 ELSE ! ISMSTR = 10
150 DO i=jft,jlt
151 signxx(i) = eigv(i,1,1)*for1(i) + eigv(i,1,2)*for2(i)
152 signyy(i) = eigv(i,2,1)*for1(i) + eigv(i,2,2)*for2(i)
153 signxy(i) = eigv(i,3,1)*for1(i) + eigv(i,3,2)*for2(i)
154 signyz(i) = sigoyz(i)+ gs(i)*depsyz(i)
155 signzx(i) = sigozx(i)+ gs(i)*depszx(i)
156 ENDDO
157C-----curveture
158C---------[F]=[F_DEF]+[1]; [B]=[F][F]^t strain-----
159 DO i=jft,jlt
160 epsxx(i)=f_def(i,6)*(two+f_def(i,6))+
161 . f_def(i,8)*f_def(i,8) +one
162 epsyy(i)=f_def(i,7)*(two+f_def(i,7))+
163 . f_def(i,5)*f_def(i,5) +one
164 epsxy(i)=f_def(i,8)+f_def(i,5)+f_def(i,6)*f_def(i,5)+
165 . f_def(i,8)*f_def(i,7)
166 ENDDO
167C principal stretch (def gradient eigenvalues)
168 DO i=jft,jlt
169 trav(i) = epsxx(i)+epsyy(i)
170 rootv(i) = sqrt((epsxx(i)-epsyy(i))*(epsxx(i)-epsyy(i))
171 . + epsxy(i)*epsxy(i)*four)
172 eb1(i) = half*(trav(i)+rootv(i))
173 eb2(i) = half*(trav(i)-rootv(i))
174 eh1(i) = half*log(eb1(i))
175 eh2(i) = half*log(eb2(i))
176 mom1(i)=b1(i)*eh1(i)+b2(i)*eh2(i)
177 mom2(i)=b1(i)*eh2(i)+b2(i)*eh1(i)
178 ENDDO
179C rot matrix (eigenvectors)
180 DO i=jft,jlt
181 IF (abs(eb2(i)-eb1(i)) < em15) THEN
182 eigv(i,1,1)=one
183 eigv(i,2,1)=zero
184 eigv(i,3,1)=zero
185 eigv(i,1,2)=zero
186 eigv(i,2,2)=one
187 eigv(i,3,2)=zero
188 ELSE
189 eigv(i,1,1) = (epsxx(i)-eb2(i)) / rootv(i)
190 eigv(i,2,1) = (epsyy(i)-eb2(i)) / rootv(i)
191 eigv(i,3,1) = (epsxy(i)) / rootv(i)
192 eigv(i,1,2) = (eb1(i)-epsxx(i)) / rootv(i)
193 eigv(i,2,2) = (eb1(i)-epsyy(i)) / rootv(i)
194 eigv(i,3,2) =-(epsxy(i)) / rootv(i)
195 ENDIF
196 ENDDO
197 DO i=jft,jlt
198 momnxx(i) = eigv(i,1,1)*mom1(i) + eigv(i,1,2)*mom2(i)
199 momnyy(i) = eigv(i,2,1)*mom1(i) + eigv(i,2,2)*mom2(i)
200 momnxy(i) = eigv(i,3,1)*mom1(i) + eigv(i,3,2)*mom2(i)
201 ENDDO
202 ENDIF ! IF (ISMSTR /= 10) THEN
203C
204 DO i=jft,jlt
205 ezz(i) = -nu(i) * (depsxx(i)+depsyy(i)) / (one-nu(i))
206 thk(i) = thk(i) * (one + ezz(i) * off(i))
207 ENDDO
208! coupling terms due to shell offset
209 IF (zshift /= zero) THEN
210 DO i=jft,jlt
211 thk08 = thk0(i)*zshift
212 signxx(i) = signxx(i) + thk08*(a1(i)*depbxx(i)+a2(i)*depbyy(i))
213 signyy(i) = signyy(i) + thk08*(a1(i)*depbyy(i)+a2(i)*depbxx(i))
214 signxy(i) = signxy(i) + thk08*g(i)*depbxy(i)
215!
216 momnxx(i) = momnxx(i) + zshift*(a1(i)*depsxx(i)+a2(i)*depsyy(i))
217 momnyy(i) = momnyy(i) + zshift*(a1(i)*depsyy(i)+a2(i)*depsxx(i))
218 momnxy(i) = momnxy(i) + zshift*g(i)*depsxy(i)
219 ENDDO
220 END IF
221C
222 RETURN
my_real
#define my_real
Definition cppsort.cpp:32