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

Go to the source code of this file.

Functions/Subroutines

subroutine main_beam3 (elbuf_str, nel, ilaw, jthe, ifail, ipm, pm, geo, tempel, off, mat, pid, ngl, time, dtime, al, npf, tf, exx, exy, exz, kxx, kyy, kzz, f1, f2, f3, m1, m2, m3, bufmat, npropg, npropmi, npropm, nummat, numgeo, sbufmat, snpc, stf, iout, istdo, nuvar, uvar, epsd, imat, for, mom, eint, ismstr, mat_param, ntable, table)

Function/Subroutine Documentation

◆ main_beam3()

subroutine main_beam3 ( type (elbuf_struct_), intent(inout) elbuf_str,
integer, intent(in) nel,
integer, intent(in) ilaw,
integer, intent(in) jthe,
integer, intent(in) ifail,
integer, dimension(npropmi,nummat), intent(in) ipm,
dimension(npropm,nummat), intent(in) pm,
dimension(npropg,numgeo), intent(in) geo,
intent(inout) tempel,
intent(inout) off,
integer, dimension(nel), intent(in) mat,
integer, dimension(nel), intent(in) pid,
integer, dimension(nel), intent(in) ngl,
intent(in) time,
intent(in) dtime,
intent(in) al,
integer, dimension(snpc), intent(in) npf,
dimension(stf), intent(in) tf,
intent(inout) exx,
intent(inout) exy,
intent(inout) exz,
intent(inout) kxx,
intent(inout) kyy,
intent(inout) kzz,
intent(inout) f1,
intent(inout) f2,
intent(inout) f3,
intent(inout) m1,
intent(inout) m2,
intent(inout) m3,
dimension(sbufmat), intent(in) bufmat,
integer, intent(in) npropg,
integer, intent(in) npropmi,
integer, intent(in) npropm,
integer, intent(in) nummat,
integer, intent(in) numgeo,
integer, intent(in) sbufmat,
integer, intent(in) snpc,
integer, intent(in) stf,
integer, intent(in) iout,
integer, intent(in) istdo,
integer, intent(in) nuvar,
intent(inout) uvar,
intent(out) epsd,
integer, intent(in) imat,
intent(inout) for,
intent(inout) mom,
intent(inout) eint,
integer, intent(in) ismstr,
type (matparam_struct_), intent(in) mat_param,
integer, intent(in) ntable,
type(ttable), dimension(ntable), intent(inout) table )

Definition at line 36 of file main_beam3.F.

