OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fail_beam3.F File Reference
#include "implicit_f.inc"
#include "comlock.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine fail_beam3 (elbuf_str, mat_param, fail, snpc, stf, nel, jthe, dpla, tempel, ngl, 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)

Function/Subroutine Documentation

◆ fail_beam3()

subroutine fail_beam3 ( type (elbuf_struct_), intent(inout) elbuf_str,
type (matparam_struct_), intent(in) mat_param,
type (fail_param_), intent(in) fail,
integer, intent(in) snpc,
integer, intent(in) stf,
integer, intent(in) nel,
integer, intent(in) jthe,
intent(in) dpla,
intent(in) tempel,
integer, dimension(nel), intent(in) ngl,
intent(inout) off,
intent(in) epsd,
integer, dimension(snpc), intent(in) npf,
intent(in) tf,
intent(in) time,
integer, intent(in) iout,
integer, intent(in) istdo,
intent(in) svm,
intent(in) pressure,
intent(in) area,
intent(in) al,
intent(inout) f1,
intent(inout) f2,
intent(inout) f3,
intent(inout) m1,
intent(inout) m2,
intent(inout) m3,
integer, intent(in) ismstr,
intent(in) epsxx,
intent(in) epsxy,
intent(in) epsxz,
intent(in) kxx,
intent(in) kyy,
intent(in) kzz,
intent(in) dtime,
integer, intent(in) ntable,
type(ttable), dimension(ntable), intent(inout) table,
intent(in) pla,
intent(in) sigy )

Definition at line 44 of file fail_beam3.F.

