OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
mulawglcpinch.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!|| mulawglcpinch ../engine/source/elements/shell/coqueba/mulawglcpinch.F
25!||--- called by ------------------------------------------------------
26!|| cmain3pinch ../engine/source/elements/shell/coqueba/cmain3pinch.F
27!||--- calls -----------------------------------------------------
28!|| sigeps01gpinch ../engine/source/materials/mat/mat001/sigeps01gpinch.F
29!|| sigeps91gpinch ../engine/source/materials/mat/mat091/sigeps91gpinch.F
30!||--- uses -----------------------------------------------------
31!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
32!|| table_mod ../engine/share/modules/table_mod.F
33!||====================================================================
34 SUBROUTINE mulawglcpinch(ELBUF_STR,
35 1 JFT ,JLT ,PM ,FOR ,MOM ,THK ,
36 2 EINT ,OFF ,GSTR ,PLA ,DIR ,SHF ,
37 3 MAT ,AREA ,EXX ,EYY ,EXY ,NEL ,
38 4 EXZ ,EYZ ,KXX ,KYY ,KXY ,DM ,
39 5 PID ,TF ,NPF ,MTN ,DT1C ,A1 ,
40 6 BUFMAT ,SSP ,RHO ,VISCMX,IOFC ,A2 ,
41 7 INDX ,NGL ,ZCFAC ,GS ,SIGY ,G ,
42 8 THK0 ,EPSP ,IPLA ,IGEO ,IPM ,TABLE ,
43 9 IR ,IS ,F_DEF ,ISMSTR,NU ,VOL0 ,
44 A KFTS ,EPINCHZZ ,EPINCHXZ ,
45 B EPINCHYZ ,FORP ,MOMP ,ALDT ,
46 C EZZAVG ,AREAPINCH)
47C-----------------------------------------------
48 USE table_mod
49 USE elbufdef_mod
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54#include "comlock.inc"
55C-----------------------------------------------
56C G l o b a l P a r a m e t e r s
57C-----------------------------------------------
58#include "mvsiz_p.inc"
59C-----------------------------------------------
60C C o m m o n B l o c k s
61C-----------------------------------------------
62#include "units_c.inc"
63#include "scr17_c.inc"
64#include "param_c.inc"
65#include "com08_c.inc"
66#include "impl1_c.inc"
67C-----------------------------------------------
68C D u m m y A r g u m e n t s
69C-----------------------------------------------
70 INTEGER JFT, JLT, MTN, IOFC, IPLA,NEL,IR,IS,ISMSTR,KFTS
71 INTEGER MAT(*), PID(*), NPF(*),NGL(*), INDX(*),IPM(NPROPMI,*)
72 my_real DM
73 my_real FOR(NEL,5), MOM(NEL,3), THK(*), EINT(JLT,2),PM(NPROPM,*),
74 . OFF(*), GSTR(NEL,8), PLA(*), DIR(*),VISCMX(*),
75 . AREA(*),TF(*),DT1C(*),
76 . EXX(*), EYY(*), EXY(*), EXZ(*), EYZ(*), EPSP(*),
77 . KXX(*), KYY(*), KXY(*),BUFMAT(*),SSP(*),RHO(*),
78 . ZCFAC(MVSIZ,2),GS(*),SIGY(*),THK0(*),SHF(*),F_DEF(MVSIZ,8),
79 . A1(MVSIZ),A2(MVSIZ),G(MVSIZ),NU(MVSIZ),VOL0(*),
80 . EPINCHZZ(MVSIZ),EPINCHXZ(MVSIZ),EPINCHYZ(MVSIZ),
81 . FORP(NEL),MOMP(NEL,2),ALDT(MVSIZ),EZZAVG(MVSIZ),AREAPINCH(MVSIZ)
82 TYPE(ttable) TABLE(*)
83 TYPE(elbuf_struct_), TARGET :: ELBUF_STR
84C-----------------------------------------------
85C L o c a l V a r i a b l e s
86C-----------------------------------------------
87 INTEGER IGTYP, I, NUVAR ,NINDX,IGEO(NPROPGI,*),
88 . MX,IOFF_DUCT(MVSIZ),ISRATE
89C REAL
90 my_real
91 . DEGMB(MVSIZ) ,DEGFX(MVSIZ) ,THKN(MVSIZ),
92 . DEPSXX(MVSIZ),DEPSYY(MVSIZ),DEPSZZ(MVSIZ),
93 . DEPSXY(MVSIZ),DEPSYZ(MVSIZ),DEPSZX(MVSIZ),
94 . sigoxx(mvsiz),sigoyy(mvsiz),sigozz(mvsiz),
95 . sigoxy(mvsiz),sigoyz(mvsiz),sigozx(mvsiz),
96 . signxx(mvsiz),signyy(mvsiz),signzz(mvsiz),
97 . signxy(mvsiz),signyz(mvsiz),signzx(mvsiz),
98 . sigvxx(mvsiz),sigvyy(mvsiz),sigvxy(mvsiz),
99 . eps_m2,eps_k2,young, visc, vol2, asrate
100 my_real
101 . depbxx(mvsiz),depbyy(mvsiz),depbxy(mvsiz),
102 . deppxz(mvsiz),deppyz(mvsiz),
103 . momoxx(mvsiz),momoyy(mvsiz),momoxy(mvsiz),
104 . momopxz(mvsiz),momopyz(mvsiz),
105 . momnxx(mvsiz),momnyy(mvsiz),momnxy(mvsiz),
106 . momnpxz(mvsiz),momnpyz(mvsiz),
107 . etse(mvsiz) ,epspl(mvsiz),epsp_loc(mvsiz)
108 my_real,
109 . DIMENSION(:) ,POINTER :: uvar
110 my_real, DIMENSION(MVSIZ) :: dt_inv
111C-------------------------------------
112 TYPE(l_bufel_) ,POINTER :: LBUF
113C-------------------------------------
114C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
115C
116 LBUF => elbuf_str%BUFLY(1)%LBUF(ir,is,1)
117C
118C pour les lois user pid(i)=pid(1)
119 igtyp = igeo(11,pid(1))
120 nuvar = elbuf_str%BUFLY(1)%NVAR_MAT
121 uvar =>elbuf_str%BUFLY(1)%MAT(ir,is,1)%VAR
122 ioff_duct(1:mvsiz) = 0
123 viscmx(1:mvsiz) = zero
124C
125 DO i=jft,jlt
126 degmb(i) = for(i,1)*exx(i)+for(i,2)*eyy(i)+for(i,3)*exy(i)
127 . + for(i,4)*eyz(i)+for(i,5)*exz(i)
128 . + half*forp(i)*epinchzz(i)
129 degfx(i) = mom(i,1)*kxx(i)+mom(i,2)*kyy(i)+mom(i,3)*kxy(i)
130 . + half*momp(i,1)*epinchxz(i)+half*momp(i,2)*epinchyz(i)
131 ENDDO
132C
133 DO i=jft,jlt
134 thkn(i) = thk(i)
135 ENDDO
136! compute the inverse of dt and save the result
137 DO i=jft,jlt
138 dt_inv(i) = dt1c(i)/max(dt1c(i)**2,em20)
139 ENDDO
140C-----------------------------------------------------------
141C EPS POINT EQUIVALENT (au sens energetique)
142C-----------------------------------------------------------
143C e = 1/t integ[1/2 E (eps_m + k z)^2 dz ]
144C e = 1/2 E eps_eq^2
145C eps_eq = sqrt[ eps_m^2 + 1/12 k^2t^2 ]
146 mx = mat(jft)
147!
148 israte = ipm(3,mx)
149 IF (israte == 1) THEN
150!#include "vectorize.inc"
151 DO i=jft,jlt
152 eps_k2 = (kxx(i)*kxx(i)+kyy(i)*kyy(i)+kxx(i)*kyy(i)
153 . + fourth*kxy(i)*kxy(i)) *thk(i)*thk(i)*one_over_9
154 eps_m2 = four_over_3*(exx(i)*exx(i)+eyy(i)*eyy(i)+exx(i)*eyy(i)
155 . + fourth*exy(i)*exy(i))
156 epsp_loc(i) = sqrt(eps_k2 + eps_m2)*dt_inv(i)
157 END DO
158 ELSE IF (israte > 1) THEN
159!#include "vectorize.inc"
160 DO i=jft,jlt
161 eps_m2 = four_over_3*(exx(i)*exx(i)+eyy(i)*eyy(i)+exx(i)*eyy(i)
162 . + fourth*exy(i)*exy(i))
163 epsp_loc(i) = sqrt(eps_m2)*dt_inv(i)
164 END DO
165 ENDIF
166!
167 IF (israte > 0) THEN ! strain rate filtering with exponential average
168 DO i=jft,jlt
169 asrate = min(one,pm(9,mx)*dt1c(i))
170 epsp(i) = asrate*epsp_loc(i) + (one-asrate)*epsp(i)
171 epspl(i) = epsp(i)
172 END DO
173 ENDIF
174C-----------------------
175C PLASTICITE GLOBALE
176C-----------------------
177 DO i=jft,jlt
178 signxx(i) = zero
179 signyy(i) = zero
180 signzz(i) = zero
181 signxy(i) = zero
182 signyz(i) = zero
183 signzx(i) = zero
184 ENDDO
185 DO i=jft,jlt
186 momnxx(i) = zero
187 momnyy(i) = zero
188 momnxy(i) = zero
189 momnpxz(i)=zero
190 momnpyz(i)=zero
191 ENDDO
192C
193 DO i=jft,jlt
194 depsxx(i) = exx(i)
195 depsyy(i) = eyy(i)
196 depszz(i) = epinchzz(i)
197 depsxy(i) = exy(i)
198 depsyz(i) = eyz(i)
199 depszx(i) = exz(i)
200 sigoxx(i) = for(i,1)
201 sigoyy(i) = for(i,2)
202 sigozz(i) = forp(i)
203 sigoxy(i) = for(i,3)
204 sigoyz(i) = for(i,4)
205 sigozx(i) = for(i,5)
206 ENDDO
207C
208 DO i=jft,jlt
209 depbxx(i) = kxx(i)
210 depbyy(i) = kyy(i)
211 depbxy(i) = kxy(i)
212 deppxz(i) = epinchxz(i)
213 deppyz(i) = epinchyz(i)
214 momoxx(i) = mom(i,1)
215 momoyy(i) = mom(i,2)
216 momoxy(i) = mom(i,3)
217 momopxz(i)= momp(i,1)
218 momopyz(i)= momp(i,2)
219 ENDDO
220C
221 IF(mtn == 1) THEN
222 CALL sigeps01gpinch(
223 1 jft ,jlt ,g ,thkn ,off ,
224 2 gs ,a1 ,a2 ,nu ,thk0 ,
225 3 nel ,ssp ,rho ,
226 4 depsxx ,depsyy ,depszz ,
227 5 depsxy ,depsyz ,depszx ,
228 6 depbxx ,depbyy ,depbxy ,
229 7 deppxz ,deppyz ,
230 8 sigoxx ,sigoyy ,sigozz ,
231 9 sigoxy ,sigoyz ,sigozx ,
232 a momoxx ,momoyy ,momoxy ,
233 b momopxz ,momopyz ,
234 c signxx ,signyy ,signzz ,
235 d signxy ,signyz ,signzx ,
236 e momnxx ,momnyy ,momnxy ,
237 f momnpxz ,momnpyz)
238 ELSEIF(mtn == 91) THEN
239 CALL sigeps91gpinch(
240 1 jft ,jlt ,nuvar ,bufmat ,rho ,
241 2 thkn ,thk0 ,nel ,ssp ,area ,
242 3 depsxx ,depsyy ,depszz ,
243 4 depsxy ,depsyz ,depszx ,
244 5 depbxx ,depbyy ,depbxy ,
245 6 deppxz ,deppyz ,
246 7 sigoxx ,sigoyy ,sigozz ,
247 8 sigoxy ,sigoyz ,sigozx ,
248 9 momoxx ,momoyy ,momoxy ,
249 a momopxz ,momopyz ,
250 b signxx ,signyy ,signzz ,
251 c signxy ,signyz ,signzx ,
252 d momnxx ,momnyy ,momnxy ,
253 e momnpxz ,momnpyz ,tt ,uvar ,dt_inv ,
254 f viscmx ,aldt ,vol0 ,ipm ,mat ,
255 g pla ,degmb ,degfx ,
256 h ngl ,ezzavg ,areapinch)
257 ENDIF
258C-----------------------
259C FORCES ET MOMENTS
260C-----------------------
261 DO i=jft,jlt
262 for(i,1) = signxx(i)
263 for(i,2) = signyy(i)
264 for(i,3) = signxy(i)
265 for(i,4) = signyz(i)
266 for(i,5) = signzx(i)
267 mom(i,1) = momnxx(i)
268 mom(i,2) = momnyy(i)
269 mom(i,3) = momnxy(i)
270 forp(i) = signzz(i)
271 momp(i,1)= momnpxz(i)
272 momp(i,2)= momnpyz(i)
273 ENDDO
274C
275 DO i=jft,jlt
276 for(i,1) =for(i,1)*off(i)
277 for(i,2) =for(i,2)*off(i)
278 for(i,3) =for(i,3)*off(i)
279 for(i,4) =for(i,4)*off(i)
280 for(i,5) =for(i,5)*off(i)
281 mom(i,1) =mom(i,1)*off(i)
282 mom(i,2) =mom(i,2)*off(i)
283 mom(i,3) =mom(i,3)*off(i)
284 forp(i) =forp(i)*off(i)
285 momp(i,1) =momp(i,1)*off(i)
286 momp(i,2) =momp(i,2)*off(i)
287 ENDDO
288C
289 DO i=jft,jlt
290 degmb(i) = degmb(i)+for(i,1)*exx(i)+for(i,2)*eyy(i)+for(i,3)*exy(i)
291 . + for(i,4)*eyz(i)+for(i,5)*exz(i)
292 . + half*forp(i)*epinchzz(i)
293 degfx(i) = degfx(i)+mom(i,1)*kxx(i)+mom(i,2)*kyy(i)+mom(i,3)*kxy(i)
294 . + half*momp(i,1)*epinchxz(i)+half*momp(i,2)*epinchyz(i)
295 vol2 = half*thk0(i)*area(i)*off(i)
296 eint(i,1) = eint(i,1) + degmb(i)*vol2
297 eint(i,2) = eint(i,2) + degfx(i)*thk0(i)*vol2
298 ENDDO
299C----------------------------
300C TEST DE RUPTURE DUCTILE
301C INDX utilise dans IELOF
302C----------------------------
303 nindx=0
304 DO i=jft,jlt
305 IF (off(i) == four_over_5 . and. ioff_duct(i) == 0) THEN
306 off(i)= zero
307 nindx=nindx+1
308 indx(nindx)=i
309 ENDIF
310 ENDDO
311 IF (nindx > 0) THEN
312 idel7nok = 1
313 IF (imconv == 1) THEN
314 DO i = 1, nindx
315#include "lockon.inc"
316 WRITE(iout, 1000) ngl(indx(i))
317 WRITE(istdo,1100) ngl(indx(i)),tt
318#include "lockoff.inc"
319 ENDDO
320 ENDIF
321 ENDIF
322 iofc = nindx
323C-----------
324 1000 FORMAT(1x,'-- RUPTURE OF SHELL ELEMENT NUMBER ',i10)
325 1100 FORMAT(1x,'-- RUPTURE OF SHELL ELEMENT :',i10,' AT TIME :',g11.4)
326C-----------
327 RETURN
328 END
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
mulawglcpinch
subroutine mulawglcpinch(elbuf_str, jft, jlt, pm, for, mom, thk, eint, off, gstr, pla, dir, shf, mat, area, exx, eyy, exy, nel, exz, eyz, kxx, kyy, kxy, dm, pid, tf, npf, mtn, dt1c, a1, bufmat, ssp, rho, viscmx, iofc, a2, indx, ngl, zcfac, gs, sigy, g, thk0, epsp, ipla, igeo, ipm, table, ir, is, f_def, ismstr, nu, vol0, kfts, epinchzz, epinchxz, epinchyz, forp, momp, aldt, ezzavg, areapinch)
Definition mulawglcpinch.F:47
table_mod
Definition table_mod.F:112
sigeps01gpinch
subroutine sigeps01gpinch(jft, jlt, g, thk, off, gs, a1, a2, nu, thk0, nel, ssp, rho, depsxx, depsyy, depszz, depsxy, depsyz, depszx, depbxx, depbyy, depbxy, deppxz, deppyz, sigoxx, sigoyy, sigozz, sigoxy, sigoyz, sigozx, momoxx, momoyy, momoxy, momopxz, momopyz, signxx, signyy, signzz, signxy, signyz, signzx, momnxx, momnyy, momnxy, momnpxz, momnpyz)
Definition sigeps01gpinch.F:44
sigeps91gpinch
subroutine sigeps91gpinch(jft, jlt, nuvar, uparam, rho0, thk, thk0, nel, ssp, area, depsxx, depsyy, depszz, depsxy, depsyz, depszx, depbxx, depbyy, depbxy, deppxz, deppyz, sigoxx, sigoyy, sigozz, sigoxy, sigoyz, sigozx, momoxx, momoyy, momoxy, momopxz, momopyz, signxx, signyy, signzz, signxy, signyz, signzx, momnxx, momnyy, momnxy, momnpxz, momnpyz, time, uvar, dt_inv, viscmx, aldt, vol0, ipm, mat, pla, degmb, degfx, ngl, ezzavg, areapinch)
Definition sigeps91gpinch.F:46
table_mod::ttable
Definition table_mod.F:125