OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sigeps62.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!|| sigeps62 ../engine/source/materials/mat/mat062/sigeps62.F
25!||--- called by ------------------------------------------------------
26!|| mulaw ../engine/source/materials/mat_share/mulaw.F90
27!|| mulaw8 ../engine/source/materials/mat_share/mulaw8.F90
28!||--- calls -----------------------------------------------------
29!|| finter ../engine/source/tools/curve/finter.F
30!|| valpvec_v ../engine/source/materials/mat/mat033/sigeps33.F
31!|| valpvecdp_v ../engine/source/materials/mat/mat033/sigeps33.F
32!||====================================================================
33 SUBROUTINE sigeps62(
34 1 NEL , NUPARAM, NUVAR , NFUNC , IFUNC , NPF ,
35 2 TF , TIME , TIMESTEP, UPARAM, RHO0 , RHO ,
36 3 VOLUME , EINT ,
37 4 EPSPXX , EPSPYY , EPSPZZ , EPSPXY, EPSPYZ, EPSPZX,
38 5 DEPSXX , DEPSYY , DEPSZZ , DEPSXY, DEPSYZ, DEPSZX,
39 6 EPSXX , EPSYY , EPSZZ , EPSXY , EPSYZ , EPSZX ,
40 7 SIGOXX , SIGOYY , SIGOZZ , SIGOXY, SIGOYZ, SIGOZX,
41 8 SIGNXX , SIGNYY , SIGNZZ , SIGNXY, SIGNYZ, SIGNZX,
42 9 SIGVXX , SIGVYY , SIGVZZ , SIGVXY, SIGVYZ, SIGVZX,
43 A SOUNDSP, VISCMAX, UVAR , OFF , ISMSTR, ET ,
44 B IHET ,OFFG , EPSTH3 , IEXPAN)
45C-----------------------------------------------
46C I M P L I C I T T Y P E S
47C-----------------------------------------------
48#include "implicit_f.inc"
49C-----------------------------------------------
50C G L O B A L P A R A M E T E R S
51C-----------------------------------------------
52#include "mvsiz_p.inc"
53C-----------------------------------------------
54C C O M M O N
55C-----------------------------------------------
56#include "scr05_c.inc"
57#include "impl1_c.inc"
58C----------------------------------------------------------------
59C I N P U T A R G U M E N T S
60C----------------------------------------------------------------
61 INTEGER NEL, NUPARAM, NUVAR,ISMSTR,IHET,IEXPAN
62
63 my_real
64 . TIME , TIMESTEP , UPARAM(NUPARAM),
65 . RHO (NEL), RHO0 (NEL), VOLUME(NEL), EINT(NEL),
66 . EPSPXX(NEL), EPSPYY(NEL), EPSPZZ(NEL),
67 . EPSPXY(NEL), EPSPYZ(NEL), EPSPZX(NEL),
68 . DEPSXX(NEL), DEPSYY(NEL), DEPSZZ(NEL),
69 . DEPSXY(NEL), DEPSYZ(NEL), DEPSZX(NEL),
70 . EPSXX (NEL), EPSYY (NEL), EPSZZ (NEL),
71 . EPSXY (NEL), EPSYZ (NEL), EPSZX (NEL),
72 . sigoxx(nel), sigoyy(nel), sigozz(nel),
73 . sigoxy(nel), sigoyz(nel), sigozx(nel), offg(nel),
74 . epsth3(nel)
75C----------------------------------------------------------------
76C O U T P U T A R G U M E N T S
77C----------------------------------------------------------------
78 my_real
79 . signxx(nel), signyy(nel), signzz(nel),
80 . signxy(nel), signyz(nel), signzx(nel),
81 . sigvxx(nel), sigvyy(nel), sigvzz(nel),
82 . sigvxy(nel), sigvyz(nel), sigvzx(nel),
83 . soundsp(nel), viscmax(nel), et(nel)
84C----------------------------------------------------------------
85C I N P U T O U T P U T A R G U M E N T S
86C----------------------------------------------------------------
87 my_real
88 . uvar(nel,nuvar), off(nel)
89C----------------------------------------------------------------
90C VARIABLES FOR FUNCTION INTERPOLATION
91C----------------------------------------------------------------
92 INTEGER NPF(*), NFUNC, IFUNC(NFUNC)
93 my_real FINTER,FINTTE,TF(*),FINT2V
94 EXTERNAL FINTER,FINTTE
95C----------------------------------------------------------------
96C L O C A L V A R I B L E S
97C----------------------------------------------------------------
98 INTEGER I,J,K,II,NORDRE,NPRONY,ITYPE,IVISC
99 my_real
100 . MU(100),AL(100),TAUX(100), GAMA(100),GAMAINF,NU,
101 . GMAX,S(MVSIZ,3),SD(MVSIZ,3),
102 . EC(MVSIZ,3),HD(MVSIZ,3,100),
103 . ev(mvsiz,3),av(mvsiz,6),dirprv(mvsiz,3,3),
104 . rv(mvsiz),p(mvsiz),rbulk,
105 . evv(mvsiz,3),ssp(mvsiz,3),fac, beta(100),
106 . dtinv,visc,fac1,
107 . e1,e2,e3,e4,e5,e6,epsp, vm,dav,visc1,sigxx(mvsiz),
108 . sigyy(mvsiz),sigzz(mvsiz),sigxy(mvsiz),sigzx(mvsiz),
109 . sigyz(mvsiz),dvisc,a11,a12,a13,a22,a21,a23,
110 . a31,a32,a33,sdg0(mvsiz,6),sdg(mvsiz,6),hg(mvsiz,6),hg0(6),h(mvsiz,3,100)
111
112 my_real
113 . a(3,3),dpra(3,3),eigen(3),amax,eti(mvsiz),efac,rv_pui,rv_pui_tab(mvsiz),
114 . cj(mvsiz),lam_al(3),rvl,lamin,cimax(mvsiz),gvis,cmax,cmax0,
115 . d11,d12,d13,d21,d22,d23,d31,d32,d33,mu2
116 my_real ,DIMENSION(NEL,3) :: ai,bi
117C----------------------------------------------------------------
118
119C SET INITIAL MATERIAL CONSTANTS
120 gmax = zero
121 nu = uparam(1)
122 nordre = nint(uparam(2))
123 nprony = nint(uparam(3))
124 gamainf = uparam(4)
125 rbulk = uparam(5)
126 dvisc = uparam(6)
127!! BETA = NU/(ONE - TWO*NU)
128 DO i= 1,nordre
129 mu(i) = uparam(6 + i)
130 al(i) = uparam(6 + nordre + i)
131 beta(i) = uparam(2*nordre + 2*nprony + 7 + i)
132C GMAX = GMAX + AL(I)*MU(I)
133 gmax = gmax + mu(i)
134 ENDDO
135!! IVISC = 0
136 ivisc = uparam(2*nordre + 2*nprony + 7 )
137 IF(nprony /= zero)THEN
138 DO i=1,nprony
139 gama(i) = uparam(6 + 2*nordre + i)
140 taux(i) = uparam(6 + 2*nordre + nprony + i)
141 ENDDO
142 ENDIF
143!! IF(IVISC == 1 .AND. ITYPE >= 2 ) IVISC = 2
144C +1
145 gmax = two*gmax
146C iniialisation des variabesl users
147 IF(time==zero)THEN
148 DO j = 1,nuvar
149 DO i = 1, nel
150 uvar(i,j) = zero
151 ENDDO
152 ENDDO
153 ENDIF
154 DO i=1,nel
155 av(i,1)=epsxx(i)
156 av(i,2)=epsyy(i)
157 av(i,3)=epszz(i)
158 av(i,4)=epsxy(i)/2
159 av(i,5)=epsyz(i)/2
160 av(i,6)=epszx(i)/2
161 ENDDO
162CEigenvalues needed to be calculated in double precision
163C for a simple precision executing
164 IF (iresp==1) THEN
165 CALL valpvecdp_v(av,evv,dirprv,nel)
166 ELSE
167 CALL valpvec_v(av,evv,dirprv,nel)
168 ENDIF
169C-ISMSTR=0-NO SMALL STRAIN OPTION:STRAINS ARE LOGARITHMIC, STRESS IS CAUCHY
170C-ISMSTR=1-SMALL STRAIN OPTION:STRAINS ARE ENGINEERING, STRESS IS CAUCHY
171C-ISMSTR=2-SMALL STRAIN OPTION:STRAINS ARE ENGINEERING, STRESS IS BIOT
172C-ISMSTR=3-NO SMALL STRAIN OPTION:STRESS IS BIOT
173 IF(ismstr==0.OR.ismstr==2.OR.ismstr==4) THEN
174 DO i=1,nel
175C ---- (STRAIN IS LOGARITHMIC)
176 ev(i,1)=exp(evv(i,1))
177 ev(i,2)=exp(evv(i,2))
178 ev(i,3)=exp(evv(i,3))
179 ENDDO
180 ELSEIF(ismstr==10.OR.ismstr==12) THEN
181 DO i =1,nel
182 IF(offg(i)<=one) THEN
183 ev(i,1)=sqrt(evv(i,1)+ one)
184 ev(i,2)=sqrt(evv(i,2)+ one)
185 ev(i,3)=sqrt(evv(i,3)+ one)
186 ELSE
187 ev(i,1)=evv(i,1)+ one
188 ev(i,2)=evv(i,2)+ one
189 ev(i,3)=evv(i,3)+ one
190 END IF
191 ENDDO
192 ELSE
193C ---- STRAIN IS ENGINEERING)
194 DO i=1,nel
195 ev(i,1)=evv(i,1)+ one
196 ev(i,2)=evv(i,2)+ one
197 ev(i,3)=evv(i,3)+ one
198 ENDDO
199 ENDIF
200C----------------
201 IF (impl_s > 0 .OR. ihet > 1) THEN
202 DO j =1,3
203 eti(1:nel)=zero
204 DO k= 1,nordre
205! ETI = ETI + MU(K)*AMAX**(AL(K)-ONE)
206 DO i=1,nel
207 amax = ev(i,j)
208 IF(amax>zero) THEN
209 IF(mu(k)/=zero) THEN
210 IF((al(k)-one)/=zero) THEN
211 eti(i) = eti(i) + mu(k)*exp((al(k)-one)*log(amax))
212 ELSE
213 eti(i) = eti(i) + mu(k)
214 ENDIF
215 ELSE
216 eti(i) = eti(i) + zero
217 ENDIF
218 ELSE
219 eti(i) = zero
220 ENDIF
221 ENDDO ! NEL
222 ENDDO ! NORDRE
223 et(1:nel)= max(eti(1:nel),et(1:nel))
224 ENDDO ! 1,3
225 efac=two
226 DO i=1,nel
227 et(i)= max(one,efac*et(i)/gmax)
228 ENDDO
229 ENDIF
230C
231C---- RV = RHO0/RHO = RELATIVE VOLUME = DET A (A = GRADIENT OF DEFORMATION)
232 rv(1:nel) = ev(1:nel,1)*ev(1:nel,2)*ev(1:nel,3)
233 ec(1:nel,1:3) = ev(1:nel,1:3)**2
234C---- RV = RHO0/RHO = RELATIVE VOLUME = DET A (A = GRADIENT OF DEFORMATION)
235C----THERM STRESS COMPUTATION-----
236 IF(iexpan > 0.AND.(ismstr==10.OR.ismstr==11.OR.ismstr==12)) rv(1:nel)= rv(1:nel)-epsth3(1:nel)
237C
238C calcul de la pression
239 p(1:nel) = zero
240 DO ii = 1,nordre
241 fac = two*mu(ii)/(al(ii))
242 IF(fac/=zero) THEN
243 IF(al(ii)/=zero) THEN
244 DO i=1,nel
245 fac1 = fac / rv(i)
246 IF(ev(i,1)>zero) p(i) = p(i) +
247 . fac1*third*( exp(al(ii)*log(ev(i,1))) )
248 IF(ev(i,2)>zero) p(i) = p(i) +
249 . fac1*third*( exp(al(ii)*log(ev(i,2))) )
250 IF(ev(i,3)>zero) p(i) = p(i) +
251 . fac1*third*( exp(al(ii)*log(ev(i,3))) )
252
253 IF(rv(i)>zero) p(i) = p(i)
254 . - fac1*exp((-beta(ii)*al(ii))*log(rv(i)))
255 ENDDO
256 ELSE
257 DO i=1,nel
258 fac1 = fac / rv(i)
259 IF(ev(i,1)/=zero) p(i) = p(i) + fac1*third
260 IF(ev(i,2)/=zero) p(i) = p(i) + fac1*third
261 IF(ev(i,3)/=zero) p(i) = p(i) + fac1*third
262 IF(rv(i)>zero) p(i) = p(i) - fac1
263 ENDDO
264 ENDIF
265 ELSE
266 p(1:nel) = zero
267 ENDIF
268 ENDDO
269C compute de piola kirchoff stress 2emme espece---
270 DO j=1,3
271 s(1:nel,j)= zero
272 DO ii= 1, nordre
273 fac = two*mu(ii)/(al(ii))
274 IF(fac/=zero) THEN
275 DO i=1,nel
276 fac1 = fac/(ev(i,j)**2)
277 IF(al(ii)/=zero) THEN
278 IF(ev(i,j)>zero) s(i,j)= s(i,j) + fac1*exp(al(ii)*log(ev(i,j)))
279 IF(rv(i)>zero) s(i,j)= s(i,j) - fac1*exp((-beta(ii)*al(ii))*log(rv(i)))
280 ELSE
281 IF(ev(i,j)==zero) THEN
282 s(i,j)= s(i,j) - fac1
283 ENDIF
284 ENDIF
285 ENDDO ! I=1,NEL
286 ENDIF ! Fac /= 0
287 ENDDO ! II=1,Nordre
288 ENDDO ! J=1,3
289CCC devaitoric piola kirchoff stress 2nd espece SD
290 DO i=1,nel
291 ! sd = rv**(dtier) * sd
292 ! sd = pui * sd
293 ! if rv > 0 --> pui = exp(dtier * ln(rv))
294 ! else pui = 0
295 IF(rv(i)>zero) THEN
296 rv_pui_tab(i) = exp( two_third*log(rv(i)) )
297 ELSE
298 rv_pui_tab(i) = zero
299 ENDIF
300 ENDDO
301 DO j=1,3
302 ssp(1:nel,j) = p(1:nel)/(ev(1:nel,j)**2)
303 sd(1:nel,j) = s(1:nel,j) - ssp(1:nel,j)*rv(1:nel)
304 sd(1:nel,j) = rv_pui_tab(1:nel)*sd(1:nel,j) ! viscosity is deviatoric
305 ENDDO
306C we consider spherical and deviatoric viscosity compared to ivisc =1
307 IF(ivisc == 2) sd(1:nel,1:3) = s(1:nel,1:3)
308
309 gvis = zero
310C Calcul de H^i_n+1
311 IF(ivisc > 0)THEN
312 gvis = gmax
313c GMAX = TWO*GMAX
314 sdg0(1:nel,1:6) = uvar(1:nel,1:6)
315 DO i=1,nel
316C
317C compute SD ---> glabal repere
318C
319 sdg(i,1) = dirprv(i,1,1)*dirprv(i,1,1)*sd(i,1)
320 . + dirprv(i,1,2)*dirprv(i,1,2)*sd(i,2)
321 . + dirprv(i,1,3)*dirprv(i,1,3)*sd(i,3)
322
323 sdg(i,2) = dirprv(i,2,2)*dirprv(i,2,2)*sd(i,2)
324 . + dirprv(i,2,3)*dirprv(i,2,3)*sd(i,3)
325 . + dirprv(i,2,1)*dirprv(i,2,1)*sd(i,1)
326
327 sdg(i,3) = dirprv(i,3,3)*dirprv(i,3,3)*sd(i,3)
328 . + dirprv(i,3,1)*dirprv(i,3,1)*sd(i,1)
329 . + dirprv(i,3,2)*dirprv(i,3,2)*sd(i,2)
330
331 sdg(i,4) = dirprv(i,1,1)*dirprv(i,2,1)*sd(i,1)
332 . + dirprv(i,1,2)*dirprv(i,2,2)*sd(i,2)
333 . + dirprv(i,1,3)*dirprv(i,2,3)*sd(i,3)
334
335 sdg(i,5) = dirprv(i,2,2)*dirprv(i,3,2)*sd(i,2)
336 . + dirprv(i,2,3)*dirprv(i,3,3)*sd(i,3)
337 . + dirprv(i,2,1)*dirprv(i,3,1)*sd(i,1)
338
339 sdg(i,6) = dirprv(i,3,3)*dirprv(i,1,3)*sd(i,3)
340 . + dirprv(i,3,1)*dirprv(i,1,1)*sd(i,1)
341 . + dirprv(i,3,2)*dirprv(i,1,2)*sd(i,2)
342
343 ENDDO
344 DO j=1,6
345 uvar(1:nel,j) = sdg(1:nel,j)
346 ENDDO
347C
348 DO ii= 1,nprony
349 fac= -timestep/taux(ii)
350 DO j= 1,6
351 DO i=1,nel
352 hg0(j) = uvar(i, 6 + (ii - 1)*6 + j)
353 hg(i,j) = exp(fac)*hg0(j) + exp(half*fac)*(sdg(i,j)-sdg0(i,j))
354 uvar(i, 6 + (ii - 1)*6 + j)= hg(i,j)
355 ENDDO
356 ENDDO
357C
358C HGlobal --- H principal
359C
360 DO i=1,nel
361 a11 = dirprv(i,1,1)*hg(i,1) + dirprv(i,2,1)*hg(i,4)
362 . + dirprv(i,3,1)*hg(i,6)
363 a12 = dirprv(i,1,2)*hg(i,1) + dirprv(i,2,2)*hg(i,4)
364 . + dirprv(i,3,2)*hg(i,6)
365 a13 = dirprv(i,1,3)*hg(i,1) + dirprv(i,2,3)*hg(i,4)
366 . + dirprv(i,3,3)*hg(i,6)
367
368!
369 a21 = dirprv(i,1,1)*hg(i,4) + dirprv(i,2,1)*hg(i,2)
370 . + dirprv(i,3,1)*hg(i,5)
371 a22 = dirprv(i,1,2)*hg(i,4) + dirprv(i,2,2)*hg(i,2)
372 . + dirprv(i,3,2)*hg(i,5)
373 a23 = dirprv(i,1,3)*hg(i,4) + dirprv(i,2,3)*hg(i,2)
374 . + dirprv(i,3,3)*hg(i,5)
375
376 a31 = dirprv(i,1,1)*hg(i,6) + dirprv(i,2,1)*hg(i,5)
377 . + dirprv(i,3,1)*hg(i,3)
378 a32 = dirprv(i,1,2)*hg(i,6) + dirprv(i,2,2)*hg(i,5)
379 . + dirprv(i,3,2)*hg(i,3)
380 a33 = dirprv(i,1,3)*hg(i,6) + dirprv(i,2,3)*hg(i,5)
381 . + dirprv(i,3,3)*hg(i,3)
382C H dans la direction H
383 h(i,1,ii) = dirprv(i,1,1)*a11 + dirprv(i,2,1)*a21
384 . + dirprv(i,3,1)*a31
385 h(i,2,ii) = dirprv(i,1,2)*a12 + dirprv(i,2,2)*a22
386 . + dirprv(i,3,2)*a32
387 h(i,3,ii) = dirprv(i,1,3)*a13 + dirprv(i,2,3)*a23
388 . + dirprv(i,3,3)*a33
389 ENDDO
390 DO j=1,3
391 DO i=1,nel
392C calcul de la partie deviatorique
393 fac = third*(h(i,1,ii)*ec(i,1) + h(i,2,ii)*ec(i,2)
394 . + h(i,3,ii)*ec(i,3))
395 hd(i,j,ii) = h(i,j,ii) - fac/max(em20,ec(i,j))
396 ENDDO
397 ENDDO
398 ENDDO
399 ENDIF
400C calcul du tenseur de piola kirchoff
401 IF(ivisc == 1) THEN ! viscosity is deviatoric
402 DO i =1,nel
403 ! FAC= GAMAINF*RV(I)**(-TWO_THIRD) = GAMAINF * PUI
404 ! if rv > 0 --> pui = exp(-dtier * ln(rv))
405 ! else pui = 0
406 IF(rv(i)>zero) THEN
407 rv_pui_tab(i) = exp( (-two_third)*log(rv(i)) )
408 ELSE
409 rv_pui_tab(i) = zero
410 ENDIF
411 ENDDO
412C
413 DO j=1,3
414 DO i =1,nel
415 fac= gamainf*rv_pui_tab(i)
416 s(i,j) = ssp(i,j)*rv(i) + fac*sd(i,j)
417 ENDDO
418 DO ii=1,nprony
419 DO i =1,nel
420 s(i,j) = s(i,j) + gama(ii)*rv_pui_tab(i)*hd(i,j,ii)
421 ENDDO
422 ENDDO
423C cauchy stress from PK2
424 DO i =1,nel
425 s(i,j) = (s(i,j)*ev(i,j)**2)/rv(i)
426 ENDDO
427 ENDDO
428 ELSEIF(ivisc == 2) THEN ! spherical + deviatoric
429 DO j=1,3
430 DO i =1,nel
431 s(i,j) = gamainf*s(i,j)
432 ENDDO
433 DO ii=1,nprony
434 DO i =1,nel
435 s(i,j) = s(i,j) + gama(ii)*h(i,j,ii)
436 ENDDO
437 ENDDO
438C cauchy stress from PK2
439 DO i =1,nel
440 s(i,j) = (s(i,j)*ev(i,j)**2)/rv(i)
441 ENDDO
442 ENDDO
443 ELSE
444C cauchy stress from PK2
445 DO j=1,3
446 DO i =1,nel
447 s(i,j) = s(i,j)*ev(i,j)**2/rv(i)
448 ENDDO
449 ENDDO
450 ENDIF
451C-----
452 cmax0 = two_third*gmax+ rbulk
453 ai(1:nel,1:3) = zero
454 bi(1:nel,1:3) = zero
455 cj(1:nel) = zero
456 DO ii = 1,nordre
457 IF(al(ii)/=zero) THEN
458 fac = one/al(ii)
459 fac1 = beta(ii)+fac
460 mu2 = two*mu(ii)
461 DO i=1,nel
462 lam_al(1:3) = exp(al(ii)*log(ev(i,1:3)))
463 rvl = fac1*exp(-beta(ii)*al(ii)*log(rv(i)))
464 cj(i) = cj(i) + mu2*rvl
465 ai(i,1:3) = ai(i,1:3) + mu2*lam_al(1:3)
466 bi(i,1:3) = bi(i,1:3) + fac*mu2*lam_al(1:3)
467 ENDDO
468 ENDIF
469 ENDDO
470C------max of Cii
471 DO i=1,nel
472 d11 = ai(i,1)-bi(i,1)+cj(i)
473 d22 = ai(i,2)-bi(i,2)+cj(i)
474 d33 = ai(i,3)-bi(i,3)+cj(i)
475C
476 fac = one/rv(i)
477 cmax = fac*max(d11,d22,d33)
478 cimax(i) = two_third*gvis + max(cmax,cmax0)
479 et(i)= one
480 ENDDO
481C TRANSFORM PRINCIPAL PIOLA KIRCHOFF STRESSES TO GLOBAL DIRECTIONS
482 DO i=1,nel
483 signxx(i) = dirprv(i,1,1)*dirprv(i,1,1)*s(i,1)
484 . + dirprv(i,1,2)*dirprv(i,1,2)*s(i,2)
485 . + dirprv(i,1,3)*dirprv(i,1,3)*s(i,3)
486
487 signyy(i) = dirprv(i,2,2)*dirprv(i,2,2)*s(i,2)
488 . + dirprv(i,2,3)*dirprv(i,2,3)*s(i,3)
489 . + dirprv(i,2,1)*dirprv(i,2,1)*s(i,1)
490
491 signzz(i) = dirprv(i,3,3)*dirprv(i,3,3)*s(i,3)
492 . + dirprv(i,3,1)*dirprv(i,3,1)*s(i,1)
493 . + dirprv(i,3,2)*dirprv(i,3,2)*s(i,2)
494
495 signxy(i) = dirprv(i,1,1)*dirprv(i,2,1)*s(i,1)
496 . + dirprv(i,1,2)*dirprv(i,2,2)*s(i,2)
497 . + dirprv(i,1,3)*dirprv(i,2,3)*s(i,3)
498
499 signyz(i) = dirprv(i,2,2)*dirprv(i,3,2)*s(i,2)
500 . + dirprv(i,2,3)*dirprv(i,3,3)*s(i,3)
501 . + dirprv(i,2,1)*dirprv(i,3,1)*s(i,1)
502
503 signzx(i) = dirprv(i,3,3)*dirprv(i,1,3)*s(i,3)
504 . + dirprv(i,3,1)*dirprv(i,1,1)*s(i,1)
505 . + dirprv(i,3,2)*dirprv(i,1,2)*s(i,2)
506C
507C calcul des contraintes visqueuse ----
508C
509* SET SOUND SPEED
510 soundsp(i)=sqrt(cimax(i)/rho(i))
511c SOUNDSP(I)=SQRT((2*GMAX/3.+ RBULK)/RHO(I))
512* SET VISCOSITY
513 viscmax(i) = zero
514 ENDDO
515 RETURN
516 END
517
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
sigeps62
subroutine sigeps62(nel, nuparam, nuvar, nfunc, ifunc, npf, tf, time, timestep, uparam, rho0, rho, volume, eint, 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 sigeps62.F:45