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

Go to the source code of this file.

Functions/Subroutines

subroutine elas24 (nel, pm, cdam, sigc, dam, ang, eps_f, crak, c44, c55, c66, de1, de2, de3, deps1, deps2, deps3, deps4, deps5, deps6, s01, s02, s03, s04, s05, s06, scal1, scal2, scal3)

Function/Subroutine Documentation

◆ elas24()

subroutine elas24 ( integer nel,
pm,
intent(inout) cdam,
intent(inout) sigc,
intent(inout) dam,
intent(in) ang,
intent(in) eps_f,
intent(inout) crak,
intent(out) c44,
intent(out) c55,
intent(out) c66,
intent(out) de1,
intent(out) de2,
intent(out) de3,
intent(inout) deps1,
intent(inout) deps2,
intent(inout) deps3,
intent(inout) deps4,
intent(inout) deps5,
intent(inout) deps6,
intent(out) s01,
intent(out) s02,
intent(out) s03,
intent(out) s04,
intent(out) s05,
intent(out) s06,
intent(out) scal1,
intent(out) scal2,
intent(out) scal3 )

Definition at line 30 of file elas24.F.

36C-----------------------------------------------
37C I m p l i c i t T y p e s
38C-----------------------------------------------
39#include "implicit_f.inc"
40C-----------------------------------------------
41C C o m m o n B l o c k s
42C-----------------------------------------------
43#include "param_c.inc"
44C-----------------------------------------------
45C D u m m y A r g u m e n t s
46C-----------------------------------------------
47 INTEGER NEL
48 my_real pm(npropm)
49 my_real, DIMENSION(NEL) ,INTENT(INOUT) :: deps1,deps2,deps3,deps4,deps5,deps6
50 my_real, DIMENSION(NEL) ,INTENT(OUT) :: s01,s02,s03,s04,s05,s06,
51 . scal1,scal2,scal3,c44,c55,c66
52 my_real, DIMENSION(NEL) ,INTENT(OUT) :: de1,de2,de3
53 my_real, DIMENSION(NEL,3) ,INTENT(IN) :: eps_f
54 my_real, DIMENSION(NEL,6) ,INTENT(IN) :: ang
55 my_real, DIMENSION(NEL,6) ,INTENT(INOUT):: sigc
56 my_real, DIMENSION(NEL,3) ,INTENT(INOUT):: dam,crak
57 my_real, DIMENSION(NEL,3,3) ,INTENT(INOUT):: cdam
58C-----------------------------------------------
59C L o c a l V a r i a b l e s
60C-----------------------------------------------
61 INTEGER I,K,NBDAMA,NBDAMA0
62 INTEGER, DIMENSION(NEL) :: DAMAI,DAMA0I
63 my_real, DIMENSION(NEL) :: de4,de5,de6
64 my_real young,nu,g,a11,a12,dsup,qq,den,depsf,hh1,hh2,hh3
65C=======================================================================
66c DEi = (1 - DAMi) coefficients d endommagement de la matrice Hook
67c-----------------------------------------------------------------------
68 young = pm(20)
69 nu = pm(21)
70 g = pm(22)
71 a11 = pm(24)
72 a12 = pm(25)
73 dsup = pm(26)
74 qq = pm(28)
75c
76c test of damaged elements
77c
78 nbdama = 0
79 nbdama0 = 0
80 DO i = 1,nel
81 IF (dam(i,1) + dam(i,2) + dam(i,3) > zero) THEN
82 nbdama = nbdama + 1
83 damai(nbdama)=i
84 ELSE
85 nbdama0 = nbdama0 + 1
86 dama0i(nbdama0)=i
87 ENDIF
88 ENDDO
89c
90c ROTATION de DEPS_IJ DANS LES DIRECTIONS D'ENDOMMAGEMENT
91c
92 IF (nbdama > 0) THEN
93 CALL rdam24(nel ,nbdama,damai ,ang ,deps1 ,
94 . deps2 ,deps3 ,deps4 ,deps5 ,deps6 )
95 ENDIF
96c-----------------------------------------------------------
97c DAMAGED ELEMENT GROUP => modified Hook matrix
98c-----------------------------------------------------------
99 DO k=1,nbdama
100 i = damai(k)
101c cumul de deformation endommagee
102 crak(i,1) = crak(i,1) + deps1(i)
103 crak(i,2) = crak(i,2) + deps2(i)
104 crak(i,3) = crak(i,3) + deps3(i)
105c
106c Damage evolution depending on EPS_F_i/EPS_i
107c
108 hh1 = max(crak(i,1),em20)
109 IF (eps_f(i,1)>zero) THEN
110 de1(i) = qq * (one- eps_f(i,1)/hh1 )
111 ELSE
112 de1(i) = zero
113 ENDIF
114 hh2 = max(crak(i,2),em20)
115 IF (eps_f(i,2)>zero) THEN
116 de2(i) = qq * (one- eps_f(i,2)/hh2 )
117 ELSE
118 de2(i) = zero
119 ENDIF
120 hh3 = max(crak(i,3),em20)
121 IF (eps_f(i,3)>zero) THEN
122 de3(i) = qq * (one - eps_f(i,3)/hh3 )
123 ELSE
124 de3(i) = zero
125 ENDIF
126c
127c DAM(I,J) SAUVEGARDE DE L'ENDOMMAGEMENT VRAI
128c
129 de1(i) = min(de1(i),dsup)
130 de2(i) = min(de2(i),dsup)
131 de3(i) = min(de3(i),dsup)
132 dam(i,1) = max(de1(i),dam(i,1))
133 dam(i,2) = max(de2(i),dam(i,2))
134 dam(i,3) = max(de3(i),dam(i,3))
135c
136 de1(i) = one - max( zero , sign(dam(i,1),crak(i,1)) )
137 de2(i) = one - max( zero , sign(dam(i,2),crak(i,2)) )
138 de3(i) = one - max( zero , sign(dam(i,3),crak(i,3)) )
139 scal1(i)= half + sign(half,de1(i)-one)
140 scal2(i)= half + sign(half,de2(i)-one)
141 scal3(i)= half + sign(half,de3(i)-one)
142 de4(i) = scal1(i)*scal2(i)
143 de5(i) = scal2(i)*scal3(i)
144 de6(i) = scal3(i)*scal1(i)
145c
146c NOUVELLE MATRICE ELASTIQUE DE HOOKE ENDOMMAGEE
147c
148 den = one - nu**2 *(de4(i) + de5(i) + de6(i)
149 . + two*nu*scal1(i)*scal2(i)*scal3(i))
150C
151 cdam(i,1,1) = young*de1(i)*(one - nu**2*de5(i))/den
152 cdam(i,2,2) = young*de2(i)*(one - nu**2*de6(i))/den
153 cdam(i,3,3) = young*de3(i)*(one - nu**2*de4(i))/den
154 cdam(i,1,2) = nu*young*de4(i) *(one + nu*scal3(i))/den
155 cdam(i,1,3) = nu*young*de6(i) *(one + nu*scal2(i))/den
156 cdam(i,2,3) = nu*young*de5(i) *(one + nu*scal1(i))/den
157 cdam(i,2,1) = cdam(i,1,2)
158 cdam(i,3,1) = cdam(i,1,3)
159 cdam(i,3,2) = cdam(i,2,3)
160c
161 c44(i) = g*de4(i)
162 c55(i) = g*de5(i)
163 c66(i) = g*de6(i)
164c
165c CONTRAINTES ELASTIQUES
166c
167 s01(i) = cdam(i,1,1)*crak(i,1) + cdam(i,1,2)*crak(i,2)
168 . + cdam(i,1,3)*crak(i,3)
169 s02(i) = cdam(i,2,1)*crak(i,1) + cdam(i,2,2)*crak(i,2)
170 . + cdam(i,2,3)*crak(i,3)
171 s03(i) = cdam(i,3,1)*crak(i,1) + cdam(i,3,2)*crak(i,2)
172 . + cdam(i,3,3)*crak(i,3)
173c formulation incrementale en cisaillement
174 s04(i) = de4(i)*sigc(i,4) + c44(i)*deps4(i)
175 s05(i) = de5(i)*sigc(i,5) + c55(i)*deps5(i)
176 s06(i) = de6(i)*sigc(i,6) + c66(i)*deps6(i)
177 ENDDO
178c------------------------------------------------------------------------
179c NON DAMAGED ELEMENT GROUP => initial Hook matrix
180c------------------------------------------------------------------------
181 DO k=1,nbdama0
182 i = dama0i(k)
183 cdam(i,1,1) = a11
184 cdam(i,2,2) = a11
185 cdam(i,3,3) = a11
186 cdam(i,1,2) = a12
187 cdam(i,2,1) = a12
188 cdam(i,1,3) = a12
189 cdam(i,3,1) = a12
190 cdam(i,2,3) = a12
191 cdam(i,3,2) = a12
192 c44(i) = g
193 c55(i) = g
194 c66(i) = g
195 ENDDO
196c--------------------
197 DO k=1,nbdama0
198 i = dama0i(k)
199 s01(i) = sigc(i,1) + cdam(i,1,1)*deps1(i)+cdam(i,1,2)*(deps2(i)+deps3(i))
200 s02(i) = sigc(i,2) + cdam(i,1,1)*deps2(i)+cdam(i,1,2)*(deps1(i)+deps3(i))
201 s03(i) = sigc(i,3) + cdam(i,1,1)*deps3(i)+cdam(i,1,2)*(deps1(i)+deps2(i))
202 s04(i) = sigc(i,4) + g*deps4(i)
203 s05(i) = sigc(i,5) + g*deps5(i)
204 s06(i) = sigc(i,6) + g*deps6(i)
205 crak(i,1) = s01(i) - nu *( s02(i) + s03(i) )
206 crak(i,2) = s02(i) - nu *( s01(i) + s03(i) )
207 crak(i,3) = s03(i) - nu *( s01(i) + s02(i) )
208 crak(i,1) = crak(i,1)/young
209 crak(i,2) = crak(i,2)/young
210 crak(i,3) = crak(i,3)/young
211 scal1(i) = one
212 scal2(i) = one
213 scal3(i) = one
214 de4(i) = one
215 de5(i) = one
216 de6(i) = one
217 de1(i) = one
218 de2(i) = one
219 de3(i) = one
220 ENDDO
221c-----------
222 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
rdam24
subroutine rdam24(nel, nbdama, damai, ang, deps1, deps2, deps3, deps4, deps5, deps6)
Definition rdam24.F:30