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

Go to the source code of this file.

Functions/Subroutines

subroutine hm_read_fail_orthenerg (fail, fail_id, irupt, lsubmodel, unitab)

Function/Subroutine Documentation

◆ hm_read_fail_orthenerg()

subroutine hm_read_fail_orthenerg ( type (fail_param_), intent(inout) fail,
integer, intent(in) fail_id,
integer, intent(in) irupt,
type (submodel_data), dimension(*), intent(in) lsubmodel,
type (unit_type_), intent(in) unitab )

Definition at line 36 of file hm_read_fail_orthenerg.F.

38C-----------------------------------------------
39c ROUTINE DESCRIPTION :
40c Read Orthenerg failure model parameters
41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE fail_param_mod
45 USE unitab_mod
46 USE message_mod
47 USE submodel_mod
48 USE hm_option_read_mod
49C-----------------------------------------------
50C I m p l i c i t T y p e s
51C-----------------------------------------------
52#include "implicit_f.inc"
53C-----------------------------------------------
54C C o m m o n B l o c k s
55C-----------------------------------------------
56#include "units_c.inc"
57C-----------------------------------------------
58C D u m m y A r g u m e n t s
59C-----------------------------------------------
60 INTEGER ,INTENT(IN) :: FAIL_ID ! failure model ID
61 INTEGER ,INTENT(IN) :: IRUPT ! failure model number
62 TYPE (UNIT_TYPE_) ,INTENT(IN) :: UNITAB ! table of input units
63 TYPE (SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(*) ! submodel table
64 TYPE (FAIL_PARAM_) ,INTENT(INOUT) :: FAIL ! failure model data structure
65C-----------------------------------------------
66C L o c a l V a r i a b l e s
67C-----------------------------------------------
68 INTEGER FAILIP,ISHAP11T,ISHAP11C,ISHAP22T,ISHAP22C,ISHAP33T,ISHAP33C,
69 . ISHAP12T,ISHAP12C,ISHAP23T,ISHAP23C,ISHAP31T,ISHAP31C,NMOD
70 my_real
71 . pthickfail,sigma_11t,g_11t,sigma_11c,g_11c,sigma_22t,g_22t,
72 . sigma_22c,g_22c,sigma_33t,g_33t,sigma_33c,g_33c,sigma_12t,g_12t,
73 . sigma_12c,g_12c,sigma_23t,g_23t,sigma_23c,g_23c,sigma_31t,g_31t,
74 . sigma_31c,g_31c
75 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
76C=======================================================================
77 is_encrypted = .false.
78 is_available = .false.
79C--------------------------------------------------
80C EXTRACT DATA (IS OPTION CRYPTED)
81C--------------------------------------------------
82 CALL hm_option_is_encrypted(is_encrypted)
83C--------------------------------------------------
84C EXTRACT DATAS (REAL VALUES)
85C--------------------------------------------------
86 ! 1st card: Element deletion parameters, strain formulation
87 CALL hm_get_floatv ('PTHICKFAIL',pthickfail ,is_available,lsubmodel,unitab)
88 CALL hm_get_intv ('NMOD' ,nmod ,is_available,lsubmodel)
89 CALL hm_get_intv ('FAILIP' ,failip ,is_available,lsubmodel)
90c
91 ! 2nd card: direction 11 failure parameters
92 CALL hm_get_floatv ('SIGMA_11T' ,sigma_11t ,is_available,lsubmodel,unitab)
93 CALL hm_get_floatv ('G_11T' ,g_11t ,is_available,lsubmodel,unitab)
94 CALL hm_get_intv ('ISHAP11T' ,ishap11t ,is_available,lsubmodel)
95 CALL hm_get_floatv ('SIGMA_11C' ,sigma_11c ,is_available,lsubmodel,unitab)
96 CALL hm_get_floatv ('G_11C' ,g_11c ,is_available,lsubmodel,unitab)
97 CALL hm_get_intv ('ISHAP11C' ,ishap11c ,is_available,lsubmodel)
98c
99 ! 3rd card: direction 22 failure parameters
100 CALL hm_get_floatv ('SIGMA_22T' ,sigma_22t ,is_available,lsubmodel,unitab)
101 CALL hm_get_floatv ('G_22T' ,g_22t ,is_available,lsubmodel,unitab)
102 CALL hm_get_intv ('ishap22t' ,ISHAP22T ,IS_AVAILABLE,LSUBMODEL)
103 CALL HM_GET_FLOATV ('sigma_22c' ,SIGMA_22C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
104 CALL HM_GET_FLOATV ('g_22c' ,G_22C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
105 CALL HM_GET_INTV ('ishap22c' ,ISHAP22C ,IS_AVAILABLE,LSUBMODEL)
106c
107 ! 4th card: direction 33 failure parameters
108 CALL HM_GET_FLOATV ('sigma_33t' ,SIGMA_33T ,IS_AVAILABLE,LSUBMODEL,UNITAB)
109 CALL HM_GET_FLOATV ('g_33t' ,G_33T ,IS_AVAILABLE,LSUBMODEL,UNITAB)
110 CALL HM_GET_INTV ('ishap33t' ,ISHAP33T ,IS_AVAILABLE,LSUBMODEL)
111 CALL HM_GET_FLOATV ('sigma_33c' ,SIGMA_33C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
112 CALL HM_GET_FLOATV ('g_33c' ,G_33C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
113 CALL HM_GET_INTV ('ishap33c' ,ISHAP33C ,IS_AVAILABLE,LSUBMODEL)
114c
115 ! 5th card: shear in plane 12 failure parameters
116 CALL HM_GET_FLOATV ('sigma_12t' ,SIGMA_12T ,IS_AVAILABLE,LSUBMODEL,UNITAB)
117 CALL HM_GET_FLOATV ('g_12t' ,G_12T ,IS_AVAILABLE,LSUBMODEL,UNITAB)
118 CALL HM_GET_INTV ('ishap12t' ,ISHAP12T ,IS_AVAILABLE,LSUBMODEL)
119 CALL HM_GET_FLOATV ('sigma_12c' ,SIGMA_12C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
120 CALL HM_GET_FLOATV ('g_12c' ,G_12C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
121 CALL HM_GET_INTV ('ishap12c' ,ISHAP12C ,IS_AVAILABLE,LSUBMODEL)
122c
123 ! 6th card: shear in plane 23 failure parameters
124 CALL HM_GET_FLOATV ('sigma_23t' ,SIGMA_23T ,IS_AVAILABLE,LSUBMODEL,UNITAB)
125 CALL HM_GET_FLOATV ('g_23t' ,G_23T ,IS_AVAILABLE,LSUBMODEL,UNITAB)
126 CALL HM_GET_INTV ('ishap23t' ,ISHAP23T ,IS_AVAILABLE,LSUBMODEL)
127 CALL HM_GET_FLOATV ('sigma_23c' ,SIGMA_23C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
128 CALL HM_GET_FLOATV ('g_23c' ,G_23C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
129 CALL HM_GET_INTV ('ishap23c' ,ISHAP23C ,IS_AVAILABLE,LSUBMODEL)
130c
131 ! 7th card: shear in plane 31 failure parameters
132 CALL HM_GET_FLOATV ('sigma_31t' ,SIGMA_31T ,IS_AVAILABLE,LSUBMODEL,UNITAB)
133 CALL HM_GET_FLOATV ('g_31t' ,G_31T ,IS_AVAILABLE,LSUBMODEL,UNITAB)
134 CALL HM_GET_INTV ('ishap31t' ,ISHAP31T ,IS_AVAILABLE,LSUBMODEL)
135 CALL HM_GET_FLOATV ('sigma_31c' ,SIGMA_31C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
136 CALL HM_GET_FLOATV ('g_31c' ,G_31C ,IS_AVAILABLE,LSUBMODEL,UNITAB)
137 CALL HM_GET_INTV ('ishap31c' ,ISHAP31C ,IS_AVAILABLE,LSUBMODEL)
138c
139 !==========================================================================
140 ! Check parameters and set defaults value
141 !==========================================================================
142 ! Percentage of failed thickness/plys/layers/ratio
143 IF (PTHICKFAIL == ZERO) PTHICKFAIL = ONE-EM06
144 PTHICKFAIL = MIN(PTHICKFAIL, ONE)
145 PTHICKFAIL = MAX(PTHICKFAIL,-ONE)
146 ! Number of failed modes prior to integration point failure
147 IF (NMOD == 0) NMOD = 1
148 NMOD = MIN(NMOD,12)
149 ! Number of failed integration points prior to solid deletion
150 IF (FAILIP == 0) FAILIP = 1
151 ! Check if value are positive
152 SIGMA_11T = ABS(SIGMA_11T)
153 SIGMA_22T = ABS(SIGMA_22T)
154 SIGMA_12T = ABS(SIGMA_12T)
155 SIGMA_11C = ABS(SIGMA_11C)
156 SIGMA_22C = ABS(SIGMA_22C)
157 SIGMA_12C = ABS(SIGMA_12C)
158 ! Critical stresses in tension
159 IF (SIGMA_11T == ZERO) SIGMA_11T = INFINITY
160 IF (SIGMA_22T == ZERO) SIGMA_22T = INFINITY
161 IF (SIGMA_33T == ZERO) SIGMA_33T = INFINITY
162 IF (SIGMA_12T == ZERO) SIGMA_12T = INFINITY
163 IF (SIGMA_23T == ZERO) SIGMA_23T = INFINITY
164 IF (SIGMA_31T == ZERO) SIGMA_31T = INFINITY
165 ! Critical stresses in compression
166 IF (SIGMA_11C == ZERO) SIGMA_11C = SIGMA_11T
167 IF (SIGMA_22C == ZERO) SIGMA_22C = SIGMA_22T
168 IF (SIGMA_33C == ZERO) SIGMA_33C = SIGMA_33T
169 IF (SIGMA_12C == ZERO) SIGMA_12C = SIGMA_12T
170 IF (SIGMA_23C == ZERO) SIGMA_23C = SIGMA_23T
171 IF (SIGMA_31C == ZERO) SIGMA_31C = SIGMA_31T
172 ! Fracture energy in tension
173 IF (G_11T == ZERO) G_11T = INFINITY
174 IF (G_22T == ZERO) G_22T = INFINITY
175 IF (G_33T == ZERO) G_33T = INFINITY
176 IF (G_12T == ZERO) G_12T = INFINITY
177 IF (G_23T == ZERO) G_23T = INFINITY
178 IF (G_31T == ZERO) G_31T = INFINITY
179 ! Fracture energy in compression
180 IF (G_11C == ZERO) G_11C = G_11T
181 IF (G_22C == ZERO) G_22C = G_22T
182 IF (G_33C == ZERO) G_33C = G_33T
183 IF (G_12C == ZERO) G_12C = G_12T
184 IF (G_23C == ZERO) G_23C = G_23T
185 IF (G_31C == ZERO) G_31C = G_31T
186 ! Softening shape flags
187 IF (ISHAP11T == ZERO) ISHAP11T = 1
188 ISHAP11T = MIN(ISHAP11T,2)
189 ISHAP11T = MAX(ISHAP11T,1)
190 IF (ISHAP11C == ZERO) ISHAP11C = 1
191 ISHAP11C = MIN(ISHAP11C,2)
192 ISHAP11C = MAX(ISHAP11C,1)
193 IF (ISHAP22T == ZERO) ISHAP22T = 1
194 ISHAP22T = MIN(ISHAP22T,2)
195 ISHAP22T = MAX(ISHAP22T,1)
196 IF (ISHAP22C == ZERO) ISHAP22C = 1
197 ISHAP22C = MIN(ISHAP22C,2)
198 ISHAP22C = MAX(ISHAP22C,1)
199 IF (ISHAP33T == ZERO) ISHAP33T = 1
200 ISHAP33T = MIN(ISHAP33T,2)
201 ISHAP33T = MAX(ISHAP33T,1)
202 IF (ISHAP33C == ZERO) ISHAP33C = 1
203 ISHAP33C = MIN(ISHAP33C,2)
204 ISHAP33C = MAX(ISHAP33C,1)
205 IF (ISHAP12T == ZERO) ISHAP12T = 1
206 ISHAP12T = MIN(ISHAP12T,2)
207 ISHAP12T = MAX(ISHAP12T,1)
208 IF (ISHAP12C == ZERO) ISHAP12C = 1
209 ISHAP12C = MIN(ISHAP12C,2)
210 ISHAP12C = MAX(ISHAP12C,1)
211 IF (ISHAP23T == ZERO) ISHAP23T = 1
212 ISHAP23T = MIN(ISHAP23T,2)
213 ISHAP23T = MAX(ISHAP23T,1)
214 IF (ISHAP23C == ZERO) ISHAP23C = 1
215 ISHAP23C = MIN(ISHAP23C,2)
216 ISHAP23C = MAX(ISHAP23C,1)
217 IF (ISHAP31T == ZERO) ISHAP31T = 1
218 ISHAP31T = MIN(ISHAP31T,2)
219 ISHAP31T = MAX(ISHAP31T,1)
220 IF (ISHAP31C == ZERO) ISHAP31C = 1
221 ISHAP31C = MIN(ISHAP31C,2)
222 ISHAP31C = MAX(ISHAP31C,1)
223c
224 !==========================================================================
225 ! Save parameters and set buffer tables size
226 !==========================================================================
227 FAIL%KEYWORD = 'orthenerg'
228 FAIL%IRUPT = IRUPT
229 FAIL%FAIL_ID = FAIL_ID
230 FAIL%NUPARAM = 39
231 FAIL%NIPARAM = 0
232 FAIL%NUVAR = 25
233 FAIL%NFUNC = 0
234 FAIL%NTABLE = 0
235 FAIL%NMOD = 0
236 FAIL%PTHK = PTHICKFAIL
237c
238 ALLOCATE (FAIL%UPARAM(FAIL%NUPARAM))
239 ALLOCATE (FAIL%IPARAM(FAIL%NIPARAM))
240 ALLOCATE (FAIL%IFUNC (FAIL%NFUNC))
241 ALLOCATE (FAIL%TABLE (FAIL%NTABLE))
242c
243 ! Failure criterion parameters
244 FAIL%UPARAM(1) = PTHICKFAIL
245 FAIL%UPARAM(2) = FAILIP
246 FAIL%UPARAM(3) = SIGMA_11T
247 FAIL%UPARAM(4) = SIGMA_11C
248 FAIL%UPARAM(5) = SIGMA_22T
249 FAIL%UPARAM(6) = SIGMA_22C
250 FAIL%UPARAM(7) = SIGMA_33T
251 FAIL%UPARAM(8) = SIGMA_33C
252 FAIL%UPARAM(9) = SIGMA_12T
253 FAIL%UPARAM(10) = SIGMA_12C
254 FAIL%UPARAM(11) = SIGMA_23T
255 FAIL%UPARAM(12) = SIGMA_23C
256 FAIL%UPARAM(13) = SIGMA_31T
257 FAIL%UPARAM(14) = SIGMA_31C
258 FAIL%UPARAM(15) = G_11T
259 FAIL%UPARAM(16) = G_11C
260 FAIL%UPARAM(17) = G_22T
261 FAIL%UPARAM(18) = G_22C
262 FAIL%UPARAM(19) = G_33T
263 FAIL%UPARAM(20) = G_33C
264 FAIL%UPARAM(21) = G_12T
265 FAIL%UPARAM(22) = G_12C
266 FAIL%UPARAM(23) = G_23T
267 FAIL%UPARAM(24) = G_23C
268 FAIL%UPARAM(25) = G_31T
269 FAIL%UPARAM(26) = G_31C
270 FAIL%UPARAM(27) = ISHAP11T
271 FAIL%UPARAM(28) = ISHAP11C
272 FAIL%UPARAM(29) = ISHAP22T
273 FAIL%UPARAM(30) = ISHAP22C
274 FAIL%UPARAM(31) = ISHAP33T
275 FAIL%UPARAM(32) = ISHAP33C
276 FAIL%UPARAM(33) = ISHAP12T
277 FAIL%UPARAM(34) = ISHAP12C
278 FAIL%UPARAM(35) = ISHAP23T
279 FAIL%UPARAM(36) = ISHAP23C
280 FAIL%UPARAM(37) = ISHAP31T
281 FAIL%UPARAM(38) = ISHAP31C
282 FAIL%UPARAM(39) = NMOD
283c
284 !==========================================================================
285 ! Print out parameters
286 !==========================================================================
287 IF (IS_ENCRYPTED) THEN
288 WRITE(IOUT,'(5x,a,//)')'confidential data'
289 ELSE
290 WRITE(IOUT,1000)
291 WRITE(IOUT,1100) PTHICKFAIL,FAILIP,NMOD
292 WRITE(IOUT,1200) SIGMA_11T,G_11T,ISHAP11T,SIGMA_11C,G_11C,ISHAP11C
293 WRITE(IOUT,1300) SIGMA_22T,G_22T,ISHAP22T,SIGMA_22C,G_22C,ISHAP22C
294 WRITE(IOUT,1400) SIGMA_33T,G_33T,ISHAP33T,SIGMA_33C,G_33C,ISHAP33C
295 WRITE(IOUT,1500) SIGMA_12T,G_12T,ISHAP12T,SIGMA_12C,G_12C,ISHAP12C
296 WRITE(IOUT,1600) SIGMA_23T,G_23T,ISHAP23T,SIGMA_23C,G_23C,ISHAP23C
297 WRITE(IOUT,1700) SIGMA_31T,G_31T,ISHAP31T,SIGMA_31C,G_31C,ISHAP31C
298 ENDIF
299c-----------------------------------------------------
3001000 FORMAT(
301 & 5X,'--------------------------------------------------'/,
302 & 5X,' energy based orthotropic failure model '/,
303 & 5X,'--------------------------------------------------'/)
3041100 FORMAT(
305 & 5X,'element deletion parameters: '/,
306 & 5X,' '/,
307 & 5X,'shell element deletion PARAMETER pthickfail. . . .=',1PG20.13,/,
308 & 5X,' > 0.0 : fraction of failed thickness ',/,
309 & 5X,' < 0.0 : fraction of failed intg. points or layers',/,
310 & 5X,'number of failed intg. pts prior to solid deletion=',I5/,
311 & 5X,'number of failed modes prior to itg. point failure=',I5/)
3121200 FORMAT(
313 & 5X,'direction 11 parameters: '/,
314 & 5X,' '/,
315 & 5x,'CRITICAL STRESS IN TENSION SIGMA_11T . . . . . . .=',1pg20.13,/,
316 & 5x,'FRACTURE ENERGY IN TENSION G_11T . . . . . . . . .=',1pg20.13,/,
317 & 5x,'STRESS SOFTENING FLAG IN TENSION ISHAP11T. . . . .=',i3/,
318 & 5x,' ISHAP11T = 1: LINEAR (DEFAULT) ',/,
319 & 5x,' ISHAP11T = 2: EXPONENTIAL ',/,
320 & 5x,'CRITICAL STRESS IN COMPRESSION SIGMA_11C . . . . .=',1pg20.13,/,
321 & 5x,'FRACTURE ENERGY IN COMPRESSION G_11C . . . . . . .=',1pg20.13,/,
322 & 5x,'STRESS SOFTENING FLAG IN COMPRESSION ISHAP11C. . .=',i3/,
323 & 5x,' ISHAP11C = 1: LINEAR (DEFAULT) ',/,
324 & 5x,' ISHAP11C = 2: EXPONENTIAL ',/)
3251300 FORMAT(
326 & 5x,'DIRECTION 22 PARAMETERS: '/,
327 & 5x,' '/,
328 & 5x,'CRITICAL STRESS IN TENSION SIGMA_22T . . . . . . .=',1pg20.13,/,
329 & 5x,'FRACTURE ENERGY IN TENSION G_22T . . . . . . . . .=',1pg20.13,/,
330 & 5x,'STRESS SOFTENING FLAG IN TENSION ISHAP22T. . . . .=',i3/,
331 & 5x,' ISHAP22T = 1: LINEAR (DEFAULT) ',/,
332 & 5x,' ISHAP22T = 2: EXPONENTIAL ',/,
333 & 5x,'CRITICAL STRESS IN COMPRESSION SIGMA_22C . . . . .=',1pg20.13,/,
334 & 5x,'FRACTURE ENERGY IN COMPRESSION G_22C . . . . . . .=',1pg20.13,/,
335 & 5x,'STRESS SOFTENING FLAG IN COMPRESSION ISHAP22C. . .=',i3/,
336 & 5x,' ISHAP22C = 1: LINEAR (DEFAULT) ',/,
337 & 5x,' ISHAP22C = 2: EXPONENTIAL ',/)
3381400 FORMAT(
339 & 5x,'DIRECTION 33 PARAMETERS: '/,
340 & 5x,' '/,
341 & 5x,'CRITICAL STRESS IN TENSION SIGMA_33T . . . . . . .=',1pg20.13,/,
342 & 5x,'FRACTURE ENERGY IN TENSION G_33T . . . . . . . . .=',1pg20.13,/,
343 & 5x,'STRESS SOFTENING FLAG IN TENSION ISHAP33T. . . . .=',i3/,
344 & 5x,' ISHAP33T = 1: LINEAR (DEFAULT) ',/,
345 & 5x,' ISHAP33T = 2: EXPONENTIAL ',/,
346 & 5x,'CRITICAL STRESS IN COMPRESSION SIGMA_33C . . . . .=',1pg20.13,/,
347 & 5x,'FRACTURE ENERGY IN COMPRESSION G_33C . . . . . . .=',1pg20.13,/,
348 & 5x,'STRESS SOFTENING FLAG IN COMPRESSION ISHAP33C. . .=',i3/,
349 & 5x,' ISHAP33C = 1: LINEAR (DEFAULT) ',/,
350 & 5x,' ISHAP33C = 2: EXPONENTIAL ',/)
3511500 FORMAT(
352 & 5x,'PLANE 12 PARAMETERS: '/,
353 & 5x,' '/,
354 & 5x,'CRITICAL STRESS IN POSITIVE SHEAR SIGMA_12T. . . .=',1pg20.13,/,
355 & 5x,'FRACTURE ENERGY IN POSITIVE SHEAR G_12T. . . . . .=',1pg20.13,/,
356 & 5x,'STRESS SOFTENING FLAG IN POSITIVE SHEAR ISHAP12T .=',i3/,
357 & 5x,' ISHAP12T = 1: LINEAR (DEFAULT) ',/,
358 & 5x,' ISHAP12T = 2: EXPONENTIAL ',/,
359 & 5x,'CRITICAL STRESS IN NEGATIVE SHEAR SIGMA_12C. . . .=',1pg20.13,/,
360 & 5x,'FRACTURE ENERGY IN NEGATIVE SHEAR G_12C. . . . . .=',1pg20.13,/,
361 & 5x,'STRESS SOFTENING FLAG IN NEGATIVE SHEAR ISHAP12C .=',i3/,
362 & 5x,' ISHAP12C = 1: LINEAR (DEFAULT) ',/,
363 & 5x,' ISHAP12C = 2: EXPONENTIAL ',/)
3641600 FORMAT(
365 & 5x,'PLANE 23 PARAMETERS: '/,
366 & 5x,' '/,
367 & 5x,'CRITICAL STRESS IN POSITIVE SHEAR SIGMA_23T. . . .=',1pg20.13,/,
368 & 5x,'FRACTURE ENERGY IN POSITIVE SHEAR G_23T. . . . . .=',1pg20.13,/,
369 & 5x,'STRESS SOFTENING FLAG IN POSITIVE SHEAR ISHAP23T .=',i3/,
370 & 5x,' ISHAP23T = 1: LINEAR (DEFAULT) ',/,
371 & 5x,' ISHAP23T = 2: EXPONENTIAL ',/,
372 & 5x,'CRITICAL STRESS IN NEGATIVE SHEAR SIGMA_23C. . . .=',1pg20.13,/,
373 & 5x,'FRACTURE ENERGY IN NEGATIVE SHEAR G_23C. . . . . .=',1pg20.13,/,
374 & 5x,'STRESS SOFTENING FLAG IN NEGATIVE SHEAR ISHAP23C .=',i3/,
375 & 5x,' ISHAP23C = 1: LINEAR (DEFAULT) ',/,
376 & 5x,' ISHAP23C = 2: EXPONENTIAL ',/)
3771700 FORMAT(
378 & 5x,'PLANE 31 PARAMETERS: '/,
379 & 5x,' '/,
380 & 5x,'CRITICAL STRESS IN POSITIVE SHEAR SIGMA_31T. . . .=',1pg20.13,/,
381 & 5x,'FRACTURE ENERGY IN POSITIVE SHEAR G_31T. . . . . .=',1pg20.13,/,
382 & 5x,'STRESS SOFTENING FLAG IN POSITIVE SHEAR ISHAP31T .=',i3/,
383 & 5x,' ISHAP31T = 1: LINEAR (DEFAULT) ',/,
384 & 5x,' ISHAP31T = 2: EXPONENTIAL ',/,
385 & 5x,'CRITICAL STRESS IN NEGATIVE SHEAR SIGMA_31C. . . .=',1pg20.13,/,
386 & 5x,'FRACTURE ENERGY IN NEGATIVE SHEAR G_31C. . . . . .=',1pg20.13,/,
387 & 5x,'STRESS SOFTENING FLAG IN NEGATIVE SHEAR ISHAP31C .=',i3/,
388 & 5x,' ISHAP31C = 1: LINEAR (DEFAULT) ',/,
389 & 5x,' ISHAP31C = 2: EXPONENTIAL ',/)
my_real
#define my_real
Definition cppsort.cpp:32
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:415
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:284
hm_option_read_mod
Definition hm_option_read_mod.F:408
message_mod
Definition message_mod.F:1249
submodel_mod
Definition submodel_mod.F:712
unitab_mod
Definition unitab_mod.F:513