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

Go to the source code of this file.

Functions/Subroutines

subroutine sigeps25cp (mat_param, jft, jlt, off, dir, shf, npt, ngl, ipt, off_old, thk0, epsd, sigy, zcfac, nel, depsxx, depsyy, depsxy, depsyz, depszx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, wplar, ioff_duct, pla, israte, asrate, epsd_pg, tsaiwu)

Function/Subroutine Documentation

◆ sigeps25cp()

subroutine sigeps25cp ( type(matparam_struct_), intent(in) mat_param,
integer jft,
integer jlt,
off,
dir,
shf,
integer npt,
integer, dimension(mvsiz) ngl,
integer ipt,
off_old,
thk0,
intent(inout) epsd,
sigy,
zcfac,
integer nel,
depsxx,
depsyy,
depsxy,
depsyz,
depszx,
sigoxx,
sigoyy,
sigoxy,
sigoyz,
sigozx,
signxx,
signyy,
signxy,
signyz,
signzx,
wplar,
integer, dimension(mvsiz) ioff_duct,
pla,
integer israte,
intent(in) asrate,
intent(in) epsd_pg,
tsaiwu )

Definition at line 32 of file sigeps25cp.F.

41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE matparam_def_mod
45C-----------------------------------------------
46C I m p l i c i t T y p e s
47C-----------------------------------------------
48#include "implicit_f.inc"
49#include "comlock.inc"
50C-----------------------------------------------
51C G l o b a l P a r a m e t e r s
52C-----------------------------------------------
53#include "mvsiz_p.inc"
54C-----------------------------------------------
55C C o m m o n B l o c k s
56C-----------------------------------------------
57#include "param_c.inc"
58C-----------------------------------------------
59C D u m m y A r g u m e n t s
60C-----------------------------------------------
61 INTEGER JFT,JLT,NPT,NEL,IPT
62 INTEGER NGL(MVSIZ),IOFF_DUCT(MVSIZ),ISRATE
63 my_real ,INTENT(IN) :: asrate
64 my_real ,DIMENSION(NEL) ,INTENT(IN) :: epsd_pg
65 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: epsd
66 my_real
67 . thk0(mvsiz),
68 . off(*),dir(*),shf(*),sigy(*),zcfac(mvsiz,2),
69 . depsxx(mvsiz),depsyy(mvsiz),depsxy(mvsiz),depsyz(mvsiz),
70 . depszx(mvsiz),sigoxx(nel),sigoyy(nel),sigoxy(nel),
71 . sigoyz(nel),sigozx(nel),signxx(nel),signyy(nel),signxy(nel),
72 . signyz(nel),signzx(nel),wplar(mvsiz),off_old(mvsiz),pla(nel),
73 . tsaiwu(nel)
74 TYPE(MATPARAM_STRUCT_) ,INTENT(IN) :: MAT_PARAM
75C-----------------------------------------------
76C L o c a l V a r i a b l e s
77C-----------------------------------------------
78 INTEGER ICC(MVSIZ),I,ILAYER,ICC_1
79 my_real
80 . yld(mvsiz),etse(mvsiz),degmb(mvsiz),cb(mvsiz),cn(mvsiz),
81 . e11(mvsiz),e22(mvsiz),nu12(mvsiz),nu21(mvsiz),g12(mvsiz),
82 . g23(mvsiz),g31(mvsiz),de(mvsiz),f1(mvsiz),f2(mvsiz),
83 . f12(mvsiz),f11(mvsiz),f22(mvsiz),f33(mvsiz),sigt1(mvsiz),
84 . sigt2(mvsiz),fmax(mvsiz),wwpla(mvsiz),wplamx(mvsiz),
85 . cc(mvsiz),epdr(mvsiz),wplaref(mvsiz),
86 . e11_1,e22_1,nu12_1,nu21_1,g12_1,g23_1,g31_1,cb_1,cn_1,
87 . cc_1,epdr_1,f1_1,f2_1,f11_1,f22_1,f33_1,f12_1,sigt1_1,sigt2_1
88C=======================================================================
89 icc_1 = mat_param%IPARAM(3)
90!
91 e11_1 = mat_param%UPARAM(1) ! PM(33)
92 e22_1 = mat_param%UPARAM(2) ! PM(34)
93 nu12_1 = mat_param%UPARAM(4) ! PM(35)
94 nu21_1 = mat_param%UPARAM(5) ! PM(36)
95 g12_1 = mat_param%UPARAM(6) ! PM(37)
96 g23_1 = mat_param%UPARAM(7) ! PM(38)
97 g31_1 = mat_param%UPARAM(8) ! PM(39)
98 cc_1 = mat_param%UPARAM(10) ! PM(50)
99 epdr_1 = mat_param%UPARAM(11) ! PM(51)
100 sigt1_1 = mat_param%UPARAM(12) ! PM(60)
101 sigt2_1 = mat_param%UPARAM(13) ! PM(61)
102
103 cb_1 = mat_param%UPARAM(27) ! PM(46)
104 cn_1 = mat_param%UPARAM(28) ! PM(47)
105
106 f1_1 = mat_param%UPARAM(21) ! PM(54)
107 f2_1 = mat_param%UPARAM(22) ! PM(55)
108 f11_1 = mat_param%UPARAM(23) ! PM(56)
109 f22_1 = mat_param%UPARAM(24) ! PM(57)
110 f33_1 = mat_param%UPARAM(25) ! PM(58)
111 f12_1 = mat_param%UPARAM(26) ! PM(59)
112!
113 DO i=jft,jlt
114 e11(i) =e11_1
115 e22(i) =e22_1
116 nu12(i) =nu12_1
117 nu21(i) =nu21_1
118 g12(i) =g12_1
119 g23(i) =g23_1
120 g31(i) =g31_1
121 de(i) = one
122 cb(i) =cb_1
123 cn(i) =cn_1
124 cc(i) =cc_1
125 epdr(i) =epdr_1
126 icc(i) =icc_1
127 f1(i) =f1_1
128 f2(i) =f2_1
129 f11(i) =f11_1
130 f22(i) =f22_1
131 f33(i) =f33_1
132 f12(i) =f12_1
133 sigt1(i) =sigt1_1
134 sigt2(i) =sigt2_1
135 yld(i) =one/sqrt(f33_1)
136 ENDDO
137!
138 DO i=jft,jlt
139 wplaref(i)= mat_param%UPARAM(9) ! PM(68)
140 wplamx(i) = mat_param%UPARAM(20) ! PM(41)
141 fmax(i) = mat_param%UPARAM(29) ! PM(49)
142 ENDDO
143!
144! CONTRAINTES PLASTIQUEMENT ADMISSIBLES
145!
146 CALL m25cplrp2(jft ,jlt ,wplaref,thk0 ,off ,etse ,
147 2 pla ,dir ,npt ,cc ,epdr ,icc ,
148 3 wwpla ,shf ,fmax ,cb ,cn ,nel ,
149 4 degmb ,f1 ,f2 ,f12 ,f11 ,f22 ,
150 5 f33 ,e11 ,e22 ,nu12 ,nu21 ,g12 ,
151 6 g23 ,g31 ,de ,epsd ,israte ,yld ,
152 7 depsxx ,depsyy ,depsxy ,depsyz ,depszx ,sigoxx ,
153 8 sigoyy ,sigoxy ,sigoyz ,sigozx ,signxx ,signyy ,
154 9 signxy ,signyz ,signzx ,tsaiwu ,asrate ,epsd_pg)
155!
156 DO i=jft,jlt
157 wplar(i) = max(wplar(i),wwpla(i))
158 END DO
159C-----------------------
160C For QEPH
161C-----------------------
162 DO i=jft,jlt
163 zcfac(i,1) = zcfac(i,1) + etse(i) / npt
164 zcfac(i,2) = min(etse(i),zcfac(i,2))
165 sigy(i) = sigy(i) + yld(i)/npt
166 ENDDO
167C----------------------------
168C TEST DE RUPTURE DUCTILE
169C---------------------------
170 DO i=jft,jlt
171 IF ( off(i) == off_old(i) .and. off_old(i) < em01 ) off(i)=zero
172 IF ( off(i) == off_old(i) .and. off_old(i) < one ) off(i)=off(i)*four_over_5
173 ENDDO
174!
175 DO i=jft,jlt
176 IF (off(i) == one .AND. wplar(i) >= wplamx(i)) THEN
177 off(i)=four_over_5
178 ioff_duct(i) = 1
179 ENDIF ! IF (OFF(I) == ONE .AND. WPLAR(I) >= WPLAMX(I))
180 ENDDO
181!
182 IF (npt == 1) THEN
183 DO i=jft,jlt
184 degmb(i)=degmb(i)*off(i)
185 ENDDO
186 ENDIF
187c-----------
188 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
m25cplrp2
subroutine m25cplrp2(jft, jlt, wplaref, thk, off, etse, wpla, dir, npt, cc, epdr, icc, wwpla, shf, fmax, cb, cn, nel, degmb, f1, f2, f12, f11, f22, f33, e11, e22, nu12, nu21, g12, g23, g31, de, epsd, israte, sigy, depsxx, depsyy, depsxy, depsyz, depszx, sigoxx, sigoyy, sigoxy, sigoyz, sigozx, signxx, signyy, signxy, signyz, signzx, tsaiwu, asrate, epsd_pg)
Definition m25cplrc2.F:37
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21