OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
czbe3.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!|| czbe3 ../engine/source/elements/shell/coquez/czbe3.F
25!||--- called by ------------------------------------------------------
26!|| czke3 ../engine/source/elements/shell/coquez/czke3.f
27!||====================================================================
28 SUBROUTINE czbe3(JFT ,JLT ,AREA ,A_I ,X13 ,
29 2 X24 ,Y13 ,Y24 ,MX13 ,MX23 ,
30 3 MX34 ,MY13 ,MY23 ,MY34 ,Z1 ,
31 4 PX1 ,PX2 ,PY1 ,PY2 ,RX ,
32 5 RY ,SX ,SY ,RX2 ,RY2 ,
33 6 SX2 ,SY2 ,RHX ,RHY ,SHX ,
34 7 SHY ,PH1 ,PH2 ,HXX ,HYY ,
35 8 HXY ,NPLAT,IPLAT)
36C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
37#include "implicit_f.inc"
38C-----------------------------------------------
39C G L O B A L P A R A M E T E R S
40C-----------------------------------------------
41#include "mvsiz_p.inc"
42C-----------------------------------------------
43C D U M M Y A R G U M E N T S
44C-----------------------------------------------
45 INTEGER JFT,JLT,NPLAT,IPLAT(*)
46 my_real
47 . X13(*),X24(*),Y13(*),Y24(*),AREA(*),Z1(*),A_I(*),
48 . MX13(*),MX23(*),MX34(*),MY13(*),MY23(*),MY34(*)
49 my_real
50 . PX1(*),PX2(*),PY1(*),PY2(*),PH1(*),PH2(*),
51 . RX(*),RY(*), SX(*),SY(*),RX2(*),RY2(*), SX2(*),SY2(*),
52 . rhx(mvsiz,4),rhy(mvsiz,4), shx(mvsiz,4),shy(mvsiz,4),hxx(*),hyy(*),hxy(*)
53C-----------------------------------------------
54C L O C A L V A R I A B L E S
55C-----------------------------------------------
56 INTEGER I,J,EP
57C-----------gama(I)=hI/4-PH(I), PH:anti-sym comme bxI------
58#include "vectorize.inc"
59 DO I=jft,jlt
60 ep=iplat(i)
61 px1(i) = y24(ep)*a_i(ep)
62 px2(i) = -y13(ep)*a_i(ep)
63 py1(i) = -x24(ep)*a_i(ep)
64 py2(i) = x13(ep)*a_i(ep)
65 ph1(i) = mx13(ep)*px1(i)+ my13(ep)*py1(i)
66 ph2(i) = mx13(ep)*px2(i)+ my13(ep)*py2(i)
67 ENDDO
68C-----------Bc = Bcr0(1+r*rI)+Bcs0(1+s*sI)------
69C-----------rhx(I)=rkxk+hkxk*SI---shx(I)=skxk+hkxk*RI-
70C----------shear----R : -1 1 1 -1; S: -1 -1 1 1--
71#include "vectorize.inc"
72 DO i=jft,jlt
73 ep=iplat(i)
74 sx(i) = mx34(ep)*a_i(ep)
75 sy(i) = my34(ep)*a_i(ep)
76 rx(i) = mx23(ep)*a_i(ep)
77 ry(i) = my23(ep)*a_i(ep)
78 rhx(i,1) = mx23(ep)-mx13(ep)
79 rhx(i,3) = mx23(ep)+mx13(ep)
80 rhy(i,1) = my23(ep)-my13(ep)
81 rhy(i,3) = my23(ep)+my13(ep)
82 shx(i,1) = mx34(ep)-mx13(ep)
83 shx(i,2) = mx34(ep)+mx13(ep)
84 shy(i,1) = my34(ep)-my13(ep)
85 shy(i,2) = my34(ep)+my13(ep)
86 ENDDO
87C
88 DO i=jft,jlt
89 sx2(i) = sx(i)*sx(i)
90 sy2(i) = sy(i)*sy(i)
91 rx2(i) = rx(i)*rx(i)
92 ry2(i) = ry(i)*ry(i)
93 rhx(i,2) = rhx(i,1)
94 rhx(i,4) = rhx(i,3)
95 rhy(i,2) = rhy(i,1)
96 rhy(i,4) = rhy(i,3)
97 shx(i,3) = shx(i,2)
98 shx(i,4) = shx(i,1)
99 shy(i,3) = shy(i,2)
100 shy(i,4) = shy(i,1)
101 hxx(i)=fivep333*(sy2(i)+ry2(i))
102 hyy(i)=fivep333*(sx2(i)+rx2(i))
103 hxy(i)=-fivep333*(sx(i)*sy(i)+rx(i)*ry(i))
104 ENDDO
105C
106 RETURN
107 END
108!||====================================================================
109!|| czber3 ../engine/source/elements/shell/coquez/czbe3.F
110!||--- called by ------------------------------------------------------
111!|| czke3 ../engine/source/elements/shell/coquez/czke3.f
112!||--- calls -----------------------------------------------------
113!|| czderirz ../engine/source/elements/shell/coquez/czdef.f
114!||====================================================================
115 SUBROUTINE czber3(JFT ,JLT ,AREA ,A_I ,X13 ,
116 1 X24 ,Y13 ,Y24 ,MX13 ,MX23 ,
117 2 MX34 ,MY13 ,MY23 ,MY34 ,Z1 ,
118 3 RX ,RY ,SX ,SY ,PRX ,
119 4 PRY ,PRXY ,PRZ ,PHKRX,PHKRY,
120 5 PHKRXY,PHERX,PHERY,PHERXY,
121 6 PHKRZ,PHERZ ,PHKX ,PHKY ,PHEX ,
122 7 PHEY ,IPLAT)
123C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
124#include "implicit_f.inc"
125#include "mvsiz_p.inc"
126C-----------------------------------------------
127C D U M M Y A R G U M E N T S
128C-----------------------------------------------
129 INTEGER JFT,JLT,IPLAT(*)
130 my_real
131 . X13(*),X24(*),Y13(*),Y24(*),AREA(*),Z1(*),A_I(*),
132 . MX13(*),MX23(*),MX34(*),MY13(*),MY23(*),MY34(*),
133 . RX(*),RY(*), SX(*),SY(*)
134 my_real
135 . PRX(4,*),PRY(4,*),PRXY(4,*),PHKRX(4,*),PHKRY(4,*),
136 . PHKRXY(4,*),PHERX(4,*),PHERY(4,*),PHERXY(4,*),PRZ(4,*),
137 . PHKRZ(4,*),PHERZ(4,*),PHKX(*) ,PHKY(*) ,PHEX(*) ,PHEY(*)
138C-----------------------------------------------
139C L O C A L V A R I A B L E S
140C-----------------------------------------------
141 INTEGER I,J,EP
142 my_real
143 . bm0rz(mvsiz,4,4),bmkrz(mvsiz,4,4),bmerz(mvsiz,4,4)
144C---------------
145 CALL czderirz(jft ,jlt ,area ,x13 ,x24 ,
146 2 y13 ,y24 ,mx13 ,mx23,mx34 ,
147 3 my13 ,my23 ,my34 ,bm0rz,bmkrz,
148 4 bmerz )
149#include "vectorize.inc"
150 DO i=jft,jlt
151 ep=iplat(i)
152 DO j=1,4
153 prx(j,i)= bm0rz(ep,1,j)*a_i(ep)
154 pry(j,i)= bm0rz(ep,2,j)*a_i(ep)
155 prxy(j,i)=bm0rz(ep,3,j)*a_i(ep)
156C
157 phkrx(j,i)= bmkrz(ep,1,j)*a_i(ep)
158 phkry(j,i)= bmkrz(ep,2,j)*a_i(ep)
159 phkrxy(j,i)=bmkrz(ep,3,j)*a_i(ep)
160C
161 pherx(j,i)= bmerz(ep,1,j)*a_i(ep)
162 phery(j,i)= bmerz(ep,2,j)*a_i(ep)
163 pherxy(j,i)=bmerz(ep,3,j)*a_i(ep)
164C
165 prz(j,i)=bm0rz(ep,4,j)*a_i(ep)
166 ENDDO
167 ENDDO
168C
169 DO j=1,4
170 DO i=jft,jlt
171 ep=iplat(i)
172 phkrz(j,i)=bmkrz(ep,4,j)*a_i(ep)
173 pherz(j,i)=bmerz(ep,4,j)*a_i(ep)
174 ENDDO
175 ENDDO
176C
177 DO i=jft,jlt
178 phex(i)= four*sy(i)
179 phey(i) =-four*sx(i)
180 phkx(i)=-four*ry(i)
181 phky(i)= four*rx(i)
182 ENDDO
183C
184 RETURN
185 END
czbe3
subroutine czbe3(jft, jlt, area, a_i, x13, x24, y13, y24, mx13, mx23, mx34, my13, my23, my34, z1, px1, px2, py1, py2, rx, ry, sx, sy, rx2, ry2, sx2, sy2, rhx, rhy, shx, shy, ph1, ph2, hxx, hyy, hxy, nplat, iplat)
Definition czbe3.F:36
czber3
subroutine czber3(jft, jlt, area, a_i, x13, x24, y13, y24, mx13, mx23, mx34, my13, my23, my34, z1, rx, ry, sx, sy, prx, pry, prxy, prz, phkrx, phkry, phkrxy, pherx, phery, pherxy, phkrz, pherz, phkx, phky, phex, phey, iplat)
Definition czbe3.F:123
czdef
subroutine czdef(jft, jlt, area, aa, v13, v24, vhi, rlxyz, vdef, vhg, x13, x24, y13, y24, mx13, mx23, mx34, my13, my23, my34, z1, dt1, off, offg, rlxyzv)
Definition czdef.F:34
czderirz
subroutine czderirz(jft, jlt, area, x13, x24, y13, y24, mx13, mx23, mx34, my13, my23, my34, bm0rz, bmkrz, bmerz)
Definition czdef.F:302
czke3
subroutine czke3(jft, jlt, nft, npt, mtn, ithk, ncycle, istrain, ipla, pm, geo, ixc, elbuf_str, bufmat, offset, indxof, etag, iddl, ndof, k_diag, k_lt, iadk, jdik, ihbe, thke, ismstr, x, ikgeo, ipm, igeo, iexpan, iparg, isubstack, stack, drape_sh4n, indx_drape, sedrape, numel_drape)
Definition czke3.F:58