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

Go to the source code of this file.

Functions/Subroutines

subroutine fail_puck_c (nel, nuparam, nuvar, uparam, uvar, time, ngl, ipg, ilay, ipt, signxx, signyy, signxy, signyz, signzx, off, foff, dmg_flag, dmg_scale, dfmax, lf_dammx, tdel, timestep)

Function/Subroutine Documentation

◆ fail_puck_c()

subroutine fail_puck_c ( integer, intent(in) nel,
integer, intent(in) nuparam,
integer, intent(in) nuvar,
intent(in) uparam,
intent(inout) uvar,
intent(in) time,
integer, dimension(nel), intent(in) ngl,
integer, intent(in) ipg,
integer, intent(in) ilay,
integer, intent(in) ipt,
intent(in) signxx,
intent(in) signyy,
intent(in) signxy,
intent(in) signyz,
intent(in) signzx,
intent(in) off,
integer, dimension(nel), intent(inout) foff,
integer, intent(out) dmg_flag,
intent(out) dmg_scale,
intent(inout) dfmax,
integer, intent(in) lf_dammx,
intent(out) tdel,
intent(in) timestep )

Definition at line 30 of file fail_puck_c.F.

36C-----------------------------------------------
37C modified Puck model ------
38C-----------------------------------------------
39C I m p l i c i t T y p e s
40C-----------------------------------------------
41#include "implicit_f.inc"
42C-----------------------------------------------
43C G l o b a l P a r a m e t e r s
44C-----------------------------------------------
45#include "units_c.inc"
46#include "comlock.inc"
47C---------+---------+---+---+--------------------------------------------
48C VAR | SIZE |TYP| RW| DEFINITION
49C---------+---------+---+---+--------------------------------------------
50C NEL | 1 | I | R | SIZE OF THE ELEMENT GROUP NEL
51C NUPARAM | 1 | I | R | SIZE OF THE USER PARAMETER ARRAY
52C UPARAM | NUPARAM | F | R | USER MATERIAL PARAMETER ARRAY
53C NUVAR | 1 | I | R | NUMBER OF USER ELEMENT VARIABLES
54C UVAR |NEL*NUVAR| F |R/W| USER ELEMENT VARIABLE ARRAY
55C---------+---------+---+---+--------------------------------------------
56C TIME | 1 | F | R | CURRENT TIME
57C---------+---------+---+---+--------------------------------------------
58C SIGNXX | NEL | F | W | NEW ELASTO PLASTIC STRESS XX
59C SIGNYY | NEL | F | W | NEW ELASTO PLASTIC STRESS YY
60C ... | | | |
61C---------+---------+---+---+--------------------------------------------
62C OFF | NEL | F | R | DELETED ELEMENT FLAG (=1. ON, =0. OFF)
63C FOFF | NEL | I |R/W| DELETED INTEGRATION POINT FLAG (=1 ON, =0 OFF)
64C DFMAX | NEL | F |R/W| MAX DAMAGE FACTOR
65C TDEL | NEL | F | W | FAILURE TIME
66C DMG_FLAG| 1 | I | W | STRESS REDUCTION FLAG DUE TO DAMAGE
67C DMG_SCALE| NEL | F | W | STRESS REDUCTION FACTOR
68C---------+---------+---+---+--------------------------------------------
69C NGL ELEMENT ID
70C IPG CURRENT GAUSS POINT (in plane)
71C ILAY CURRENT LAYER
72C IPT CURRENT INTEGRATION POINT IN THE LAYER
73C---------+---------+---+---+--------------------------------------------
74C I N P U T A r g u m e n t s
75C-----------------------------------------------
76 INTEGER ,INTENT(IN) :: NEL,NUPARAM,NUVAR,IPG,ILAY,IPT,LF_DAMMX
77 INTEGER ,DIMENSION(NEL) ,INTENT(IN) :: NGL
78 my_real ,INTENT(IN) :: time,timestep
79 my_real ,DIMENSION(NEL) ,INTENT(IN) :: off,
80 . signxx,signyy,signxy,signyz,signzx
81 my_real,DIMENSION(NUPARAM) ,INTENT(IN) :: uparam
82C-----------------------------------------------
83C I N P U T O U T P U T A r g u m e n t s
84C-----------------------------------------------
85 INTEGER ,INTENT(OUT) ::DMG_FLAG
86 INTEGER ,DIMENSION(NEL) ,INTENT(INOUT) :: FOFF
87 my_real ,DIMENSION(NEL,LF_DAMMX),INTENT(INOUT) :: dfmax
88 my_real ,DIMENSION(NEL) ,INTENT(OUT) :: tdel,dmg_scale
89 my_real ,DIMENSION(NEL,NUVAR) ,INTENT(INOUT) :: uvar
90C-----------------------------------------------
91C L o c a l V a r i a b l e s
92C-----------------------------------------------
93 INTEGER :: I,J,NINDX,FSMOOTH
94 INTEGER ,DIMENSION(NEL) :: INDX
95 my_real :: sigt1,sigt2, sigc1,sigc2,fsig12,f1,fa,fb,fc,
96 . pn12,pp12,pn22,fac,tmax,dammx,fcut,asrate,
97 . sxx(nel),syy(nel),sxy(nel)
98C=======================================================================
99 sigt1 = uparam(1)
100 sigt2 = uparam(2)
101 sigc1 = uparam(3)
102 sigc2 = uparam(4)
103 fsig12 = uparam(5)
104 pp12 = uparam(6)
105 pn12 = uparam(7)
106 pn22 = uparam(8)
107 tmax = uparam(9)
108 fcut = uparam(12)
109 IF (fcut > zero) THEN
110 fsmooth = 1
111 asrate = two*pi*fcut*timestep
112 asrate = asrate/(one+asrate)
113 ELSE
114 fsmooth = 0
115 ENDIF
116c
117 dmg_flag = 1
118 nindx = 0
119 IF (fsmooth > 0) THEN
120 DO i = 1,nel
121 sxx(i) = asrate*signxx(i) + (one - asrate)*uvar(i,2)
122 syy(i) = asrate*signyy(i) + (one - asrate)*uvar(i,3)
123 sxy(i) = asrate*signxy(i) + (one - asrate)*uvar(i,4)
124 uvar(i,2) = sxx(i)
125 uvar(i,3) = syy(i)
126 uvar(i,4) = sxy(i)
127 ENDDO
128 ELSE
129 DO i = 1,nel
130 sxx(i) = signxx(i)
131 syy(i) = signyy(i)
132 sxy(i) = signxy(i)
133 ENDDO
134 ENDIF
135c-------------------------------
136 DO i=1,nel
137 IF (off(i) == one .and. foff(i) == 1) THEN
138 f1 = zero
139 fa = zero
140 fb = zero
141 fc = zero
142C
143 IF (uvar(i,1) == zero) THEN
144c
145c fiber criteria
146c
147 IF (sxx(i) >= zero) THEN
148 f1 = sxx(i)/sigt1
149 dfmax(i,2) = max(dfmax(i,2),f1)
150 dfmax(i,2) = min(dfmax(i,2),one)
151 ELSE
152 f1 = -sxx(i)/sigc1
153 dfmax(i,3) = max(dfmax(i,3),f1)
154 dfmax(i,3) = min(dfmax(i,3),one)
155 ENDIF
156c
157c matrix criteria
158c
159 IF (syy(i) >= zero) THEN
160 fac = one - pp12*sigt2/fsig12
161 fac = fac*syy(i)/sigt2
162 fa = sqrt((sxy(i)/fsig12)**2 + fac*fac)
163 . + pp12*syy(i)/fsig12
164 dfmax(i,4) = max(dfmax(i,4),fa)
165 dfmax(i,4) = min(dfmax(i,4),one)
166 ELSE
167 fac = half/(one + pn22)/fsig12
168 fc = (sxy(i)*fac)**2 + (syy(i)/sigc2)**2
169 fc =-fc*sigc2/syy(i)
170 dfmax(i,6) = max(dfmax(i,6),fc)
171 dfmax(i,6) = min(dfmax(i,6),one)
172 ENDIF
173 fb = sqrt(sxy(i)**2 + (pn12*syy(i))**2 )
174 . + pn12*syy(i)
175 fb = fb/fsig12
176 dfmax(i,5) = max(dfmax(i,5),fb)
177 dfmax(i,5) = min(dfmax(i,5),one)
178C
179 dammx = min(one,max(f1,fa,fb,fc))
180 dfmax(i,1) = min(one,dammx)
181 IF (dammx >= one) THEN
182 nindx = nindx+1
183 indx(nindx) = i
184 uvar(i,1) = time
185 ENDIF
186 ENDIF
187c
188 IF (uvar(i,1) > zero) THEN
189 dmg_scale(i) = exp(-(time - uvar(i,1))/tmax)
190 IF (dmg_scale(i) < em02) THEN
191 foff(i) = 0
192 tdel(i) = time
193 dmg_scale(i) = zero
194 ENDIF
195 ENDIF
196 ENDIF
197 ENDDO
198c------------------------
199c print
200c------------------------
201 IF (nindx > 0) THEN
202 DO j=1,nindx
203 i = indx(j)
204#include "lockon.inc"
205 WRITE(iout, 2000) ngl(i),ipg,ilay,ipt
206 WRITE(istdo,2100) ngl(i),ipg,ilay,ipt,time
207#include "lockoff.inc"
208 END DO
209 END IF
210c------------------------
211 2000 FORMAT(1x,'FAILURE (PUCK) OF SHELL ELEMENT ',i10,1x,',GAUSS PT',
212 . i2,1x,',LAYER',i3,1x,',INTEGRATION PT',i3)
213 2100 FORMAT(1x,'FAILURE (PUCK) OF SHELL ELEMENT ',i10,1x,',GAUSS PT',
214 . i2,1x,',LAYER',i3,1x,',INTEGRATION PT',i3,1x,'AT TIME :',1pe12.4)
215c------------------------
216 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21