OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cbafint_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!|| cbafint_ply ../engine/source/properties/composite_options/stack/cbafint_ply.F
25!||--- called by ------------------------------------------------------
26!|| cbaforc3 ../engine/source/elements/shell/coqueba/cbaforc3.F
27!||--- uses -----------------------------------------------------
28!|| element_mod ../common_source/modules/elements/element_mod.F90
29!||====================================================================
30 SUBROUTINE cbafint_ply(JFT,JLT,NPT,NG,NPLAT,IPLAT,CDET,THKLY,TH12,
31 . VOL,FF0,BM,BC,TC,SIG_IPLY,VNI,AREA,
32 . VF,VFI,IXC)
33 use element_mod , only : nixc
34C--------------------------------------------------------------------------------------------------
35C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
36C CALCUL DES FORCES INTERNES AUX NOEUDS
37C ENTREES : NEL,MNT,NPT,NP,NIXC,IXC,HZ,OFF,
38C CDET,EPAIS,BUFCON,BM,,,BC
39C SORTIES : VF
40C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
41C--------------------------------------------------------------------------------------------------
42C-----------------------------------------------
43C I M P L I C I T T Y P E S
44C-----------------------------------------------
45#include "implicit_f.inc"
46#include "mvsiz_p.inc"
47C-----------------------------------------------
48C D U M M Y A R G U M E N T S
49C-----------------------------------------------
50 INTEGER JFT,JLT,NPT,NPLAT,IPLAT(*),NG,IXC(NIXC,*)
51 my_real
52c . FF(5,MVSIZ,NPT),CDET(*),VOL(MVSIZ,NPT),
53 . ff0(mvsiz,5,npt),cdet(*),vol(*),
54 . bm(mvsiz,36),bc(mvsiz,40),
55 . thkly(mvsiz,*),vf(mvsiz,12,npt),th12(*),
56 . tc(mvsiz,2,2),sig_iply(mvsiz,3,npt),
57 . area(mvsiz),vni(4,4),
58 . vfi(mvsiz,12,npt)
59C-----------------------------------------------
60C L O C A L V A R I A B L E S
61C-----------------------------------------------
62 INTEGER EP,I,J,JSUP,JINF,IPLY,NPLAT0
63 my_real
64 . c1,c2,bc1,bc3,bc5,bcx,bcy,fl1,fl2,fl3,
65 . ff(mvsiz,5,npt)
66
67C--------------------------------------------------------------------------------------------------$
68!! NPLAT0 = JLT
69 nplat0 = nplat
70 DO j=1,npt
71#include "vectorize.inc"
72 DO i=jft,nplat0
73 ep=iplat(i)
74 ff(ep,1,j)=ff0(ep,1,j)
75 ff(ep,2,j)=ff0(ep,2,j)
76 ff(ep,3,j)=ff0(ep,3,j)
77 ff(ep,5,j)=ff0(ep,4,j)
78 ff(ep,4,j)=ff0(ep,5,j)
79C
80cc C2=VOL(EP)/NPT
81!! C2=VOL(EP)
82 c2=thkly(ep,j)*cdet(ep)
83 c1=th12(ep)*cdet(ep)
84 bc1=bc(ep,1)*ff(ep,4,j)
85 bc3=bc(ep,3)*ff(ep,4,j)
86 bc5=bc(ep,5)*ff(ep,4,j)
87C----------NOEUD 1 ------------
88 vf(ep,1,j) = vf(ep,1,j) + c2*(bm(ep,1)*ff(ep,1,j))
89 vf(ep,2,j) = vf(ep,2,j) + c2*(bm(ep,5)*ff(ep,2,j))
90 vf(ep,3,j) = vf(ep,3,j) + c2*(bc1+bc(ep,2)*ff(ep,5,j))
91C complete integration of the shear
92cc VF(1,EP,J) = VF(1,EP,J) + C2*(BM(EP,5)*FF(3,EP,J))
93cc VF(2,EP,J) = VF(2,EP,J) + C2*(BM(EP,1)*FF(3,EP,J))
94
95C
96 vf(ep,7,j) = vf(ep,7,j) + c2*(bm(ep,3)*ff(ep,1,j))
97 vf(ep,8,j) = vf(ep,8,j) + c2*(bm(ep,7)*ff(ep,2,j))
98C complete integration of the shear
99cc VF(7,EP,J) = VF(7,EP,J) + C2*(BM(7,EP)*FF(3,EP,J))
100cc VF(8,EP,J) = VF(8,EP,J) + C2*(BM(3,EP)*FF(3,EP,J))
101C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
102 vf(ep,9,j) = vf(ep,9,j) + c2*(bc(ep,13)*ff(ep,4,j)
103 . + bc(ep,14)*ff(ep,5,j))
104C----------NOEUD 2 ------------
105 vf(ep,4,j) = vf(ep,4,j) + c2*(bm(ep,2)*ff(ep,1,j))
106 vf(ep,5,j) = vf(ep,5,j) + c2*(bm(ep,6)*ff(ep,2,j))
107 vf(ep,6,j) = vf(ep,6,j) + c2*(bc(ep,7)*ff(ep,4,j)
108 . + bc(ep,8)*ff(ep,5,j))
109C complete integration of the shear
110cc VF(4,EP,J) = VF(4,EP,J) + C2*(BM(EP,6)*FF(3,EP,J))
111cc VF(5,EP,J) = VF(5,EP,J) + C2*(BM(EP,2)*FF(3,EP,J))
112C
113 vf(ep,10,j) =-vf(ep,7,j)
114 vf(ep,11,j) =-vf(ep,8,j)
115 vf(ep,12,j) =-vf(ep,9,j)
116 ENDDO
117#include "vectorize.inc"
118 DO i=nplat0+1,jlt
119 ep=iplat(i)
120C
121 ff(ep,1,j)=ff0(ep,1,j)
122 ff(ep,2,j)=ff0(ep,2,j)
123 ff(ep,3,j)=ff0(ep,3,j)
124 ff(ep,5,j)=ff0(ep,4,j)
125 ff(ep,4,j)=ff0(ep,5,j)
126C
127!! C2=VOL(EP)
128 c2=thkly(ep,j)*cdet(ep)
129 c1=th12(ep)*cdet(ep)
130 bcx=tc(ep,1,1)*ff(ep,4,j)+tc(ep,1,2)*ff(ep,5,j)
131 bcy=tc(ep,2,1)*ff(ep,4,j)+tc(ep,2,2)*ff(ep,5,j)
132C J=1 K0=0 K1=0 K2=0
133 vf(ep,1,j) = vf(ep,1,j) + c2*(bm(ep,1)*ff(ep,1,j)+
134 1 bm(ep,2)*ff(ep,2,j)+
135 3 bc(ep,1)*bcx+
136 4 bc(ep,2)*bcy )
137c 5 +C1*(BMF(EP,1)*MM(1,EP)+
138c 6 BMF(EP,2)*MM(2,EP)+
139c 7 BMF(EP,3)*MM(3,EP))
140
141 vf(ep,2,j) = vf(ep,2,j )+ c2*(bm(ep,4)*ff(ep,1,j)+
142 2 bm(ep,5)*ff(ep,2,j)+
143 4 bc(ep,3)*bcx+
144 5 bc(ep,4)*bcy )
145c 6 +C1*(BMF(EP,4,EP)*MM(1,EP)+
146c 7 BMF(EP,5,EP)*MM(2,EP)+
147c 8 BMF(EP,6,EP)*MM(3,EP))
148
149C I1=3 I2=3
150 vf(ep,3,j)=vf(ep,3,j)+ c2*(bm(ep,7)*ff(ep,1,j)+
151 2 bm(ep,8)*ff(ep,2,j)+
152 4 bc(ep,5)*bcx+
153 5 bc(ep,6)*bcy )
154c 6 +C1*(BMF(EP,7)*MM(1,EP)+
155c 7 BMF(EP,8)*MM(2,EP)+
156c 8 BMF(EP,9)*MM(3,EP))
157C J=2 K0=2 K1=3 K2=5 I1=4 I2=6
158 vf(ep,4,j) = vf(ep,4,j)+c2*(bm(ep,10)*ff(ep,1,j)+
159 2 bm(ep,11)*ff(ep,2,j)+
160 4 bc(ep,11)*bcx+
161 5 bc(ep,12)*bcy )
162c 6 +C1*(BMF(EP,10)*MM(1,EP)+
163c 7 BMF(EP,11)*MM(2,EP)+
164c 8 BMF(EP,12)*MM(3,EP))
165
166C I1=5 I2=7
167 vf(ep,5,j)=vf(ep,5,j)+c2*(bm(ep,13)*ff(ep,1,j)+
168 2 bm(ep,14)*ff(ep,2,j)+
169 4 bc(ep,13)*bcx+
170 5 bc(ep,14)*bcy )
171c 6 +C1*(BMF(EP,13)*MM(1,EP)+
172c 7 BMF(EP,14)*MM(2,EP)+
173c 8 BMF(EP,15)*MM(3,EP))
174C I1=6 I2=8
175 vf(ep,6,j)=vf(ep,6,j)+c2*(bm(ep,16)*ff(ep,1,j)+
176 2 bm(ep,17)*ff(ep,2,j)+
177 4 bc(ep,15)*bcx+
178 5 bc(ep,16)*bcy )
179c 6 +C1*(BMF(EP,16)*MM(1,EP)+
180c 7 BMF(EP,17)*MM(2,EP)+
181c 8 BMF(EP,18)*MM(3,EP))
182C J=3 K0=4 K1=6 K2=10 I1=7 I2=11
183 vf(ep,7,j)=vf(ep,7,j)+c2*(bm(ep,19)*ff(ep,1,j)+
184 2 bm(ep,20)*ff(ep,2,j)+
185 4 bc(ep,21)*bcx+
186 5 bc(ep,22)*bcy )
187c 6 +C1*(BMF(EP,19)*MM(1,EP)+
188c 7 BMF(EP,20)*MM(2,EP)+
189c 8 BMF(EP,21)*MM(3,EP))
190C I1=8 I2=12
191 vf(ep,8,j)=vf(ep,8,j)+c2*(bm(ep,22)*ff(ep,1,j)+
192 2 bm(ep,23)*ff(ep,2,j)+
193 4 bc(ep,23)*bcx+
194 5 bc(ep,24)*bcy )
195c 6 +C1*(BMF(EP,22)*MM(1,EP)+
196c 7 BMF(EP,23)*MM(2,EP)+
197c 8 BMF(EP,24)*MM(3,EP))
198C I1=9 I2=13
199 vf(ep,9,j)=vf(ep,9,j)+c2*(bm(ep,25)*ff(ep,1,j)+
200 2 bm(ep,26)*ff(ep,2,j)+
201 4 bc(ep,25)*bcx+
202 5 bc(ep,26)*bcy )
203cc 6 +C1*(BMF(EP,25)*MM(1,EP)+
204c 7 BMF(EP,26)*MM(2,EP)+
205c 8 BMF(EP,27)*MM(3,EP))
206C J=4 K0=6 K1=9 K2=15 I=1 I1=10 I2=16
207 vf(ep,10,j)=vf(ep,10,j)+c2*(bm(ep,28)*ff(ep,1,j)+
208 2 bm(ep,29)*ff(ep,2,j)+
209 4 bc(ep,31)*bcx+
210 5 bc(ep,32)*bcy )
211c 6 +C1*(BMF(EP,28)*MM(1,EP)+
212c 7 BMF(EP,29)*MM(2,EP)+
213c 8 BMF(EP,30)*MM(3,EP))
214C I1=11 I2=17
215 vf(ep,11,j)=vf(ep,11,j)+c2*(bm(ep,31)*ff(ep,1,j)+
216 2 bm(ep,32)*ff(ep,2,j)+
217 4 bc(ep,33)*bcx+
218 5 bc(ep,34)*bcy )
219c 6 +C1*(BMF(EP,31)*MM(1,EP)+
220c 7 BMF(EP,32)*MM(2,EP)+
221c 8 BMF(EP,33)*MM(3,EP))
222C I1=12 I2=18
223 vf(ep,12,j)=vf(ep,12,j)+c2*(bm(ep,34)*ff(ep,1,j)+
224 2 bm(ep,35)*ff(ep,2,j)+
225 4 bc(ep,35)*bcx+
226 5 bc(ep,36)*bcy )
227c 6 +C1*(BMF(EP,34)*MM(1,EP)+
228c 7 BMF(EP,35)*MM(2,EP)+
229c 8 BMF(EP,36)*MM(3,EP))
230C
231 ENDDO
232 ENDDO
233C
234C interface
235 DO iply = 1 , npt - 1
236 jinf = iply
237 jsup = iply + 1
238 DO ep=jft,jlt
239C
240 fl1 = cdet(ep)*sig_iply(ep,2,iply)
241 fl2 = cdet(ep)*sig_iply(ep,1,iply)
242 fl3 = cdet(ep)*sig_iply(ep,3,iply)
243C----------NOEUD 1 ------------
244 vfi(ep,1,jsup) = vfi(ep,1,jsup) + fl1*vni(1,ng)
245 vfi(ep,2,jsup) = vfi(ep,2,jsup) + fl2*vni(1,ng)
246 vfi(ep,3,jsup) = vfi(ep,3,jsup) + fl3*vni(1,ng)
247C ----------NOEUD 2 ------------
248 vfi(ep,4,jsup) = vfi(ep,4,jsup) + fl1*vni(2,ng)
249 vfi(ep,5,jsup) = vfi(ep,5,jsup) + fl2*vni(2,ng)
250 vfi(ep,6,jsup) = vfi(ep,6,jsup) + fl3*vni(2,ng)
251C----------NOEUD 3 ------------
252 vfi(ep,7,jsup) = vfi(ep,7,jsup) + fl1*vni(3,ng)
253 vfi(ep,8,jsup) = vfi(ep,8,jsup) + fl2*vni(3,ng)
254 vfi(ep,9,jsup) = vfi(ep,9,jsup) + fl3*vni(3,ng)
255C----------NOEUD 4 ------------
256 vfi(ep,10,jsup) = vfi(ep,10,jsup) + fl1*vni(4,ng)
257 vfi(ep,11,jsup) = vfi(ep,11,jsup) + fl2*vni(4,ng)
258 vfi(ep,12,jsup) = vfi(ep,12,jsup) + fl3*vni(4,ng)
259cc -----------ply inf
260C----------NOEUD 1 ------------
261 vfi(ep,1,jinf) = vfi(ep,1,jinf) - fl1*vni(1,ng)
262 vfi(ep,2,jinf) = vfi(ep,2,jinf) - fl2*vni(1,ng)
263 vfi(ep,3,jinf) = vfi(ep,3,jinf) - fl3*vni(1,ng)
264C ----------NOEUD 2 ------------
265 vfi(ep,4,jinf) = vfi(ep,4,jinf) - fl1*vni(2,ng)
266 vfi(ep,5,jinf) = vfi(ep,5,jinf) - fl2*vni(2,ng)
267 vfi(ep,6,jinf) = vfi(ep,6,jinf) - fl3*vni(2,ng)
268C----------NOEUD 3 ------------
269 vfi(ep,7,jinf) = vfi(ep,7,jinf) - fl1*vni(3,ng)
270 vfi(ep,8,jinf) = vfi(ep,8,jinf) - fl2*vni(3,ng)
271 vfi(ep,9,jinf) = vfi(ep,9,jinf) - fl3*vni(3,ng)
272C----------NOEUD 4 ------------
273 vfi(ep,10,jinf) = vfi(ep,10,jinf) - fl1*vni(4,ng)
274 vfi(ep,11,jinf) = vfi(ep,11,jinf) - fl2*vni(4,ng)
275 vfi(ep,12,jinf) = vfi(ep,12,jinf) - fl3*vni(4,ng)
276 ENDDO
277 ENDDO
278 RETURN
279 END
cbafint_ply
subroutine cbafint_ply(jft, jlt, npt, ng, nplat, iplat, cdet, thkly, th12, vol, ff0, bm, bc, tc, sig_iply, vni, area, vf, vfi, ixc)
Definition cbafint_ply.F:33
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567