OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fsigcini.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!|| fsigcini ../starter/source/constraints/fxbody/fsigcini.F
25!||--- called by ------------------------------------------------------
26!|| fxbsini ../starter/source/constraints/fxbody/fxbsini.F
27!||--- calls -----------------------------------------------------
28!|| ccoefi ../starter/source/constraints/fxbody/fsigcini.F
29!|| ccurvi ../starter/source/constraints/fxbody/fsigcini.F
30!|| cdefli ../starter/source/constraints/fxbody/fsigcini.f
31!|| ceveci ../starter/source/elements/shell/coque/ceveci.F
32!|| cm1inif ../starter/source/elements/shell/coque/cm1inif.f
33!|| cpxpyi ../starter/source/elements/shell/coque/cepsini.F
34!||====================================================================
35 SUBROUTINE fsigcini(FXBELM, IPARG , X , PM, IXC ,
36 . GEO , FXBMOD, FXBSIG, R , NELC)
37C-----------------------------------------------
38C I m p l i c i t T y p e s
39C-----------------------------------------------
40#include "implicit_f.inc"
41C-----------------------------------------------
42C G l o b a l P a r a m e t e r s
43C-----------------------------------------------
44#include "mvsiz_p.inc"
45C-----------------------------------------------
46C C o m m o n B l o c k s
47C-----------------------------------------------
48#include "param_c.inc"
49C-----------------------------------------------
50C D u m m y A r g u m e n t s
51C-----------------------------------------------
52 INTEGER FXBELM(*), IPARG(NPARG,*), IXC(NIXC,*), NELC
53 my_real
54 . fxbsig(*), x(3,*), pm(npropm), fxbmod(*),
55 . geo(npropg,*), r(3,*)
56C-----------------------------------------------
57C L o c a l V a r i a b l e s
58C-----------------------------------------------
59 INTEGER IG,OFFSET,FIRST,LAST,NFT,I,NG,IEL,
60 . n1,n2,n3,n4,ism,ithk,npt,nfs
61 INTEGER MAT(MVSIZ), PROP(MVSIZ)
62 INTEGER II,J
63 my_real
64 . d11, d12, d13, d21, d22, d23, d31, d32, d33, d41, d42, d43,
65 . dr11, dr12, dr13, dr21, dr22, dr23, dr31, dr32, dr33,
66 . dr41, dr42, dr43
67 my_real vl(3,4,mvsiz), vrl(3,4,mvsiz),
68 . px1g(mvsiz), px2g(mvsiz), py1g(mvsiz), py2g(mvsiz),
69 . px1(mvsiz) , px2(mvsiz) , py1(mvsiz) , py2(mvsiz),
70 . x2s(mvsiz), y2s(mvsiz), x3s(mvsiz),
71 . y3s(mvsiz), x4s(mvsiz), y4s(mvsiz),
72 . x1(mvsiz) , x2(mvsiz) , x3(mvsiz) , x4(mvsiz) ,
73 . y1(mvsiz) , y2(mvsiz) , y3(mvsiz) , y4(mvsiz) ,
74 . z1(mvsiz) , z2(mvsiz) , z3(mvsiz) , z4(mvsiz),
75 . e1x(mvsiz), e1y(mvsiz), e1z(mvsiz),
76 . e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
77 . e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),
78 . exx(mvsiz),eyy(mvsiz),exy(mvsiz),eyz(mvsiz),ezx(mvsiz),
79 . kxx(mvsiz),kyy(mvsiz),kxy(mvsiz),
80 . gstrbid(8,mvsiz), for(5,mvsiz), mom(3,mvsiz),
81 . eint(2,mvsiz) , thk(mvsiz), area(mvsiz),
82 . nu(mvsiz), g(mvsiz), a1(mvsiz), a2(mvsiz), gs(mvsiz)
83C=======================================================================
84 first=1
85 DO ig=1,nelc,mvsiz
86 offset=ig-1
87 last=min(mvsiz,nelc-offset)
88 nft=offset*10
89 nfs=offset*10
90 DO i=1,last
91 ng=fxbelm(nft+10*(i-1)+1)
92 iel=iparg(3,ng)+fxbelm(nft+10*(i-1)+2)
93 mat(i)=ixc(1,iel)
94 prop(i)=ixc(6,iel)
95 thk(i)=geo(1,prop(i))
96 x1(i)=0.
97 y1(i)=0.
98 z1(i)=0.
99 x2(i)=x(1,ixc(3,iel))-x(1,ixc(2,iel))
100 y2(i)=x(2,ixc(3,iel))-x(2,ixc(2,iel))
101 z2(i)=x(3,ixc(3,iel))-x(3,ixc(2,iel))
102 x3(i)=x(1,ixc(4,iel))-x(1,ixc(2,iel))
103 y3(i)=x(2,ixc(4,iel))-x(2,ixc(2,iel))
104 z3(i)=x(3,ixc(4,iel))-x(3,ixc(2,iel))
105 x4(i)=x(1,ixc(5,iel))-x(1,ixc(2,iel))
106 y4(i)=x(2,ixc(5,iel))-x(2,ixc(2,iel))
107 z4(i)=x(3,ixc(5,iel))-x(3,ixc(2,iel))
108 n1=fxbelm(nft+10*(i-1)+3)
109 n2=fxbelm(nft+10*(i-1)+4)
110 n3=fxbelm(nft+10*(i-1)+5)
111 n4=fxbelm(nft+10*(i-1)+6)
112 d11=fxbmod(6*(n1-1)+1)
113 d12=fxbmod(6*(n1-1)+2)
114 d13=fxbmod(6*(n1-1)+3)
115 d21=fxbmod(6*(n2-1)+1)
116 d22=fxbmod(6*(n2-1)+2)
117 d23=fxbmod(6*(n2-1)+3)
118 d31=fxbmod(6*(n3-1)+1)
119 d32=fxbmod(6*(n3-1)+2)
120 d33=fxbmod(6*(n3-1)+3)
121 d41=fxbmod(6*(n4-1)+1)
122 d42=fxbmod(6*(n4-1)+2)
123 d43=fxbmod(6*(n4-1)+3)
124 vl(1,1,i)=r(1,1)*d11+r(1,2)*d12+r(1,3)*d13
125 vl(2,1,i)=r(2,1)*d11+r(2,2)*d12+r(2,3)*d13
126 vl(3,1,i)=r(3,1)*d11+r(3,2)*d12+r(3,3)*d13
127 vl(1,2,i)=r(1,1)*d21+r(1,2)*d22+r(1,3)*d23
128 vl(2,2,i)=r(2,1)*d21+r(2,2)*d22+r(2,3)*d23
129 vl(3,2,i)=r(3,1)*d21+r(3,2)*d22+r(3,3)*d23
130 vl(1,3,i)=r(1,1)*d31+r(1,2)*d32+r(1,3)*d33
131 vl(2,3,i)=r(2,1)*d31+r(2,2)*d32+r(2,3)*d33
132 vl(3,3,i)=r(3,1)*d31+r(3,2)*d32+r(3,3)*d33
133 vl(1,4,i)=r(1,1)*d41+r(1,2)*d42+r(1,3)*d43
134 vl(2,4,i)=r(2,1)*d41+r(2,2)*d42+r(2,3)*d43
135 vl(3,4,i)=r(3,1)*d41+r(3,2)*d42+r(3,3)*d43
136 dr11=fxbmod(6*(n1-1)+4)
137 dr12=fxbmod(6*(n1-1)+5)
138 dr13=fxbmod(6*(n1-1)+6)
139 dr21=fxbmod(6*(n2-1)+4)
140 dr22=fxbmod(6*(n2-1)+5)
141 dr23=fxbmod(6*(n2-1)+6)
142 dr31=fxbmod(6*(n3-1)+4)
143 dr32=fxbmod(6*(n3-1)+5)
144 dr33=fxbmod(6*(n3-1)+6)
145 dr41=fxbmod(6*(n4-1)+4)
146 dr42=fxbmod(6*(n4-1)+5)
147 dr43=fxbmod(6*(n4-1)+6)
148 vrl(1,1,i)=r(1,1)*dr11+r(1,2)*dr12+r(1,3)*dr13
149 vrl(2,1,i)=r(2,1)*dr11+r(2,2)*dr12+r(2,3)*dr13
150 vrl(3,1,i)=r(3,1)*dr11+r(3,2)*dr12+r(3,3)*dr13
151 vrl(1,2,i)=r(1,1)*dr21+r(1,2)*dr22+r(1,3)*dr23
152 vrl(2,2,i)=r(2,1)*dr21+r(2,2)*dr22+r(2,3)*dr23
153 vrl(3,2,i)=r(3,1)*dr21+r(3,2)*dr22+r(3,3)*dr23
154 vrl(1,3,i)=r(1,1)*dr31+r(1,2)*dr32+r(1,3)*dr33
155 vrl(2,3,i)=r(2,1)*dr31+r(2,2)*dr32+r(2,3)*dr33
156 vrl(3,3,i)=r(3,1)*dr31+r(3,2)*dr32+r(3,3)*dr33
157 vrl(1,4,i)=r(1,1)*dr41+r(1,2)*dr42+r(1,3)*dr43
158 vrl(2,4,i)=r(2,1)*dr41+r(2,2)*dr42+r(2,3)*dr43
159 vrl(3,4,i)=r(3,1)*dr41+r(3,2)*dr42+r(3,3)*dr43
160 gstrbid(1,i)=zero
161 gstrbid(2,i)=zero
162 gstrbid(3,i)=zero
163 for(1,i)=zero
164 for(2,i)=zero
165 for(3,i)=zero
166 for(4,i)=zero
167 for(5,i)=zero
168 mom(1,i)=zero
169 mom(2,i)=zero
170 mom(3,i)=zero
171 px1g(i)=zero
172 px2g(i)=zero
173 py1g(i)=zero
174 py2g(i)=zero
175 ism=1
176 ithk=0
177 ENDDO
178C
179 CALL ceveci(first, last , area,
180 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
181 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
182 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z )
183 npt=1
184 CALL ccoefi(last, pm , geo, nu ,
185 . g , a1 , a2 , gs , thk,
186 . mat , prop, npt, area)
187 CALL cpxpyi(first,last, ism,
188 . px1g ,px2g ,py1g ,py2g ,area ,
189 . px1 ,px2 ,py1 ,py2 ,
190 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
191 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
192 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z ,
193 . x2s ,y2s ,x3s ,y3s ,x4s ,y4s )
194 CALL cdefli(last ,vl ,gstrbid,
195 . px1 ,px2 ,py1 ,py2 ,
196 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z ,
197 . exx ,eyy ,exy ,eyz ,ezx ,area )
198 CALL ccurvi(last ,vrl ,
199 . px1 ,px2 ,py1 ,py2 ,
200 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z ,
201 . eyz ,ezx ,kxx ,kyy ,kxy ,area )
202 ithk=0
203 CALL cm1inif(first,last ,for ,mom ,ithk ,
204 . thk ,eint ,nu ,g ,a1 ,
205 . a2 ,gs ,kxx ,kyy ,kxy ,
206 . exx ,eyy ,exy ,eyz ,ezx ,
207 . area )
208C
209 DO i=1,last
210 fxbsig(nfs+10*(i-1)+1)=for(1,i)
211 fxbsig(nfs+10*(i-1)+2)=for(2,i)
212 fxbsig(nfs+10*(i-1)+3)=for(3,i)
213 fxbsig(nfs+10*(i-1)+4)=for(4,i)
214 fxbsig(nfs+10*(i-1)+5)=for(5,i)
215 fxbsig(nfs+10*(i-1)+6)=mom(1,i)
216 fxbsig(nfs+10*(i-1)+7)=mom(2,i)
217 fxbsig(nfs+10*(i-1)+8)=mom(3,i)
218 fxbsig(nfs+10*(i-1)+9)=eint(1,i)
219 fxbsig(nfs+10*(i-1)+10)=eint(2,i)
220 ENDDO
221 ENDDO
222C-----------
223 RETURN
224 END SUBROUTINE fsigcini
225!||====================================================================
226!|| ccoefi ../starter/source/constraints/fxbody/fsigcini.F
227!||--- called by ------------------------------------------------------
228!|| fsigcini ../starter/source/constraints/fxbody/fsigcini.F
229!|| fsigtini ../starter/source/constraints/fxbody/fsigtini.F
230!||====================================================================
231 SUBROUTINE ccoefi(NEL , PM , GEO, NU ,
232 . G , A1 , A2 , GS , THK,
233 . MAT , PROP, NPT, AREA )
234C-----------------------------------------------
235C I m p l i c i t T y p e s
236C-----------------------------------------------
237#include "implicit_f.inc"
238C-----------------------------------------------
239C G l o b a l P a r a m e t e r s
240C-----------------------------------------------
241#include "mvsiz_p.inc"
242C-----------------------------------------------
243C C o m m o n B l o c k s
244C-----------------------------------------------
245#include "param_c.inc"
246C-----------------------------------------------
247C D u m m y A r g u m e n t s
248C-----------------------------------------------
249 INTEGER :: NEL, MAT(*), PROP(*), NPT
250 my_real AREA(MVSIZ),
251 . PM(NPROPM,*), GEO(NPROPG,*), NU(*), G(*), A1(*),
252 . a2(*), gs(*), thk(*)
253C-----------------------------------------------
254C L o c a l V a r i a b l e s
255C-----------------------------------------------
256 INTEGER I,IM,IP,ISH
257 my_real THK02,FAC1,FSH,SHF
258C=======================================================================
259 DO i=1,nel
260 im=mat(i)
261 ip=prop(i)
262 nu(i)=pm(21,im)
263 g(i)=pm(22,im)
264 a1(i)=pm(24,im)
265 a2(i)=pm(25,im)
266 IF (npt > 1) THEN
267 thk02=thk(i)*thk(i)
268 fac1=two*(one+nu(i))*thk02
269 ish=nint(geo(37,ip))
270 fsh=geo(38,ip)
271 shf=fsh*(one-ish+ish*fac1/(fsh*area(i)+fac1))
272 ELSE
273 shf=zero
274 ENDIF
275 gs(i)=g(i)*shf
276 ENDDO
277C
278 RETURN
279 END SUBROUTINE ccoefi
280!||====================================================================
281!|| cdefli ../starter/source/constraints/fxbody/fsigcini.F
282!||--- called by ------------------------------------------------------
283!|| fsigcini ../starter/source/constraints/fxbody/fsigcini.F
284!||====================================================================
285 SUBROUTINE cdefli (NEL ,VL ,GSTR ,
286 . PX1 ,PX2 ,PY1 ,PY2 ,
287 . E1X ,E2X ,E3X ,E1Y ,E2Y ,E3Y ,E1Z ,E2Z ,E3Z ,
288 . EXX ,EYY ,EXY ,EYZ ,EZX ,AREA )
289C-----------------------------------------------
290C I m p l i c i t T y p e s
291C-----------------------------------------------
292#include "implicit_f.inc"
293C-----------------------------------------------
294C G l o b a l P a r a m e t e r s
295C-----------------------------------------------
296#include "mvsiz_p.inc"
297C-----------------------------------------------
298C D u m m y A r g u m e n t s
299C-----------------------------------------------
300 INTEGER :: NEL
301 my_real VL(3,4,*),GSTR(8,*),PX1(*) ,PX2(*) ,PY1(*) ,PY2(*)
302 my_real EXX(MVSIZ),EYY(MVSIZ),EXY(MVSIZ),EZX(MVSIZ),EYZ(MVSIZ),
303 . E1X(MVSIZ) , E1Y(MVSIZ) , E1Z(MVSIZ) ,
304 . E2X(MVSIZ) , E2Y(MVSIZ) , E2Z(MVSIZ) ,
305 . E3X(MVSIZ) , E3Y(MVSIZ) , E3Z(MVSIZ) ,AREA(MVSIZ)
306C-----------------------------------------------
307C L o c a l V a r i a b l e s
308C-----------------------------------------------
309 INTEGER I
310 my_real FAC1
311 my_real
312 . VX1(MVSIZ) , VX2(MVSIZ) , VX3(MVSIZ) , VX4(MVSIZ) ,
313 . VY1(MVSIZ) , VY2(MVSIZ) , VY3(MVSIZ) , VY4(MVSIZ) ,
314 . VZ1(MVSIZ) , VZ2(MVSIZ) , VZ3(MVSIZ) , VZ4(MVSIZ) ,
315 . vx13(mvsiz), vx24(mvsiz), vy13(mvsiz), vy24(mvsiz),
316 . vz13(mvsiz), vz24(mvsiz)
317C=======================================================================
318 DO i=1,nel
319 vx1(i)=e1x(i)*vl(1,1,i)+e1y(i)*vl(2,1,i)+e1z(i)*vl(3,1,i)
320 vx2(i)=e1x(i)*vl(1,2,i)+e1y(i)*vl(2,2,i)+e1z(i)*vl(3,2,i)
321 vx3(i)=e1x(i)*vl(1,3,i)+e1y(i)*vl(2,3,i)+e1z(i)*vl(3,3,i)
322 vx4(i)=e1x(i)*vl(1,4,i)+e1y(i)*vl(2,4,i)+e1z(i)*vl(3,4,i)
323C
324 vy4(i)=e2x(i)*vl(1,4,i)+e2y(i)*vl(2,4,i)+e2z(i)*vl(3,4,i)
325 vy3(i)=e2x(i)*vl(1,3,i)+e2y(i)*vl(2,3,i)+e2z(i)*vl(3,3,i)
326 vy2(i)=e2x(i)*vl(1,2,i)+e2y(i)*vl(2,2,i)+e2z(i)*vl(3,2,i)
327 vy1(i)=e2x(i)*vl(1,1,i)+e2y(i)*vl(2,1,i)+e2z(i)*vl(3,1,i)
328C
329 vz1(i)=e3x(i)*vl(1,1,i)+e3y(i)*vl(2,1,i)+e3z(i)*vl(3,1,i)
330 vz2(i)=e3x(i)*vl(1,2,i)+e3y(i)*vl(2,2,i)+e3z(i)*vl(3,2,i)
331 vz3(i)=e3x(i)*vl(1,3,i)+e3y(i)*vl(2,3,i)+e3z(i)*vl(3,3,i)
332 vz4(i)=e3x(i)*vl(1,4,i)+e3y(i)*vl(2,4,i)+e3z(i)*vl(3,4,i)
333 ENDDO
334C
335 DO i=1,nel
336 vz13(i)=vz1(i)-vz3(i)
337 vz24(i)=vz2(i)-vz4(i)
338C
339 vx13(i)=vx1(i)-vx3(i)
340 vx24(i)=vx2(i)-vx4(i)
341C
342 exx(i)=px1(i)*vx13(i)+px2(i)*vx24(i)
343 exy(i)=py1(i)*vx13(i)+py2(i)*vx24(i)
344C
345 vy13(i)=vy1(i)-vy3(i)
346 vy24(i)=vy2(i)-vy4(i)
347C
348 exy(i)=exy(i)+px1(i)*vy13(i)+px2(i)*vy24(i)
349 eyy(i)=py1(i)*vy13(i)+py2(i)*vy24(i)
350 eyz(i)=py1(i)*vz13(i)+py2(i)*vz24(i)
351 ezx(i)=px1(i)*vz13(i)+px2(i)*vz24(i)
352 ENDDO
353C
354 DO i=1,nel
355 fac1 = one/area(i)
356 exx(i)=exx(i)*fac1
357 eyy(i)=eyy(i)*fac1
358 exy(i)=exy(i)*fac1
359 ENDDO
360C
361 DO i=1,nel
362 gstr(1,i)=gstr(1,i)+exx(i)
363 gstr(2,i)=gstr(2,i)+eyy(i)
364 gstr(3,i)=gstr(3,i)+exy(i)
365 ENDDO
366 RETURN
367 END SUBROUTINE cdefli
368!||====================================================================
369!|| ccurvi ../starter/source/constraints/fxbody/fsigcini.F
370!||--- called by ------------------------------------------------------
371!|| fsigcini ../starter/source/constraints/fxbody/fsigcini.F
372!||====================================================================
373 SUBROUTINE ccurvi(NEL ,VRL ,
374 . PX1 ,PX2 ,PY1 ,PY2 ,
375 . E1X ,E2X ,E3X ,E1Y ,E2Y ,E3Y ,E1Z ,E2Z ,E3Z ,
376 . EYZ ,EZX ,KXX ,KYY ,KXY ,AREA )
377C-----------------------------------------------
378C I m p l i c i t T y p e s
379C-----------------------------------------------
380#include "implicit_f.inc"
381C-----------------------------------------------
382C G l o b a l P a r a m e t e r s
383C-----------------------------------------------
384#include "mvsiz_p.inc"
385C-----------------------------------------------
386C D u m m y A r g u m e n t s
387C-----------------------------------------------
388 INTEGER :: NEL
389 my_real VRL(3,4,*),PX1(*),PX2(*),PY1(*),PY2(*)
390 my_real E1X(MVSIZ), E1Y(MVSIZ), E1Z(MVSIZ), E2X(MVSIZ),
391 . E2Y(MVSIZ), E2Z(MVSIZ), E3X(MVSIZ), E3Y(MVSIZ),
392 . E3Z(MVSIZ), KXX(MVSIZ), KYY(MVSIZ), KXY(MVSIZ),
393 . AREA(MVSIZ), EZX(MVSIZ), EYZ(MVSIZ)
394C-----------------------------------------------
395C L o c a l V a r i a b l e s
396C-----------------------------------------------
397 INTEGER I
398 my_real
399 . rx1(mvsiz), rx2(mvsiz), rx3(mvsiz), rx4(mvsiz),
400 . ry1(mvsiz), ry2(mvsiz), ry3(mvsiz), ry4(mvsiz)
401 my_real
402 . rxavta,rx13ta,rx24ta,ry13ta,ryavta,ry24ta,fac1
403C=======================================================================
404 DO i=1,nel
405 rx1(i)=e1x(i)*vrl(1,1,i)+e1y(i)*vrl(2,1,i)+e1z(i)*vrl(3,1,i)
406 rx2(i)=e1x(i)*vrl(1,2,i)+e1y(i)*vrl(2,2,i)+e1z(i)*vrl(3,2,i)
407 rx3(i)=e1x(i)*vrl(1,3,i)+e1y(i)*vrl(2,3,i)+e1z(i)*vrl(3,3,i)
408 rx4(i)=e1x(i)*vrl(1,4,i)+e1y(i)*vrl(2,4,i)+e1z(i)*vrl(3,4,i)
409 ry1(i)=e2x(i)*vrl(1,1,i)+e2y(i)*vrl(2,1,i)+e2z(i)*vrl(3,1,i)
410 ry2(i)=e2x(i)*vrl(1,2,i)+e2y(i)*vrl(2,2,i)+e2z(i)*vrl(3,2,i)
411 ry3(i)=e2x(i)*vrl(1,3,i)+e2y(i)*vrl(2,3,i)+e2z(i)*vrl(3,3,i)
412 ry4(i)=e2x(i)*vrl(1,4,i)+e2y(i)*vrl(2,4,i)+e2z(i)*vrl(3,4,i)
413 ENDDO
414C
415 DO i=1,nel
416 rx13ta =rx1(i)-rx3(i)
417 rxavta =rx1(i)+rx2(i)+rx3(i)+rx4(i)
418 rx24ta =rx2(i)-rx4(i)
419C
420 kyy(i)=-py1(i)*rx13ta-py2(i)*rx24ta
421 kxy(i)= px1(i)*rx13ta+px2(i)*rx24ta
422C
423 ry13ta = ry1(i)-ry3(i)
424 ryavta = ry1(i)+ry2(i)+ry3(i)+ry4(i)
425 ry24ta = ry2(i)-ry4(i)
426C
427 kxx(i)= px1(i)*ry13ta+px2(i)*ry24ta
428 kxy(i)= py1(i)*ry13ta+py2(i)*ry24ta
429 + -kxy(i)
430C
431 ezx(i)=ezx(i)+ryavta*(fourth*area(i))
432 eyz(i)=eyz(i)-rxavta*(fourth*area(i))
433 ENDDO
434C
435 DO i=1,nel
436 fac1 = one/area(i)
437 ezx(i)=ezx(i)*fac1
438 eyz(i)=eyz(i)*fac1
439 kxx(i)=kxx(i)*fac1
440 kyy(i)=kyy(i)*fac1
441 kxy(i)=kxy(i)*fac1
442 ENDDO
443C
444 RETURN
445 END SUBROUTINE ccurvi
446
cpxpyi
subroutine cpxpyi(jft, jlt, ismstr, area, px1g, px2g, py1g, py2g, px1, px2, py1, py2, x1g, x2g, x3g, x4g, y1g, y2g, y3g, y4g, z1g, z2g, z3g, z4g, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, x2l, y2l, x3l, y3l, x4l, y4l)
Definition cepsini.F:350
ceveci
subroutine ceveci(jft, jlt, area, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition ceveci.F:37
cm1inif
subroutine cm1inif(jft, jlt, for, mom, ithk, thk, eint, nu, g, a1, a2, gs, kxx, kyy, kxy, exx, eyy, exy, eyz, exz, area)
Definition cm1inif.F:34
my_real
#define my_real
Definition cppsort.cpp:32
cdefli
subroutine cdefli(nel, vl, gstr, px1, px2, py1, py2, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, exx, eyy, exy, eyz, ezx, area)
Definition fsigcini.F:289
ccurvi
subroutine ccurvi(nel, vrl, px1, px2, py1, py2, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, eyz, ezx, kxx, kyy, kxy, area)
Definition fsigcini.F:377
ccoefi
subroutine ccoefi(nel, pm, geo, nu, g, a1, a2, gs, thk, mat, prop, npt, area)
Definition fsigcini.F:234
fsigcini
subroutine fsigcini(fxbelm, iparg, x, pm, ixc, geo, fxbmod, fxbsig, r, nelc)
Definition fsigcini.F:37
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
min
#define min(a, b)
Definition macros.h:20
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
starter
program starter
Definition starter.F:39