OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cupdt_ply.F File Reference
#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "parit_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine cupdt_ply (jft, jlt, nvc, offg, off, iadc, iel, inod, ixc, ms, in, ms_ply, zi_ply, istack, posly, fly11, fly12, fly13, fly14, fly21, fly22, fly23, fly24, fly31, fly32, fly33, fly34, fac, sti, msz2, nft, npt)

Function/Subroutine Documentation

◆ cupdt_ply()

subroutine cupdt_ply ( integer jft,
integer jlt,
integer nvc,
offg,
off,
integer, dimension(4,*) iadc,
integer, dimension(*) iel,
integer, dimension(*) inod,
integer, dimension(nixc,mvsiz) ixc,
ms,
in,
ms_ply,
zi_ply,
integer, dimension(mvsiz,npt) istack,
posly,
fly11,
fly12,
fly13,
fly14,
fly21,
fly22,
fly23,
fly24,
fly31,
fly32,
fly33,
fly34,
fac,
sti,
msz2,
integer, intent(in) nft,
integer, intent(in) npt )

Definition at line 30 of file cupdt_ply.F.

39C-----------------------------------------------
40C M o d u l e s
41C-----------------------------------------------
42 USE plyxfem_mod
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47#include "comlock.inc"
48C-----------------------------------------------
49C G l o b a l P a r a m e t e r s
50C-----------------------------------------------
51#include "mvsiz_p.inc"
52#include "parit_c.inc"
53C-----------------------------------------------
54C C o m m o n B l o c k s
55C-----------------------------------------------
56#include "param_c.inc"
57C-----------------------------------------------
58C D u m m y A r g u m e n t s
59C-----------------------------------------------
60 INTEGER, INTENT(IN) :: NFT
61 INTEGER, INTENT(IN) :: NPT
62 INTEGER JFT, JLT, NVC
63 INTEGER IADC(4,*),IXC(NIXC,MVSIZ),ISTACK(MVSIZ,NPT),IEL(*),
64 . INOD(*)
65C REAL
66 my_real
67 . offg(*), off(*), in(*), ms(*),ms_ply(nplyxfe,*),
68 . posly(mvsiz,npt),zi_ply(nplyxfe,*),
69 . fly11(mvsiz, npt), fly21(mvsiz, npt), fly31(mvsiz, npt),
70 . fly12(mvsiz, npt), fly22(mvsiz, npt), fly32(mvsiz, npt),
71 . fly13(mvsiz, npt), fly23(mvsiz, npt), fly33(mvsiz, npt),
72 . fly14(mvsiz, npt), fly24(mvsiz, npt), fly34(mvsiz, npt),
73 . pth(mvsiz,npt),fac(mvsiz,2),sti(mvsiz,*),msz2(*)
74C-----------------------------------------------
75C L o c a l V a r i a b l e s
76C-----------------------------------------------
77 INTEGER NVC1, NVC2, NVC3, NVC4, I, J, MX,IOFF,IP,I1,I2,I3,I4,N,
78 . IN1,IN2,IN3,IN4
79 my_real
80 . off_l,
81 . fply11,fply21,fply31,fply12,fply22,fply32,fply13,fply23,
82 . fply33,fply14,fply24,fply34,fply,mply,
83 . mply11,mply21,mply31,mply12,mply22,mply32,mply13,mply23,
84 . mply33,mply14,mply24,mply34
85C-----------------------------------------------
86C
87C cumul de l'energie des elements deletes AU moment du delete
88C
89 off_l = zero
90 DO i=jft,jlt
91 IF(off(i)<one)offg(i) = off(i)
92 off_l = min(off_l,offg(i))
93 ENDDO
94 IF(off_l<zero)THEN
95 DO i=jft,jlt
96 IF(offg(i)<zero)THEN
97 DO j=1,npt
98 fly11(i,j)=zero
99 fly21(i,j)=zero
100 fly31(i,j)=zero
101 fly12(i,j)=zero
102 fly22(i,j)=zero
103 fly32(i,j)=zero
104 fly13(i,j)=zero
105 fly23(i,j)=zero
106 fly33(i,j)=zero
107 fly14(i,j)=zero
108 fly24(i,j)=zero
109 fly34(i,j)=zero
110 ENDDO
111 ENDIF
112 ENDDO
113 ENDIF
114C
115 IF(iparit == 0) THEN
116 nvc1= nvc/8
117 nvc2=(nvc-nvc1*8)/4
118 nvc3=(nvc-nvc1*8-nvc2*4)/2
119 nvc4=(nvc-nvc1*8-nvc2*4-nvc3*2)
120C to do
121C
122 DO i=jft,jlt
123
124 i1 = inod(ixc(2,i))
125 i2 = inod(ixc(3,i))
126 i3 = inod(ixc(4,i))
127 i4 = inod(ixc(5,i))
128C
129 in1 = ixc(2,i)
130 in2 = ixc(3,i)
131 in3 = ixc(4,i)
132 in4 = ixc(5,i)
133 n = iel(i + nft)
134 DO j=1,npt
135 ip = istack(i,j)
136C PRINT *,J,IP
137 ply(ip)%A(1,i1)= ply(ip)%A(1,i1) - (fly11(i,j) )
138 ply(ip)%A(2,i1)= ply(ip)%A(2,i1) - (fly21(i,j) )
139 ply(ip)%A(3,i1)= ply(ip)%A(3,i1) - (fly31(i,j) )
140 ply(ip)%A(4,i1)= ply(ip)%A(4,i1) + sti(i,j)*fac(i,1)
141c
142 ply(ip)%A(1,i2)= ply(ip)%A(1,i2) - (fly12(i,j) )
143 ply(ip)%A(2,i2)= ply(ip)%A(2,i2) - (fly22(i,j) )
144 ply(ip)%A(3,i2)= ply(ip)%A(3,i2) - (fly32(i,j) )
145 ply(ip)%A(4,i2)= ply(ip)%A(4,i2) + sti(i,j)*fac(i,2)
146
147c
148 ply(ip)%A(1,i3)= ply(ip)%A(1,i3) - (fly13(i,j) )
149 ply(ip)%A(2,i3)= ply(ip)%A(2,i3) - (fly23(i,j) )
150 ply(ip)%A(3,i3)= ply(ip)%A(3,i3) - (fly33(i,j) )
151 ply(ip)%A(4,i3)= ply(ip)%A(4,i3) + sti(i,j)*fac(i,1)
152c
153 ply(ip)%A(1,i4)= ply(ip)%A(1,i4) - (fly14(i,j) )
154 ply(ip)%A(2,i4)= ply(ip)%A(2,i4) - (fly24(i,j) )
155 ply(ip)%A(3,i4)= ply(ip)%A(3,i4) - (fly34(i,j) )
156 ply(ip)%A(4,i4)= ply(ip)%A(4,i4) + sti(i,j)*fac(i,2)
157C
158 ENDDO
159 ENDDO
160
161 ELSE
162 DO i=jft,jlt
163
164 i1 = inod(ixc(2,i))
165 i2 = inod(ixc(3,i))
166 i3 = inod(ixc(4,i))
167 i4 = inod(ixc(5,i))
168C
169 in1 = ixc(2,i)
170 in2 = ixc(3,i)
171 in3 = ixc(4,i)
172 in4 = ixc(5,i)
173C
174 n = iel(i + nft)
175 DO j=1,npt
176 ip = istack(i,j)
177 plysky(ip)%FSKY(1,iadc(1,n))= -(fly11(i,j) )
178 plysky(ip)%FSKY(2,iadc(1,n))= -(fly21(i,j) )
179 plysky(ip)%FSKY(3,iadc(1,n))= -(fly31(i,j) )
180C stifness by ply
181 plysky(ip)%FSKY(4,iadc(1,n))= sti(i,j)*fac(i,1)
182C
183 plysky(ip)%FSKY(1,iadc(2,n))= -(fly12(i,j) )
184 plysky(ip)%FSKY(2,iadc(2,n))= -(fly22(i,j) )
185 plysky(ip)%FSKY(3,iadc(2,n))= -(fly32(i,j) )
186C stifness by ply
187 plysky(ip)%FSKY(4,iadc(2,n))= sti(i,j)*fac(i,2)
188c
189 plysky(ip)%FSKY(1,iadc(3,n))= -(fly13(i,j))
190 plysky(ip)%FSKY(2,iadc(3,n))= -(fly23(i,j) )
191 plysky(ip)%FSKY(3,iadc(3,n))= -(fly33(i,j) )
192C stifness by ply
193 plysky(ip)%FSKY(4,iadc(3,n))= sti(i,j)*fac(i,1)
194
195 plysky(ip)%FSKY(1,iadc(4,n))= -(fly14(i,j) )
196 plysky(ip)%FSKY(2,iadc(4,n))= -(fly24(i,j) )
197 plysky(ip)%FSKY(3,iadc(4,n))= -(fly34(i,j) )
198C stifness by ply
199 plysky(ip)%FSKY(4,iadc(4,n))= sti(i,j)*fac(i,2)
200C
201 ENDDO
202 ENDDO
203C
204 ENDIF
205C
206 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
plyxfem_mod
Definition plyxfem_mod.F:31
plyxfem_mod::ply
type(ply_data), dimension(:), allocatable ply
Definition plyxfem_mod.F:91
plyxfem_mod::plysky
type(ply_data), dimension(:), allocatable plysky
Definition plyxfem_mod.F:91