54C-----------------------------------------------
55C M o d u l e s
56C-----------------------------------------------
57 USE mat_elem_mod
58 USE elbufdef_mod
59 USE fail_visual_b_mod
60 USE fail_gene1_b_mod
61 USE fail_inievo_b_mod
62 USE fail_tab2_b_mod
63C-----------------------------------------------
64C I m p l i c i t T y p e s
65C-----------------------------------------------
66#include "implicit_f.inc"
67C-----------------------------------------------
68C C o m m o n B l o c k s
69C-----------------------------------------------
70#include "comlock.inc"
71C-----------------------------------------------
72C D u m m y A r g u m e n t s
73C-----------------------------------------------
74 INTEGER ,INTENT(IN) :: NEL ! size of element group
75 INTEGER ,INTENT(IN) :: JTHE ! thermal dependency flag
76 INTEGER ,INTENT(IN) :: SNPC
77 INTEGER ,INTENT(IN) :: STF
78 INTEGER ,INTENT(IN) :: IOUT ! output file unit
79 INTEGER ,INTENT(IN) :: ISTDO ! output file unit
80 INTEGER ,INTENT(IN) :: ISMSTR
81 INTEGER ,DIMENSION(SNPC) ,INTENT(IN) :: NPF
82 INTEGER ,DIMENSION(NEL) ,INTENT(IN) :: NGL ! table of element identifiers
83 my_real ,INTENT(IN) :: time
84 my_real ,INTENT(IN) :: area
85 my_real ,DIMENSION(NEL) ,INTENT(IN) :: dpla
86 my_real ,DIMENSION(NEL) ,INTENT(IN) :: al
87 my_real ,DIMENSION(NEL) ,INTENT(IN) :: svm
88 my_real ,DIMENSION(NEL) ,INTENT(IN) :: pressure
89 my_real ,DIMENSION(NEL) ,INTENT(IN) :: epsd
90 my_real ,DIMENSION(NEL) ,INTENT(IN) :: epsxx,epsxy,epsxz
91 my_real ,DIMENSION(NEL) ,INTENT(IN) :: kxx,kyy,kzz
92 my_real ,DIMENSION(NEL) ,INTENT(IN) :: tempel
93 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: f1,f2,f3
94 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: m1,m2,m3
95
96 my_real ,DIMENSION(STF) ,INTENT(IN) :: tf
97 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: off
98 my_real ,DIMENSION(NEL) ,INTENT(IN) :: pla ! the plastic strain of the current element
99 my_real ,INTENT(IN) :: dtime
100 INTEGER ,INTENT(IN) :: NTABLE
101 TYPE(TTABLE), DIMENSION(NTABLE), INTENT(INOUT) :: TABLE ! TABLE DATA
102 my_real ,DIMENSION(NEL) ,INTENT(IN) :: sigy
103C
104 TYPE (ELBUF_STRUCT_) ,INTENT(INOUT) :: ELBUF_STR
105 TYPE (MATPARAM_STRUCT_) ,INTENT(IN) :: MAT_PARAM
106 TYPE (FAIL_PARAM_) ,INTENT(IN) :: FAIL
107C-----------------------------------------------
108C L o c a l V a r i a b l e s
109C-----------------------------------------------
110
111 INTEGER :: I,IFL,NFUNC,NPARAM,IRUPT
112 my_real :: t0, tm
113 my_real ,DIMENSION(NEL) :: tstar
114
115C
116
117C=======================================================================
118c to avoid compilation error with unused arguments
119c they will be necessary for next development step
120c-----------------------------------------------------
121 ifl = 1 ! only one failure model for beams
122C--------------------------------------
123 nparam = fail%NUPARAM
124 nfunc = fail%NFUNC
125c-----------------------------------------
126 irupt = elbuf_str%GBUF%FAIL(1)%ILAWF
127c-------------------------------------
128c progressive element erosion
129c
130 DO i = 1,nel
131
132 IF (off(i) < em01) off(i) = zero
133 IF (off(i) < one ) off(i) = off(i)*four_over_5
134 ENDDO
135
136c-----------------------------------------
137 SELECT CASE (irupt)
138c------------------------------------
139 CASE (1) ! Johnson-Cook
140 ! Tstar computation for Jhonson-Cook failure : T* = (T-T0)/(TM-T0)
141 IF (jthe > 0 .or. elbuf_str%gbuf%g_temp > 0) THEN
142 t0 = mat_param%THERM%TREF
143 tm = mat_param%THERM%TMELT
144 tstar(1:nel) = max(zero,(tempel(1:nel)-t0)) / max((tm-t0),em20)
145 ELSE
146 tstar(1:nel) = zero
147 ENDIF
148
149 CALL fail_johnson_b(nel ,ngl ,nparam ,fail%UPARAM,
150 . time ,tstar ,svm ,pressure ,
151 . dpla ,epsd ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX,
152 . elbuf_str%GBUF%FAIL(1)%TDEL ,iout ,istdo )
153c---------------
154 CASE (10) ! Tension Strain failure model
155 IF (jthe > 0 .or. elbuf_str%gbuf%g_temp > 0) THEN
156 t0 = mat_param%THERM%TREF
157 tm = mat_param%THERM%TMELT
158 tstar(1:nel) = max(zero,(tempel(1:nel)-t0)) / max((tm-t0),em20)
159 ELSE
160 tstar(1:nel) = zero
161 ENDIF
162
163 CALL fail_tensstrain_b(
164 . nel ,ngl ,nparam ,fail%UPARAM,
165 . time ,epsd ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX,
166 . elbuf_str%GBUF%FAIL(1)%TDEL,iout ,istdo ,fail%IFUNC ,
167 . epsxx ,al ,tstar ,
168 . snpc ,npf ,stf ,
169 . tf ,ismstr ,elbuf_str%GBUF%FAIL(1)%NVAR, elbuf_str%GBUF%FAIL(1)%VAR,
170 . elbuf_str%GBUF%DMGSCL)
171
172c---------------
173 CASE (11) ! Energy failure model
174 CALL fail_energy_b(nel ,ngl ,nparam ,fail%UPARAM,
175 . time ,epsd ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX,
176 . elbuf_str%GBUF%FAIL(1)%TDEL ,iout ,istdo ,fail%IFUNC ,
177 . area ,f1 ,f2 ,f3 ,
178 . m1 ,m2 ,m3 ,snpc ,npf ,stf ,
179 . tf ,epsxx ,epsxy ,epsxz ,kxx ,kyy ,
180 . kzz ,elbuf_str%GBUF%FAIL(1)%NVAR,elbuf_str%GBUF%FAIL(1)%VAR,
181 . elbuf_str%GBUF%DMGSCL)
182c---------------
183
184 CASE (30) ! BIQUAD
185 CALL fail_biquad_b(nel ,ngl ,nparam ,fail%UPARAM,
186 . time ,svm ,pressure ,
187 . dpla ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX,
188 . elbuf_str%GBUF%FAIL(1)%TDEL ,iout ,istdo ,fail%IFUNC ,
189 . nfunc,snpc ,npf ,stf ,tf,
190 . elbuf_str%GBUF%FAIL(1)%NVAR, elbuf_str%GBUF%FAIL(1)%VAR,al,
191 . elbuf_str%GBUF%DMGSCL)
192c---------------
193 CASE (36) ! VISUAL
194 CALL fail_visual_b(nel ,ngl ,nparam ,fail%UPARAM,
195 . time ,elbuf_str%GBUF%FAIL(1)%DAMMX,
196 . iout ,istdo ,
197 . area ,f1 ,
198 . epsxx ,epsxy ,epsxz ,
199 . elbuf_str%GBUF%FAIL(1)%NVAR,
200 . elbuf_str%GBUF%FAIL(1)%VAR,ismstr, dtime)
201
202 CASE (39) ! Gene1 failure model
203 CALL fail_gene1_b(
204 . nel ,fail%NUPARAM ,elbuf_str%GBUF%FAIL(1)%NVAR ,fail%NFUNC,
205 . fail%IFUNC ,npf ,tf ,
206 . time ,dtime ,fail%UPARAM ,
207 . ngl ,elbuf_str%GBUF%DT ,epsd ,elbuf_str%GBUF%FAIL(1)%VAR ,
208 . off ,
209 . area ,f1 ,
210 . epsxx ,epsxy ,epsxz ,
211 . tempel ,
212 . elbuf_str%GBUF%FAIL(1)%DAMMX ,al ,table,
213 . fail%NTABLE,
214 . fail%TABLE ,elbuf_str%GBUF%FAIL(1)%LF_DAMMX ,fail%NIPARAM ,fail%IPARAM,
215 . snpc ,stf ,ntable)
216
217 CASE (42) ! INIEVO failure model
218 CALL fail_inievo_b(
219 . nel ,fail%NUPARAM ,elbuf_str%GBUF%FAIL(1)%NVAR ,
220 . table ,fail%NTABLE ,fail%TABLE ,time ,fail%UPARAM ,
221 . ngl ,al ,dpla ,epsd , elbuf_str%GBUF%FAIL(1)%VAR,
222 . f1 ,area ,
223 . pla ,sigy ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX ,
224 . elbuf_str%GBUF%FAIL(1)%TDEL ,
225 . elbuf_str%GBUF%FAIL(1)%DAMINI,
226 . ntable,elbuf_str%GBUF%DMGSCL)
227
228 CASE (41) ! TAB2 failure model
229 CALL fail_tab2_b(
230 . nel ,fail%NUPARAM ,elbuf_str%GBUF%FAIL(1)%NVAR ,fail%NFUNC ,fail%IFUNC ,
231 . npf ,table ,tf ,time ,fail%UPARAM ,
232 . ngl ,al ,dpla ,epsd ,elbuf_str%GBUF%FAIL(1)%VAR,
233 . f1 ,area ,
234 . tempel ,off ,elbuf_str%GBUF%FAIL(1)%DAMMX ,elbuf_str%GBUF%FAIL(1)%TDEL ,
235 . fail%NTABLE ,fail%TABLE,
236 . snpc ,stf ,ntable ,elbuf_str%GBUF%DMGSCL)
237c-------------
238 END SELECT
239c-----------
240 DO i= 1,nel
241 IF (off(i) == four_over_5) THEN
242#include "lockon.inc"
243 WRITE(iout, 1000) ngl(i),time
244 WRITE(istdo,1000) ngl(i),time
245#include "lockoff.inc"
246 END IF
247 END DO
248c------------------
249 1000 FORMAT(1x,'DELETED BEAM ELEMENT in 1',i10,1x,'AT TIME :',1pe12.4)
250
251
252 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
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