OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
c3coor3_crk.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com08_c.inc"
#include "com_xfem1.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine c3coor3_crk (jft, jlt, offg, off, sigy, dt1c, thke, vl1, vl2, vl3, vrl1, vrl2, vrl3, x1, x2, x3, y1, y2, y3, z1, z2, z3, ilev, iel_crk, iadc_crk, nft)
subroutine c3coor3_crk2 (jft, jlt, nft, ilev, iel_crk, x, v, vr, vl1, vl2, vl3, vrl1, vrl2, vrl3, x1, x2, x3, y1, y2, y3, z1, z2, z3, iadtg_crk)

Function/Subroutine Documentation

◆ c3coor3_crk()

subroutine c3coor3_crk ( integer jft,
integer jlt,
offg,
off,
sigy,
dt1c,
thke,
vl1,
vl2,
vl3,
vrl1,
vrl2,
vrl3,
real(kind=8), dimension(mvsiz), intent(out) x1,
real(kind=8), dimension(mvsiz), intent(out) x2,
real(kind=8), dimension(mvsiz), intent(out) x3,
real(kind=8), dimension(mvsiz), intent(out) y1,
real(kind=8), dimension(mvsiz), intent(out) y2,
real(kind=8), dimension(mvsiz), intent(out) y3,
real(kind=8), dimension(mvsiz), intent(out) z1,
real(kind=8), dimension(mvsiz), intent(out) z2,
real(kind=8), dimension(mvsiz), intent(out) z3,
integer ilev,
integer, dimension(*) iel_crk,
integer, dimension(3,*) iadc_crk,
integer nft )

Definition at line 30 of file c3coor3_crk.F.

