OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
desacti.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!|| desacti ../engine/source/elements/desacti.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- calls -----------------------------------------------------
28!|| convecoff ../engine/source/constraints/thermic/convecoff.F
29!|| eloff ../engine/source/elements/eloff.F
30!|| radiatoff ../engine/source/constraints/thermic/radiatoff.F
31!||--- uses -----------------------------------------------------
32!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
33!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
34!|| groupdef_mod ../common_source/modules/groupdef_mod.F
35!|| sensor_mod ../common_source/modules/sensor_mod.F90
36!||====================================================================
37 SUBROUTINE desacti(IXS ,IXQ ,IXC ,IXP ,IXT ,
38 . IXR ,IXTG ,IPARG ,IACTIV ,
39 . NSENSOR ,SENSOR_TAB,FSKY ,X ,ELBUF_TAB,
40 . IBCV ,FCONV ,IBCR ,FRADIA ,IGROUPS ,
41 . FACTIV ,TEMP ,MCP ,PM ,MCP_OFF ,
42 . IGRBRIC ,IGRQUAD ,IGRSH4N ,IGRSH3N ,IGRTRUSS ,
43 . IGRBEAM ,IGRSPRING,GLOB_THERM)
44C-----------------------------------------------
45C M o d u l e s
46C-----------------------------------------------
47 USE elbufdef_mod
48 USE groupdef_mod
49 USE sensor_mod
50 use glob_therm_mod
51C-----------------------------------------------
52C I m p l i c i t T y p e s
53C-----------------------------------------------
54#include "implicit_f.inc"
55C-----------------------------------------------
56C C o m m o n B l o c k s
57C-----------------------------------------------
58#include "param_c.inc"
59#include "com01_c.inc"
60#include "com04_c.inc"
61#include "com08_c.inc"
62#include "parit_c.inc"
63C-----------------------------------------------
64C D u m m y A r g u m e n t s
65C-----------------------------------------------
66 INTEGER ,INTENT(IN) :: NSENSOR
67 INTEGER IACTIV(LACTIV,*),IPARG(NPARG,*),
68 . IXS(NIXS,*), IXQ(NIXQ,*), IXC(NIXC,*),
69 . IXT(NIXT,*),IXP(NIXP,*), IXR(NIXR,*), IXTG(NIXTG,*)
70 INTEGER IBCV(*), IBCR(*), IGROUPS(*)
71 my_real FSKY(*), FCONV(*), FRADIA(*)
72 my_real x(3,*), factiv(lractiv,*), temp(*), mcp(*), pm(npropm,*)
73 my_real mcp_off(*)
74 TYPE(elbuf_struct_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
75 TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) ,INTENT(IN) :: SENSOR_TAB
76 type (glob_therm_) ,intent(inout) :: glob_therm
77C-----------------------------------------------
78 TYPE (GROUP_) , DIMENSION(NGRBRIC) :: IGRBRIC
79 TYPE (GROUP_) , DIMENSION(NGRQUAD) :: IGRQUAD
80 TYPE (GROUP_) , DIMENSION(NGRSHEL) :: IGRSH4N
81 TYPE (GROUP_) , DIMENSION(NGRSH3N) :: IGRSH3N
82 TYPE (GROUP_) , DIMENSION(NGRTRUS) :: IGRTRUSS
83 TYPE (GROUP_) , DIMENSION(NGRBEAM) :: IGRBEAM
84 TYPE (group_) , DIMENSION(NGRSPRI) :: igrspring
85C-----------------------------------------------
86C L o c a l V a r i a b l e s
87C-----------------------------------------------
88 INTEGER I,N,IFLAG,ISENS,IGSH,IGSH3,IGBR,IGQU,IGBM,IGTR,IGSP,OFFCHANGE
89 INTEGER ITHERM_FE
90 my_real startt, stopt
91C======================================================================|
92 offchange = 0
93 itherm_fe = glob_therm%ITHERM_FE
94!
95 DO n = 1,nactiv
96 isens = iactiv(2,n)
97 iflag = iactiv(1,n)
98 IF (iactiv(10,n) == 1) THEN
99 IF (isens > 0) THEN
100C----------------------------------------------------------------------
101C DESACTIVATION DES ELEMENTS
102C----------------------------------------------------------------------
103 IF (iflag == 1 .AND. tt > sensor_tab(isens)%TSTART) THEN
104 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
105 . ixr ,ixtg ,iparg ,
106 . iactiv ,tt ,iflag ,n ,elbuf_tab,
107 . x ,temp ,mcp ,pm ,igroups ,
108 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
109 . igrtruss,igrbeam ,igrspring,itherm_fe)
110
111 iactiv(1,n) = 0
112 offchange=1
113 IF (iparit /= 0) THEN
114 DO i=1,8*lsky
115 fsky(i)=zero
116 ENDDO
117 ENDIF
118C----------------------------------------------------------------------
119C ACTIVATION DES ELEMENTS
120C----------------------------------------------------------------------
121 ELSEIF (iflag == 0 .AND. tt <= sensor_tab(isens)%TSTART) THEN
122 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
123 . ixr ,ixtg ,iparg ,
124 . iactiv ,tt ,iflag ,n ,elbuf_tab,
125 . x ,temp ,mcp ,pm ,igroups ,
126 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
127 . igrtruss,igrbeam ,igrspring,itherm_fe)
128
129 iactiv(1,n) = 1
130 offchange=1
131 ENDIF
132 ENDIF
133
134 ELSEIF(iactiv(10,n) == 2) THEN
135C----------------------------------------------------------------------
136C DESACTIVATION DES ELEMENTS AT TIME ZERO
137C----------------------------------------------------------------------
138 IF(tt == zero) THEN
139 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
140 . ixr ,ixtg ,iparg ,
141 . iactiv ,tt ,1 ,n ,elbuf_tab,
142 . x ,temp ,mcp ,pm ,igroups ,
143 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
144 . igrtruss,igrbeam ,igrspring,itherm_fe)
145
146 iactiv(1,n) = 0
147 offchange=1
148 ENDIF
149C----------------------------------------------------------------------
150C ACTIVATION DES ELEMENTS
151C----------------------------------------------------------------------
152 startt = factiv(1,n)
153 stopt = factiv(2,n)
154 iflag = iactiv(1,n)
155 IF(itherm_fe > 0) THEN
156 startt = startt / glob_therm%THEACCFACT
157 stopt = stopt / glob_therm%THEACCFACT
158 ENDIF
159 IF(iflag == 0 .AND. tt >= startt .AND. tt < stopt) THEN
160 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
161 . ixr ,ixtg ,iparg ,
162 . iactiv ,tt ,iflag ,n ,elbuf_tab,
163 . x ,temp ,mcp ,pm ,igroups ,
164 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
165 . igrtruss,igrbeam ,igrspring,itherm_fe)
166 iactiv(1,n) = 1
167 offchange=1
168 ENDIF
169C----------------------------------------------------------------------
170C DESACTIVATION DES ELEMENTS
171C----------------------------------------------------------------------
172 IF(iflag == 1 .AND. tt > stopt) THEN
173 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
174 . ixr ,ixtg ,iparg ,
175 . iactiv ,tt ,1 ,n ,elbuf_tab,
176 . x ,temp ,mcp ,pm ,igroups ,
177 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
178 . igrtruss,igrbeam ,igrspring,itherm_fe)
179 iactiv(1,n) = 0
180 offchange=1
181 ENDIF
182 ENDIF
183 ENDDO ! N=1,NACTIV
184C
185 IF(itherm_fe > 0 .AND. offchange == 1) THEN
186 IF(glob_therm%NUMCONV > 0) CALL convecoff(ibcv, fconv, iparg, igroups, elbuf_tab,glob_therm)
187 IF(glob_therm%NUMRADIA > 0) CALL radiatoff(ibcr, fradia,iparg, igroups, ixs, elbuf_tab,glob_therm)
188 ENDIF
189
190 IF(itherm_fe > 0 .AND. tt == zero) THEN
191C Final call to set MCP_OFF(i) = 0 if all nodes
192 iflag=-1
193 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
194 . ixr ,ixtg ,iparg ,
195 . iactiv ,tt ,iflag ,n ,elbuf_tab,
196 . x ,temp ,mcp ,pm ,igroups ,
197 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
198 . igrtruss,igrbeam ,igrspring,itherm_fe)
199 ENDIF
200
201
202C
203 RETURN
204 END
205
convecoff
subroutine convecoff(ibcv, fconv, iparg, igroups, elbuf_tab, glob_therm)
Definition convecoff.F:32
desacti
subroutine desacti(ixs, ixq, ixc, ixp, ixt, ixr, ixtg, iparg, iactiv, nsensor, sensor_tab, fsky, x, elbuf_tab, ibcv, fconv, ibcr, fradia, igroups, factiv, temp, mcp, pm, mcp_off, igrbric, igrquad, igrsh4n, igrsh3n, igrtruss, igrbeam, igrspring, glob_therm)
Definition desacti.F:44
eloff
subroutine eloff(ixs, ixq, ixc, ixp, ixt, ixr, ixtg, iparg, iactiv, time, iflag, nn, elbuf_tab, x, temp, mcp, pm, igroups, mcp_off, igrbric, igrquad, igrsh4n, igrsh3n, igrtruss, igrbeam, igrspring, itherm_fe)
Definition eloff.F:42
groupdef_mod
Definition groupdef_mod.F:662
radiatoff
subroutine radiatoff(ibcr, fradia, iparg, igroups, ixs, elbuf_tab, glob_therm)
Definition radiatoff.F:32
groupdef_mod::group_
Definition groupdef_mod.F:700