OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
check_pthickfail.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!|| check_pthickfail ../starter/source/materials/fail/check_pthickfail.F
25!||--- called by ------------------------------------------------------
26!|| lectur ../starter/source/starter/lectur.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!||--- uses -----------------------------------------------------
30!|| message_mod ../starter/share/message_module/message_mod.F
31!|| stack_mod ../starter/share/modules1/stack_mod.F
32!||====================================================================
33 SUBROUTINE check_pthickfail(ELBUF_STR,MAT_PARAM,IPARG ,GEO ,
34 . IPM ,STACK ,IGEO ,
35 . NUMMAT ,NUMGEO ,NGROUP ,NPARG ,
36 . NPROPMI ,NPROPGI ,NPROPG )
37C-----------------------------------------------
38C M o d u l e s
39C-----------------------------------------------
40 USE mat_elem_mod
41 USE stack_mod
42 USE message_mod
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C D u m m y A r g u m e n t s
49C-----------------------------------------------
50 INTEGER, INTENT(IN) :: IPARG(NPARG,NGROUP),NPROPG,
51 . IPM(NPROPMI,NUMMAT),IGEO(NPROPGI,NUMGEO),
52 . NUMMAT,NUMGEO,NGROUP,NPARG,NPROPMI,NPROPGI
53 my_real, INTENT(IN) :: geo(npropg,numgeo)
54 TYPE (ELBUF_STRUCT_) ,DIMENSION(NGROUP) :: ELBUF_STR
55 TYPE (STACK_PLY) :: STACK
56 TYPE (MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(INOUT) :: MAT_PARAM
57C-----------------------------------------------
58C L o c a l V a r i a b l e s
59C-----------------------------------------------
60 INTEGER :: NG,ITY,IGTYP,PID,NLAY,IFL,NFAIL,NPTTOT,IL,
61 . IMAT,ISUBSTACK,IPTHKLY,PROP_ID,MAT_ID
62 my_real :: P_THICKG,P_THICKFAIL,PTHKLY
63c
64 !===============================================================================
65 ! CHECK PTHICKFAIL VALUE FOR SHELL PROPERTIES
66 !===============================================================================
67 ! Loop over groups
68 DO ng = 1,ngroup
69c
70 ! Group parameters
71 igtyp = iparg(38,ng) ! -> Property type
72 pid = iparg(62,ng) ! -> Property internal ID
73 isubstack = iparg(71,ng) ! -> Substack number
74c
75 ! Check group property
76 IF ((igtyp == 1).OR.(igtyp == 9).OR.(igtyp == 10).OR.(igtyp == 11).OR.
77 . (igtyp == 16).OR.(igtyp == 17).OR.(igtyp == 51).OR.(igtyp == 52)) THEN
78c
79 ! Property parameters
80 p_thickg = geo(42,pid) ! -> Global Pthickfail parameter of the property
81 prop_id = igeo(1,pid) ! -> Property user ID
82 nlay = elbuf_str(ng)%NLAY ! -> Number of layers
83 ipthkly = 1 + 4*nlay ! -> Address of PTHKFAIL for PLYs
84 npttot = 0
85 DO il = 1,nlay
86 npttot = npttot + elbuf_str(ng)%BUFLY(il)%NPTT
87 ENDDO
88c
89 ! Loop over layers
90 DO il = 1,nlay
91c
92 ! Layer parameters
93 nfail = elbuf_str(ng)%BUFLY(il)%NFAIL ! -> Number of failure criteria
94 imat = elbuf_str(ng)%BUFLY(il)%IMAT ! -> Material internal ID
95 IF ((igtyp == 51).OR.(igtyp == 52)) THEN
96 pthkly = stack%GEO(ipthkly + il,isubstack) ! -> Pthickfail parameter for ply
97 ENDIF
98 mat_id = ipm(1,imat) ! -> Material user ID
99c
100 ! Loop over failure criteria
101 DO ifl = 1,nfail
102c
103 ! Recover criterion PTHKFAIL value
104 p_thickfail = mat_param(imat)%FAIL(ifl)%PTHK
105c
106 ! Check conflicts between failure criteria and property value
107 ! -> For 1 layer properties, check between global and failure criterion value
108 IF (nlay == 1) THEN
109 IF (((p_thickfail > zero).AND.(p_thickg < zero)).OR.
110 . ((p_thickfail < zero).AND.(p_thickg > zero))) THEN
111 CALL ancmsg(msgid=3029,
112 . msgtype=msgwarning,
113 . anmode=aninfo_blind_1,
114 . i1=prop_id,
115 . r1=p_thickg,
116 . i2=mat_id,
117 . r2=p_thickfail,
118 . prmod=msg_cumu)
119 ENDIF
120 ! -> For multiple layer properties with 1 intg. point, check for failure criteria value
121 ELSEIF (nlay == npttot) THEN
122 IF (p_thickfail /= zero) THEN
123 CALL ancmsg(msgid=3030,
124 . msgtype=msgwarning,
125 . anmode=aninfo_blind_1,
126 . i1=prop_id,
127 . i2=igtyp,
128 . prmod=msg_cumu)
129 ENDIF
130 ! -> For multiple layer properties with several intg. points, check between ply and failure criterion value
131 ELSE
132 IF (((p_thickfail > zero).AND.(pthkly < zero)).OR.
133 . ((p_thickfail < zero).AND.(pthkly > zero))) THEN
134 CALL ancmsg(msgid=3031,
135 . msgtype=msgwarning,
136 . anmode=aninfo_blind_1,
137 . i1=prop_id,
138 . i2=il,
139 . r1=pthkly,
140 . i3=mat_id,
141 . r2=p_thickfail,
142 . prmod=msg_cumu)
143 ENDIF
144 ENDIF
145 ENDDO
146 ENDDO
147 ENDIF
148 ENDDO
149c
150 ! Printing out warning messages
151 CALL ancmsg(msgid=3029,
152 . msgtype=msgwarning,
153 . anmode=aninfo_blind_1,
154 . prmod=msg_print)
155 CALL ancmsg(msgid=3030,
156 . msgtype=msgwarning,
157 . anmode=aninfo_blind_1,
158 . prmod=msg_print)
159 CALL ancmsg(msgid=3031,
160 . msgtype=msgwarning,
161 . anmode=aninfo_blind_1,
162 . prmod=msg_print)
163 !===============================================================================
164c
165 END
check_pthickfail
subroutine check_pthickfail(elbuf_str, mat_param, iparg, geo, ipm, stack, igeo, nummat, numgeo, ngroup, nparg, npropmi, npropgi, npropg)
Definition check_pthickfail.F:37
my_real
#define my_real
Definition cppsort.cpp:32
message_mod
Definition message_mod.F:1249
stack_mod
Definition stack_mod.F:72
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889