OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fail_beam18.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_beam18 ../engine/source/elements/beam/fail_beam18.F
25!||--- called by ------------------------------------------------------
26!|| main_beam18 ../engine/source/elements/beam/main_beam18.f
27!||--- calls -----------------------------------------------------
28!|| fail_biquad_ib ../engine/source/materials/fail/biquad/fail_biquad_ib.f
29!|| fail_cockroft_ib ../engine/source/materials/fail/cockroft_latham/fail_cockroft_ib.F
30!|| fail_energy_ib ../engine/source/materials/fail/energy/fail_energy_ib.F
31!|| fail_gene1_ib ../engine/source/materials/fail/gene1/fail_gene1_ib.F90
32!|| fail_inievo_ib ../engine/source/materials/fail/inievo/fail_inievo_ib.F90
33!|| fail_johnson_ib ../engine/source/materials/fail/johnson_cook/fail_johnson_ib.F
34!|| fail_tab2_ib ../engine/source/materials/fail/tabulated/fail_tab2_ib.F90
35!|| fail_tensstrain_ib ../engine/source/materials/fail/tensstrain/fail_tensstrain_ib.F
36!|| fail_visual_ib ../engine/source/materials/fail/visual/fail_visual_ib.F90
37!||--- uses -----------------------------------------------------
38!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
39!|| fail_gene1_ib_mod ../engine/source/materials/fail/gene1/fail_gene1_ib.F90
40!|| fail_inievo_ib_mod ../engine/source/materials/fail/inievo/fail_inievo_ib.F90
41!|| fail_tab2_ib_mod ../engine/source/materials/fail/tabulated/fail_tab2_ib.F90
42!|| fail_visual_ib_mod ../engine/source/materials/fail/visual/fail_visual_ib.F90
43!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
44!||====================================================================
45 SUBROUTINE fail_beam18(ELBUF_STR,FAIL ,NUMMAT ,NUMGEO ,
46 . NPROPM ,NPROPG ,SNPC ,STF ,
47 . NEL ,NPT ,IMAT ,IPROP ,JTHE ,
48 . TEMPEL ,NGL ,PM ,GEO ,
49 . OFF ,EPSD ,NPF ,TF ,
50 . DPLA ,EINT ,TIME ,IOUT ,ISTDO ,
51 . AL ,ISMSTR ,EXX ,EXY ,EXZ ,
52 . KXX ,KYY ,KZZ ,DTIME ,
53 . NTABLE ,TABLE ,SIGY )
54C-----------------------------------------------
55C M o d u l e s
56C-----------------------------------------------
57 USE mat_elem_mod
58 USE elbufdef_mod
59 USE fail_visual_ib_mod
60 USE fail_gene1_ib_mod
61 USE fail_inievo_ib_mod
62 USE fail_tab2_ib_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) :: IMAT ! material law number
76 INTEGER ,INTENT(IN) :: IPROP ! beam property number
77 INTEGER ,INTENT(IN) :: NPT ! number of integration points in beam section
78 INTEGER ,INTENT(IN) :: JTHE ! thermal dependency flag
79 INTEGER ,INTENT(IN) :: NUMMAT ! number of defined materials
80 INTEGER ,INTENT(IN) :: NUMGEO ! number of defined properties
81 INTEGER ,INTENT(IN) :: NPROPM ! size of real material parameter table
82 INTEGER ,INTENT(IN) :: NPROPG ! size of real property parameter table
83 INTEGER ,INTENT(IN) :: SNPC
84 INTEGER ,INTENT(IN) :: STF
85 INTEGER ,INTENT(IN) :: NTABLE
86 TYPE(ttable), DIMENSION(NTABLE), INTENT(INOUT) :: TABLE ! TABLE DATA
87 INTEGER ,INTENT(IN) :: IOUT ! output file unit
88 INTEGER ,INTENT(IN) :: ISTDO ! output file unit
89 INTEGER ,INTENT(IN) :: ISMSTR
90 INTEGER ,DIMENSION(SNPC) ,INTENT(IN) :: NPF
91 INTEGER ,DIMENSION(NEL) ,INTENT(IN) :: NGL ! table of element identifiers
92 my_real ,INTENT(IN) :: time
93 my_real ,DIMENSION(NPROPM ,NUMMAT) ,INTENT(IN) :: pm
94 my_real ,DIMENSION(NPROPG ,NUMGEO) ,INTENT(IN) :: geo
95 my_real ,DIMENSION(NEL) ,INTENT(IN) :: epsd
96 my_real ,DIMENSION(NEL) ,INTENT(IN) :: al
97 my_real ,DIMENSION(NEL) ,INTENT(IN) :: tempel
98 my_real ,DIMENSION(STF) ,INTENT(IN) :: tf
99 my_real ,DIMENSION(NEL) ,INTENT(IN) :: exx,exy,exz,kxx,kyy,kzz
100 my_real ,DIMENSION(NEL,NPT) ,INTENT(IN) :: dpla
101 my_real ,DIMENSION(NEL,2) ,INTENT(IN) :: eint
102 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: off
103 my_real ,INTENT(IN) :: dtime
104 my_real ,DIMENSION(NEL,NPT) ,INTENT(IN) :: sigy
105C
106 TYPE (ELBUF_STRUCT_) ,INTENT(INOUT) :: ELBUF_STR
107 TYPE (FAIL_PARAM_) ,INTENT(IN) :: FAIL
108 TARGET :: elbuf_str
109C-----------------------------------------------
110C L o c a l V a r i a b l e s
111C-----------------------------------------------
112 INTEGER :: I,IFL,IPT,NFUNC,IPY,IPZ,IPA,NVARF,NPARAM,IRUPT
113 INTEGER :: II(3)
114 INTEGER ,DIMENSION(NEL) :: COUNT
115 my_real :: t0,tm
116 my_real ,DIMENSION(NEL) :: shfact,epsxx,epsxy,epsxz,
117 . signxx,signxy,signxz,ypt,zpt,apt,tstar,depsxx,depsxy,depsxz
118 my_real ,DIMENSION(NEL) :: dpla_ipt,pla_ipt,sigy_ipt
119 my_real :: bidon
120C
121 TYPE(l_bufel_) ,POINTER :: LBUF
122 TYPE(buf_fail_),POINTER :: FBUF
123 my_real, DIMENSION(:), POINTER :: uvarf,dfmax,tdel
124 INTEGER, DIMENSION(:), POINTER :: FOFF
125C=======================================================================
126c to avoid compilation error with unused arguments
127c they will be necessary for next development step
128 bidon = eint(1,1)
129c-----------------------------------------------------
130 ipy = 200
131 ipz = 300
132 ipa = 400
133 ifl = 1 ! only one failure model for beams
134 shfact = five_over_6
135C--------------------------------------
136 DO i=1,3
137 ii(i) = nel*(i-1)
138 ENDDO
139 count(1:nel) = 0
140 nparam = fail%NUPARAM
141 nfunc = fail%NFUNC
142C---------------------------------------
143C START LOOP OVER INTEGRATION POINTS
144C---------------------------------------
145 DO ipt = 1,npt
146 lbuf => elbuf_str%BUFLY(1)%LBUF(1,1,ipt)
147 fbuf => elbuf_str%BUFLY(1)%FAIL(1,1,ipt)
148 uvarf => fbuf%FLOC(ifl)%VAR
149 nvarf = fbuf%FLOC(ifl)%NVAR
150 irupt = fbuf%FLOC(ifl)%ILAWF
151 dfmax => fbuf%FLOC(ifl)%DAMMX
152 tdel => fbuf%FLOC(ifl)%TDEL
153 foff => fbuf%FLOC(ifl)%OFF
154 pla_ipt = elbuf_str%BUFLY(1)%LBUF(1,1,ipt)%PLA
155 dpla_ipt = dpla(:,ipt)
156 sigy_ipt = sigy(:,ipt)
157
158c
159C--- Coordinates of integration points
160 DO i=1,nel
161 ypt(i) = geo(ipy+ipt,iprop)
162 zpt(i) = geo(ipz+ipt,iprop)
163 apt(i) = geo(ipa+ipt,iprop)
164 ENDDO
165C
166 DO i=1,nel
167 signxx(i) = lbuf%SIG(ii(1)+i)
168 signxy(i) = lbuf%SIG(ii(2)+i)
169 signxz(i) = lbuf%SIG(ii(3)+i)
170 ENDDO
171c--- Total strain
172 DO i= 1,nel
173 epsxx(i) = lbuf%STRA(ii(1)+i)
174 epsxy(i) = lbuf%STRA(ii(2)+i)
175 epsxz(i) = lbuf%STRA(ii(3)+i)
176 END DO
177c--- Incremental strain
178 DO i = 1,nel
179 depsxx(i) = exx(i) - ypt(i)*kzz(i) + zpt(i)*kyy(i)
180 depsxy(i) = exy(i) + zpt(i)*kxx(i)
181 depsxz(i) = exz(i) - ypt(i)*kxx(i)
182 depsxy(i) = depsxy(i) / shfact(i)
183 depsxz(i) = depsxz(i) / shfact(i)
184 ENDDO
185c
186c------------------------------------
187 SELECT CASE (irupt)
188c------------------------------------
189c
190c---------------
191 CASE (1) ! Johnson-Cook
192 ! Tstar computation for Jhonson-Cook failure : T* = (T-T0)/(TM-T0)
193 IF (jthe > 0) THEN
194 t0 = pm(79, imat)
195 tm = pm(80, imat)
196 DO i=1,nel
197 tstar(i) = max(zero,(tempel(i)-t0)/(tm-t0))
198 ENDDO
199 ELSE
200 tstar(1:nel) = zero
201 ENDIF
202 CALL fail_johnson_ib(
203 . nel ,ngl ,ipt ,nparam ,fail%UPARAM,
204 . time ,tstar ,signxx ,signxy ,signxz ,
205 . dpla_ipt ,epsd ,off ,foff ,dfmax ,
206 . tdel ,iout ,istdo ,elbuf_str%GBUF%UELR ,npt)
207c---------------
208 CASE (10) ! Tension Strain failure model
209 IF (jthe > 0) THEN
210 t0 = pm(79, imat)
211 tm = pm(80, imat)
212 DO i=1,nel
213 tstar(i) = max(zero,(tempel(i)-t0)/(tm-t0))
214 ENDDO
215 ELSE
216 tstar(1:nel) = zero
217 ENDIF
218
219 CALL fail_tensstrain_ib(
220 . nel ,ngl ,nparam ,fail%UPARAM ,
221 . time ,epsd ,off ,dfmax,
222 . tdel ,iout ,istdo ,fail%IFUNC ,
223 . epsxx ,al ,tstar ,lbuf%DMGSCL ,
224 . snpc ,npf ,stf ,uvarf ,nvarf,
225 . tf ,ipt ,foff ,ismstr ,elbuf_str%GBUF%UELR ,npt)
226c---------------
227 CASE (11) ! Energy failure model
228 CALL fail_energy_ib(
229 . nel ,ngl ,nparam ,fail%UPARAM ,
230 . time ,epsd ,off ,dfmax,
231 . tdel ,iout ,istdo ,fail%IFUNC ,
232 . lbuf%DMGSCL,uvarf ,nvarf ,
233 . snpc ,npf ,stf ,
234 . tf ,ipt ,foff ,
235 . signxx ,signxy ,signxz ,
236 . depsxx ,depsxy ,depsxz ,elbuf_str%GBUF%UELR ,npt)
237
238c---------------
239 CASE (30) ! BIQUAD
240 CALL fail_biquad_ib (
241 . nel ,ngl ,nparam ,fail%UPARAM,
242 . time ,off ,dfmax ,tdel ,
243 . iout ,istdo ,nfunc ,fail%IFUNC,lbuf%DMGSCL,
244 . uvarf ,nvarf ,snpc ,npf ,
245 . stf ,tf ,ipt ,foff ,
246 . signxx ,signxy ,signxz ,dpla_ipt ,al,elbuf_str%GBUF%UELR ,npt)
247
248c---------------
249 CASE (34) ! cockroft failure model
250
251 CALL fail_cockroft_ib(
252 . nel ,ngl ,nparam ,fail%UPARAM ,
253 . time ,dpla_ipt ,off ,dfmax,
254 . tdel ,iout ,istdo ,epsxx ,
255 . ipt ,signxx ,signxy ,signxz ,
256 . nvarf ,uvarf ,foff ,elbuf_str%GBUF%UELR ,npt)
257
258C--------------
259 CASE (36) ! visual failure model
260 CALL fail_visual_ib(
261 . nel ,ngl ,nparam ,fail%UPARAM,
262 . time ,dfmax,
263 . iout ,istdo ,
264 . signxx ,signxy ,signxz ,
265 . epsxx ,epsxy ,epsxz ,
266 . nvarf ,
267 . uvarf ,ismstr ,dtime , ipt)
268c-------------
269 CASE (39) ! GENE1 failure model
270 CALL fail_gene1_ib(
271 . nel ,fail%NUPARAM ,fbuf%FLOC(ifl)%NVAR ,fail%NFUNC,
272 . fail%IFUNC ,npf ,tf ,
273 . time ,dtime ,fail%UPARAM ,ipt ,
274 . ngl ,elbuf_str%GBUF%DT,epsd ,fbuf%FLOC(ifl)%VAR,
275 . off ,epsxx ,
276 . epsxy ,epsxz ,
277 . signxx ,signxy ,
278 . signxz ,tempel ,
279 . fbuf%FLOC(ifl)%DAMMX , al, table,
280 . fail%NTABLE ,
281 . fail%TABLE , fbuf%FLOC(ifl)%LF_DAMMX, fail%NIPARAM ,fail%IPARAM,
282 . snpc ,stf ,ntable, foff,elbuf_str%GBUF%UELR ,npt, tdel)
283
284
285c-------------
286 CASE (42) ! Inievo failure model
287
288 CALL fail_inievo_ib(
289 . nel ,fail%NUPARAM ,fbuf%FLOC(ifl)%NVAR ,
290 . table ,fail%NTABLE ,fail%TABLE ,time ,fail%UPARAM ,
291 . ngl ,al ,dpla_ipt ,epsd , fbuf%FLOC(ifl)%VAR,
292 . signxx ,signxy ,signxz ,
293 . pla_ipt ,sigy_ipt ,foff ,fbuf%FLOC(ifl)%DAMMX ,
294 . tdel ,lbuf%DMGSCL ,elbuf_str%GBUF%UELR ,ipt ,npt,
295 . fbuf%FLOC(ifl)%DAMINI ,
296 . ntable ,off)
297
298c-------------
299 CASE (41) ! TAB2 failure model
300
301 CALL fail_tab2_ib(
302 . nel ,fail%NUPARAM ,fbuf%FLOC(ifl)%NVAR ,fail%NFUNC ,fail%IFUNC ,
303 . npf ,table ,tf ,time ,fail%UPARAM ,
304 . ngl ,al ,dpla_ipt ,epsd ,fbuf%FLOC(ifl)%VAR,
305 . signxx ,signxy ,signxz ,
306 . tempel ,off ,fbuf%FLOC(ifl)%DAMMX ,tdel ,lbuf%DMGSCL,
307 . ipt, npt, foff, fail%NTABLE ,fail%TABLE,
308 . elbuf_str%GBUF%UELR,
309 . snpc ,stf ,ntable)
310 DO i= 1,nel
311 IF (foff(i) == 0) count(i) = count(i) + 1
312 END DO
313
314c-------------------------------------
315 END SELECT
316
317C-------------------------------------
318 ENDDO ! IPT = 1,NPT : END LOOP OVER INTEGRATION POINTS
319c------------------
320 RETURN
321 END SUBROUTINE fail_beam18
my_real
#define my_real
Definition cppsort.cpp:32
fail_beam18
subroutine fail_beam18(elbuf_str, fail, nummat, numgeo, npropm, npropg, snpc, stf, nel, npt, imat, iprop, jthe, tempel, ngl, pm, geo, off, epsd, npf, tf, dpla, eint, time, iout, istdo, al, ismstr, exx, exy, exz, kxx, kyy, kzz, dtime, ntable, table, sigy)
Definition fail_beam18.F:54
fail_biquad_ib
subroutine fail_biquad_ib(nel, ngl, nuparam, uparam, time, off, dfmax, tdel, iout, istdo, nfunc, ifunc, damscl, uvar, nuvar, snpc, npf, stf, tf, ipt, foff, signxx, signxy, signxz, dpla, al, uelr, npg)
Definition fail_biquad_ib.F:37
fail_cockroft_ib
subroutine fail_cockroft_ib(nel, ngl, nuparam, uparam, time, dpla, off, dfmax, tdel, iout, istdo, epsxx, ipt, signxx, signxy, signxz, nvarf, uvar, foff, uelr, npg)
Definition fail_cockroft_ib.F:34
fail_energy_ib
subroutine fail_energy_ib(nel, ngl, nuparam, uparam, time, epsd, off, dfmax, tdel, iout, istdo, ifunc, damscl, uvar, nvar, snpc, npf, stf, tf, ipt, foff, signxx, signxy, signxz, depsxx, depsxy, depsxz, uelr, npg)
Definition fail_energy_ib.F:39
fail_johnson_ib
subroutine fail_johnson_ib(nel, ngl, ipt, nuparam, uparam, time, tstar, signxx, signxy, signzx, dpla, epsd, off, foff, dfmax, tdel, iout, istdo, uelr, npg)
Definition fail_johnson_ib.F:33
fail_tensstrain_ib
subroutine fail_tensstrain_ib(nel, ngl, nuparam, uparam, time, epsd, off, dfmax, tdel, iout, istdo, ifunc, epsxx, al, tstar, damscl, snpc, npf, stf, uvar, nvarf, tf, ipt, foff, ismstr, uelr, npg)
Definition fail_tensstrain_ib.F:37
max
#define max(a, b)
Definition macros.h:21
main_beam18
subroutine main_beam18(elbuf_str, nel, npt, mtn, imat, pid, ngl, pm, ipm, geo, off, for, mom, eint, al, epsd, bufmat, npf, tf, exx, exy, exz, kxx, kyy, kzz, f1, f2, f3, m1, m2, m3, jthe, tempel, ifail, sbufmat, snpc, stf, nummat, numgeo, iout, istdo, npropmi, npropm, npropg, time, dtime, idel7nok, isigi, imconv, ismstr, mat_param, ntable, table)
Definition main_beam18.F:48
bidon
subroutine bidon
Definition machine.F:41