OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
elas24.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!|| elas24 ../engine/source/materials/mat/mat024/elas24.F
25!||--- called by ------------------------------------------------------
26!|| conc24 ../engine/source/materials/mat/mat024/conc24.F
27!||--- calls -----------------------------------------------------
28!|| rdam24 ../engine/source/materials/mat/mat024/rdam24.F
29!||====================================================================
30 SUBROUTINE elas24(NEL ,PM ,
31 . CDAM ,SIGC ,DAM ,ANG ,EPS_F ,CRAK ,
32 . C44 ,C55 ,C66 ,DE1 ,DE2 ,DE3 ,
33 . DEPS1 ,DEPS2 ,DEPS3 ,DEPS4 ,DEPS5 ,DEPS6 ,
34 . S01 ,S02 ,S03 ,S04 ,S05 ,S06 ,
35 . SCAL1 ,SCAL2 ,SCAL3 )
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
223 END
elas24
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)
Definition elas24.F:36
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