48C-----------------------------------------------
49C M o d u l e s
50C-----------------------------------------------
51 USE mat_elem_mod
52 USE elbufdef_mod
53C-----------------------------------------------
54C I m p l i c i t T y p e s
55C-----------------------------------------------
56#include "implicit_f.inc"
57C-----------------------------------------------
58C D u m m y A r g u m e n t s
59C-----------------------------------------------
60 INTEGER ,INTENT(IN) :: NEL
61 INTEGER ,INTENT(IN) :: ILAW
62 INTEGER ,INTENT(IN) :: IMAT
63 INTEGER ,INTENT(IN) :: JTHE
64 INTEGER ,INTENT(IN) :: IFAIL
65 INTEGER ,INTENT(IN) :: SBUFMAT
66 INTEGER ,INTENT(IN) :: SNPC
67 INTEGER ,INTENT(IN) :: STF
68 INTEGER ,INTENT(IN) :: NUMMAT
69 INTEGER ,INTENT(IN) :: NUMGEO
70 INTEGER ,INTENT(IN) :: NPROPMI
71 INTEGER ,INTENT(IN) :: NPROPM
72 INTEGER ,INTENT(IN) :: NPROPG
73 INTEGER ,INTENT(IN) :: IOUT
74 INTEGER ,INTENT(IN) :: ISTDO
75 INTEGER ,INTENT(IN) :: NUVAR
76 INTEGER ,INTENT(IN) :: ISMSTR
77 INTEGER ,INTENT(IN) :: NTABLE
78 TYPE(TTABLE), DIMENSION(NTABLE), INTENT(INOUT) :: TABLE ! TABLE DATA
79 INTEGER ,DIMENSION(SNPC) ,INTENT(IN) :: NPF
80 INTEGER ,DIMENSION(NPROPMI,NUMMAT) ,INTENT(IN) :: IPM
81 INTEGER ,DIMENSION(NEL) ,INTENT(IN) :: MAT
82 INTEGER ,DIMENSION(NEL) ,INTENT(IN) :: PID
83 INTEGER ,DIMENSION(NEL) ,INTENT(IN) :: NGL
84 my_real ,INTENT(IN) :: time
85 my_real ,INTENT(IN) :: dtime
86 my_real ,INTENT(IN) :: pm(npropm,nummat)
87 my_real ,INTENT(IN) :: geo(npropg,numgeo)
88 my_real ,INTENT(IN) :: tf(stf)
89 my_real ,INTENT(IN) :: bufmat(sbufmat)
90 my_real ,DIMENSION(NEL) ,INTENT(IN) :: al
91 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: tempel
92 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: exx,exy,exz
93 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: kxx,kyy,kzz
94 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: off
95 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: f1,f2,f3
96 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: m1,m2,m3
97 my_real ,DIMENSION(NEL,2),INTENT(INOUT) :: eint
98 my_real ,DIMENSION(NEL,3),INTENT(INOUT) :: for,mom
99 my_real ,DIMENSION(NEL) ,INTENT(OUT) :: epsd
100 my_real ,DIMENSION(NEL,NUVAR) ,INTENT(INOUT) :: uvar
101 TYPE (ELBUF_STRUCT_) ,INTENT(INOUT) :: ELBUF_STR
102 TYPE (MATPARAM_STRUCT_) ,INTENT(IN) :: MAT_PARAM
103C-----------------------------------------------
104C L o c a l V a r i a b l e s
105C-----------------------------------------------
106 INTEGER :: I,IPID,IPLA,ISRATE,NUPARAM,NFUNC,IADBUF,IFUNC(100)
107 my_real :: area,iyy,izz,ixx,asrate,epsdi,fact
108 my_real ,DIMENSION(NEL) :: dpla,svm,pressure,degmb,degfx,degsh,sigy,
109 . svm0
110c=======================================================================
111 ipid = pid(1)
112c-------------------
113c STRAIN RATE
114c-------------------
115 israte = ipm(3,imat)
116 asrate = min(one, pm(9,imat)*dtime)
117c
118 ! Geometric properties of the beam element
119 area = geo(1 ,ipid)
120 ixx = geo(4 ,ipid)
121 iyy = geo(2 ,ipid)
122 izz = geo(18,ipid)
123c
124 DO i = 1,nel
125 epsdi = half*(exx(i)**2) + (half*exy(i))**2 + (half*exz(i))**2
126 epsdi = al(i)*sqrt(three*epsdi)/three_half
127 sigy(i) = 1.e30
128 IF (israte > 0) THEN
129 epsd(i)= asrate*epsdi + (one - asrate)*epsd(i)
130 ELSE
131 epsd(i)= epsdi
132 ENDIF
133 ENDDO
134c-------------------
135c STRAIN
136c-------------------
137 DO i=1,nel
138 ! Strain increment
139 exx(i) = exx(i) * dtime
140 exy(i) = exy(i) * dtime
141 exz(i) = exz(i) * dtime
142 kxx(i) = kxx(i) * dtime
143 kyy(i) = kyy(i) * dtime
144 kzz(i) = kzz(i) * dtime
145 ! Initialization of internal energy increment
146 degmb(i) = for(i,1)*exx(i)
147 degsh(i) = for(i,2)*exy(i)+for(i,3)*exz(i)
148 degfx(i) = mom(i,1)*kxx(i)+mom(i,2)*kyy(i)+mom(i,3)*kzz(i)
149 ENDDO
150c
151 ! Compute undamaged forces and internal energy increment
152 IF (elbuf_str%GBUF%G_DMGSCL > 0) THEN
153 DO i = 1,nel
154 for(i,1) = for(i,1)/max(em20,elbuf_str%GBUF%DMGSCL(i))
155 for(i,2) = for(i,2)/max(em20,elbuf_str%GBUF%DMGSCL(i))
156 for(i,3) = for(i,3)/max(em20,elbuf_str%GBUF%DMGSCL(i))
157 mom(i,1) = mom(i,1)/max(em20,elbuf_str%GBUF%DMGSCL(i))
158 mom(i,2) = mom(i,2)/max(em20,elbuf_str%GBUF%DMGSCL(i))
159 mom(i,3) = mom(i,3)/max(em20,elbuf_str%GBUF%DMGSCL(i))
160 ENDDO
161 ENDIF
162c
163 ipla = elbuf_str%GBUF%G_PLA
164 ! Cumulated plastic strain increment
165 IF (ipla > 0) THEN
166 DO i = 1,nel
167 dpla(i) = elbuf_str%GBUF%PLA(i)
168 svm0(i) = for(i,1)*for(i,1) + three * area *
169 . ( mom(i,1)*mom(i,1) / max(ixx,em20)
170 . + mom(i,2)*mom(i,2) / max(iyy,em20)
171 . + mom(i,3)*mom(i,3) / max(izz,em20))
172 svm0(i) = sqrt(svm0(i)) / area
173 ENDDO
174 ENDIF
175c
176c---------------------------
177c material models
178c---------------------------
179 SELECT CASE(ilaw)
180c
181 CASE (1)
182 CALL m1lawp(
183 . pm ,for ,mom ,geo,
184 . off ,exx ,exy ,exz ,kxx,
185 . kyy ,kzz ,al ,f1 ,f2 ,
186 . f3 ,m1 ,m2 ,m3 ,nel,
187 . mat ,pid )
188c
189 CASE (2) ! Johnson-Cook
190 CALL m2lawp (mat_param ,
191 . for ,mom ,eint ,geo ,
192 . off ,elbuf_str%GBUF%PLA ,exx ,exy ,exz ,
193 . kxx ,kyy ,kzz ,al ,f1 ,
194 . f2 ,f3 ,m1 ,m2 ,m3 ,
195 . nel ,pid ,ngl ,
196 . nuvar ,uvar ,sigy )
197c
198 CASE (44) ! Cowper-Symonds
199 iadbuf = ipm(7 ,imat)
200 nuparam = ipm(9 ,imat)
201 nfunc = ipm(10,imat)
202 DO i=1,nfunc
203 ifunc(i) = ipm(10+i,imat)
204 ENDDO
205 CALL sigeps44p(
206 . nel ,ngl ,mat ,pid ,nuparam ,bufmat(iadbuf),
207 . geo ,off ,elbuf_str%GBUF%PLA ,al ,
208 . exx ,exy ,exz ,kxx ,kyy ,kzz ,
209 . f1 ,f2 ,f3 ,m1 ,m2 ,m3 ,
210 . for ,mom ,pm ,nuvar ,uvar ,nfunc ,
211 . ifunc ,tf ,npf ,sigy )
212c
213 END SELECT
214c
215c---------------------------
216c Plastic Work
217c---------------------------
218 IF (ipla > 0 .OR. ifail > 0) THEN
219 DO i = 1,nel
220 dpla(i) = elbuf_str%GBUF%PLA(i) - dpla(i)
221 ENDDO
222
223 DO i = 1,nel
224 svm(i) = f1(i)*f1(i) + three * area *
225 . ( m1(i)*m1(i) / max(ixx,em20)
226 . + m2(i)*m2(i) / max(iyy,em20)
227 . + m3(i)*m3(i) / max(izz,em20) )
228 svm(i) = sqrt(svm(i)) / area
229 pressure(i) = third * f1(i) / area
230 ENDDO
231
232 IF (ipla > 0) THEN
233 DO i = 1,nel
234 IF (elbuf_str%GBUF%G_WPLA > 0) elbuf_str%GBUF%WPLA(i) = elbuf_str%GBUF%WPLA(i) +
235 . half*(svm(i)+svm0(i))*dpla(i)*area*al(i)
236 ENDDO
237 ENDIF
238
239 ENDIF
240c
241c---------------------------
242c failure models
243c---------------------------
244 IF (ifail > 0) THEN
245c
246 CALL fail_beam3(elbuf_str ,mat_param,mat_param%FAIL(1),
247 . snpc ,stf ,nel ,jthe ,dpla ,
248 . tempel ,ngl ,off ,epsd ,npf ,
249 . tf ,time ,iout ,istdo ,
250 . svm ,pressure,area ,al ,
251 . f1 ,f2 ,f3 ,m1 ,m2 ,
252 . m3 ,ismstr ,exx ,exy ,exz ,
253 . kxx ,kyy ,kzz ,dtime ,
254 . ntable ,table ,elbuf_str%GBUF%PLA , sigy )
255c
256 ! Compute damaged forces
257 IF (elbuf_str%GBUF%G_DMGSCL > 0) THEN
258 DO i = 1,nel
259 f1(i) = f1(i)*elbuf_str%GBUF%DMGSCL(i)
260 f2(i) = f2(i)*elbuf_str%GBUF%DMGSCL(i)
261 f3(i) = f3(i)*elbuf_str%GBUF%DMGSCL(i)
262 m1(i) = m1(i)*elbuf_str%GBUF%DMGSCL(i)
263 m2(i) = m2(i)*elbuf_str%GBUF%DMGSCL(i)
264 m3(i) = m3(i)*elbuf_str%GBUF%DMGSCL(i)
265 ENDDO
266 ENDIF
267c
268 END IF
269C-------------------------------------
270 DO i=1,nel
271 for(i,1)=f1(i)*off(i)
272 for(i,2)=for(i,2)*off(i)
273 for(i,3)=for(i,3)*off(i)
274 mom(i,1)=m1(i)*off(i)
275 mom(i,2)=m2(i)*off(i)
276 mom(i,3)=m3(i)*off(i)
277 ENDDO
278C
279 DO i=1,nel
280 f1(i) = for(i,1)
281 f2(i) = for(i,2)
282 f3(i) = for(i,3)
283 m1(i) = mom(i,1)
284 m2(i) = mom(i,2)
285 m3(i) = mom(i,3)
286 ENDDO
287c
288 ! Update internal energy
289 DO i=1,nel
290 degmb(i) = degmb(i) + for(i,1)*exx(i)
291 degsh(i) = degsh(i) + for(i,2)*exy(i) + for(i,3)*exz(i)
292 degfx(i) = degfx(i) + mom(i,1)*kxx(i) + mom(i,2)*kyy(i) + mom(i,3)*kzz(i)
293 fact = half*off(i)*al(i)
294 eint(i,1) = eint(i,1) + (degsh(i)+degmb(i))*fact
295 eint(i,2) = eint(i,2) + fact*degfx(i)
296 ENDDO
297C-----------------------------------------------
298 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
fail_beam3
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)
Definition fail_beam3.F:54
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
m1lawp
subroutine m1lawp(pm, for, mom, geo, off, exx, exy, exz, kxx, kyy, kzz, al, f1, f2, f3, m1, m2, m3, nel, mat, pid)
Definition m1lawp.F:34
m2lawp
subroutine m2lawp(mat_param, for, mom, eint, geo, off, pla, exx, exy, exz, kxx, kyy, kzz, al, fa1, fa2, fa3, ma1, ma2, ma3, nel, pid, ngl, nuvar, uvar, sigy)
Definition m2lawp.F:37
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
sigeps44p
subroutine sigeps44p(nel, ngl, mat, pid, nuparam, uparam, geo, off, pla, al, exx, exy, exz, kxx, kyy, kzz, fa1, fa2, fa3, ma1, ma2, ma3, for, mom, pm, nuvar, uvar, nfunc, ifunc, tf, npf, sigy)
Definition sigeps44p.F:37