OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cbadef_ply.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!|| cbadef_ply ../engine/source/properties/composite_options/stack/cbadef_ply.f
25!||--- called by ------------------------------------------------------
26!|| cbaforc3 ../engine/source/elements/shell/coqueba/cbaforc3.F
27!||====================================================================
28 SUBROUTINE cbadef_ply(JFT,JLT,NG, NPT,NPLAT,IPLAT,VQ,
29 1 VXYZ,VETA,VKSI,BM,BC,TC,DT1C,
30 2 EXX, EYY , EYZ, EZX)
31C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
32C CALCUL DES DEFORMATIONS GENERALISEES AUX POINTS DE GAUSS ET [B]
33C ENTREES : NEL,NG,VQN,VXYZ,VASM,VASTN
34C SORTIES : VDEF,BM (MEMBRANE),BMF(COUPLAGE MEM-FLEXION),BF (FLEXION),BC(C.T)
35C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
36C-----------------------------------------------
37C I M P L I C I T T Y P E S
38C-----------------------------------------------
39#include "implicit_f.inc"
40C-----------------------------------------------
41C G L O B A L P A R A M E T E R S
42C-----------------------------------------------
43#include "mvsiz_p.inc"
44C-----------------------------------------------
45C D U M M Y A R G U M E N T S
46C-----------------------------------------------
47 INTEGER NG,JFT,JLT,NPLAT,IPLAT(*),NPT
48C PARAMETER (NPG = 4)
49 my_real
50 . vxyz(mvsiz,12,npt),vksi(4,4),veta(4,4),
51 . bm(mvsiz,36),bc(mvsiz,40)
52 my_real
53 . vq(mvsiz,3,3,4),dt1c(*),
54 . exx(mvsiz,npt), eyy(mvsiz,npt) , eyz(mvsiz,npt),
55 . ezx(mvsiz,npt)
56 my_real
57 . tc(mvsiz,2,2)
58C-----------------------------------------------
59C L O C A L V A R I A B L E S
60C-----------------------------------------------
61 INTEGER PT,PT00,PT0,I,J,EP,IUN,NG1,K,NPLAT0
62 MY_REAL
63 . thk,detj,det,
64 . tfn(3,2),bcx,bcy,
65 . v1(2),v2(2),rv1,rv2,
66 . c1,c2,vt1,vt2,bc1,bc2,vb1,
67 . vjf1(2,3),vjf(3,3),tbi(2,2),tbc(2,2)
68 my_real
69 . a_1,c11,c12,c21,c22,cc,beta1,ksi1,ksiy1,beta2,ksi2,ksiy2
70c PARAMETER (PG=.577350269189626)
71c PARAMETER (PG1=-.577350269189626)
72C--------------------------
73C INITIALISATION
74C--------------------------
75 DATA iun/1/
76C
77 nplat0 = nplat
78 DO j= 1, npt
79#include "vectorize.inc"
80 DO i=jft,nplat0
81 ep=iplat(i)
82C
83 exx(ep,j) =
84 . bm(ep,1)*vxyz(ep,1,j) + bm(ep,2)*vxyz(ep,4 , j)
85 . + bm(ep,3)*vxyz(ep,7 , j)
86 eyy(ep,j)=
87 . bm(ep,5)*vxyz(ep,2,j) + bm(ep,6)*vxyz(ep,5 , j)
88 . + bm(ep,7)*vxyz(ep,8 , j)
89C
90!! EYZ(EP,J)=
91 ezx(ep,j)=
92 . bc(ep,1)*vxyz(ep,3,j) + bc(ep,7 )*vxyz(ep,6 , j)
93 . + bc(ep,13)*vxyz(ep,9 , j)
94!! EZX(EP,J)=
95 eyz(ep,j)=
96 . bc(ep,2)*vxyz(ep,3,j) + bc(ep,8 )*vxyz(ep,6 , j)
97 . + bc(ep,14)*vxyz(ep,9 , j)
98 exx(ep,j) = exx(ep,j)*dt1c(ep)
99 eyy(ep,j) = eyy(ep,j)*dt1c(ep)
100 eyz(ep,j) = eyz(ep,j)*dt1c(ep)
101 ezx(ep,j) = ezx(ep,j)*dt1c(ep)
102 ENDDO
103#include "vectorize.inc"
104 DO 150 i=nplat0+1,jlt
105 ep=iplat(i)
106C----------VT1=<T1>*VXYZ AU CHAQUE NOEUD---------
107 vt1=
108 . vq(ep,1,1,ng)*vxyz(ep,1, j) + vq(ep,2,1,ng)*vxyz(ep,2, j)
109 . + vq(ep,3,1,ng)*vxyz(ep,3, j)
110C----------VT2=<T2>*VXYZ AU CHAQUE NOEUD---------
111 vt2=
112 . vq(ep,1,2,ng)*vxyz(ep,1, j) + vq(ep,2,2,ng)*vxyz(ep,2, j)
113 . + vq(ep,3,2,ng)*vxyz(ep,3, j)
114C--------TERMES DE [B0]=<T>CI-------
115 c1=vksi(1,ng)*tc(ep,1,1) + veta(1,ng)*tc(ep,2,1)
116 c2=vksi(1,ng)*tc(ep,1,2) + veta(1,ng)*tc(ep,2,2)
117C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
118c BC1=VKSI(1,NG)*TBC(1,1)+VETA(1,NG)*TBC(2,1)
119c BC2=VKSI(1,NG)*TBC(1,2)+VETA(1,NG)*TBC(2,2)
120
121C---------DEF MEMBRANE------
122 exx(ep ,j)= c1*vt1
123 eyy(ep ,j)= c2*vt2
124C
125C---------FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
126C----------[B01]---------
127c VDEF(6,EP)= BC1*VT1
128c VDEF(7,EP)= BC2*VT2
129c VDEF(8,EP)= BC1*VT2+BC2*VT1
130C--------J=2---II=(J-1)*2 =2 KK = 3*(J-1) =3
131C----------VT1=<T1>*VXYZ AU CHAQUE NOEUD---------
132 vt1 = vq(ep,1,1,ng)*vxyz(ep,4,j)+vq(ep,2,1,ng)*vxyz(ep,5,j)
133 1 +vq(ep,3,1,ng)*vxyz(ep,6,j)
134C----------VT2=<T2>*VXYZ AU CHAQUE NOEUD---------
135 vt2 = vq(ep,1,2,ng)*vxyz(ep,4,j)+vq(ep,2,2,ng)*vxyz(ep,5,j)
136 1 +vq(ep,3,2,ng)*vxyz(ep,6,j)
137C--------- [B0], [B01]---------
138C--------TERMES DE [B0]=<T>CI-------
139 c1=vksi(2,ng)*tc(ep,1,1)+veta(2,ng)*tc(ep,2,1)
140 c2=vksi(2,ng)*tc(ep,1,2)+veta(2,ng)*tc(ep,2,2)
141C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
142c BC1=VKSI(2,NG)*TBC(1,1)+VETA(2,NG)*TBC(2,1)
143c BC2=VKSI(2,NG)*TBC(1,2)+VETA(2,NG)*TBC(2,2)
144C
145C I = 1; JJ = 3*KK =9
146C---------DEF MEMBRANE------
147 exx(ep,j) = exx(ep,j) + c1*vt1
148 eyy(ep,j) = eyy(ep,j) + c2*vt2
149C---------DEF DE FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
150C----------[B01]---------
151c VDEF(6,EP)= VDEF(6,EP)+ BC1*VT1
152cc VDEF(7,EP)= VDEF(7,EP)+ BC2*VT2
153c VDEF(8,EP)= VDEF(8,EP)+ BC1*VT2+BC2*VT1
154C---------- TERMES ASSOCIEES AU BETA ----------
155C
156C--------J=3---II=(J-1)*2 =4 KK = 3*(J-1)=6
157C----------VT1=<T1>*VXYZ AU CHAQUE NOEUD---------
158 vt1=
159 . vq(ep,1,1,ng)*vxyz(ep,7,j) + vq(ep,2,1,ng)*vxyz(ep,8,j)
160 1 + vq(ep,3,1,ng)*vxyz(ep,9,j)
161C----------VT2=<T2>*VXYZ AU CHAQUE NOEUD---------
162 vt2=
163 . vq(ep,1,2,ng)*vxyz(ep,7,j) + vq(ep,2,2,ng)*vxyz(ep,8,j)
164 1 + vq(ep,3,2,ng)*vxyz(ep,9,j)
165C--------- [B0], [B01]---------
166C--------TERMES DE [B0]=<T>CI-------
167 c1=vksi(3,ng)*tc(ep,1,1)+veta(3,ng)*tc(ep,2,1)
168 c2=vksi(3,ng)*tc(ep,1,2)+veta(3,ng)*tc(ep,2,2)
169C--------TERMES ASSOCIEES AUX U,V,W DE [B1]=<T>BCI-------
170c BC1=VKSI(3,NG)*TBC(1,1)+VETA(3,NG)*TBC(2,1)
171c BC2=VKSI(3,NG)*TBC(1,2)+VETA(3,NG)*TBC(2,2)
172C I = 1;JJ = 3*KK=18
173C---------:EPS-X:EPS-Y:EPS-XY
174C---------DEF MEMBRANE------
175 exx(ep,j) = exx(ep,j) + c1*vt1
176 eyy(ep,j) = eyy(ep,j) + c2*vt2
177C
178C---------DEF DE FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
179C----------[B01]---------
180c VDEF(6,EP)= VDEF(6,EP)+ BC1*VT1
181c VDEF(7,EP)= VDEF(7,EP)+ BC2*VT2
182c VDEF(8,EP)= VDEF(8,EP)+ BC1*VT2+BC2*VT1
183C---------- TERMES ASSOCIEES AU BETA ----------
184C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
185C----------VT1=<T1>*VXYZ AU CHAQUE NOEUD---------
186 vt1=
187 . vq(ep,1,1,ng)*vxyz(ep,10,j) + vq(ep,2,1,ng)*vxyz(ep,11,j)
188 1 + vq(ep,3,1,ng)*vxyz(ep,12,j)
189C----------VT2=<T2>*VXYZ AU CHAQUE NOEUD---------
190 vt2=
191 . vq(ep,1,2,ng)*vxyz(ep,10,j) + vq(ep,2,2,ng)*vxyz(ep,11,j)
192 1 + vq(ep,3,2,ng)*vxyz(ep,12,j)
193C--------- [B0], [B01]---------
194C--------TERMES DE [B0]=<T>CI-------
195 c1=vksi(4,ng)*tc(ep,1,1)+veta(4,ng)*tc(ep,2,1)
196 c2=vksi(4,ng)*tc(ep,1,2)+veta(4,ng)*tc(ep,2,2)
197C
198C---------DEF MEMBRANE------
199 exx(ep,j)= exx(ep,j)+ c1*vt1
200 eyy(ep,j)= eyy(ep,j)+ c2*vt2
201C---------DEF DE FLEXION ASSOCIEE AUX U,V,W ([BMF]=2H[B0]+[B01])------
202C----------SEULEMENT LA PARTIE DE [B01]---------
203c VDEF(6,EP)= VDEF(6,EP)+ BC1*VT1
204c VDEF(7,EP)= VDEF(7,EP)+ BC2*VT2
205c VDEF(8,EP)= VDEF(8,EP)+ BC1*VT2+BC2*VT1
206C---------- TERMES ASSOCIEES AU BETA ----------
207C---------PRODUIT SCALAIRE <T(PG)>*<-T^2 T^1>----------
208C--------------------------
209C DEFORMATION OUT-PLANE (C.T)
210C--------------------------
211C [BC](2,5*NPG),VDEF(4,5)
212C--------------------------
213C
214 bcx= bc(ep,1)*vxyz(ep,1,j) + bc(ep,3)*vxyz(ep,2,j)
215 1 + bc(ep,5)*vxyz(ep,3,j)
216 bcy= bc(ep,2)*vxyz(ep,1,j) + bc(ep,4)*vxyz(ep,2,j)
217 1 + bc(ep,6)*vxyz(ep,3,j)
218C------- ---J=2 II=(J-1)*3 =3 JJ=(J-1)*2 =2 KK=10*(J-1) =10
219 bcx=bcx
220 1 + bc(ep,11)*vxyz(ep,4,j) + bc(ep,13)*vxyz(ep,5,j)
221 2 + bc(ep,15)*vxyz(ep,6,j)
222 bcy=bcy
223 1 + bc(ep,12)*vxyz(ep,4,j) + bc(ep,14)*vxyz(ep,5,j)
224 2 + bc(ep,16)*vxyz(ep,6,j)
225C------- ---J=3-- II=(J-1)*3=6 JJ=(J-1)*2=4 KK=10*(J-1)=20
226 bcx=bcx
227 1 + bc(ep,21)*vxyz(ep,7,j) + bc(ep,23)*vxyz(ep,8,j)
228 2 + bc(ep,25)*vxyz(ep,9,j)
229 bcy=bcy
230 1 +bc(ep,22)*vxyz(ep,7,j) + bc(ep,24)*vxyz(ep,8,j)
231 2 +bc(ep,26)*vxyz(ep,9,j)
232C------- ---J=4-- II=(J-1)*3 =9 JJ=(J-1)*2 =6 KK=10*(J-1)=30
233 bcx=bcx
234 1 +bc(ep,31)*vxyz(ep,10,j) + bc(ep,33)*vxyz(ep,11,j)
235 2 +bc(ep,35)*vxyz(ep,12,j)
236 bcy=bcy
237 1 +bc(ep,32)*vxyz(ep,10,j) + bc(ep,34)*vxyz(ep,11,j)
238 2 +bc(ep,36)*vxyz(ep,12,j)
239
240 ezx(ep,j)= tc(ep,1,1)*bcx + tc(ep,2,1)*bcy
241 eyz(ep,j)= tc(ep,1,2)*bcx + tc(ep,2,2)*bcy
242C
243 exx(ep,j) = exx(ep,j)*dt1c(ep)
244 eyy(ep,j) = eyy(ep,j)*dt1c(ep)
245 eyz(ep,j) = eyz(ep,j)*dt1c(ep)
246 ezx(ep,j) = ezx(ep,j)*dt1c(ep)
247C
248 150 CONTINUE
249 ENDDO
250 RETURN
251 END
252C
253!||====================================================================
254!|| cbadefsh_ply ../engine/source/properties/composite_options/stack/cbadef_ply.F
255!||--- called by ------------------------------------------------------
256!|| cbaforc3 ../engine/source/elements/shell/coqueba/cbaforc3.F
257!||====================================================================
258 SUBROUTINE cbadefsh_ply(JFT,JLT,NPT ,NPLAT, IPLAT,X13 ,X24 ,
259 . Y13 ,Y24, VXYZ,DT1C, EXY )
260C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
261C CALCUL 'membrane shear traitement' use only PARTIE CONSTANTE
262C-----------------------------------------------
263C I M P L I C I T T Y P E S
264C-----------------------------------------------
265#include "implicit_f.inc"
266C-----------------------------------------------
267C G L O B A L P A R A M E T E R S
268C-----------------------------------------------
269#include "mvsiz_p.inc"
270C-----------------------------------------------
271C D U M M Y A R G U M E N T S
272C-----------------------------------------------
273 INTEGER JFT ,JLT,NPLAT,IPLAT(*),NPT
274 my_real
275 . X13(*) ,X24(*) ,Y13(*) ,Y24(*) ,
276 . exy(mvsiz,npt) ,vxyz(mvsiz,12,npt),dt1c(*)
277C-----------------------------------------------
278C L O C A L V A R I A B L E S
279C-----------------------------------------------
280 INTEGER I ,J,EP,NPLAT0
281C
282!! NPLAT0= JLT
283 nplat0= nplat
284 DO j=1,npt
285 DO ep=jft,nplat0
286 i=iplat(ep)
287 exy(i,j)= y24(i)*vxyz(i,2,j) - y13(i)*vxyz(i,5,j)
288 . -x24(i)*vxyz(i,1,j) + x13(i)*vxyz(i,4,j)
289 exy(i,j) = exy(i,j)*dt1c(i)
290 ENDDO
291 DO ep=nplat0+1,jlt
292 i=iplat(ep)
293 exy(i,j) = y24(i)*( vxyz(i,2,j) - vxyz(i,8 ,j))
294 1 + y13(i)*(-vxyz(i,5,j) + vxyz(i,11 ,j))
295 2 - x24(i)*( vxyz(i,1,j) - vxyz(i,7 ,j))
296 3 + x13(i)*( vxyz(i,4,j) - vxyz(i,10 ,j))
297 exy(i,j) = exy(i,j)*dt1c(i)
298 ENDDO
299 ENDDO
300 RETURN
301 END
302
cbadefsh_ply
subroutine cbadefsh_ply(jft, jlt, npt, nplat, iplat, x13, x24, y13, y24, vxyz, dt1c, exy)
Definition cbadef_ply.F:260
cbadef_ply
subroutine cbadef_ply(jft, jlt, ng, npt, nplat, iplat, vq, vxyz, veta, vksi, bm, bc, tc, dt1c, exx, eyy, eyz, ezx)
Definition cbadef_ply.F:31
my_real
#define my_real
Definition cppsort.cpp:32