OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
czke3.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com04_c.inc"
#include "param_c.inc"
#include "impl1_c.inc"

Go to the source code of this file.

Functions/Subroutines

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)
subroutine czfkeu3 (jft, jlt, ke11, ke22, ke33, ke44, ke12, ke13, ke14, ke23, ke24, ke34, ui, ri, fi, mi)
subroutine czfikijuj (keij, uj, rj, fi, mi, itran)

Function/Subroutine Documentation

◆ czfikijuj()

subroutine czfikijuj ( keij,
uj,
rj,
fi,
mi,
integer itran )

Definition at line 608 of file czke3.F.

610C-----------------------------------------------
611C I M P L I C I T T Y P E S
612C-----------------------------------------------
613#include "implicit_f.inc"
614C-----------------------------------------------
615C D U M M Y A R G U M E N T S
616C-----------------------------------------------
617 INTEGER ITRAN
618 my_real
619 . uj(3),rj(3),fi(3),mi(3)
620 my_real
621 . keij(6,6)
622C-----------------------------------------------
623C L O C A L V A R I A B L E S
624C-----------------------------------------------
625 INTEGER
626 . I, J
627C
628 DO i=1,3
629 fi(i)=zero
630 mi(i)=zero
631 END DO
632C
633 IF (itran == 0) THEN
634 DO i=1,3
635 DO j=1,3
636 fi(i)=fi(i)+keij(i,j)*uj(j)+keij(i,j+3)*rj(j)
637 mi(i)=mi(i)+keij(i+3,j)*uj(j)+keij(i+3,j+3)*rj(j)
638 END DO
639 END DO
640 ELSE
641 DO i=1,3
642 DO j=1,3
643 fi(i)=fi(i)+keij(j,i)*uj(j)+keij(j+3,i)*rj(j)
644 mi(i)=mi(i)+keij(j,i+3)*uj(j)+keij(j+3,i+3)*rj(j)
645 END DO
646 END DO
647 END IF
648 RETURN
my_real
#define my_real
Definition cppsort.cpp:32

◆ czfkeu3()

subroutine czfkeu3 ( integer jft,
integer jlt,
ke11,
ke22,
ke33,
ke44,
ke12,
ke13,
ke14,
ke23,
ke24,
ke34,
ui,
ri,
fi,
mi )

Definition at line 427 of file czke3.F.

