OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_fail_tsaiwu.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!|| hm_read_fail_tsaiwu ../starter/source/materials/fail/tsaiwu/hm_read_fail_tsaiwu.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_fail ../starter/source/materials/fail/hm_read_fail.F
27!||--- calls -----------------------------------------------------
28!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
29!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
30!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
31!||--- uses -----------------------------------------------------
32!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
33!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
34!|| message_mod ../starter/share/message_module/message_mod.F
35!|| submodel_mod ../starter/share/modules1/submodel_mod.F
36!||====================================================================
37 SUBROUTINE hm_read_fail_tsaiwu(
38 . FAIL ,MAT_ID ,FAIL_ID ,IRUPT ,
39 . LSUBMODEL,UNITAB ,FAIL_TAG )
40C-----------------------------------------------
41c ROUTINE DESCRIPTION :
42c Read Tsai-Wu failure model parameters
43C-----------------------------------------------
44C M o d u l e s
45C-----------------------------------------------
46 USE fail_param_mod
47 USE unitab_mod
48 USE message_mod
49 USE submodel_mod
50 USE hm_option_read_mod
51 USE elbuftag_mod
52C-----------------------------------------------
53C I m p l i c i t T y p e s
54C-----------------------------------------------
55#include "implicit_f.inc"
56C-----------------------------------------------
57C C o m m o n B l o c k s
58C-----------------------------------------------
59#include "units_c.inc"
60C-----------------------------------------------
61C D u m m y A r g u m e n t s
62C-----------------------------------------------
63 INTEGER ,INTENT(IN) :: FAIL_ID !< failure model ID
64 INTEGER ,INTENT(IN) :: MAT_ID !< material law ID
65 INTEGER ,INTENT(IN) :: IRUPT !< failure model number
66 TYPE (unit_type_) ,INTENT(IN) :: unitab !< table of input units
67 TYPE (SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(*) !< submodel table
68 TYPE (FAIL_PARAM_) ,INTENT(INOUT) :: FAIL !< failure model data structure
69 TYPE (FAIL_TAG_) ,INTENT(INOUT) :: FAIL_TAG !< failure model tag for buffer allocation
70C-----------------------------------------------
71C L o c a l V a r i a b l e s
72C-----------------------------------------------
73 INTEGER :: IFAIL_SH,IFAIL_SO,NUVAR
74 my_real :: SIGT1,SIGT2,SIGT12,SIGC1,SIGC2,ALPHA,TMAX,FCUT,PTHKF,
75 . f1,f2,f11,f22,f66,f12
76 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
77C--------------------------------------------------
78C S o u r c e L i n e s
79C--------------------------------------------------
80 is_encrypted = .false.
81 is_available = .false.
82
83 CALL hm_option_is_encrypted(is_encrypted)
84
85 ! 1st line of the card
86 CALL hm_get_floatv ('SIGMA_1T' ,sigt1 ,is_available,lsubmodel,unitab)
87 CALL hm_get_floatv ('SIGMA_2T' ,sigt2 ,is_available,lsubmodel,unitab)
88 CALL hm_get_floatv ('sigma_12' ,SIGT12 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
89 CALL HM_GET_FLOATV ('sigma_1c' ,SIGC1 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
90 CALL HM_GET_FLOATV ('sigma_2c' ,SIGC2 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
91
92 ! 2nd line of the card
93 CALL HM_GET_FLOATV ('alpha' ,ALPHA ,IS_AVAILABLE,LSUBMODEL,UNITAB)
94 CALL HM_GET_FLOATV ('tau_max' ,TMAX ,IS_AVAILABLE,LSUBMODEL,UNITAB)
95 CALL HM_GET_FLOATV ('fcut' ,FCUT ,IS_AVAILABLE,LSUBMODEL,UNITAB)
96 CALL HM_GET_INTV ('ifail_sh' ,IFAIL_SH ,IS_AVAILABLE,LSUBMODEL)
97 CALL HM_GET_INTV ('ifail_so' ,IFAIL_SO ,IS_AVAILABLE,LSUBMODEL)
98
99 ! Check parameters and set defaults value
100 IF (SIGT1 == ZERO) SIGT1 = INFINITY
101 IF (SIGT2 == ZERO) SIGT2 = INFINITY
102 IF (SIGC1 == ZERO) SIGC1 = INFINITY
103 IF (SIGC2 == ZERO) SIGC2 = INFINITY
104 IF (TMAX == ZERO) TMAX = INFINITY
105 IF (SIGT12 == ZERO) SIGT12 = INFINITY
106 IFAIL_SH = MIN(IFAIL_SH,2)
107 IFAIL_SH = MAX(0,IFAIL_SH)
108 IFAIL_SO = MIN(IFAIL_SO,2)
109 IFAIL_SO = MAX(0,IFAIL_SO)
110 FCUT = MAX(ZERO,FCUT)
111c
112 ! Shell element deletion flag
113 IF (IFAIL_SH == 1) THEN
114 PTHKF = EM06
115 ELSEIF (IFAIL_SH == 2) THEN
116 PTHKF = ONE
117 ENDIF
118c
119 ! Compute Tsai-Wu criterion parameters
120 F1 = ONE/SIGT1-(ONE/SIGC1)
121 F2 = ONE/SIGT2-(ONE/SIGC2)
122 F11 = ONE/(SIGT1*SIGC1)
123 F22 = ONE/(SIGT2*SIGC2)
124 F66 = ONE/(SIGT12**2)
125 F12 = -ALPHA*HALF*SQRT(F11*F22)
126c---------------------------
127 IF (FCUT > ZERO) THEN
128 NUVAR = 7
129 ELSE
130 NUVAR = 1
131 ENDIF
132c---------------------------
133 FAIL%KEYWORD = 'tsai-wu'
134 FAIL%IRUPT = IRUPT
135 FAIL%FAIL_ID = FAIL_ID
136 FAIL%NUPARAM = 10
137 FAIL%NIPARAM = 0
138 FAIL%NUVAR = NUVAR
139 FAIL%NFUNC = 0
140 FAIL%NTABLE = 0
141 FAIL%NMOD = 1
142 FAIL%PTHK = PTHKF
143c
144 ALLOCATE (FAIL%UPARAM(FAIL%NUPARAM))
145 ALLOCATE (FAIL%IPARAM(FAIL%NIPARAM))
146 ALLOCATE (FAIL%IFUNC (FAIL%NFUNC))
147 ALLOCATE (FAIL%TABLE (FAIL%NTABLE))
148c
149 ! Modes of failure
150 FAIL_TAG%LF_DAMMX = FAIL_TAG%LF_DAMMX + FAIL%NMOD
151 ALLOCATE (FAIL%MODE(FAIL%NMOD))
152 FAIL%MODE(1) = "Reserve factor"
153c
154 FAIL%UPARAM(1) = F1
155 FAIL%UPARAM(2) = F2
156 FAIL%UPARAM(3) = F11
157 FAIL%UPARAM(4) = F22
158 FAIL%UPARAM(5) = F66
159 FAIL%UPARAM(6) = F12
160 FAIL%UPARAM(7) = TMAX
161 FAIL%UPARAM(8) = FCUT
162 FAIL%UPARAM(9) = IFAIL_SH
163 FAIL%UPARAM(10) = IFAIL_SO
164c---------------------------
165C
166 ! Print parameter in 0.out file
167 IF (IS_ENCRYPTED) THEN
168 WRITE(IOUT,'(5x,a,//)')'confidential data'
169 ELSE
170 WRITE(IOUT,1000)
171 WRITE(IOUT,1100) SIGT1,SIGT2,SIGC1,SIGC2,SIGT12,ALPHA
172.OR. IF (IFAIL_SH > 0 IFAIL_SO > 0) THEN
173 WRITE(IOUT,1200) TMAX
174 ENDIF
175 IF (FCUT > ZERO) THEN
176 WRITE(IOUT,1300) FCUT
177 ENDIF
178 WRITE(IOUT,1400) IFAIL_SH
179 WRITE(IOUT,1500) IFAIL_SO
180 ENDIF
181c-----------
182 1000 FORMAT(
183 & 5X,' ---------------------------------------------------- ',/
184 & 5X,' failure criterion : tsai-wu ',/,
185 & 5X,' ---------------------------------------------------- ',/)
186 1100 FORMAT(
187 & 5X,'tensile strength in direction 1 sigma_1t . . . . . . . =',1PG20.13/
188 & 5X,'tensile strength in direction 2 sigma_2t . . . . . . . =',1PG20.13/
189 & 5X,'compressive strength in direction 1 sigma_1c . . . . . =',1PG20.13/
190 & 5X,'compressive strength in direction 2 sigma_2c . . . . . =',1PG20.13/
191 & 5X,'shear strength sigma_12 . . . . . . . . . . . . . . . . =',1PG20.13/
192 & 5X,'interaction PARAMETER alpha . . . . . . . . . . . . . . =',1PG20.13/)
193 1200 FORMAT(
194 & 5X,'stress softening activated: ',/,
195 & 5x,'RELAXATION TIME TAU_MAX . . . . . . . . . . . . . . . . =',1pg20.13/)
196 1300 FORMAT(
197 & 5x,'STRESS TENSOR FILTERING ACTIVATED: ',/,
198 & 5x,'FILTERING FREQUENCY FCUT . . . . . . . . . . . . . . . =',1pg20.13/)
199 1400 FORMAT(
200 & 5x,'SHELL DELETION FLAG IFAIL_SH . . . . . . . . . . . . . =',i10,/
201 & 5x,' = 0: SHELL NEVER DELETED AND NO STRESS SOFTENING ',/,
202 & 5x,' = 1: shell deleted IF damage is reached for one layer ',/,
203 & 5X,' = 2: shell deleted IF damage is reached for all layers',/)
204 1500 FORMAT(
205 & 5X,'solid deletion flag ifail_so . . . . . . . . . . . . . =',I10,/
206 & 5X,' = 0: solid never deleted and no stress softening ',/,
207 & 5X,' = 1: solid deleted IF damage is reached for one integration point ',/,
208 & 5X,' = 2: solid deleted IF damage is reached for all integration points',/)
209c-----------
210 RETURN
211 END
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:288
hm_read_fail_tsaiwu
subroutine hm_read_fail_tsaiwu(fail, mat_id, fail_id, irupt, lsubmodel, unitab, fail_tag)
Definition hm_read_fail_tsaiwu.F:40
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
elbuftag_mod
Definition elbuftag_mod.F:230
hm_option_read_mod
Definition hm_option_read_mod.F:412
message_mod
Definition message_mod.F:1257
submodel_mod
Definition submodel_mod.F:717
unitab_mod
Definition unitab_mod.F:517
unitab_mod::unit_type_
Definition unitab_mod.F:524