36C-----------------------------------------------
37 USE crackxfem_mod
38C-----------------------------------------------
39C I m p l i c i t T y p e s
40C-----------------------------------------------
41#include "implicit_f.inc"
42C-----------------------------------------------
43C G l o b a l P a r a m e t e r s
44C-----------------------------------------------
45#include "mvsiz_p.inc"
46C-----------------------------------------------
47C C o m m o n B l o c k s
48C-----------------------------------------------
49#include "com08_c.inc"
50#include "com_xfem1.inc"
51C-----------------------------------------------
52C D u m m y A r g u m e n t s
53C-----------------------------------------------
54 INTEGER JFT, JLT
55 INTEGER ILEV,IEL_CRK(*),IADC_CRK(3,*),NFT
56 my_real offg(*), off(*),dt1c(*),thke(*),
57 . vl1(mvsiz,3),vl2(mvsiz,3),vl3(mvsiz,3),
58 . vrl1(mvsiz,3),vrl2(mvsiz,3),vrl3(mvsiz,3),sigy(*)
59! SP issue :
60 REAL(kind=8), dimension(mvsiz), INTENT(out) ::x1,x2,x3
61 REAL(kind=8), dimension(mvsiz), INTENT(out) ::y1,y2,y3
62 REAL(kind=8), dimension(mvsiz), INTENT(out) ::z1,z2,z3
63C-----------------------------------------------
64C L o c a l V a r i a b l e s
65C-----------------------------------------------
66 INTEGER I,J,NC1,NC2,NC3,ELCRK,ELCRK0,IADC1,IADC2,IADC3
67 my_real off_l
68C=======================================================================
69C get elem from activ level: (superposed element)
70C-----------------------------------------------
71 DO i=jft,jlt
72 elcrk0 = iel_crk(i+nft)
73 elcrk = elcrk0 + ecrkxfec
74 iadc1 = iadc_crk(1,elcrk0)
75 iadc2 = iadc_crk(2,elcrk0)
76 iadc3 = iadc_crk(3,elcrk0)
77C--------------
78C COORDINATES:
79C--------------
80C node 1:
81 x1(i) = crkavx(ilev)%X(1,iadc1)
82 y1(i) = crkavx(ilev)%X(2,iadc1)
83 z1(i) = crkavx(ilev)%X(3,iadc1)
84C node 2:
85 x2(i) = crkavx(ilev)%X(1,iadc2)
86 y2(i) = crkavx(ilev)%X(2,iadc2)
87 z2(i) = crkavx(ilev)%X(3,iadc2)
88C node 3:
89 x3(i) = crkavx(ilev)%X(1,iadc3)
90 y3(i) = crkavx(ilev)%X(2,iadc3)
91 z3(i) = crkavx(ilev)%X(3,iadc3)
92C--------------
93C--------------
94C VELOCITIES :
95C--------------
96C node 1:
97 vl1(i,1) = crkavx(ilev)%V(1,iadc1)
98 vl1(i,2) = crkavx(ilev)%V(2,iadc1)
99 vl1(i,3) = crkavx(ilev)%V(3,iadc1)
100 vrl1(i,1) = crkavx(ilev)%VR(1,iadc1)
101 vrl1(i,2) = crkavx(ilev)%VR(2,iadc1)
102 vrl1(i,3) = crkavx(ilev)%VR(3,iadc1)
103C node 2:
104 vl2(i,1) = crkavx(ilev)%V(1,iadc2)
105 vl2(i,2) = crkavx(ilev)%V(2,iadc2)
106 vl2(i,3) = crkavx(ilev)%V(3,iadc2)
107 vrl2(i,1) = crkavx(ilev)%VR(1,iadc2)
108 vrl2(i,2) = crkavx(ilev)%VR(2,iadc2)
109 vrl2(i,3) = crkavx(ilev)%VR(3,iadc2)
110C node 3:
111 vl3(i,1) = crkavx(ilev)%V(1,iadc3)
112 vl3(i,2) = crkavx(ilev)%V(2,iadc3)
113 vl3(i,3) = crkavx(ilev)%V(3,iadc3)
114 vrl3(i,1) = crkavx(ilev)%VR(1,iadc3)
115 vrl3(i,2) = crkavx(ilev)%VR(2,iadc3)
116 vrl3(i,3) = crkavx(ilev)%VR(3,iadc3)
117 END DO
118C-----------------------------------------------
119 off_l = zero
120 DO 40 i=jft,jlt
121 dt1c(i) = dt1
122 off(i) = min(one,abs(offg(i)))
123 off_l = min(off_l,offg(i))
124 sigy(i) = ep30
125 40 CONTINUE
126C---
127 IF(off_l<zero)THEN
128 DO i=jft,jlt
129 IF(offg(i)<zero)THEN
130 vl1(i,1)=zero
131 vl1(i,2)=zero
132 vl1(i,3)=zero
133 vl2(i,1)=zero
134 vl2(i,2)=zero
135 vl2(i,3)=zero
136 vl3(i,1)=zero
137 vl3(i,2)=zero
138 vl3(i,3)=zero
139 vrl1(i,1)=zero
140 vrl1(i,2)=zero
141 vrl1(i,3)=zero
142 vrl2(i,1)=zero
143 vrl2(i,2)=zero
144 vrl2(i,3)=zero
145 vrl3(i,1)=zero
146 vrl3(i,2)=zero
147 vrl3(i,3)=zero
148 ENDIF
149 ENDDO
150 ENDIF
151C---
152 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
crackxfem_mod
Definition crackxfem_mod.F:121
crackxfem_mod::crkavx
type(xfem_avx_), dimension(:), allocatable crkavx
Definition crackxfem_mod.F:163

◆ c3coor3_crk2()

subroutine c3coor3_crk2 ( integer jft,
integer jlt,
integer nft,
integer ilev,
integer, dimension(*) iel_crk,
x,
v,
vr,
vl1,
vl2,
vl3,
vrl1,
vrl2,
vrl3,
x1,
x2,
x3,
y1,
y2,
y3,
z1,
z2,
z3,
integer, dimension(3,*) iadtg_crk )

Definition at line 161 of file c3coor3_crk.F.