433C-----------------------------------------------
434C I M P L I C I T T Y P E S
435C-----------------------------------------------
436#include "implicit_f.inc"
437C-----------------------------------------------
438C G L O B A L P A R A M E T E R S
439C-----------------------------------------------
440#include "mvsiz_p.inc"
441C-----------------------------------------------
442C D U M M Y A R G U M E N T S
443C-----------------------------------------------
444 INTEGER JFT ,JLT
445 my_real
446 . ui(3,4,mvsiz),ri(3,4,mvsiz),fi(3,4,mvsiz),mi(3,4,mvsiz)
447 my_real
448 . ke11(36,mvsiz),ke22(36,mvsiz),ke33(36,mvsiz),ke44(36,mvsiz),
449 . ke12(36,mvsiz),ke13(36,mvsiz),ke14(36,mvsiz),ke23(36,mvsiz),
450 . ke24(36,mvsiz),ke34(36,mvsiz)
451C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
452C L O C A L V A R I A B L E S
453C--------------------------------
454 INTEGER
455 . I, J ,EP,ITRAN0,ITRAN1
456 my_real
457 . fj(3),mj(3),uj(3),rj(3)
458C
459 itran0=0
460 itran1=1
461 DO i=1,3
462 DO j=1,4
463 DO ep=jft ,jlt
464 fi(i,j,ep)=zero
465 mi(i,j,ep)=zero
466 END DO
467 END DO
468 END DO
469C
470 DO ep=jft ,jlt
471 j=1
472 DO i=1,3
473 uj(i)=ui(i,j,ep)
474 rj(i)=ri(i,j,ep)
475 END DO
476 print *,'KE(1,4),KE(1,5),KE(1,6)='
477 print *,ke11(19,1),ke11(25,1),ke11(31,1)
478 print *,ke11(19,1)*rj(1),ke11(25,1)*rj(2),ke11(31,1)*rj(3)
479 CALL czfikijuj (ke11(1,ep),uj ,rj ,fj ,mj ,
480 . itran0 )
481 DO i=1,3
482 fi(i,1,ep)=fi(i,1,ep)+fj(i)
483 mi(i,1,ep)=mi(i,1,ep)+mj(i)
484 END DO
485 print *,'KE12(1,4),KE(1,5),KE(1,6)='
486 print *,ke12(19,1),ke12(25,1),ke12(31,1)
487 print *,ke12(19,1)*rj(1),ke12(25,1)*rj(2),ke12(31,1)*rj(3)
488 CALL czfikijuj (ke12(1,ep),uj ,rj ,fj ,mj ,
489 . itran1 )
490 DO i=1,3
491 fi(i,2,ep)=fi(i,2,ep)+fj(i)
492 mi(i,2,ep)=mi(i,2,ep)+mj(i)
493 END DO
494 print *,'KE(1,4),KE(1,5),KE(1,6)='
495 print *,ke13(19,1),ke13(25,1),ke13(31,1)
496 print *,ke13(19,1)*rj(1),ke13(25,1)*rj(2),ke13(31,1)*rj(3)
497 CALL czfikijuj (ke13(1,ep),uj ,rj ,fj ,mj ,
498 . itran1 )
499 DO i=1,3
500 fi(i,3,ep)=fi(i,3,ep)+fj(i)
501 mi(i,3,ep)=mi(i,3,ep)+mj(i)
502 END DO
503 print *,'KE14(1,4),KE(1,5),KE(1,6)='
504 print *,ke14(19,1),ke14(25,1),ke14(31,1)
505 print *,ke14(19,1)*rj(1),ke14(25,1)*rj(2),ke14(31,1)*rj(3)
506 CALL czfikijuj (ke14(1,ep),uj ,rj ,fj ,mj ,
507 . itran1 )
508 DO i=1,3
509 fi(i,4,ep)=fi(i,4,ep)+fj(i)
510 mi(i,4,ep)=mi(i,4,ep)+mj(i)
511 END DO
512 j=2
513 DO i=1,3
514 uj(i)=ui(i,j,ep)
515 rj(i)=ri(i,j,ep)
516 END DO
517 CALL czfikijuj (ke12(1,ep),uj ,rj ,fj ,mj ,
518 . itran0 )
519 DO i=1,3
520 fi(i,1,ep)=fi(i,1,ep)+fj(i)
521 mi(i,1,ep)=mi(i,1,ep)+mj(i)
522 END DO
523 CALL czfikijuj (ke22(1,ep),uj ,rj ,fj ,mj ,
524 . itran0 )
525 DO i=1,3
526 fi(i,2,ep)=fi(i,2,ep)+fj(i)
527 mi(i,2,ep)=mi(i,2,ep)+mj(i)
528 END DO
529 CALL czfikijuj (ke23(1,ep),uj ,rj ,fj ,mj ,
530 . itran1 )
531 DO i=1,3
532 fi(i,3,ep)=fi(i,3,ep)+fj(i)
533 mi(i,3,ep)=mi(i,3,ep)+mj(i)
534 END DO
535 CALL czfikijuj (ke24(1,ep),uj ,rj ,fj ,mj ,
536 . itran1 )
537 DO i=1,3
538 fi(i,4,ep)=fi(i,4,ep)+fj(i)
539 mi(i,4,ep)=mi(i,4,ep)+mj(i)
540 END DO
541 j=3
542 DO i=1,3
543 uj(i)=ui(i,j,ep)
544 rj(i)=ri(i,j,ep)
545 END DO
546 CALL czfikijuj (ke13(1,ep),uj ,rj ,fj ,mj ,
547 . itran0 )
548 DO i=1,3
549 fi(i,1,ep)=fi(i,1,ep)+fj(i)
550 mi(i,1,ep)=mi(i,1,ep)+mj(i)
551 END DO
552 CALL czfikijuj (ke23(1,ep),uj ,rj ,fj ,mj ,
553 . itran0 )
554 DO i=1,3
555 fi(i,2,ep)=fi(i,2,ep)+fj(i)
556 mi(i,2,ep)=mi(i,2,ep)+mj(i)
557 END DO
558 CALL czfikijuj (ke33(1,ep),uj ,rj ,fj ,mj ,
559 . itran0 )
560 DO i=1,3
561 fi(i,3,ep)=fi(i,3,ep)+fj(i)
562 mi(i,3,ep)=mi(i,3,ep)+mj(i)
563 END DO
564 CALL czfikijuj (ke34(1,ep),uj ,rj ,fj ,mj ,
565 . itran1 )
566 DO i=1,3
567 fi(i,4,ep)=fi(i,4,ep)+fj(i)
568 mi(i,4,ep)=mi(i,4,ep)+mj(i)
569 END DO
570 j=4
571 DO i=1,3
572 uj(i)=ui(i,j,ep)
573 rj(i)=ri(i,j,ep)
574 END DO
575 CALL czfikijuj (ke14(1,ep),uj ,rj ,fj ,mj ,
576 . itran0 )
577 DO i=1,3
578 fi(i,1,ep)=fi(i,1,ep)+fj(i)
579 mi(i,1,ep)=mi(i,1,ep)+mj(i)
580 END DO
581 CALL czfikijuj (ke24(1,ep),uj ,rj ,fj ,mj ,
582 . itran0 )
583 DO i=1,3
584 fi(i,2,ep)=fi(i,2,ep)+fj(i)
585 mi(i,2,ep)=mi(i,2,ep)+mj(i)
586 END DO
587 CALL czfikijuj (ke34(1,ep),uj ,rj ,fj ,mj ,
588 . itran0 )
589 DO i=1,3
590 fi(i,3,ep)=fi(i,3,ep)+fj(i)
591 mi(i,3,ep)=mi(i,3,ep)+mj(i)
592 END DO
593 CALL czfikijuj (ke44(1,ep),uj ,rj ,fj ,mj ,
594 . itran0 )
595 DO i=1,3
596 fi(i,4,ep)=fi(i,4,ep)+fj(i)
597 mi(i,4,ep)=mi(i,4,ep)+mj(i)
598 END DO
599 END DO
600C
601 RETURN
czfikijuj
subroutine czfikijuj(keij, uj, rj, fi, mi, itran)
Definition czke3.F:610

