OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sigeps94.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!|| sigeps94 ../engine/source/materials/mat/mat094/sigeps94.f
25!||--- called by ------------------------------------------------------
26!|| mulaw ../engine/source/materials/mat_share/mulaw.F90
27!||--- calls -----------------------------------------------------
28!|| finter ../engine/source/tools/curve/finter.F
29!|| valpvec_v ../engine/source/materials/mat/mat033/sigeps33.F
30!|| valpvecdp_v ../engine/source/materials/mat/mat033/sigeps33.F
31!||====================================================================
32 SUBROUTINE sigeps94(
33 1 NEL , NUPARAM, NUVAR , NFUNC , IFUNC , NPF ,
34 2 TF , TIME , TIMESTEP, UPARAM, RHO0 , RHO ,
35 3 VOLUME , EINT , NGL ,
36 4 EPSPXX , EPSPYY , EPSPZZ , EPSPXY, EPSPYZ, EPSPZX,
37 5 DEPSXX , DEPSYY , DEPSZZ , DEPSXY, DEPSYZ, DEPSZX,
38 6 EPSXX , EPSYY , EPSZZ , EPSXY , EPSYZ , EPSZX ,
39 7 SIGOXX , SIGOYY , SIGOZZ , SIGOXY, SIGOYZ, SIGOZX,
40 8 SIGNXX , SIGNYY , SIGNZZ , SIGNXY, SIGNYZ, SIGNZX,
41 9 SIGVXX , SIGVYY , SIGVZZ , SIGVXY, SIGVYZ, SIGVZX,
42 A SOUNDSP, VISCMAX, UVAR , OFF , ISMSTR, ET ,
43 B IHET ,OFFG , EPSTH3, IEXPAN)
44C-----------------------------------------------
45C I M P L I C I T T Y P E S
46C-----------------------------------------------
47#include "implicit_f.inc"
48C-----------------------------------------------
49C G L O B A L P A R A M E T E R S
50C-----------------------------------------------
51#include "mvsiz_p.inc"
52C-----------------------------------------------
53C C O M M O N
54C-----------------------------------------------
55#include "scr05_c.inc"
56#include "impl1_c.inc"
57C----------------------------------------------------------------
58C I N P U T A R G U M E N T S
59C----------------------------------------------------------------
60 INTEGER NEL, NUPARAM, NUVAR,ISMSTR,NGL(*),IHET,IEXPAN
61 my_real
62 . TIME , TIMESTEP , UPARAM(NUPARAM),
63 . RHO (NEL), RHO0 (NEL), VOLUME(NEL), EINT(NEL),
64 . EPSPXX(NEL), EPSPYY(NEL), EPSPZZ(NEL),
65 . EPSPXY(NEL), EPSPYZ(NEL), EPSPZX(NEL),
66 . DEPSXX(NEL), DEPSYY(NEL), DEPSZZ(NEL),
67 . DEPSXY(NEL), DEPSYZ(NEL), DEPSZX(NEL),
68 . EPSXX (NEL), EPSYY (NEL), EPSZZ (NEL),
69 . EPSXY (NEL), EPSYZ (NEL), EPSZX (NEL),
70 . SIGOXX(NEL), SIGOYY(NEL), SIGOZZ(NEL),
71 . sigoxy(nel), sigoyz(nel), sigozx(nel),offg(nel),
72 . epsth3(nel)
73C----------------------------------------------------------------
74C O U T P U T A R G U M E N T S
75C----------------------------------------------------------------
76 my_real
77 . signxx(nel), signyy(nel), signzz(nel),
78 . signxy(nel), signyz(nel), signzx(nel),
79 . sigvxx(nel), sigvyy(nel), sigvzz(nel),
80 . sigvxy(nel), sigvyz(nel), sigvzx(nel),
81 . soundsp(nel), viscmax(nel), et(nel)
82C----------------------------------------------------------------
83C I N P U T O U T P U T A R G U M E N T S
84C----------------------------------------------------------------
85 my_real
86 . uvar(nel,nuvar), off(nel)
87C----------------------------------------------------------------
88C VARIABLES FOR FUNCTION INTERPOLATION
89C----------------------------------------------------------------
90 INTEGER NPF(*), NFUNC, IFUNC(NFUNC)
91 my_real FINTER,FINTTE,TF(*),FINT2V
92 EXTERNAL FINTER,FINTTE
93C----------------------------------------------------------------
94C L O C A L V A R I B L E S
95C----------------------------------------------------------------
96 INTEGER I,J,II
97
98 my_real
99 . G,RBULK,AA,BB,CC,P,TRACE,C(5),
100 . T1(MVSIZ), T2(MVSIZ),T3(MVSIZ),AV(MVSIZ,6),EV(MVSIZ,3),
101 . evv(mvsiz,3),rv(mvsiz),rvd,l1di1lam1,l2di1lam2,l3di1lam3,
102 . dirprv(mvsiz,3,3),evd(3),c0(3),d(3),invr,eti,et1,et2,et3,
103 . ev2(3),clam(3),lam_2(3),lam_4(3),amax,clp
104 my_real ,DIMENSION(NEL,3) :: evm,cii
105 my_real ,DIMENSION(NEL) :: gtmax,rkmax,trac3
106C----------------------------------------------------------------
107C material parameters
108 g = uparam(1)
109 rbulk = uparam(2)
110 c0(1) = uparam(4)
111 c0(2) = uparam(5)
112 c0(3) = uparam(6)
113C
114 d(1) = uparam(7)
115 d(2) = uparam(8)
116 d(3) = uparam(9)
117
118
119C iniialisation des variabesl users
120
121 IF(time==zero)THEN
122 DO j = 1,nuvar
123 DO i = 1, nel
124 uvar(i,j) = zero
125 ENDDO
126 ENDDO
127 ENDIF
128C
129 DO i=1,nel
130 av(i,1)=epsxx(i)
131 av(i,2)=epsyy(i)
132 av(i,3)=epszz(i)
133 av(i,4)=epsxy(i)/2
134 av(i,5)=epsyz(i)/2
135 av(i,6)=epszx(i)/2
136 ENDDO
137CEigenvalues needed to be calculated in double precision
138C for a simple precision executing
139 IF (iresp==1) THEN
140 CALL valpvecdp_v(av,evv,dirprv,nel)
141 ELSE
142 CALL valpvec_v(av,evv,dirprv,nel)
143 ENDIF
144C-ISMSTR=0-NO SMALL STRAIN OPTION:STRAINS ARE LOGARITHMIC, STRESS IS CAUCHY
145C-ISMSTR=1-SMALL STRAIN OPTION:STRAINS ARE ENGINEERING, STRESS IS CAUCHY
146C-ISMSTR=2-SMALL STRAIN OPTION:STRAINS ARE ENGINEERING, STRESS IS BIOT
147C-ISMSTR=3-NO SMALL STRAIN OPTION:STRESS IS BIOT
148 IF(ismstr==0.OR.ismstr==2.OR.ismstr==4) THEN
149 DO i=1,nel
150C ---- (STRAIN IS LOGARITHMIC)
151 ev(i,1)=exp(evv(i,1))
152 ev(i,2)=exp(evv(i,2))
153 ev(i,3)=exp(evv(i,3))
154 ENDDO
155 ELSEIF(ismstr==10) THEN
156 DO i =1,nel
157 IF(offg(i)<=one) THEN
158 ev(i,1)=sqrt(evv(i,1)+ one)
159 ev(i,2)=sqrt(evv(i,2)+ one)
160 ev(i,3)=sqrt(evv(i,3)+ one)
161 ELSE
162 ev(i,1)=evv(i,1)+ one
163 ev(i,2)=evv(i,2)+ one
164 ev(i,3)=evv(i,3)+ one
165 END IF
166 ENDDO
167 ELSE
168C ---- STRAIN IS ENGINEERING)
169 DO i=1,nel
170 ev(i,1)=evv(i,1)+ one
171 ev(i,2)=evv(i,2)+ one
172 ev(i,3)=evv(i,3)+ one
173 ENDDO
174 ENDIF
175
176 IF (impl_s > 0 .OR. ihet > 1) THEN
177 DO i=1,nel
178C
179 ev2(1) = ev(i,1)*ev(i,1)
180 ev2(2) = ev(i,2)*ev(i,2)
181 ev2(3) = ev(i,3)*ev(i,3)
182C
183 trace = ev2(1) + ev2(2) + ev2(3)
184 aa = (trace - three)
185 bb = aa*aa
186 cc = (c0(1) + two*c0(2)*aa + three*c0(3)*bb)
187 et1 = ev2(1)*cc + ev2(1)**2*(two*c0(2) + six*c0(3)*aa)
188 et2 = ev2(2)*cc + ev2(2)**2*(two*c0(2) + six*c0(3)*aa)
189 et3 = ev2(3)*cc + ev2(3)**2*(two*c0(2) + six*c0(3)*aa)
190 et(i)= max(et1,et2,et3)
191 et(i) = four*et(i)
192 et(i)= max(one,et(i)/max(em20,g))
193 ENDDO
194 ENDIF
195C----------------
196C
197C
198 DO i = 1,nel
199C---- RV = RHO0/RHO = RELATIVE VOLUME = DET A (A = GRADIENT OF DEFORMATION)
200 rv(i) = ev(i,1)*ev(i,2)*ev(i,3)
201 ENDDO
202C----THERM STRESS COMPUTATION-----
203 IF(iexpan > 0.AND.(ismstr==10.OR.ismstr==11.OR.ismstr==12)) THEN
204 DO i=1,nel
205 rv(i) = rv(i) -epsth3(i)
206 ENDDO
207 ENDIF
208C----------------
209 DO i = 1,nel
210 invr = one/(max(rv(i), em20))
211
212! RVD = RV**(-THIRD)
213 ! if rv > 0 --> rvd = exp(-third * ln(rv))
214 ! else --> rvd = 0
215 rvd = zero
216 IF(rv(i) > zero) rvd = exp((-third)*log(rv(i)))
217
218 evd(1) = ev(i,1)*rvd
219 evd(2) = ev(i,2)*rvd
220 evd(3) = ev(i,3)*rvd
221 evm(i,1:3) = evd(1:3)
222C
223 trace = evd(1)*evd(1) + evd(2)*evd(2) + evd(3)*evd(3)
224 l1di1lam1 = two*(evd(1)**2 - third*trace)
225 l2di1lam2 = two*(evd(2)**2 - third*trace)
226 l3di1lam3 = two*(evd(3)**2 - third*trace)
227
228 aa = (trace - three)
229 bb = aa*aa
230 cc = (c0(1) + two*c0(2)*aa + three*c0(3)*bb)*invr
231 trac3(i) = aa
232C
233
234 t1(i) = l1di1lam1*cc
235 t2(i) = l2di1lam2*cc
236 t3(i) = l3di1lam3*cc
237C RBilk = 2/D
238 p = rbulk*(rv(i) - one) + four*d(2)*(rv(i) - one)**3 + six*d(3)*(rv(i) - one)**5
239
240 t1(i) = t1(i) + p
241 t2(i) = t2(i) + p
242 t3(i) = t3(i) + p
243 ENDDO
244C
245c----compute GT_MAX
246 gtmax(1:nel) = g
247 rkmax(1:nel) = rbulk
248 cii(1:nel,1:3) = zero
249 DO ii = 1,3
250 clp = four*ii*c0(ii)
251 DO i=1,nel
252 IF (abs(trac3(i))<em10) cycle
253 lam_2(1:3) = evm(i,1:3)*evm(i,1:3)
254 lam_4(1:3) = lam_2(1:3)*lam_2(1:3)
255 aa = one_over_9*ii*trac3(i)**ii
256 bb = zero
257 cc = zero
258 IF (ii>1) bb =third*(3-ii)*trac3(i)**(ii-1)
259 IF (ii>2) cc =(ii-1)*trac3(i)**(ii-2)
260 cii(i,1:3) = cii(i,1:3) +clp*(aa+bb*lam_2(1:3) +cc*lam_4(1:3))
261 ENDDO
262 ENDDO
263 DO i=1,nel
264 amax= max(cii(i,1),cii(i,2),cii(i,3))
265C---------reduce old result change
266 eti = max(one,amax*0.81)
267 gtmax(i) = g*eti
268C---- see if many QA return for HEPH
269 rkmax(i) = rbulk+twelve*d(2)*(rv(i) - one)**2 + 30*d(3)*(rv(i) - one)**4
270 rkmax(i) = max(rbulk,rkmax(i))
271 et(i)= max(one,amax)
272 ENDDO
273C cauchy to glabale
274 DO i=1,nel
275 signxx(i) = dirprv(i,1,1)*dirprv(i,1,1)*t1(i)
276 . + dirprv(i,1,2)*dirprv(i,1,2)*t2(i)
277 . + dirprv(i,1,3)*dirprv(i,1,3)*t3(i)
278
279 signyy(i) = dirprv(i,2,2)*dirprv(i,2,2)*t2(i)
280 . + dirprv(i,2,3)*dirprv(i,2,3)*t3(i)
281 . + dirprv(i,2,1)*dirprv(i,2,1)*t1(i)
282
283 signzz(i) = dirprv(i,3,3)*dirprv(i,3,3)*t3(i)
284 . + dirprv(i,3,1)*dirprv(i,3,1)*t1(i)
285 . + dirprv(i,3,2)*dirprv(i,3,2)*t2(i)
286
287 signxy(i) = dirprv(i,1,1)*dirprv(i,2,1)*t1(i)
288 . + dirprv(i,1,2)*dirprv(i,2,2)*t2(i)
289 . + dirprv(i,1,3)*dirprv(i,2,3)*t3(i)
290
291 signyz(i) = dirprv(i,2,2)*dirprv(i,3,2)*t2(i)
292 . + dirprv(i,2,3)*dirprv(i,3,3)*t3(i)
293 . + dirprv(i,2,1)*dirprv(i,3,1)*t1(i)
294
295 signzx(i) = dirprv(i,3,3)*dirprv(i,1,3)*t3(i)
296 . + dirprv(i,3,1)*dirprv(i,1,1)*t1(i)
297 . + dirprv(i,3,2)*dirprv(i,1,2)*t2(i)
298C
299* SET SOUND SPEED
300 soundsp(i)=sqrt((four_over_3*gtmax(i) + rkmax(i))/rho(i))
301* SET VISCOSITY
302 viscmax(i) = zero
303 ENDDO
304C
305 RETURN
306 END
307
max
#define max(a, b)
Definition macros.h:21
valpvecdp_v
subroutine valpvecdp_v(sig, val, vec, nel)
Definition sigeps33.F:2231
valpvec_v
subroutine valpvec_v(sig, val, vec, nel)
Definition sigeps33.F:1902
sigeps94
subroutine sigeps94(nel, nuparam, nuvar, nfunc, ifunc, npf, tf, time, timestep, uparam, rho0, rho, volume, eint, ngl, epspxx, epspyy, epspzz, epspxy, epspyz, epspzx, depsxx, depsyy, depszz, depsxy, depsyz, depszx, epsxx, epsyy, epszz, epsxy, epsyz, epszx, sigoxx, sigoyy, sigozz, sigoxy, sigoyz, sigozx, signxx, signyy, signzz, signxy, signyz, signzx, sigvxx, sigvyy, sigvzz, sigvxy, sigvyz, sigvzx, soundsp, viscmax, uvar, off, ismstr, et, ihet, offg, epsth3, iexpan)
Definition sigeps94.F:44