OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
dtnodarayl.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!|| dtnodarayl ../engine/source/time_step/dtnodarayl.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- calls -----------------------------------------------------
28!|| damping_vref_compute_dampa ../engine/source/assembly/damping_vref_compute_dampa.F90
29!||--- uses -----------------------------------------------------
30!|| damping_vref_compute_dampa_mod ../engine/source/assembly/damping_vref_compute_dampa.F90
31!|| groupdef_mod ../common_source/modules/groupdef_mod.F
32!||====================================================================
33 SUBROUTINE dtnodarayl(MS ,IN ,STIFN ,STIFR ,DT2T ,
34 1 IGRNOD ,DAMPR)
35C-----------------------------------------------
36C M o d u l e s
37C-----------------------------------------------
38 USE groupdef_mod
39 use damping_vref_compute_dampa_mod
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C C o m m o n B l o c k s
46C-----------------------------------------------
47#include "com01_c.inc"
48#include "com04_c.inc"
49#include "com06_c.inc"
50#include "com08_c.inc"
51#include "scr18_c.inc"
52#include "param_c.inc"
53#include "stati_c.inc"
54#include "statr_c.inc"
55C-----------------------------------------------------------------
56C D u m m y A r g u m e n t s
57C-----------------------------------------------
58C REAL
59 my_real , INTENT(IN) :: dt2t
60 my_real , DIMENSION(NUMNOD), INTENT(IN) ::ms
61 my_real , DIMENSION(NUMNOD*IRODDL), INTENT(IN) ::in
62 my_real , DIMENSION(NUMNOD), INTENT(INOUT) ::stifn
63 my_real , DIMENSION(NUMNOD*IRODDL), INTENT(INOUT) ::stifr
64 my_real , DIMENSION(NRDAMP,NDAMP), INTENT(INOUT) ::dampr
65C-----------------------------------------------
66 TYPE (GROUP_) , DIMENSION(NGRNOD), INTENT(IN) :: IGRNOD
67C-----------------------------------------------
68C L o c a l V a r i a b l e s
69C-----------------------------------------------
70 INTEGER N,I,J,ISK,IGR,ND,K,L,LL,ITYPE,FL_VREL,ID_RBY,FL_FREQ_RANGE
71 my_real factb,dampbr,dampar,d_tstart,d_tstop,dampt,
72 . dt2n(numnod),bb,coeff,dtn,dtn2,dampai,dampbi,dt0,dampa3,
73 . alpha05,beta1,dtmax,fac,dtmi1,dtmi2,bbi,faci,damp_a(3)
74C----------------------------------------------------------
75 dtmi1 = dtmin1(11)
76 dtmax = ep04
77 dt0 = dtmax
78 faci = zero
79 IF (dt1>zero) THEN
80 dt0=min(dt0,dt1)
81 dtmi1 = max(dtmi1,dt0)
82 ELSE
83 IF (dtmi1>zero) dt0=min(dt0,dtmi1)
84 END IF
85 dtmi2 = dtmi1/dtfac1(11)
86c--
87 DO i=1,ndamp
88C--
89 fl_vrel = 0
90 fl_freq_range = 0
91 itype = nint(dampr(21,i))
92 SELECT CASE (itype)
93 CASE(2)
94 fl_vrel = 1
95 CASE(3)
96 fl_freq_range = 1
97 END SELECT
98! FL_FREQ_RANGE = NINT(DAMPR(31,I))
99 IF ((dampr(19,i)/=0).OR.(fl_freq_range > 0)) cycle
100 d_tstart = dampr(17,i)
101 d_tstop = dampr(18,i)
102 IF ((tt<d_tstart).OR.(tt>d_tstop)) cycle
103 igr = nint(dampr(2,i))
104 dtn = dtmax
105 DO j=1,igrnod(igr)%NENTITY
106 nd=igrnod(igr)%ENTITY(j)
107 dt2n(nd) = dtmax
108 IF(stifn(nd)>em20 .AND. ms(nd)>zero)THEN
109 dt2n(nd) = sqrt(two*ms(nd)/stifn(nd))
110 dt2n(nd) = min(dtmax,dt2n(nd))
111 dtn = min(dt2n(nd),dtn)
112 END IF
113 END DO
114 factb = one
115 dampai = max(dampr(3,i),dampr(5,i),dampr(7,i))
116 dampai = max(dampai,dampr(9,i),dampr(11,i),dampr(13,i))
117 dampbi = max(dampr(4,i),dampr(6,i),dampr(8,i))
118 dampbi = max(dampbi,dampr(10,i),dampr(12,i),dampr(14,i))
119 dampt = dt0*factb
120 dampbi = min(dampbi,dampt)
121C--
122C-- /DAMP/VREL - recompute damping parameters at current time
123! FL_VREL = NINT(DAMPR(21,I))
124 IF (fl_vrel==1) THEN
125 call damping_vref_compute_dampa(i,ndamp,nrdamp,dampr,dt1,tt,damp_a)
126 dampai = max(damp_a(1),damp_a(2),damp_a(3))
127 dampbi = zero
128 ENDIF
129C--
130 IF (dtmi2 > zero) THEN
131 bbi = dampbi+half*dampai*dtmi2*dtmi2
132 dtn2= sqrt(dtmi2*dtmi2+two*bbi*dtmi2)
133 bbi = dampbi/dtn2+half*dampai*dtn2
134 faci = sqrt(bbi*bbi+one) -bbi
135 END IF
136 DO j=1,igrnod(igr)%NENTITY
137 nd=igrnod(igr)%ENTITY(j)
138 IF(stifn(nd)>em20 .AND. ms(nd)>zero)THEN
139 bb=dampbi/dt2n(nd)+half*dampai*dt2n(nd)
140 fac = sqrt(bb*bb+one) -bb
141 IF (fac*dt2n(nd) < dtmi2) THEN
142 fac = max(fac,faci)
143 END IF
144 coeff =one/fac**2
145 stifn(nd) =stifn(nd)*coeff
146 END IF
147 END DO
148 IF (iroddl/=0) THEN
149 DO j=1,igrnod(igr)%NENTITY
150 nd=igrnod(igr)%ENTITY(j)
151 IF(stifr(nd)>em20 .AND. in(nd)>zero)THEN
152 dt2n(nd) = sqrt(two*in(nd)/stifr(nd))
153 dt2n(nd) = min(dtmax,dt2n(nd))
154 bb=dampbi/dt2n(nd)+half*dampai*dt2n(nd)
155 fac = sqrt(bb*bb+one) -bb
156 IF (fac*dt2n(nd) < dtmi2) THEN
157 fac = max(fac,faci)
158 END IF
159 coeff =one/fac**2
160 stifr(nd) =stifr(nd)*coeff
161 END IF
162 END DO
163 END IF
164 END DO !I=1,NDAMP
165C----------global damping
166 IF(idamp>0.OR.istat==3)THEN
167 IF (istat==3 .AND.tt > zero) THEN
168 dampa3 = two*betate/(one + betate * dt12)
169 ELSE
170 dampa3 = zero
171 END IF
172 dtn2 = dtmax
173C------- still factor HALF* in another program
174 alpha05 = half*(dampa+dampa3)
175 beta1 = min(dampb,dt0)
176 IF (dtmi2 > zero) THEN
177 bbi = beta1+alpha05*dtmi2*dtmi2
178 dtn2= sqrt(dtmi2*dtmi2+two*bbi*dtmi2)
179 bbi = beta1/dtn2+alpha05*dtn2
180 faci = sqrt(bbi*bbi+one) -bbi
181 END IF
182 DO nd=1,numnod
183 coeff =one
184 IF(stifn(nd)>em20 .AND. ms(nd)>zero)THEN
185 dt2n(nd) = sqrt(two*ms(nd)/stifn(nd))
186 dt2n(nd) = min(dtmax,dt2n(nd))
187 bb = beta1/dt2n(nd) + alpha05*dt2n(nd)
188 fac = sqrt(bb*bb+one) -bb
189 IF (fac*dt2n(nd) < dtmi2) THEN
190 fac = max(fac,faci)
191 END IF
192 coeff =one/fac**2
193 stifn(nd) =stifn(nd)*coeff
194 END IF
195 END DO
196 IF (iroddl/=0) THEN
197 DO nd=1,numnod
198 coeff =one
199 IF(stifr(nd)>em20 .AND. in(nd)>zero)THEN
200 dt2n(nd) = sqrt(two*in(nd)/stifr(nd))
201 dt2n(nd) = min(dtmax,dt2n(nd))
202 bb = beta1/dt2n(nd) + alpha05*dt2n(nd)
203 fac = sqrt(bb*bb+one) -bb
204 IF (fac*dt2n(nd) < dtmi2) THEN
205 fac = max(fac,faci)
206 END IF
207 coeff =one/fac**2
208 stifr(nd) =stifr(nd)*coeff
209 END IF
210 END DO
211 END IF
212 END IF
213C
214 RETURN
215 END SUBROUTINE dtnodarayl
my_real
#define my_real
Definition cppsort.cpp:32
dtnodarayl
subroutine dtnodarayl(ms, in, stifn, stifr, dt2t, igrnod, dampr)
Definition dtnodarayl.F:35
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
groupdef_mod
Definition groupdef_mod.F:662