OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
czke3.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!|| czke3 ../engine/source/elements/shell/coquez/czke3.F
25!||--- called by ------------------------------------------------------
26!|| imp_glob_k ../engine/source/implicit/imp_glob_k.F
27!|| imp_glob_k0 ../engine/source/implicit/imp_glob_k.F
28!||--- calls -----------------------------------------------------
29!|| assem_c4 ../engine/source/implicit/assem_c4.F
30!|| c4eoff ../engine/source/elements/shell/coque/c4eoff.F
31!|| cmatc3 ../engine/source/elements/shell/coqueba/cmatc3.F
32!|| cmatch3 ../engine/source/elements/shell/coqueba/cmatc3.F
33!|| cmatip3 ../engine/source/elements/shell/coqueba/cmatc3.F
34!|| czbe3 ../engine/source/elements/shell/coquez/czbe3.F
35!|| czber3 ../engine/source/elements/shell/coquez/czbe3.f
36!|| czcoork3 ../engine/source/elements/shell/coquez/czcoork3.F
37!|| czlkec3 ../engine/source/elements/shell/coquez/czlkec3.F
38!|| czlkecg3 ../engine/source/elements/shell/coquez/czlkecg3.F
39!|| czlkecr3 ../engine/source/elements/shell/coquez/czlkec3.F
40!|| czlkect3 ../engine/source/elements/shell/coquez/czlkect3.F
41!|| czlken3 ../engine/source/elements/shell/coquez/czlken3.F
42!|| czlkenr3 ../engine/source/elements/shell/coquez/czlken3.F
43!|| czsumg3 ../engine/source/elements/shell/coquez/czsumg3.F
44!||--- uses -----------------------------------------------------
45!|| drape_mod ../engine/share/modules/drape_mod.F
46!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
47!|| stack_mod ../engine/share/modules/stack_mod.F
48!||====================================================================
49 SUBROUTINE czke3 (
50 1 JFT ,JLT ,NFT ,NPT ,MTN ,
51 2 ITHK ,NCYCLE ,ISTRAIN,IPLA ,PM ,
52 3 GEO ,IXC ,ELBUF_STR ,BUFMAT ,
53 4 OFFSET ,INDXOF ,ETAG , IDDL ,NDOF ,
54 5 K_DIAG ,K_LT , IADK ,JDIK ,
55 6 IHBE ,THKE ,ISMSTR ,X ,IKGEO ,
56 7 IPM ,IGEO ,IEXPAN ,IPARG ,ISUBSTACK,
57 8 STACK ,DRAPE_SH4N , INDX_DRAPE,SEDRAPE,NUMEL_DRAPE)
58C-----------------------------------------------
59C M o d u l e s
60C-----------------------------------------------
61 USE elbufdef_mod
62 USE stack_mod
63 USE drape_mod
64C-----------------------------------------------
65C I M P L I C I T T Y P E S
66C-----------------------------------------------
67#include "implicit_f.inc"
68C-----------------------------------------------
69C G L O B A L P A R A M E T E R S
70C-----------------------------------------------
71#include "mvsiz_p.inc"
72C-----------------------------------------------
73C C O M M O N B L O C K S
74C-----------------------------------------------
75#include "com04_c.inc"
76#include "param_c.inc"
77#include "impl1_c.inc"
78C-----------------------------------------------
79C D U M M Y A R G U M E N T S
80C-----------------------------------------------
81 INTEGER JFT ,JLT ,NFT ,NPT ,
82 . MTN ,ITHK ,NCYCLE,ISUBSTACK,
83 . ISTRAIN ,IPLA ,OFFSET,IHBE ,ISMSTR,IKGEO,IEXPAN
84 INTEGER , INTENT(IN) :: SEDRAPE,NUMEL_DRAPE
85 INTEGER IXC(NIXC,*),IGEO(NPROPGI,*),IPM(*),IPARG(*)
86 INTEGER INDXOF(MVSIZ),
87 . ETAG(*),IDDL(*) ,NDOF(*) ,IADK(*) ,JDIK(*)
88 INTEGER, DIMENSION(SEDRAPE) :: INDX_DRAPE
89
90C REAL OU REAL*8
91 my_real
92 . pm(npropm,*),geo(npropg,*),bufmat(*),x(3,*),thke(*),
93 . off(mvsiz),k_diag(*) ,k_lt(*)
94 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
95 TYPE (STACK_PLY) :: STACK
96 TYPE (DRAPE_) :: DRAPE_SH4N(NUMELC_DRAPE)
97C=======================================================================
98c FUNCTION: [K] stiffness Matrix of QEPH element
99c
100c Note:
101c ARGUMENTS: (I: input, O: output, IO: input * output, W: workspace)
102c
103c TYPE NAME FUNCTION
104c I JFT,JLT,NFT - local element id limit and first id(global)of this element group
105c NEL=JLT-JFT+1
106c I NPT,MTN - num. of integrating point in thickness and material type id
107c I ITHK - flag of thickness updating (if >0)
108c I NCYCLE - cycle(increment) number
109c I ISTRAIN - total strain output flag
110c I IPLA - radial return plasticity compute option
111c I PM ,GEO - Material and geometrical property data
112c I IXC(NIXC,NEL) - connectivity and mid,pid integer data
113c I BUFMAT() - internal material data
114c I INDXOF(NEL) - deleted element flag (not used in this subroutine)
115c I ETAG(NEL) - activating element flag for Eigenvalue analysis
116c I IDDL(NUMNOD) - DOF id of node N =IDDL(N)+1,NDOF
117c I NDOF(NUMNOD) - Num of DOF (nodal)
118c IO K_DIAG(NDDL) - Diagnale terms of assembled [K] NDDL: total number of model DOF
119c IO K_LT(NNZK) - terms of strick triagular of assembled [K] NNZK: number of no-zero terms
120c I IADK(NDDL),JDIK(NNZK) - Indice integer tables of Compress format of [K]
121c I IHBE - Shell formulation flag (Ishell)
122c I THKE - initial thickness
123c I ISMSTR - small strain flag
124c I IKGEO , - geometrical stiffness matrix including (if >0)
125c I X(3,NUMNOD) co-ordinate
126c I IPM ,IGEO - Material and geometrical property integer data
127c I IEXPAN - thermo flag used in materials
128c I IPARG(NG) element group data
129C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
130C L O C A L V A R I A B L E S
131C--------------------------------
132 INTEGER
133 . I, J,J1,J2, NEL, NPLAT,IPLAT(MVSIZ), NLAY,L_DIRA,L_DIRB,
134 . IREP,IBID,EP,IDRIL,IBID1
135 INTEGER MAT(MVSIZ), PID(MVSIZ), NGL(MVSIZ),IORTH,IGTYP,IUN
136 my_real
137 . X13(MVSIZ), X24(MVSIZ), Y13(MVSIZ), Y24(MVSIZ),
138 . MX13(MVSIZ), MX23(MVSIZ), MX34(MVSIZ),
139 . MY13(MVSIZ), MY23(MVSIZ), MY34(MVSIZ), Z1(MVSIZ),
140 . px1(mvsiz), px2(mvsiz), py1(mvsiz),py2(mvsiz),
141 . sx(mvsiz), sy(mvsiz), rx(mvsiz),ry(mvsiz),
142 . sx2(mvsiz), sy2(mvsiz), rx2(mvsiz),ry2(mvsiz),
143 . rhx(mvsiz,4),rhy(mvsiz,4),shx(mvsiz,4),shy(mvsiz,4),
144 . ph1(mvsiz),ph2(mvsiz),hxx(mvsiz),hyy(mvsiz),hxy(mvsiz)
145 my_real
146 . vq(mvsiz,9),area(mvsiz), vqn(mvsiz,12),thk0(mvsiz),vol(mvsiz),
147 . a_i(mvsiz), thk2(mvsiz),hm(mvsiz,4),hf(mvsiz,4),hc(mvsiz,2),
148 . hz(mvsiz),dhz(mvsiz),hmor(mvsiz,2),hfor(mvsiz,2),
149 . gs(mvsiz),hmfor(mvsiz,6)
150 my_real
151 . corelv(mvsiz,2,4)
152 my_real
153 . k11(9,mvsiz),k12(9,mvsiz),k13(9,mvsiz),k14(9,mvsiz),
154 . k22(9,mvsiz),k23(9,mvsiz),k24(9,mvsiz),k33(9,mvsiz),
155 . m11(9,mvsiz),m12(9,mvsiz),m13(9,mvsiz),m14(9,mvsiz),
156 . m22(9,mvsiz),m23(9,mvsiz),m24(9,mvsiz),m33(9,mvsiz),
157 . mf11(9,mvsiz),mf12(9,mvsiz),mf13(9,mvsiz),mf14(9,mvsiz),
158 . mf22(9,mvsiz),mf23(9,mvsiz),mf24(9,mvsiz),mf33(9,mvsiz),
159 . fm12(9,mvsiz),fm13(9,mvsiz),fm14(9,mvsiz),
160 . fm23(9,mvsiz),fm24(9,mvsiz),fm34(9,mvsiz),
161 . k34(9,mvsiz),k44(9,mvsiz),m34(9,mvsiz),m44(9,mvsiz),
162 . mf34(9,mvsiz),mf44(9,mvsiz)
163 my_real
164 . prx(4,mvsiz),pry(4,mvsiz),prxy(4,mvsiz),phkrx(4,mvsiz),
165 . phkry(4,mvsiz),phkrxy(4,mvsiz),pherx(4,mvsiz),phery(4,mvsiz),
166 . pherxy(4,mvsiz),prz(4,mvsiz),phkrz(4,mvsiz),pherz(4,mvsiz),
167 . phkx(mvsiz),phky(mvsiz),phex(mvsiz),phey(mvsiz)
168 my_real
169 . ke11(36,mvsiz),ke22(36,mvsiz),ke33(36,mvsiz),ke44(36,mvsiz),
170 . ke12(36,mvsiz),ke13(36,mvsiz),ke14(36,mvsiz),ke23(36,mvsiz),
171 . ke24(36,mvsiz),ke34(36,mvsiz)
172C-----------------------------------------------
173 my_real,
174 . DIMENSION(:) ,POINTER :: dir_a, dir_b
175 my_real,
176 . ALLOCATABLE, DIMENSION(:), TARGET :: dira,dirb
177 TYPE(g_bufel_) ,POINTER :: GBUF
178C------------|---------|------------------------------------------------------
179C--------------------------
180C INITIALISATION
181C--------------------------
182C OPEN(UNIT=17,FILE='DEBZ.TMP',STATUS='UNKNOWN',FORM='FORMATTED')
183c
184 gbuf => elbuf_str%GBUF
185 nel=jlt-jft+1
186 idril = iparg(41)
187C
188 igtyp = igeo(11,ixc(6,1))
189 irep = igeo(6 ,ixc(6,1))
190 nlay = elbuf_str%NLAY
191 l_dira = elbuf_str%BUFLY(1)%LY_DIRA
192 l_dirb = elbuf_str%BUFLY(1)%LY_DIRB
193 ALLOCATE(dira(nlay*nel*l_dira))
194 ALLOCATE(dirb(nlay*nel*l_dirb))
195 dira = zero
196 dirb = zero
197 dir_a => dira(1:nlay*nel*l_dira)
198 dir_b => dirb(1:nlay*nel*l_dirb)
199 IF (irep == 0) THEN
200 DO j=1,nlay
201 j1 = 1+(j-1)*l_dira*nel
202 j2 = j*l_dira*nel
203 dira(j1:j2) = elbuf_str%BUFLY(j)%DIRA(1:nel*l_dira)
204 ENDDO
205 ENDIF
206C
207 CALL czcoork3(jft ,jlt ,x ,ixc ,pm ,
208 1 gbuf%OFF,area,a_i,vqn ,vq ,
209 2 x13 ,x24 ,y13 ,y24 ,mx13,
210 3 mx23,mx34 ,my13 ,my23 ,my34,
211 4 z1 , geo ,
212 5 elbuf_str,gbuf%SMSTR,nlay,
213 6 irep,npt,ismstr,
214 7 dir_a,dir_b,pid,mat,ngl,nplat,iplat ,
215 8 corelv,off,thke,nel)
216 IF (ikproj>0.OR.idril>0) THEN
217 DO i=1,9
218 DO ep=jft,jlt
219 m11(i,ep) =zero
220 m22(i,ep) =zero
221 m33(i,ep) =zero
222 m44(i,ep) =zero
223 m12(i,ep) =zero
224 m13(i,ep) =zero
225 m14(i,ep) =zero
226 m23(i,ep) =zero
227 m24(i,ep) =zero
228 m34(i,ep) =zero
229 mf11(i,ep) =zero
230 mf22(i,ep) =zero
231 mf33(i,ep) =zero
232 mf44(i,ep) =zero
233 mf12(i,ep) =zero
234 mf13(i,ep) =zero
235 mf14(i,ep) =zero
236 mf23(i,ep) =zero
237 mf24(i,ep) =zero
238 mf34(i,ep) =zero
239 fm12(i,ep) =zero
240 fm13(i,ep) =zero
241 fm14(i,ep) =zero
242 fm23(i,ep) =zero
243 fm24(i,ep) =zero
244 fm34(i,ep) =zero
245 ENDDO
246 ENDDO
247 ENDIF
248 IF (irep>0) THEN
249 CALL cmatc3(jft ,jlt ,pm ,mat ,geo ,
250 1 pid ,area ,thk0 ,thk2 ,gbuf%THK ,
251 2 thke ,vol ,mtn ,npt ,ithk ,
252 3 hm ,hf ,hc ,hz ,igtyp ,
253 4 iorth ,hmor ,hfor ,dir_a ,igeo ,
254 5 idril ,ihbe ,hmfor ,gs ,isubstack,
255 6 stack ,elbuf_str ,nlay ,drape_sh4n ,nft ,
256 7 nel ,indx_drape,sedrape,numel_drape)
257 ELSE
258 CALL cmatc3(jft ,jlt ,pm ,mat ,geo ,
259 1 pid ,area ,thk0 ,thk2 ,gbuf%THK ,
260 2 thke ,vol ,mtn ,npt ,ithk ,
261 3 hm ,hf ,hc ,hz ,igtyp ,
262 4 iorth ,hmor ,hfor ,dira ,igeo ,
263 5 idril ,ihbe ,hmfor ,gs ,isubstack,
264 6 stack ,elbuf_str ,nlay ,drape_sh4n ,nft ,
265 7 nel ,indx_drape,sedrape,numel_drape)
266 ENDIF
267C-----------------------------------------------
268C IF [KT]
269C-----------------------------------------------
270 iun = 1
271 CALL cmatip3(jft ,jlt ,pm ,mat ,pid ,
272 1 mtn ,npt ,hm ,hf ,iorth ,
273 2 hmor ,hfor ,hmfor ,iun )
274C
275 IF (iorth >0 .AND.ikproj<=0 .AND.idril==0 ) THEN
276 DO i=1,9
277 DO ep=jft,jlt
278 mf11(i,ep) =zero
279 mf22(i,ep) =zero
280 mf33(i,ep) =zero
281 mf44(i,ep) =zero
282 mf12(i,ep) =zero
283 mf13(i,ep) =zero
284 mf14(i,ep) =zero
285 mf23(i,ep) =zero
286 mf24(i,ep) =zero
287 mf34(i,ep) =zero
288 fm12(i,ep) =zero
289 fm13(i,ep) =zero
290 fm14(i,ep) =zero
291 fm23(i,ep) =zero
292 fm24(i,ep) =zero
293 fm34(i,ep) =zero
294 ENDDO
295 ENDDO
296 ENDIF
297C-----------------------------------------------
298C MATRICE [B]---index changed from here JFT-NPLAT (plat els)+,JLT(warped)
299C-----------------------------------------------
300 CALL czbe3(jft ,jlt ,area ,a_i ,x13 ,
301 2 x24 ,y13 ,y24 ,mx13 ,mx23 ,
302 3 mx34 ,my13 ,my23 ,my34 ,z1 ,
303 4 px1 ,px2 ,py1 ,py2 ,rx ,
304 5 ry ,sx ,sy ,rx2 ,ry2 ,
305 6 sx2 ,sy2 ,rhx ,rhy ,shx ,
306 7 shy ,ph1 ,ph2 ,hxx ,hyy ,
307 8 hxy ,nplat,iplat)
308C----------------------------------
309C SOUS-MATRICE DE RIGIDITE [K]
310C----------------------------------
311C--------------------------
312C 1. PARTIE CONSTANTE
313C--------------------------
314 CALL czlkec3(jft ,jlt ,vol ,thk0 ,thk2 ,
315 2 hm ,hf ,hz ,a_i ,z1 ,
316 3 px1 ,px2 ,py1 ,py2 ,nplat,
317 4 iplat,dhz ,
318 4 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
319 5 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
320 6 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
321 7 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
322 8 iorth,hmor,hfor,hmfor)
323C--------------------------
324C 2. Cisaillement Transversale (const+hourglass):
325C--------------------------
326 CALL czlkect3(jft ,jlt ,vol ,hc ,rx ,
327 4 ry ,sx ,sy ,rx2 ,ry2 ,
328 5 sx2 ,sy2 ,rhx ,rhy ,shx ,
329 6 shy ,gs ,nplat ,iplat,
330 9 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
331 a m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
332 b mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
333 c mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34)
334 IF (idril>0) THEN
335 CALL czber3(jft ,jlt ,area ,a_i ,x13 ,
336 1 x24 ,y13 ,y24 ,mx13 ,mx23 ,
337 2 mx34 ,my13 ,my23 ,my34 ,z1 ,
338 3 rx ,ry ,sx ,sy ,prx ,
339 4 pry ,prxy ,prz ,phkrx,phkry,
340 5 phkrxy,pherx,phery,pherxy,
341 6 phkrz,pherz ,phkx ,phky ,phex ,
342 7 phey ,iplat)
343 CALL czlkecr3(jft ,jlt ,vol ,thk0 ,thk2 ,
344 2 hm ,hf ,hz ,a_i ,z1 ,
345 3 px1 ,px2 ,py1 ,py2 ,
346 6 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
347 7 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
348 8 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
349 9 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
350 a iorth,hmor,hfor ,iplat,dhz ,
351 4 prx ,pry ,prxy ,prz ,hmfor,nplat)
352 ENDIF
353C--------------------------
354C 3. PARTIE HOURGLASS
355C--------------------------
356C--------------modif here!!! with IKT>0---add ET in HM,HF
357 IF ( iorth >0 .OR. mtn == 27) THEN
358c IF (MTN==19.OR.MTN==15.OR.MTN==25) THEN
359C-----------keep elastic-isotropic for hourglass part ------------
360 CALL cmatch3(jft ,jlt ,pm ,mat ,geo ,
361 1 pid ,mtn ,idril ,igeo ,hm ,
362 2 hf ,hz )
363 ENDIF
364 CALL czlken3(jft ,jlt ,vol ,thk0 ,thk2 ,
365 2 hm ,hz ,a_i ,px1 ,px2 ,
366 3 py1 ,py2 ,hxx ,hyy ,hxy ,
367 4 ph1 ,ph2 ,z1 ,nplat,iplat,dhz ,
368 5 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
369 6 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
370 7 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
371 8 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
372 9 idril )
373 IF (idril>0) THEN
374 CALL czlkenr3(jft ,jlt ,vol ,thk0 ,thk2 ,
375 2 hm ,hz ,a_i ,px1 ,px2 ,
376 3 py1 ,py2 ,hxx ,hyy ,hxy ,
377 4 ph1 ,ph2 ,z1 ,nplat,iplat,dhz ,
378 5 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
379 6 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
380 7 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
381 8 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
382 9 phkrx,phkry,phkrxy,pherx,phery,pherxy,
383 a phkrz,pherz,phkx ,phky ,phex ,phey )
384 ENDIF
385 IF (ikgeo ==1)
386 . CALL czlkecg3(jft ,jlt ,vol ,thk0 ,thk2 ,
387 1 px1 ,px2 ,py1 ,py2 ,rx ,
388 2 ry ,sx ,sy ,rx2 ,ry2 ,
389 3 sx2 ,sy2 ,rhx ,rhy ,shx ,
390 4 shy ,nplat ,iplat,gbuf%FOR,gbuf%MOM,
391 5 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
392 6 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
393 7 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
394 8 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
395 9 idril,iorth ,nel)
396C--------------------------
397C ASSEMBLE+LOCAL->GLOBAL
398C--------------------------
399 CALL czsumg3(
400 1 jft ,jlt ,vqn ,vq ,nplat,
401 2 iplat ,
402 3 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
403 4 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
404 5 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
405 6 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
406 7 ke11,ke22,ke33,ke44,ke12,ke13,ke14,ke23,
407 8 ke24,ke34,corelv,z1 ,idril ,iorth)
408C
409 IF (neig>0) CALL c4eoff(
410 1 jft, jlt, ixc, etag, off)
411C
412 CALL assem_c4(
413 1 ixc ,nel ,iddl ,ndof ,k_diag ,
414 2 k_lt ,iadk ,jdik ,ke11 ,ke12 ,
415 3 ke13 ,ke14 ,ke22 ,ke23 ,ke24 ,
416 5 ke33 ,ke34 ,ke44 ,off )
417C
418 RETURN
419 END
420!||====================================================================
421!|| czfkeu3 ../engine/source/elements/shell/coquez/czke3.F
422!||--- calls -----------------------------------------------------
423!|| czfikijuj ../engine/source/elements/shell/coquez/czke3.F
424!||====================================================================
425 SUBROUTINE czfkeu3 (
426 1 JFT ,JLT ,
427 8 KE11 ,KE22 ,KE33 ,KE44 ,
428 8 KE12 ,KE13 ,KE14 ,KE23 ,
429 8 KE24 ,KE34 ,UI ,RI ,
430 9 FI ,MI )
431C-----------------------------------------------
432C I M P L I C I T T Y P E S
433C-----------------------------------------------
434#include "implicit_f.inc"
435C-----------------------------------------------
436C G L O B A L P A R A M E T E R S
437C-----------------------------------------------
438#include "mvsiz_p.inc"
439C-----------------------------------------------
440C D U M M Y A R G U M E N T S
441C-----------------------------------------------
442 INTEGER JFT ,JLT
443 my_real
444 . UI(3,4,MVSIZ),RI(3,4,MVSIZ),FI(3,4,MVSIZ),MI(3,4,MVSIZ)
445 my_real
446 . KE11(36,MVSIZ),KE22(36,MVSIZ),KE33(36,MVSIZ),KE44(36,MVSIZ),
447 . KE12(36,MVSIZ),KE13(36,MVSIZ),KE14(36,MVSIZ),KE23(36,MVSIZ),
448 . KE24(36,MVSIZ),KE34(36,MVSIZ)
449C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
450C L O C A L V A R I A B L E S
451C--------------------------------
452 INTEGER
453 . I, J, M ,EP,ITRAN0,ITRAN1
454 my_real
455 . fj(3),mj(3),uj(3),rj(3)
456C
457 itran0=0
458 itran1=1
459 DO i=1,3
460 DO j=1,4
461 DO ep=jft ,jlt
462 fi(i,j,ep)=zero
463 mi(i,j,ep)=zero
464 END DO
465 END DO
466 END DO
467C
468 DO ep=jft ,jlt
469 j=1
470 DO i=1,3
471 uj(i)=ui(i,j,ep)
472 rj(i)=ri(i,j,ep)
473 END DO
474 print *,'KE(1,4),KE(1,5),KE(1,6)='
475 print *,ke11(19,1),ke11(25,1),ke11(31,1)
476 print *,ke11(19,1)*rj(1),ke11(25,1)*rj(2),ke11(31,1)*rj(3)
477 CALL czfikijuj (ke11(1,ep),uj ,rj ,fj ,mj ,
478 . itran0 )
479 DO i=1,3
480 fi(i,1,ep)=fi(i,1,ep)+fj(i)
481 mi(i,1,ep)=mi(i,1,ep)+mj(i)
482 END DO
483 print *,'KE12(1,4),KE(1,5),KE(1,6)='
484 print *,ke12(19,1),ke12(25,1),ke12(31,1)
485 print *,ke12(19,1)*rj(1),ke12(25,1)*rj(2),ke12(31,1)*rj(3)
486 CALL czfikijuj (ke12(1,ep),uj ,rj ,fj ,mj ,
487 . itran1 )
488 DO i=1,3
489 fi(i,2,ep)=fi(i,2,ep)+fj(i)
490 mi(i,2,ep)=mi(i,2,ep)+mj(i)
491 END DO
492 print *,'KE(1,4),KE(1,5),KE(1,6)='
493 print *,ke13(19,1),ke13(25,1),ke13(31,1)
494 print *,ke13(19,1)*rj(1),ke13(25,1)*rj(2),ke13(31,1)*rj(3)
495 CALL czfikijuj (ke13(1,ep),uj ,rj ,fj ,mj ,
496 . itran1 )
497 DO i=1,3
498 fi(i,3,ep)=fi(i,3,ep)+fj(i)
499 mi(i,3,ep)=mi(i,3,ep)+mj(i)
500 END DO
501 print *,'KE14(1,4),KE(1,5),KE(1,6)='
502 print *,ke14(19,1),ke14(25,1),ke14(31,1)
503 print *,ke14(19,1)*rj(1),ke14(25,1)*rj(2),ke14(31,1)*rj(3)
504 CALL czfikijuj (ke14(1,ep),uj ,rj ,fj ,mj ,
505 . itran1 )
506 DO i=1,3
507 fi(i,4,ep)=fi(i,4,ep)+fj(i)
508 mi(i,4,ep)=mi(i,4,ep)+mj(i)
509 END DO
510 j=2
511 DO i=1,3
512 uj(i)=ui(i,j,ep)
513 rj(i)=ri(i,j,ep)
514 END DO
515 CALL czfikijuj (ke12(1,ep),uj ,rj ,fj ,mj ,
516 . itran0 )
517 DO i=1,3
518 fi(i,1,ep)=fi(i,1,ep)+fj(i)
519 mi(i,1,ep)=mi(i,1,ep)+mj(i)
520 END DO
521 CALL czfikijuj (ke22(1,ep),uj ,rj ,fj ,mj ,
522 . itran0 )
523 DO i=1,3
524 fi(i,2,ep)=fi(i,2,ep)+fj(i)
525 mi(i,2,ep)=mi(i,2,ep)+mj(i)
526 END DO
527 CALL czfikijuj (ke23(1,ep),uj ,rj ,fj ,mj ,
528 . itran1 )
529 DO i=1,3
530 fi(i,3,ep)=fi(i,3,ep)+fj(i)
531 mi(i,3,ep)=mi(i,3,ep)+mj(i)
532 END DO
533 CALL czfikijuj (ke24(1,ep),uj ,rj ,fj ,mj ,
534 . itran1 )
535 DO i=1,3
536 fi(i,4,ep)=fi(i,4,ep)+fj(i)
537 mi(i,4,ep)=mi(i,4,ep)+mj(i)
538 END DO
539 j=3
540 DO i=1,3
541 uj(i)=ui(i,j,ep)
542 rj(i)=ri(i,j,ep)
543 END DO
544 CALL czfikijuj (ke13(1,ep),uj ,rj ,fj ,mj ,
545 . itran0 )
546 DO i=1,3
547 fi(i,1,ep)=fi(i,1,ep)+fj(i)
548 mi(i,1,ep)=mi(i,1,ep)+mj(i)
549 END DO
550 CALL czfikijuj (ke23(1,ep),uj ,rj ,fj ,mj ,
551 . itran0 )
552 DO i=1,3
553 fi(i,2,ep)=fi(i,2,ep)+fj(i)
554 mi(i,2,ep)=mi(i,2,ep)+mj(i)
555 END DO
556 CALL czfikijuj (ke33(1,ep),uj ,rj ,fj ,mj ,
557 . itran0 )
558 DO i=1,3
559 fi(i,3,ep)=fi(i,3,ep)+fj(i)
560 mi(i,3,ep)=mi(i,3,ep)+mj(i)
561 END DO
562 CALL czfikijuj (ke34(1,ep),uj ,rj ,fj ,mj ,
563 . itran1 )
564 DO i=1,3
565 fi(i,4,ep)=fi(i,4,ep)+fj(i)
566 mi(i,4,ep)=mi(i,4,ep)+mj(i)
567 END DO
568 j=4
569 DO i=1,3
570 uj(i)=ui(i,j,ep)
571 rj(i)=ri(i,j,ep)
572 END DO
573 CALL czfikijuj (ke14(1,ep),uj ,rj ,fj ,mj ,
574 . itran0 )
575 DO i=1,3
576 fi(i,1,ep)=fi(i,1,ep)+fj(i)
577 mi(i,1,ep)=mi(i,1,ep)+mj(i)
578 END DO
579 CALL czfikijuj (ke24(1,ep),uj ,rj ,fj ,mj ,
580 . itran0 )
581 DO i=1,3
582 fi(i,2,ep)=fi(i,2,ep)+fj(i)
583 mi(i,2,ep)=mi(i,2,ep)+mj(i)
584 END DO
585 CALL czfikijuj (ke34(1,ep),uj ,rj ,fj ,mj ,
586 . itran0 )
587 DO i=1,3
588 fi(i,3,ep)=fi(i,3,ep)+fj(i)
589 mi(i,3,ep)=mi(i,3,ep)+mj(i)
590 END DO
591 CALL czfikijuj (ke44(1,ep),uj ,rj ,fj ,mj ,
592 . itran0 )
593 DO i=1,3
594 fi(i,4,ep)=fi(i,4,ep)+fj(i)
595 mi(i,4,ep)=mi(i,4,ep)+mj(i)
596 END DO
597 END DO
598C
599 RETURN
600 END
601!||====================================================================
602!|| czfikijuj ../engine/source/elements/shell/coquez/czke3.f
603!||--- called by ------------------------------------------------------
604!|| czfkeu3 ../engine/source/elements/shell/coquez/czke3.F
605!||====================================================================
606 SUBROUTINE czfikijuj (KEIJ ,UJ ,RJ ,FI ,MI ,
607 . ITRAN )
608C-----------------------------------------------
609C I M P L I C I T T Y P E S
610C-----------------------------------------------
611#include "implicit_f.inc"
612C-----------------------------------------------
613C D U M M Y A R G U M E N T S
614C-----------------------------------------------
615 INTEGER ITRAN
616 my_real
617 . UJ(3),RJ(3),FI(3),MI(3)
618 my_real
619 . KEIJ(6,6)
620C-----------------------------------------------
621C L O C A L V A R I A B L E S
622C-----------------------------------------------
623 INTEGER
624 . I, J, M
625C
626 DO I=1,3
627 fi(i)=zero
628 mi(i)=zero
629 END DO
630C
631 IF (itran == 0) THEN
632 DO i=1,3
633 DO j=1,3
634 fi(i)=fi(i)+keij(i,j)*uj(j)+keij(i,j+3)*rj(j)
635 mi(i)=mi(i)+keij(i+3,j)*uj(j)+keij(i+3,j+3)*rj(j)
636 END DO
637 END DO
638 ELSE
639 DO i=1,3
640 DO j=1,3
641 fi(i)=fi(i)+keij(j,i)*uj(j)+keij(j+3,i)*rj(j)
642 mi(i)=mi(i)+keij(j,i+3)*uj(j)+keij(j+3,i+3)*rj(j)
643 END DO
644 END DO
645 END IF
646 RETURN
647 END
assem_c4
subroutine assem_c4(ixc, nel, iddl, ndof, k_diag, k_lt, iadk, jdik, kc11, kc12, kc13, kc14, kc22, kc23, kc24, kc33, kc34, kc44, off)
Definition assem_c4.F:37
c4eoff
subroutine c4eoff(jft, jlt, ixc, etag, off)
Definition c4eoff.F:30
cmatip3
subroutine cmatip3(jft, jlt, pm, mat, pid, mtn, npt, hm, hf, iorth, hmor, hfor, hmfor, ipg)
Definition cmatc3.F:582
cmatch3
subroutine cmatch3(jft, jlt, pm, mat, geo, pid, mtn, idril, igeo, hm, hf, hz)
Definition cmatc3.F:821
cmatc3
subroutine cmatc3(jft, jlt, pm, mat, geo, pid, area, thk0, thk02, thk, thke, volg, mtn, npt, ithk, hm, hf, hc, hz, igtyp, iorth, hmor, hfor, dir, igeo, idril, ihbe, hmfor, gs, isubstack, stack, elbuf_str, nlay, drape, nft, nel, indx_drape, sedrape, numel_drape)
Definition cmatc3.F:47
my_real
#define my_real
Definition cppsort.cpp:32
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
czcoork3
subroutine czcoork3(jft, jlt, x, ixc, pm, offg, area, area_i, vqn, vq, x13, x24, y13, y24, mx13, mx23, mx34, my13, my23, my34, z1, geo, elbuf_str, smstr, nlay, irep, npt, ismstr, dir_a, dir_b, pid, mat, ngl, nplat, iplat, corelv, off, thk, nel)
Definition czcoork3.F:42
czfkeu3
subroutine czfkeu3(jft, jlt, ke11, ke22, ke33, ke44, ke12, ke13, ke14, ke23, ke24, ke34, ui, ri, fi, mi)
Definition czke3.F:431
czfikijuj
subroutine czfikijuj(keij, uj, rj, fi, mi, itran)
Definition czke3.F:608
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
czlkec3
subroutine czlkec3(jft, jlt, vol, thk0, thk2, hm, hf, hz, a_i, z1, px1, px2, py1, py2, nplat, iplat, dhz, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34, iorth, hmor, hfor, hmfor)
Definition czlkec3.F:37
czlkecr3
subroutine czlkecr3(jft, jlt, vol, thk0, thk2, hm, hf, hz, a_i, z1, px1, px2, py1, py2, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34, iorth, hmor, hfor, iplat, dhz, prx, pry, prxy, prz, hmfor, nplat)
Definition czlkec3.F:645
czlkecg3
subroutine czlkecg3(jft, jlt, vol, thk0, thk2, px1, px2, py1, py2, rx, ry, sx, sy, rx2, ry2, sx2, sy2, rhx, rhy, shx, shy, nplat, iplat, for, mom, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34, idril, iorth, nel)
Definition czlkecg3.F:38
czlkect3
subroutine czlkect3(jft, jlt, vol, hc, rx, ry, sx, sy, rx2, ry2, sx2, sy2, rhx, rhy, shx, shy, gs, nplat, iplat, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34)
Definition czlkect3.F:36
czlken3
subroutine czlken3(jft, jlt, vol, thk0, thk2, hm, hz, a_i, px1, px2, py1, py2, hxx, hyy, hxy, ph1, ph2, z1, nplat, iplat, dhz, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34, idril)
Definition czlken3.F:37
czlkenr3
subroutine czlkenr3(jft, jlt, vol, thk0, thk2, hm, hz, a_i, px1, px2, py1, py2, hxx, hyy, hxy, ph1, ph2, z1, nplat, iplat, dhz, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34, phkrx, phkry, phkrxy, pherx, phery, pherxy, phkrz, pherz, phkx, phky, phex, phey)
Definition czlken3.F:287
czsumg3
subroutine czsumg3(jft, jlt, vqn, vq, nplat, iplat, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34, ke11, ke22, ke33, ke44, ke12, ke13, ke14, ke23, ke24, ke34, corelv, z1, idril, iorth)
Definition czsumg3.F:48
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
drape_mod
Definition drape_mod.F:54
stack_mod
Definition stack_mod.F:72