167C-----------------------------------------------
168 USE crackxfem_mod
169C-----------------------------------------------
170C I m p l i c i t T y p e s
171C-----------------------------------------------
172#include "implicit_f.inc"
173C-----------------------------------------------
174C G l o b a l P a r a m e t e r s
175C-----------------------------------------------
176#include "mvsiz_p.inc"
177C-----------------------------------------------
178C C o m m o n B l o c k s
179C-----------------------------------------------
180#include "com_xfem1.inc"
181C-----------------------------------------------
182C D u m m y A r g u m e n t s
183C-----------------------------------------------
184 INTEGER JFT,JLT,NFT,ILEV
185 INTEGER IEL_CRK(*),IADTG_CRK(3,*)
186C REAL
187 my_real
188 . x(3,*),v(3,*),vr(3,*),vl1(mvsiz,3),vl2(mvsiz,3),vl3(mvsiz,3),
189 . vrl1(mvsiz,3),vrl2(mvsiz,3),vrl3(mvsiz,3)
190 my_real
191 . x1(mvsiz), x2(mvsiz), x3(mvsiz),
192 . y1(mvsiz), y2(mvsiz), y3(mvsiz),
193 . z1(mvsiz), z2(mvsiz), z3(mvsiz)
194C-----------------------------------------------
195C L o c a l V a r i a b l e s
196C-----------------------------------------------
197 INTEGER I,ELCRK,IADC1,IADC2,IADC3,ELCRK0
198C-----------------------------------------------
199 DO i=jft,jlt
200 elcrk0 = iel_crk(i+nft)
201 elcrk = elcrk0 + ecrkxfec
202 iadc1 = iadtg_crk(1,elcrk0)
203 iadc2 = iadtg_crk(2,elcrk0)
204 iadc3 = iadtg_crk(3,elcrk0)
205C--------------
206C COORDINATES:
207C--------------
208C node 1:
209 x1(i) = crkavx(ilev)%X(1,iadc1)
210 y1(i) = crkavx(ilev)%X(2,iadc1)
211 z1(i) = crkavx(ilev)%X(3,iadc1)
212C node 2:
213 x2(i) = crkavx(ilev)%X(1,iadc2)
214 y2(i) = crkavx(ilev)%X(2,iadc2)
215 z2(i) = crkavx(ilev)%X(3,iadc2)
216C node 3:
217 x3(i) = crkavx(ilev)%X(1,iadc3)
218 y3(i) = crkavx(ilev)%X(2,iadc3)
219 z3(i) = crkavx(ilev)%X(3,iadc3)
220C--------------
221C VELOCITIES :
222C--------------
223C node 1:
224 vl1(i,1) = crkavx(ilev)%V(1,iadc1)
225 vl1(i,2) = crkavx(ilev)%V(2,iadc1)
226 vl1(i,3) = crkavx(ilev)%V(3,iadc1)
227 vrl1(i,1) = crkavx(ilev)%VR(1,iadc1)
228 vrl1(i,2) = crkavx(ilev)%VR(2,iadc1)
229 vrl1(i,3) = crkavx(ilev)%VR(3,iadc1)
230C node 2:
231 vl2(i,1) = crkavx(ilev)%V(1,iadc2)
232 vl2(i,2) = crkavx(ilev)%V(2,iadc2)
233 vl2(i,3) = crkavx(ilev)%V(3,iadc2)
234 vrl2(i,1) = crkavx(ilev)%VR(1,iadc2)
235 vrl2(i,2) = crkavx(ilev)%VR(2,iadc2)
236 vrl2(i,3) = crkavx(ilev)%VR(3,iadc2)
237C node 3:
238 vl3(i,1) = crkavx(ilev)%V(1,iadc3)
239 vl3(i,2) = crkavx(ilev)%V(2,iadc3)
240 vl3(i,3) = crkavx(ilev)%V(3,iadc3)
241 vrl3(i,1) = crkavx(ilev)%VR(1,iadc3)
242 vrl3(i,2) = crkavx(ilev)%VR(2,iadc3)
243 vrl3(i,3) = crkavx(ilev)%VR(3,iadc3)
244 END DO
245C-----------------------------------------------
246C
247 RETURN