OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fail_beam3.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!|| fail_beam3 ../engine/source/elements/beam/fail_beam3.F
25!||--- called by ------------------------------------------------------
26!|| main_beam3 ../engine/source/elements/beam/main_beam3.F
27!||--- calls -----------------------------------------------------
28!|| fail_biquad_b ../engine/source/materials/fail/biquad/fail_biquad_b.F
29!|| fail_energy_b ../engine/source/materials/fail/energy/fail_energy_b.F
30!|| fail_gene1_b ../engine/source/materials/fail/gene1/fail_gene1_b.F90
31!|| fail_inievo_b ../engine/source/materials/fail/inievo/fail_inievo_b.F90
32!|| fail_johnson_b ../engine/source/materials/fail/johnson_cook/fail_johnson_b.F
33!|| fail_tab2_b ../engine/source/materials/fail/tabulated/fail_tab2_b.F90
34!|| fail_tensstrain_b ../engine/source/materials/fail/tensstrain/fail_tensstrain_b.F
35!|| fail_visual_b ../engine/source/materials/fail/visual/fail_visual_b.F90
36!||--- uses -----------------------------------------------------
37!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
38!|| fail_gene1_b_mod ../engine/source/materials/fail/gene1/fail_gene1_b.F90
39!|| fail_inievo_b_mod ../engine/source/materials/fail/inievo/fail_inievo_b.F90
40!|| fail_tab2_b_mod ../engine/source/materials/fail/tabulated/fail_tab2_b.F90
41!|| fail_visual_b_mod ../engine/source/materials/fail/visual/fail_visual_b.F90
42!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
43!||====================================================================
44 SUBROUTINE fail_beam3(ELBUF_STR ,FAIL ,NUMMAT ,
45 . NPROPM ,SNPC ,STF ,
46 . NEL ,IMAT ,JTHE ,DPLA ,
47 . TEMPEL ,NGL ,PM ,
48 . OFF ,EPSD ,NPF ,TF ,
49 . TIME ,IOUT ,ISTDO ,
50 . SVM ,PRESSURE,AREA ,AL ,
51 . F1 ,F2 ,F3 ,M1 ,M2 ,
52 . M3 ,ISMSTR ,EPSXX ,EPSXY ,EPSXZ ,
53 . KXX ,KYY ,KZZ ,DTIME ,
54 . NTABLE ,TABLE ,PLA ,SIGY )
55C-----------------------------------------------
56C M o d u l e s
57C-----------------------------------------------
58 USE mat_elem_mod
59 USE elbufdef_mod
60 USE fail_visual_b_mod
61 USE fail_gene1_b_mod
62 USE fail_inievo_b_mod
63 USE fail_tab2_b_mod
64C-----------------------------------------------
65C I m p l i c i t T y p e s
66C-----------------------------------------------
67#include "implicit_f.inc"
68C-----------------------------------------------
69C C o m m o n B l o c k s
70C-----------------------------------------------
71#include "comlock.inc"
72C-----------------------------------------------
73C D u m m y A r g u m e n t s
74C-----------------------------------------------
75 INTEGER ,INTENT(IN) :: NEL ! size of element group
76 INTEGER ,INTENT(IN) :: IMAT ! material law number
77 INTEGER ,INTENT(IN) :: JTHE ! thermal dependency flag
78 INTEGER ,INTENT(IN) :: NUMMAT ! number of defined materials
79 INTEGER ,INTENT(IN) :: NPROPM ! size of real material parameter table
80 INTEGER ,INTENT(IN) :: SNPC
81 INTEGER ,INTENT(IN) :: STF
82 INTEGER ,INTENT(IN) :: IOUT ! output file unit
83 INTEGER ,INTENT(IN) :: ISTDO ! output file unit
84 INTEGER ,INTENT(IN) :: ISMSTR
85 INTEGER ,DIMENSION(SNPC) ,INTENT(IN) :: NPF
86 INTEGER ,DIMENSION(NEL) ,INTENT(IN) :: NGL ! table of element identifiers
87 my_real ,INTENT(IN) :: time
88 my_real ,INTENT(IN) :: area
89 my_real ,DIMENSION(NPROPM ,NUMMAT) ,INTENT(IN) :: pm
90 my_real ,DIMENSION(NEL) ,INTENT(IN) :: dpla
91 my_real ,DIMENSION(NEL) ,INTENT(IN) :: al
92 my_real ,DIMENSION(NEL) ,INTENT(IN) :: svm
93 my_real ,DIMENSION(NEL) ,INTENT(IN) :: pressure
94 my_real ,DIMENSION(NEL) ,INTENT(IN) :: epsd
95 my_real ,DIMENSION(NEL) ,INTENT(IN) :: epsxx,epsxy,epsxz
96 my_real ,DIMENSION(NEL) ,INTENT(IN) :: kxx,kyy,kzz
97 my_real ,DIMENSION(NEL) ,INTENT(IN) :: tempel
98 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: f1,f2,f3
99 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: m1,m2,m3
100
101 my_real ,DIMENSION(STF) ,INTENT(IN) :: tf
102 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: off
103 my_real ,DIMENSION(NEL) ,INTENT(IN) :: pla ! the plastic strain of the current element
104 my_real ,INTENT(IN) :: dtime
105 INTEGER ,INTENT(IN) :: NTABLE
106 TYPE(ttable), DIMENSION(NTABLE), INTENT(INOUT) :: TABLE ! TABLE DATA
107 my_real ,DIMENSION(NEL) ,INTENT(IN) :: sigy
108C
109 TYPE (ELBUF_STRUCT_), INTENT(INOUT) :: ELBUF_STR
110 TYPE (FAIL_PARAM_) ,INTENT(IN) :: FAIL
111C-----------------------------------------------
112C L o c a l V a r i a b l e s
113C-----------------------------------------------
114 INTEGER :: IFUNC(100)
115 INTEGER :: I,IFL,NFUNC,NPARAM,IRUPT
116 my_real :: t0,tm,g_module
117 my_real ,DIMENSION(NEL) :: tstar
118 my_real :: bidon
119C
120 my_real, DIMENSION(:) ,POINTER :: uparamf
121C=======================================================================
122c to avoid compilation error with unused arguments
123c they will be necessary for next development step
124c-----------------------------------------------------
125 ifl = 1 ! only one failure model for beams
126C--------------------------------------
127 nparam = fail%NUPARAM
128 nfunc = fail%NFUNC
129c-----------------------------------------
130 irupt = elbuf_str%GBUF%FAIL(1)%ILAWF
131c-------------------------------------
132c progressive element erosion
133c
134 DO i = 1,nel
135
136 IF (off(i) < em01) off(i) = zero
137 IF (off(i) < one ) off(i) = off(i)*four_over_5
138 ENDDO
139
140c-----------------------------------------
141 SELECT CASE (irupt)
142c------------------------------------
143 CASE (1) ! Johnson-Cook
144 ! Tstar computation for Jhonson-Cook failure : T* = (T-T0)/(TM-T0)
145 IF (jthe > 0) THEN
146 t0 = pm(79, imat)
147 tm = pm(80, imat)
148 DO i=1,nel
149 tstar(i) = max(zero,(tempel(i)-t0)/(tm-t0))
150 ENDDO
151 ELSE
152 tstar(1:nel) = zero
153 ENDIF
154 CALL fail_johnson_b(nel ,ngl ,nparam ,fail%UPARAM,
155 . time ,tstar ,svm ,pressure ,
156 . dpla ,epsd ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX,
157 . elbuf_str%GBUF%FAIL(1)%TDEL ,iout ,istdo )
158c---------------
159 CASE (10) ! Tension Strain failure model
160 IF (jthe > 0) THEN
161 t0 = pm(79, imat)
162 tm = pm(80, imat)
163 DO i=1,nel
164 tstar(i) = max(zero,(tempel(i)-t0)/(tm-t0))
165 ENDDO
166 ELSE
167 tstar(1:nel) = zero
168 ENDIF
169
170 CALL fail_tensstrain_b(
171 . nel ,ngl ,nparam ,fail%UPARAM,
172 . time ,epsd ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX,
173 . elbuf_str%GBUF%FAIL(1)%TDEL,iout ,istdo ,fail%IFUNC ,
174 . epsxx ,al ,tstar ,
175 . snpc ,npf ,stf ,
176 . tf ,ismstr ,elbuf_str%GBUF%FAIL(1)%NVAR, elbuf_str%GBUF%FAIL(1)%VAR,
177 . elbuf_str%GBUF%DMGSCL)
178
179c---------------
180 CASE (11) ! Energy failure model
181 CALL fail_energy_b(nel ,ngl ,nparam ,fail%UPARAM,
182 . time ,epsd ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX,
183 . elbuf_str%GBUF%FAIL(1)%TDEL ,iout ,istdo ,fail%IFUNC ,
184 . area ,f1 ,f2 ,f3 ,
185 . m1 ,m2 ,m3 ,snpc ,npf ,stf ,
186 . tf ,epsxx ,epsxy ,epsxz ,kxx ,kyy ,
187 . kzz ,elbuf_str%GBUF%FAIL(1)%NVAR,elbuf_str%GBUF%FAIL(1)%VAR,
188 . elbuf_str%GBUF%DMGSCL)
189c---------------
190
191 CASE (30) ! BIQUAD
192 CALL fail_biquad_b(nel ,ngl ,nparam ,fail%UPARAM,
193 . time ,svm ,pressure ,
194 . dpla ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX,
195 . elbuf_str%GBUF%FAIL(1)%TDEL ,iout ,istdo ,fail%IFUNC ,
196 . nfunc,snpc ,npf ,stf ,tf,
197 . elbuf_str%GBUF%FAIL(1)%NVAR, elbuf_str%GBUF%FAIL(1)%VAR,al,
198 . elbuf_str%GBUF%DMGSCL)
199c---------------
200 CASE (36) ! VISUAL
201 CALL fail_visual_b(nel ,ngl ,nparam ,fail%UPARAM,
202 . time ,elbuf_str%GBUF%FAIL(1)%DAMMX,
203 . iout ,istdo ,
204 . area ,f1 ,
205 . epsxx ,epsxy ,epsxz ,
206 . elbuf_str%GBUF%FAIL(1)%NVAR,
207 . elbuf_str%GBUF%FAIL(1)%VAR,ismstr, dtime)
208
209 CASE (39) ! Gene1 failure model
210 CALL fail_gene1_b(
211 . nel ,fail%NUPARAM ,elbuf_str%GBUF%FAIL(1)%NVAR ,fail%NFUNC,
212 . fail%IFUNC ,npf ,tf ,
213 . time ,dtime ,fail%UPARAM ,
214 . ngl ,elbuf_str%GBUF%DT ,epsd ,elbuf_str%GBUF%FAIL(1)%VAR ,
215 . off ,
216 . area ,f1 ,
217 . epsxx ,epsxy ,epsxz ,
218 . tempel ,
219 . elbuf_str%GBUF%FAIL(1)%DAMMX ,al ,table,
220 . fail%NTABLE,
221 . fail%TABLE ,elbuf_str%GBUF%FAIL(1)%LF_DAMMX ,fail%NIPARAM ,fail%IPARAM,
222 . snpc ,stf ,ntable)
223
224 CASE (42) ! INIEVO failure model
225 CALL fail_inievo_b(
226 . nel ,fail%NUPARAM ,elbuf_str%GBUF%FAIL(1)%NVAR ,
227 . table ,fail%NTABLE ,fail%TABLE ,time ,fail%UPARAM ,
228 . ngl ,al ,dpla ,epsd , elbuf_str%GBUF%FAIL(1)%VAR,
229 . f1 ,area ,
230 . pla ,sigy ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX ,
231 . elbuf_str%GBUF%FAIL(1)%TDEL ,
232 . elbuf_str%GBUF%FAIL(1)%DAMINI,
233 . ntable,elbuf_str%GBUF%DMGSCL)
234
235 CASE (41) ! TAB2 failure model
236 CALL fail_tab2_b(
237 . nel ,fail%NUPARAM ,elbuf_str%GBUF%FAIL(1)%NVAR ,fail%NFUNC ,fail%IFUNC ,
238 . npf ,table ,tf ,time ,fail%UPARAM ,
239 . ngl ,al ,dpla ,epsd ,elbuf_str%GBUF%FAIL(1)%VAR,
240 . f1 ,area ,
241 . tempel ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX ,elbuf_str%GBUF%FAIL(1)%TDEL ,
242 . fail%NTABLE ,fail%TABLE,
243 . snpc ,stf ,ntable ,elbuf_str%GBUF%DMGSCL)
244c-------------
245 END SELECT
246c-----------
247 DO i= 1,nel
248 IF (off(i) == four_over_5) THEN
249#include "lockon.inc"
250 WRITE(iout, 1000) ngl(i),time
251 WRITE(istdo,1000) ngl(i),time
252#include "lockoff.inc"
253 END IF
254 END DO
255c------------------
256 1000 FORMAT(1x,'DELETED BEAM ELEMENT in 1',i10,1x,'AT TIME :',1pe12.4)
257
258
259 RETURN
260 END SUBROUTINE fail_beam3
my_real
#define my_real
Definition cppsort.cpp:32
fail_beam3
subroutine fail_beam3(elbuf_str, fail, nummat, npropm, snpc, stf, nel, imat, jthe, dpla, tempel, ngl, pm, off, epsd, npf, tf, time, iout, istdo, svm, pressure, area, al, f1, f2, f3, m1, m2, m3, ismstr, epsxx, epsxy, epsxz, kxx, kyy, kzz, dtime, ntable, table, pla, sigy)
Definition fail_beam3.F:55
fail_biquad_b
subroutine fail_biquad_b(nel, ngl, nuparam, uparam, time, svm, pressure, dpla, off, dfmax, tdel, iout, istdo, ifunc, nfunc, snpc, npf, stf, tf, nuvar, uvar, al, dmgscl)
Definition fail_biquad_b.F:38
fail_energy_b
subroutine fail_energy_b(nel, ngl, nuparam, uparam, time, epsd, off, dfmax, tdel, iout, istdo, ifunc, area, f1, f2, f3, m1, m2, m3, snpc, npf, stf, tf, epsxx, epsxy, epsxz, kxx, kyy, kzz, nvar, uvar, dmgscl)
Definition fail_energy_b.F:40
fail_johnson_b
subroutine fail_johnson_b(nel, ngl, nuparam, uparam, time, tstar, svm, pressure, dpla, epsd, off, dfmax, tdel, iout, istdo)
Definition fail_johnson_b.F:33
fail_tensstrain_b
subroutine fail_tensstrain_b(nel, ngl, nuparam, uparam, time, epsd, off, dfmax, tdel, iout, istdo, ifunc, epsxx, al, tstar, snpc, npf, stf, tf, ismstr, nuvar, uvar, dmgscl)
Definition fail_tensstrain_b.F:38
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
max
#define max(a, b)
Definition macros.h:21
bidon
subroutine bidon
Definition machine.F:41