◆ czke3()

subroutine czke3 ( integer jft,
integer jlt,
integer nft,
integer npt,
integer mtn,
integer ithk,
integer ncycle,
integer istrain,
integer ipla,
pm,
geo,
integer, dimension(nixc,*) ixc,
type (elbuf_struct_), target elbuf_str,
bufmat,
integer offset,
integer, dimension(mvsiz) indxof,
integer, dimension(*) etag,
integer, dimension(*) iddl,
integer, dimension(*) ndof,
k_diag,
k_lt,
integer, dimension(*) iadk,
integer, dimension(*) jdik,
integer ihbe,
thke,
integer ismstr,
x,
integer ikgeo,
integer, dimension(*) ipm,
integer, dimension(npropgi,*) igeo,
integer iexpan,
integer, dimension(*) iparg,
integer isubstack,
type (stack_ply) stack,
type (drape_), dimension(numelc_drape) drape_sh4n,
integer, dimension(sedrape) indx_drape,
integer, intent(in) sedrape,
integer, intent(in) numel_drape )

Definition at line 50 of file czke3.F.

59C-----------------------------------------------
60C M o d u l e s
61C-----------------------------------------------
62 USE elbufdef_mod
63 USE stack_mod
64 USE drape_mod
65 use element_mod , only : nixc
66C-----------------------------------------------
67C I M P L I C I T T Y P E S
68C-----------------------------------------------
69#include "implicit_f.inc"
70C-----------------------------------------------
71C G L O B A L P A R A M E T E R S
72C-----------------------------------------------
73#include "mvsiz_p.inc"
74C-----------------------------------------------
75C C O M M O N B L O C K S
76C-----------------------------------------------
77#include "com04_c.inc"
78#include "param_c.inc"
79#include "impl1_c.inc"
80C-----------------------------------------------
81C D U M M Y A R G U M E N T S
82C-----------------------------------------------
83 INTEGER JFT ,JLT ,NFT ,NPT ,
84 . MTN ,ITHK ,NCYCLE,ISUBSTACK,
85 . ISTRAIN ,IPLA ,OFFSET,IHBE ,ISMSTR,IKGEO,IEXPAN
86 INTEGER , INTENT(IN) :: SEDRAPE,NUMEL_DRAPE
87 INTEGER IXC(NIXC,*),IGEO(NPROPGI,*),IPM(*),IPARG(*)
88 INTEGER INDXOF(MVSIZ),
89 . ETAG(*),IDDL(*) ,NDOF(*) ,IADK(*) ,JDIK(*)
90 INTEGER, DIMENSION(SEDRAPE) :: INDX_DRAPE
91
92C real or real*8
93 my_real
94 . pm(npropm,*),geo(npropg,*),bufmat(*),x(3,*),thke(*),
95 . off(mvsiz),k_diag(*) ,k_lt(*)
96 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
97 TYPE (STACK_PLY) :: STACK
98 TYPE (DRAPE_) :: DRAPE_SH4N(NUMELC_DRAPE)
99C=======================================================================
100c FUNCTION: [K] stiffness Matrix of QEPH element
101c
102c Note:
103c ARGUMENTS: (I: input, O: output, IO: input * output, W: workspace)
104c
105c TYPE NAME FUNCTION
106c I JFT,JLT,NFT - local element id limit and first id(global)of this element group
107c NEL=JLT-JFT+1
108c I NPT,MTN - num. of integrating point in thickness and material type id
109c I ITHK - flag of thickness updating (if >0)
110c I NCYCLE - cycle(increment) number
111c I ISTRAIN - total strain output flag
112c I IPLA - radial return plasticity compute option
113c I PM ,GEO - Material and geometrical property data
114c I IXC(NIXC,NEL) - connectivity and mid,pid integer data
115c I BUFMAT() - internal material data
116c I INDXOF(NEL) - deleted element flag (not used in this subroutine)
117c I ETAG(NEL) - activating element flag for Eigenvalue analysis
118c I IDDL(NUMNOD) - DOF id of node N =IDDL(N)+1,NDOF
119c I NDOF(NUMNOD) - Num of DOF (nodal)
120c IO K_DIAG(NDDL) - Diagnale terms of assembled [K] NDDL: total number of model DOF
121c IO K_LT(NNZK) - terms of strick triagular of assembled [K] NNZK: number of no-zero terms
122c I IADK(NDDL),JDIK(NNZK) - Indice integer tables of Compress format of [K]
123c I IHBE - Shell formulation flag (Ishell)
124c I THKE - initial thickness
125c I ISMSTR - small strain flag
126c I IKGEO , - geometrical stiffness matrix including (if >0)
127c I X(3,NUMNOD) co-ordinate
128c I IPM ,IGEO - Material and geometrical property integer data
129c I IEXPAN - thermo flag used in materials
130c I IPARG(NG) element group data
131C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
132C L O C A L V A R I A B L E S
133C--------------------------------
134 INTEGER
135 . I, J,J1,J2, NEL, NPLAT,IPLAT(MVSIZ), NLAY,L_DIRA,L_DIRB,
136 . IREP,EP,IDRIL
137 INTEGER MAT(MVSIZ), PID(MVSIZ), NGL(MVSIZ),IORTH,IGTYP,IUN
138 my_real
139 . x13(mvsiz), x24(mvsiz), y13(mvsiz), y24(mvsiz),
140 . mx13(mvsiz), mx23(mvsiz), mx34(mvsiz),
141 . my13(mvsiz), my23(mvsiz), my34(mvsiz), z1(mvsiz),
142 . px1(mvsiz), px2(mvsiz), py1(mvsiz),py2(mvsiz),
143 . sx(mvsiz), sy(mvsiz), rx(mvsiz),ry(mvsiz),
144 . sx2(mvsiz), sy2(mvsiz), rx2(mvsiz),ry2(mvsiz),
145 . rhx(mvsiz,4),rhy(mvsiz,4),shx(mvsiz,4),shy(mvsiz,4),
146 . ph1(mvsiz),ph2(mvsiz),hxx(mvsiz),hyy(mvsiz),hxy(mvsiz)
147 my_real
148 . vq(mvsiz,9),area(mvsiz), vqn(mvsiz,12),thk0(mvsiz),vol(mvsiz),
149 . a_i(mvsiz), thk2(mvsiz),hm(mvsiz,4),hf(mvsiz,4),hc(mvsiz,2),
150 . hz(mvsiz),dhz(mvsiz),hmor(mvsiz,2),hfor(mvsiz,2),
151 . gs(mvsiz),hmfor(mvsiz,6)
152 my_real
153 . corelv(mvsiz,2,4)
154 my_real
155 . k11(9,mvsiz),k12(9,mvsiz),k13(9,mvsiz),k14(9,mvsiz),
156 . k22(9,mvsiz),k23(9,mvsiz),k24(9,mvsiz),k33(9,mvsiz),
157 . m11(9,mvsiz),m12(9,mvsiz),m13(9,mvsiz),m14(9,mvsiz),
158 . m22(9,mvsiz),m23(9,mvsiz),m24(9,mvsiz),m33(9,mvsiz),
159 . mf11(9,mvsiz),mf12(9,mvsiz),mf13(9,mvsiz),mf14(9,mvsiz),
160 . mf22(9,mvsiz),mf23(9,mvsiz),mf24(9,mvsiz),mf33(9,mvsiz),
161 . fm12(9,mvsiz),fm13(9,mvsiz),fm14(9,mvsiz),
162 . fm23(9,mvsiz),fm24(9,mvsiz),fm34(9,mvsiz),
163 . k34(9,mvsiz),k44(9,mvsiz),m34(9,mvsiz),m44(9,mvsiz),
164 . mf34(9,mvsiz),mf44(9,mvsiz)
165 my_real
166 . prx(4,mvsiz),pry(4,mvsiz),prxy(4,mvsiz),phkrx(4,mvsiz),
167 . phkry(4,mvsiz),phkrxy(4,mvsiz),pherx(4,mvsiz),phery(4,mvsiz),
168 . pherxy(4,mvsiz),prz(4,mvsiz),phkrz(4,mvsiz),pherz(4,mvsiz),
169 . phkx(mvsiz),phky(mvsiz),phex(mvsiz),phey(mvsiz)
170 my_real
171 . ke11(36,mvsiz),ke22(36,mvsiz),ke33(36,mvsiz),ke44(36,mvsiz),
172 . ke12(36,mvsiz),ke13(36,mvsiz),ke14(36,mvsiz),ke23(36,mvsiz),
173 . ke24(36,mvsiz),ke34(36,mvsiz)
174C-----------------------------------------------
175 my_real,
176 . DIMENSION(:) ,POINTER :: dir_a, dir_b
177 my_real,
178 . ALLOCATABLE, DIMENSION(:), TARGET :: dira,dirb
179 TYPE(G_BUFEL_) ,POINTER :: GBUF
180C------------|---------|------------------------------------------------------
181C--------------------------
182C INITIALISATION
183C--------------------------
184C OPEN(UNIT=17,FILE='DEBZ.TMP',STATUS='UNKNOWN',FORM='FORMATTED')
185c
186 gbuf => elbuf_str%GBUF
187 nel=jlt-jft+1
188 idril = iparg(41)
189C
190 igtyp = igeo(11,ixc(6,1))
191 irep = igeo(6 ,ixc(6,1))
192 nlay = elbuf_str%NLAY
193 l_dira = elbuf_str%BUFLY(1)%LY_DIRA
194 l_dirb = elbuf_str%BUFLY(1)%LY_DIRB
195 ALLOCATE(dira(nlay*nel*l_dira))
196 ALLOCATE(dirb(nlay*nel*l_dirb))
197 dira = zero
198 dirb = zero
199 dir_a => dira(1:nlay*nel*l_dira)
200 dir_b => dirb(1:nlay*nel*l_dirb)
201 IF (irep == 0) THEN
202 DO j=1,nlay
203 j1 = 1+(j-1)*l_dira*nel
204 j2 = j*l_dira*nel
205 dira(j1:j2) = elbuf_str%BUFLY(j)%DIRA(1:nel*l_dira)
206 ENDDO
207 ENDIF
208C
209 CALL czcoork3(jft ,jlt ,x ,ixc ,pm ,
210 1 gbuf%OFF,area,a_i,vqn ,vq ,
211 2 x13 ,x24 ,y13 ,y24 ,mx13,
212 3 mx23,mx34 ,my13 ,my23 ,my34,
213 4 z1 , geo ,
214 5 elbuf_str,gbuf%SMSTR,nlay,
215 6 irep,npt,ismstr,
216 7 dir_a,dir_b,pid,mat,ngl,nplat,iplat ,
217 8 corelv,off,thke,nel)
218 IF (ikproj>0.OR.idril>0) THEN
219 DO i=1,9
220 DO ep=jft,jlt
221 m11(i,ep) =zero
222 m22(i,ep) =zero
223 m33(i,ep) =zero
224 m44(i,ep) =zero
225 m12(i,ep) =zero
226 m13(i,ep) =zero
227 m14(i,ep) =zero
228 m23(i,ep) =zero
229 m24(i,ep) =zero
230 m34(i,ep) =zero
231 mf11(i,ep) =zero
232 mf22(i,ep) =zero
233 mf33(i,ep) =zero
234 mf44(i,ep) =zero
235 mf12(i,ep) =zero
236 mf13(i,ep) =zero
237 mf14(i,ep) =zero
238 mf23(i,ep) =zero
239 mf24(i,ep) =zero
240 mf34(i,ep) =zero
241 fm12(i,ep) =zero
242 fm13(i,ep) =zero
243 fm14(i,ep) =zero
244 fm23(i,ep) =zero
245 fm24(i,ep) =zero
246 fm34(i,ep) =zero
247 ENDDO
248 ENDDO
249 ENDIF
250 IF (irep>0) THEN
251 CALL cmatc3(jft ,jlt ,pm ,mat ,geo ,
252 1 pid ,area ,thk0 ,thk2 ,gbuf%THK ,
253 2 thke ,vol ,mtn ,npt ,ithk ,
254 3 hm ,hf ,hc ,hz ,igtyp ,
255 4 iorth ,hmor ,hfor ,dir_a ,igeo ,
256 5 idril ,ihbe ,hmfor ,gs ,isubstack,
257 6 stack ,elbuf_str ,nlay ,drape_sh4n ,nft ,
258 7 nel ,indx_drape,sedrape,numel_drape)
259 ELSE
260 CALL cmatc3(jft ,jlt ,pm ,mat ,geo ,
261 1 pid ,area ,thk0 ,thk2 ,gbuf%THK ,
262 2 thke ,vol ,mtn ,npt ,ithk ,
263 3 hm ,hf ,hc ,hz ,igtyp ,
264 4 iorth ,hmor ,hfor ,dira ,igeo ,
265 5 idril ,ihbe ,hmfor ,gs ,isubstack,
266 6 stack ,elbuf_str ,nlay ,drape_sh4n ,nft ,
267 7 nel ,indx_drape,sedrape,numel_drape)
268 ENDIF
269C-----------------------------------------------
270C IF [KT]
271C-----------------------------------------------
272 iun = 1
273 CALL cmatip3(jft ,jlt ,pm ,mat ,pid ,
274 1 mtn ,npt ,hm ,hf ,iorth ,
275 2 hmor ,hfor ,hmfor ,iun )
276C
277 IF (iorth >0 .AND.ikproj<=0 .AND.idril==0 ) THEN
278 DO i=1,9
279 DO ep=jft,jlt
280 mf11(i,ep) =zero
281 mf22(i,ep) =zero
282 mf33(i,ep) =zero
283 mf44(i,ep) =zero
284 mf12(i,ep) =zero
285 mf13(i,ep) =zero
286 mf14(i,ep) =zero
287 mf23(i,ep) =zero
288 mf24(i,ep) =zero
289 mf34(i,ep) =zero
290 fm12(i,ep) =zero
291 fm13(i,ep) =zero
292 fm14(i,ep) =zero
293 fm23(i,ep) =zero
294 fm24(i,ep) =zero
295 fm34(i,ep) =zero
296 ENDDO
297 ENDDO
298 ENDIF
299C-----------------------------------------------
300C MATRICE [B]---index changed from here JFT-NPLAT (plat els)+,JLT(warped)
301C-----------------------------------------------
302 CALL czbe3(jft ,jlt ,area ,a_i ,x13 ,
303 2 x24 ,y13 ,y24 ,mx13 ,mx23 ,
304 3 mx34 ,my13 ,my23 ,my34 ,z1 ,
305 4 px1 ,px2 ,py1 ,py2 ,rx ,
306 5 ry ,sx ,sy ,rx2 ,ry2 ,
307 6 sx2 ,sy2 ,rhx ,rhy ,shx ,
308 7 shy ,ph1 ,ph2 ,hxx ,hyy ,
309 8 hxy ,nplat,iplat)
310C----------------------------------
311C SOUS-MATRICE DE RIGIDITE [K]
312C----------------------------------
313C--------------------------
314C 1. PARTIE CONSTANTE
315C--------------------------
316 CALL czlkec3(jft ,jlt ,vol ,thk0 ,thk2 ,
317 2 hm ,hf ,hz ,a_i ,z1 ,
318 3 px1 ,px2 ,py1 ,py2 ,nplat,
319 4 iplat,dhz ,
320 4 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
321 5 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
322 6 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
323 7 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
324 8 iorth,hmor,hfor,hmfor)
325C--------------------------
326C 2. Transversal shear (Const+Hourglass):
327C--------------------------
328 CALL czlkect3(jft ,jlt ,vol ,hc ,rx ,
329 4 ry ,sx ,sy ,rx2 ,ry2 ,
330 5 sx2 ,sy2 ,rhx ,rhy ,shx ,
331 6 shy ,gs ,nplat ,iplat,
332 9 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
333 a m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
334 b mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
335 c mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34)
336 IF (idril>0) THEN
337 CALL czber3(jft ,jlt ,area ,a_i ,x13 ,
338 1 x24 ,y13 ,y24 ,mx13 ,mx23 ,
339 2 mx34 ,my13 ,my23 ,my34 ,z1 ,
340 3 rx ,ry ,sx ,sy ,prx ,
341 4 pry ,prxy ,prz ,phkrx,phkry,
342 5 phkrxy,pherx,phery,pherxy,
343 6 phkrz,pherz ,phkx ,phky ,phex ,
344 7 phey ,iplat)
345 CALL czlkecr3(jft ,jlt ,vol ,thk0 ,thk2 ,
346 2 hm ,hf ,hz ,a_i ,z1 ,
347 3 px1 ,px2 ,py1 ,py2 ,
348 6 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
349 7 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
350 8 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
351 9 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
352 a iorth,hmor,hfor ,iplat,dhz ,
353 4 prx ,pry ,prxy ,prz ,hmfor,nplat)
354 ENDIF
355C--------------------------
356C 3. PARTIE HOURGLASS
357C--------------------------
358C--------------modif here!!! with IKT>0---add ET in HM,HF
359 IF ( iorth >0 .OR. mtn == 27) THEN
360c IF (MTN==19.OR.MTN==15.OR.MTN==25) THEN
361C-----------keep elastic-isotropic for hourglass part ------------
362 CALL cmatch3(jft ,jlt ,pm ,mat ,geo ,
363 1 pid ,mtn ,idril ,igeo ,hm ,
364 2 hf ,hz )
365 ENDIF
366 CALL czlken3(jft ,jlt ,vol ,thk0 ,thk2 ,
367 2 hm ,hz ,a_i ,px1 ,px2 ,
368 3 py1 ,py2 ,hxx ,hyy ,hxy ,
369 4 ph1 ,ph2 ,z1 ,nplat,iplat,dhz ,
370 5 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
371 6 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
372 7 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
373 8 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
374 9 idril )
375 IF (idril>0) THEN
376 CALL czlkenr3(jft ,jlt ,vol ,thk0 ,thk2 ,
377 2 hm ,hz ,a_i ,px1 ,px2 ,
378 3 py1 ,py2 ,hxx ,hyy ,hxy ,
379 4 ph1 ,ph2 ,z1 ,nplat,iplat,dhz ,
380 5 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
381 6 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
382 7 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
383 8 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
384 9 phkrx,phkry,phkrxy,pherx,phery,pherxy,
385 a phkrz,pherz,phkx ,phky ,phex ,phey )
386 ENDIF
387 IF (ikgeo ==1)
388 . CALL czlkecg3(jft ,jlt ,vol ,thk0 ,thk2 ,
389 1 px1 ,px2 ,py1 ,py2 ,rx ,
390 2 ry ,sx ,sy ,rx2 ,ry2 ,
391 3 sx2 ,sy2 ,rhx ,rhy ,shx ,
392 4 shy ,nplat ,iplat,gbuf%FOR,gbuf%MOM,
393 5 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
394 6 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
395 7 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
396 8 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
397 9 idril,iorth ,nel)
398C--------------------------
399C ASSEMBLE+LOCAL->GLOBAL
400C--------------------------
401 CALL czsumg3(
402 1 jft ,jlt ,vqn ,vq ,nplat,
403 2 iplat ,
404 3 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
405 4 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
406 5 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
407 6 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
408 7 ke11,ke22,ke33,ke44,ke12,ke13,ke14,ke23,
409 8 ke24,ke34,corelv,z1 ,idril ,iorth)
410C
411 IF (neig>0) CALL c4eoff(
412 1 jft, jlt, ixc, etag, off)
413C
414 CALL assem_c4(
415 1 ixc ,nel ,iddl ,ndof ,k_diag ,
416 2 k_lt ,iadk ,jdik ,ke11 ,ke12 ,
417 3 ke13 ,ke14 ,ke22 ,ke23 ,ke24 ,
418 5 ke33 ,ke34 ,ke44 ,off )
419C
420 RETURN
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:39
c4eoff
subroutine c4eoff(jft, jlt, ixc, etag, off)
Definition c4eoff.F:32
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
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:43
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:75