OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
tshgeodel3.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/.
23C-----------------------------------------------
24!||====================================================================
25!|| tshgeodel3 ../engine/source/elements/thickshell/solidec/tshgeodel3.F
26!||--- called by ------------------------------------------------------
27!|| s16forc3 ../engine/source/elements/thickshell/solide16/s16forc3.F
28!|| s6cforc3 ../engine/source/elements/thickshell/solide6c/s6cforc3.F
29!|| s8cforc3 ../engine/source/elements/thickshell/solide8c/s8cforc3.F
30!|| scforc3 ../engine/source/elements/thickshell/solidec/scforc3.F
31!||--- uses -----------------------------------------------------
32!|| dt_mod ../engine/source/modules/dt_mod.F
33!|| message_mod ../engine/share/message_module/message_mod.F
34!||====================================================================
35 SUBROUTINE tshgeodel3(NGL, OFFG ,VOLG ,AREA, VOLG0 ,L_SH,GEO ,NNOD,DT,NEL)
36C-----------------------------------------------
37C M o d u l e s
38C-----------------------------------------------
39 USE message_mod
40 USE dt_mod
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45#include "mvsiz_p.inc"
46#include "comlock.inc"
47C-----------------------------------------------
48C C o m m o n B l o c k s
49C-----------------------------------------------
50#include "units_c.inc"
51#include "param_c.inc"
52#include "scr17_c.inc"
53C-----------------------------------------------
54C D u m m y A r g u m e n t s
55C-----------------------------------------------
56 INTEGER ,INTENT(IN) :: NNOD,NEL
57 INTEGER ,DIMENSION(MVSIZ) ,INTENT(IN) :: NGL
58 my_real, DIMENSION(NPROPG) , INTENT(IN) :: geo
59 my_real, DIMENSION(NEL) , INTENT(IN) :: volg0
60 my_real, DIMENSION(MVSIZ) , INTENT(IN) :: volg,area,l_sh ! L_SH=L_MAX**2
61 my_real, DIMENSION(NEL) , INTENT(INOUT) :: offg
62 TYPE(dt_), INTENT(INOUT) :: DT
63C-----------------------------------------------
64C L o c a l V a r i a b l e s
65C-----------------------------------------------
66 INTEGER I,J,ICOUNT,LIST(MVSIZ)
67 my_real
68 . tet_colla(mvsiz),def_v(mvsiz),aps_r(mvsiz),
69 . c_min,vdefmin,vdefmax,aspmax,colmin,c_max,gdel
70C---------------------------------------------
71 vdefmin = geo(190)
72 vdefmax = geo(191)
73 aspmax = geo(192)
74 colmin = geo(193)
75C-----
76 gdel=vdefmin+vdefmax+aspmax+colmin
77 IF(gdel==zero.AND.dt%IDEL_BRICK==0) RETURN
78C-----
79 IF(dt%BRICK_DEL_DEFV_MAX>zero .AND. (dt%BRICK_DEL_DEFV_MAX<vdefmax
80 . .OR.vdefmax==zero) ) vdefmax = dt%BRICK_DEL_DEFV_MAX
81 IF(dt%BRICK_DEL_ASP_MAX>zero .AND. ( dt%BRICK_DEL_ASP_MAX<aspmax
82 . .OR.aspmax==zero) ) aspmax = dt%BRICK_DEL_ASP_MAX
83 c_min = one
84 c_max = one
85C--- penta6
86 IF (nnod ==6) THEN
87 c_max = fourth*sqr3
88 c_min = one/c_max
89 END IF
90C
91 IF ((vdefmin+vdefmax) >zero) THEN
92 DO i=1,nel
93 IF (offg(i)==zero) cycle
94 def_v(i) = volg(i)/volg0(i)
95 END DO
96 END IF
97 IF (vdefmin >zero) THEN
98 icount = 0
99 DO i=1,nel
100 IF (offg(i)==zero) cycle
101 IF (def_v(i) < vdefmin) THEN
102 offg(i) = zero
103 icount=icount+1
104 list(icount)=i
105 END IF
106 END DO
107 IF (icount > 0) THEN
108#include "lockon.inc"
109 WRITE(iout,1001)
110 WRITE(istdo,1001)
111 DO j=1,icount
112 i=list(j)
113 WRITE(iout,'(6X,I10,1PG20.13)')ngl(i),def_v(i)
114 WRITE(istdo,'(6X,I10,1PG20.13)')ngl(i),def_v(i)
115 ENDDO
116 WRITE(iout,*)
117 WRITE(istdo,*)
118#include "lockoff.inc"
119 idel7nok = 1
120 END IF !(ICOUNT > 0) THEN
121 END IF
122C
123 IF (vdefmax >zero) THEN
124 icount = 0
125 DO i=1,nel
126 IF (offg(i)==zero) cycle
127 IF (def_v(i) > vdefmax) THEN
128 offg(i) = zero
129 icount=icount+1
130 list(icount)=i
131 END IF
132 END DO
133 IF (icount > 0) THEN
134#include "lockon.inc"
135 WRITE(iout,1002)
136 WRITE(istdo,1002)
137 DO j=1,icount
138 i=list(j)
139 WRITE(iout,'(6X,I10,1PG20.13)')ngl(i),def_v(i)
140 WRITE(istdo,'(6X,I10,1PG20.13)')ngl(i),def_v(i)
141 ENDDO
142 WRITE(iout,*)
143 WRITE(istdo,*)
144#include "lockoff.inc"
145 idel7nok = 1
146 END IF !(ICOUNT > 0) THEN
147 END IF
148C--- Hexa : APS_r=lmax/lmin lmax~sqrt(L_SH); lmin~A/lmax ->
149C--- APS_r~lmax^2/A
150 IF (aspmax>zero) THEN
151 icount = 0
152 DO i=1,nel
153 IF (offg(i)==zero) cycle
154 aps_r(i) = c_max*l_sh(i)/area(i)
155 IF (aps_r(i) > aspmax) THEN
156 offg(i) = zero
157 icount=icount+1
158 list(icount)=i
159 END IF
160 END DO
161 IF (icount > 0) THEN
162#include "lockon.inc"
163 WRITE(iout,1003)
164 WRITE(istdo,1003)
165 DO j=1,icount
166 i=list(j)
167 WRITE(iout,'(6X,I10,1PG20.13)')ngl(i),aps_r(i)
168 WRITE(istdo,'(6X,I10,1PG20.13)')ngl(i),aps_r(i)
169 ENDDO
170 WRITE(iout,*)
171 WRITE(istdo,*)
172#include "lockoff.inc"
173 idel7nok = 1
174 END IF !(ICOUNT > 0) THEN
175 END IF
176C
177C--- TET_c =1/ASP_r
178 IF (colmin>zero) THEN
179 icount = 0
180 DO i=1,nel
181 IF (offg(i)==zero) cycle
182 tet_colla(i) = c_min *area(i)/l_sh(i)
183 IF (tet_colla(i) < colmin) THEN
184 offg(i) = zero
185 icount=icount+1
186 list(icount)=i
187 END IF
188 END DO
189 IF (icount > 0) THEN
190#include "lockon.inc"
191 WRITE(iout,1004)
192 WRITE(istdo,1004)
193 DO j=1,icount
194 i=list(j)
195 WRITE(iout,'(6X,I10,1PG20.13)')ngl(i),tet_colla(i)
196 WRITE(istdo,'(6X,I10,1PG20.13)')ngl(i),tet_colla(i)
197 ENDDO
198 WRITE(iout,*)
199 WRITE(istdo,*)
200#include "lockoff.inc"
201 idel7nok = 1
202 END IF !(ICOUNT > 0) THEN
203 END IF
204C
205 RETURN
206 1001 FORMAT(/
207 . 'THE FOLLOWING THICK-SHELL ELEMENTS WILL BE DELETED BY MINIMUM VOLUME RATIO' /)
208 1002 FORMAT(/
209 . 'THE FOLLOWING THICK-SHELL ELEMENTS WILL BE DELETED BY MAXIMUM VOLUME RATIO' /)
210 1003 FORMAT(/
211 . 'THE FOLLOWING THICK-SHELL ELEMENTS WILL BE DELETED BY MAXIMUM ASPECT RATIO' /)
212 1004 FORMAT(/
213 . 'THE FOLLOWING THICK-SHELL ELEMENTS WILL BE DELETED BY MINIMUM COLLAPSE RATIO' /)
214 END
my_real
#define my_real
Definition cppsort.cpp:32
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
dt_mod
Definition dt_mod.F:80
message_mod
Definition message_mod.F:1257
dt_mod::dt_
Definition dt_mod.F:86
tshgeodel3
subroutine tshgeodel3(ngl, offg, volg, area, volg0, l_sh, geo, nnod, dt, nel)
Definition tshgeodel3.F:36