OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
crk_vitesse.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!|| crk_vitesse ../engine/source/elements/xfem/crk_vitesse.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- uses -----------------------------------------------------
28!|| crackxfem_mod ../engine/share/modules/crackxfem_mod.F
29!||====================================================================
30 SUBROUTINE crk_vitesse(
31 . ADDCNE_CRK,INOD_CRK ,NODLEVXF ,NODFT ,NODLT ,
32 . X ,V ,VR ,A ,AR ,
33 . ITAB )
34C-----------------------------------------------
35 USE crackxfem_mod
36C-----------------------------------------------
37C I m p l i c i t T y p e s
38C-----------------------------------------------
39#include "implicit_f.inc"
40C-----------------------------------------------
41C C o m m o n B l o c k s
42C-----------------------------------------------
43#include "com04_c.inc"
44#include "com08_c.inc"
45#include "com_xfem1.inc"
46C-----------------------------------------------------------------
47C D u m m y A r g u m e n t s
48C-----------------------------------------------
49 INTEGER NODFT,NODLT
50 INTEGER ADDCNE_CRK(*),INOD_CRK(*),NODLEVXF(*),ITAB(NUMNOD)
51 my_real X(3,*),V(3,*),VR(3,*),A(3,*),AR(3,*)
52C-----------------------------------------------
53C L o c a l V a r i a b l e s
54C-----------------------------------------------
55 INTEGER I,J,K,KK,IL,N,IAD,IAD1,IAD2,NN,NC,NCT,EN0,EN,EN1,IFI,ILAY,
56 . IXEL,ILEV,NLEV,N1,N2,NX1
57 my_real X1,X2,X3,V1,V2,V3,A1,A2,A3,AR1,AR2,AR3,DX,DY,DZ
58C=======================================================================
59 DO n = nodft,nodlt
60 nn = inod_crk(n)
61 IF (nn <= 0) cycle
62 nct = addcne_crk(nn)-1
63 nc = addcne_crk(nn+1)-addcne_crk(nn)
64 nlev = nodlevxf(nn)
65c---
66 DO ilev = 1,nlev
67 ixel = mod(ilev-1, nxel) + 1
68 ilay = (ilev-ixel)/nxel + 1
69 DO kk = nct+1, nct+nc
70 en0 = crklvset(ilev)%ENR0(2,kk) ! enr initial du debut de cycle
71 en = crklvset(ilev)%ENR0(1,kk) ! enr mise a jour dans le cycle
72 ifi = xfem_phantom(ilay)%IFI(kk)
73c---
74 IF (en0 > 0) THEN ! free node
75 a1 = crkavx(ilev)%A(1,kk)
76 a2 = crkavx(ilev)%A(2,kk)
77 a3 = crkavx(ilev)%A(3,kk)
78 ar1 = crkavx(ilev)%AR(1,kk)
79 ar2 = crkavx(ilev)%AR(2,kk)
80 ar3 = crkavx(ilev)%AR(3,kk)
81c new velocity
82 crkavx(ilev)%V(1,kk) = crkavx(ilev)%V(1,kk) + dt12*a1
83 crkavx(ilev)%V(2,kk) = crkavx(ilev)%V(2,kk) + dt12*a2
84 crkavx(ilev)%V(3,kk) = crkavx(ilev)%V(3,kk) + dt12*a3
85 crkavx(ilev)%VR(1,kk) = crkavx(ilev)%VR(1,kk)+ dt12*ar1
86 crkavx(ilev)%VR(2,kk) = crkavx(ilev)%VR(2,kk)+ dt12*ar2
87 crkavx(ilev)%VR(3,kk) = crkavx(ilev)%VR(3,kk)+ dt12*ar3
88 crkavx(ilev)%A(1,kk) = zero
89 crkavx(ilev)%A(2,kk) = zero
90 crkavx(ilev)%A(3,kk) = zero
91 crkavx(ilev)%AR(1,kk) = zero
92 crkavx(ilev)%AR(2,kk) = zero
93 crkavx(ilev)%AR(3,kk) = zero
94c new coordinates
95 dx = dt2 * crkavx(ilev)%V(1,kk)
96 dy = dt2 * crkavx(ilev)%V(2,kk)
97 dz = dt2 * crkavx(ilev)%V(3,kk)
98 crkavx(ilev)%X(1,kk) = crkavx(ilev)%X(1,kk) + dx
99 crkavx(ilev)%X(2,kk) = crkavx(ilev)%X(2,kk) + dy
100 crkavx(ilev)%X(3,kk) = crkavx(ilev)%X(3,kk) + dz
101 ELSEIF (en0 <= 0) THEN
102 crkavx(ilev)%X(1,kk) = x(1,n)
103 crkavx(ilev)%X(2,kk) = x(2,n)
104 crkavx(ilev)%X(3,kk) = x(3,n)
105 IF (en > 0) THEN ! noeud relache => vit noeud std instead of zero
106 crkavx(ilev)%V(1,kk) = v(1,n)
107 crkavx(ilev)%V(2,kk) = v(2,n)
108 crkavx(ilev)%V(3,kk) = v(3,n)
109 crkavx(ilev)%VR(1,kk) = vr(1,n)
110 crkavx(ilev)%VR(2,kk) = vr(2,n)
111 crkavx(ilev)%VR(3,kk) = vr(3,n)
112 ELSEIF (ifi /= 0) THEN ! fissure statique => noeud standard
113 crkavx(ilev)%V(1,kk) = v(1,n)
114 crkavx(ilev)%V(2,kk) = v(2,n)
115 crkavx(ilev)%V(3,kk) = v(3,n)
116 crkavx(ilev)%VR(1,kk) = vr(1,n)
117 crkavx(ilev)%VR(2,kk) = vr(2,n)
118 crkavx(ilev)%VR(3,kk) = vr(3,n)
119 ENDIF
120 END IF
121c---
122 ENDDO ! DO KK = NCT+1, NCT+NC
123 ENDDO ! DO ilev= nlev,1,-1
124 ENDDO ! N = NODFT,NODLT
125c-----------
126 RETURN
127 END
128c
129!||====================================================================
130!|| crk_zero_accel ../engine/source/elements/xfem/crk_vitesse.F
131!||--- called by ------------------------------------------------------
132!|| resol ../engine/source/engine/resol.F
133!||--- uses -----------------------------------------------------
134!|| crackxfem_mod ../engine/share/modules/crackxfem_mod.F
135!||====================================================================
136 SUBROUTINE crk_zero_accel(ADDCNE_CRK,INOD_CRK,NODFT,NODLT,NODLEVXF)
137C-----------------------------------------------
138 USE crackxfem_mod
139C-----------------------------------------------
140C I m p l i c i t T y p e s
141C-----------------------------------------------
142#include "implicit_f.inc"
143C-----------------------------------------------
144C C o m m o n B l o c k s
145C-----------------------------------------------
146C-----------------------------------------------------------------
147C D u m m y A r g u m e n t s
148C-----------------------------------------------
149 INTEGER ADDCNE_CRK(*),INOD_CRK(*),NODFT,NODLT,NODLEVXF(*)
150C-----------------------------------------------
151C L o c a l V a r i a b l e s
152C-----------------------------------------------
153 INTEGER I,KK,N,NN,NCT,NC,ILEV,NLEV
154C=======================================================================
155 DO N = nodft,nodlt
156 nn = inod_crk(n)
157 IF (nn > 0) THEN
158C---
159 nct = addcne_crk(nn)-1
160 nc = addcne_crk(nn+1)-addcne_crk(nn)
161 nlev = nodlevxf(nn)
162 DO kk = nct+1, nct+nc
163 DO ilev=1,nlev
164 crkavx(ilev)%A(1,kk) = zero
165 crkavx(ilev)%A(2,kk) = zero
166 crkavx(ilev)%A(3,kk) = zero
167 crkavx(ilev)%AR(1,kk) = zero
168 crkavx(ilev)%AR(2,kk) = zero
169 crkavx(ilev)%AR(3,kk) = zero
170 END DO
171 END DO
172 ENDIF
173 ENDDO
174C-----------
175 RETURN
176 END
177!||====================================================================
178!|| crk_zero_fsky ../engine/source/elements/xfem/crk_vitesse.F
179!||--- called by ------------------------------------------------------
180!|| resol ../engine/source/engine/resol.F
181!||--- uses -----------------------------------------------------
182!|| crackxfem_mod ../engine/share/modules/crackxfem_mod.F
183!||====================================================================
184 SUBROUTINE crk_zero_fsky(CRKSKY,ADDCNE_CRK,INOD_CRK,NODFT,NODLT,
185 . NODLEVXF)
186C-----------------------------------------------
187 USE crackxfem_mod
188C-----------------------------------------------
189C I m p l i c i t T y p e s
190C-----------------------------------------------
191#include "implicit_f.inc"
192C-----------------------------------------------
193C C o m m o n B l o c k s
194C-----------------------------------------------
195C-----------------------------------------------------------------
196C D u m m y A r g u m e n t s
197C-----------------------------------------------
198 INTEGER ADDCNE_CRK(*),INOD_CRK(*),NODFT,NODLT,NODLEVXF(*)
199 TYPE(XFEM_SKY_) , DIMENSION(*) :: CRKSKY
200C-----------------------------------------------
201C L o c a l V a r i a b l e s
202C-----------------------------------------------
203 INTEGER I,KK,N,NN,NCT,NC,ILEV,NLEV
204C=======================================================================
205 DO N = nodft,nodlt
206 nn = inod_crk(n)
207 IF (nn > 0) THEN
208 nct = addcne_crk(nn)-1
209 nc = addcne_crk(nn+1)-addcne_crk(nn)
210 nlev = nodlevxf(nn)
211 DO kk = nct+1, nct+nc
212 DO ilev=1,nlev
213 crksky(ilev)%FSKY(1,kk) = zero
214 crksky(ilev)%FSKY(2,kk) = zero
215 crksky(ilev)%FSKY(3,kk) = zero
216 crksky(ilev)%FSKY(4,kk) = zero
217 crksky(ilev)%FSKY(5,kk) = zero
218 crksky(ilev)%FSKY(6,kk) = zero
219 END DO
220 END DO
221 ENDIF
222 ENDDO
223C-----------
224 RETURN
225 END
crk_zero_accel
subroutine crk_zero_accel(addcne_crk, inod_crk, nodft, nodlt, nodlevxf)
Definition crk_vitesse.F:137
crk_zero_fsky
subroutine crk_zero_fsky(crksky, addcne_crk, inod_crk, nodft, nodlt, nodlevxf)
Definition crk_vitesse.F:186
crk_vitesse
subroutine crk_vitesse(addcne_crk, inod_crk, nodlevxf, nodft, nodlt, x, v, vr, a, ar, itab)
Definition crk_vitesse.F:34
crackxfem_mod
Definition crackxfem_mod.F:121
crackxfem_mod::xfem_phantom
type(xfem_phantom_), dimension(:), allocatable xfem_phantom
Definition crackxfem_mod.F:139
crackxfem_mod::crkavx
type(xfem_avx_), dimension(:), allocatable crkavx
Definition crackxfem_mod.F:163
crackxfem_mod::crklvset
type(xfem_lvset_), dimension(:), allocatable crklvset
Definition crackxfem_mod.F:146