OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cbake3.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!|| cbake3 ../engine/source/elements/shell/coqueba/cbake3.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!|| cbabe3 ../engine/source/elements/shell/coqueba/cbabe3.F
32!|| cbabec3 ../engine/source/elements/shell/coqueba/cbabe3.F
33!|| cbaber3 ../engine/source/elements/shell/coqueba/cbabe3.F
34!|| cbacoork ../engine/source/elements/shell/coqueba/cbacoork.F
35!|| cbaderirz ../engine/source/elements/shell/coqueba/cbadef.F
36!|| cbaini3 ../engine/source/elements/shell/coqueba/cbake3.f
37!|| cbalke3 ../engine/source/elements/shell/coqueba/cbalke3.F
38!|| cbalkec3 ../engine/source/elements/shell/coqueba/cbalke3.F
39!|| cbalkerz ../engine/source/elements/shell/coqueba/cbalke3.F
40!|| cbasumg3 ../engine/source/elements/shell/coqueba/cbasumg3.F
41!|| cmatc3 ../engine/source/elements/shell/coqueba/cmatc3.F
42!|| cmatip3 ../engine/source/elements/shell/coqueba/cmatc3.F
43!||--- uses -----------------------------------------------------
44!|| drape_mod ../engine/share/modules/drape_mod.F
45!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
46!|| stack_mod ../engine/share/modules/stack_mod.F
47!||====================================================================
48 SUBROUTINE cbake3 (
49 1 JFT ,JLT ,NFT ,NPT ,MTN ,
50 2 ITHK ,NCYCLE ,
51 3 ISTRAIN,IPLA ,PM ,GEO ,IXC ,
52 4 ELBUF_STR ,BUFMAT ,OFFSET ,INDXOF ,
53 5 ETAG , IDDL ,NDOF ,K_DIAG ,K_LT , IADK ,JDIK ,
54 6 IHBE ,THKE ,ISMSTR ,X ,IKGEO ,
55 7 IPM ,IGEO ,IEXPAN ,IPARG ,ISUBSTACK ,STACK ,
56 8 DRAPE_SH4N ,INDX_DRAPE,SEDRAPE,NUMEL_DRAPE)
57C-----------------------------------------------
58C M o d u l e s
59C-----------------------------------------------
60 USE elbufdef_mod
61 USE stack_mod
62 USE drape_mod
63C-----------------------------------------------
64C I M P L I C I T T Y P E S
65C-----------------------------------------------
66#include "implicit_f.inc"
67C-----------------------------------------------
68C G L O B A L P A R A M E T E R S
69C-----------------------------------------------
70#include "mvsiz_p.inc"
71C-----------------------------------------------
72C C O M M O N B L O C K S
73C-----------------------------------------------
74#include "com04_c.inc"
75#include "param_c.inc"
76C-----------------------------------------------
77C D U M M Y A R G U M E N T S
78C-----------------------------------------------
79 INTEGER JFT ,JLT ,NFT ,NPT ,MTN ,ITHK ,
80 . NCYCLE,ISTRAIN ,IPLA ,OFFSET,IHBE ,ISMSTR,IKGEO,IEXPAN
81 INTEGER , INTENT(IN) :: SEDRAPE,NUMEL_DRAPE
82 INTEGER IXC(NIXC,*) ,IGEO(NPROPGI,*),IPM(*),IPARG(*)
83 INTEGER INDXOF(MVSIZ),ISUBSTACK,
84 . ETAG(*),IDDL(*) ,NDOF(*) ,IADK(*) ,JDIK(*)
85 INTEGER, DIMENSION(SEDRAPE) :: INDX_DRAPE
86C REAL OU REAL*8
87 my_real
88 . pm(npropm,*),geo(npropg,*),bufmat(*), x(3,*),thke(*)
89 my_real
90 . ke11(36,mvsiz),ke22(36,mvsiz),ke33(36,mvsiz),ke44(36,mvsiz),
91 . ke12(36,mvsiz),ke13(36,mvsiz),ke14(36,mvsiz),ke23(36,mvsiz),
92 . ke24(36,mvsiz),ke34(36,mvsiz),off(mvsiz),k_diag(*) ,k_lt(*)
93 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
94 TYPE (STACK_PLY) :: STACK
95 TYPE (DRAPE_) :: DRAPE_SH4N(NUMELC_DRAPE)
96C-----------------------------------------------
97C L O C A L V A R I A B L E S
98C-----------------------------------------------
99 INTEGER
100 . NPLAT,NLAY,IPLAT(MVSIZ)
101 INTEGER
102 . i, j,j1,j2, ir, is, nel, iun, mx,l_dira ,l_dirb ,
103 . ep,ng,npg,nnod,irep,pt1,pt2,pt3,lenf,lenm,nptr,npts,
104 . pt0,ptf,ptm,pte,ptep,pts,pptf,pptm,ppte,pptep,ppts
105 INTEGER MAT(MVSIZ), PID(MVSIZ), NGL(MVSIZ)
106 INTEGER MPT,IORTH,IBID,IDRIL
107 parameter(npg = 4)
108 parameter(nnod = 4)
109 my_real
110 . vcore(mvsiz,3*nnod),
111 . vqn(mvsiz,9*nnod),vqg(mvsiz,9*nnod),vnrm(mvsiz,3*nnod),
112 . bm(mvsiz,9*nnod),bmf(mvsiz,9*nnod),bf(mvsiz,6*nnod),
113 . bc(mvsiz,10*nnod),vq(mvsiz,9),vjfi(mvsiz,6,4),
114 . tc(mvsiz,4),jac(mvsiz,npg),hx(mvsiz,npg),hy(mvsiz,npg),
115 . veta(4,npg),vksi(4,npg),bzz(mvsiz,2*nnod)
116 my_real
117 . vastn(mvsiz,4*nnod),area(mvsiz),
118 . cdet(mvsiz),thk2(mvsiz)
119 INTEGER
120 . NEL8,NEL5,NEL3,NPTM,IGTYP,PTMAT,NBM_S,NBDIR,NB16A
121 MY_REAL
122 . SIGY(MVSIZ),VOL0(MVSIZ),THK0(MVSIZ),
123 . X13(MVSIZ) ,Y13(MVSIZ), X24(MVSIZ) ,HZ(MVSIZ),
124 . volg(mvsiz),y24(mvsiz),hm(mvsiz,4),hf(mvsiz,4),hc(mvsiz,2),
125 . hmor(mvsiz,2),hfor(mvsiz,2),hmfor(mvsiz,6),gs(mvsiz)
126 my_real
127 . k11(9,mvsiz),k12(9,mvsiz),k13(9,mvsiz),k14(9,mvsiz),
128 . k22(9,mvsiz),k23(9,mvsiz),k24(9,mvsiz),k33(9,mvsiz),
129 . m11(9,mvsiz),m12(9,mvsiz),m13(9,mvsiz),m14(9,mvsiz),
130 . m22(9,mvsiz),m23(9,mvsiz),m24(9,mvsiz),m33(9,mvsiz),
131 . mf11(9,mvsiz),mf12(9,mvsiz),mf13(9,mvsiz),mf14(9,mvsiz),
132 . mf22(9,mvsiz),mf23(9,mvsiz),mf24(9,mvsiz),mf33(9,mvsiz),
133 . fm12(9,mvsiz),fm13(9,mvsiz),fm14(9,mvsiz),
134 . fm23(9,mvsiz),fm24(9,mvsiz),fm34(9,mvsiz),
135 . k34(9,mvsiz),k44(9,mvsiz),m34(9,mvsiz),m44(9,mvsiz),
136 . mf34(9,mvsiz),mf44(9,mvsiz),
137 . bm0rz(mvsiz,4,4),bmkrz(mvsiz,4,4),bmerz(mvsiz,4,4),
138 . bmrz(mvsiz,3,4),brz(mvsiz,4,4)
139C-----------------------------------------------
140 my_real,
141 . DIMENSION(:) ,POINTER :: dir_a, dir_b
142 my_real,
143 . ALLOCATABLE, DIMENSION(:), TARGET :: dira,dirb
144 TYPE(g_bufel_) ,POINTER :: GBUF
145C-----------------------------------------------
146C INITIALISATION
147C--------------------------
148 gbuf => elbuf_str%GBUF
149 iun = 1
150 nel=jlt-jft+iun
151 IF (mtn==1) npt=0
152 mpt=iabs(npt)
153 idril = iparg(41)
154C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
155 nel3 = nel*3
156 nel5 = nel*5
157 nel8 = nel*8
158 nptm = max(1,mpt)
159 nlay = elbuf_str%NLAY
160C
161 igtyp = igeo(11,ixc(6,1))
162 irep = igeo(6 ,ixc(6,1))
163 l_dira = elbuf_str%BUFLY(1)%LY_DIRA
164 l_dirb = elbuf_str%BUFLY(1)%LY_DIRB
165 ALLOCATE(dira(nlay*nel*l_dira))
166 ALLOCATE(dirb(nlay*nel*l_dirb))
167 dira = zero
168 dirb = zero
169 dir_a => dira(1:nlay*nel*l_dira)
170 dir_b => dirb(1:nlay*nel*l_dirb)
171 IF (irep == 0) THEN
172 DO j=1,nlay
173 j1 = 1+(j-1)*l_dira*nel
174 j2 = j*l_dira*nel
175 dira(j1:j2) = elbuf_str%BUFLY(j)%DIRA(1:nel*l_dira)
176 ENDDO
177 ENDIF
178C
179 CALL cbacoork(jft,jlt,x,ixc,pm,gbuf%OFF,
180 1 geo,area,vcore,jac,hx,hy,
181 2 vqn,vqg,vq,vjfi,vnrm,vastn,nplat,iplat,
182 3 x13 ,x24 ,y13,y24,
183 4 elbuf_str,nlay, gbuf%SMSTR,
184 5 irep,npt,ismstr,dir_a,dir_b ,
185 6 pid ,mat,ngl,off,idril,nel)
186 CALL cbaini3(jft,jlt,vksi,veta,
187 1 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
188 2 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
189 3 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
190 4 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34)
191C
192 IF (irep>0) THEN
193 CALL cmatc3(jft ,jlt ,pm ,mat ,geo ,
194 1 pid ,area ,thk0 ,thk2 ,gbuf%THK,
195 2 thke ,volg ,mtn ,npt ,ithk ,
196 3 hm ,hf ,hc ,hz ,igtyp ,
197 4 iorth ,hmor ,hfor ,dir_a ,igeo ,
198 5 idril ,ihbe ,hmfor ,gs ,isubstack,
199 6 stack ,elbuf_str,nlay ,drape_sh4n ,nft ,
200 7 nel ,indx_drape,sedrape,numel_drape)
201 ELSE
202 CALL cmatc3(jft ,jlt ,pm ,mat ,geo ,
203 1 pid ,area ,thk0 ,thk2 ,gbuf%THK,
204 2 thke ,volg ,mtn ,npt ,ithk ,
205 3 hm ,hf ,hc ,hz ,igtyp ,
206 4 iorth ,hmor ,hfor ,dira ,igeo ,
207 5 idril ,ihbe ,hmfor ,gs ,isubstack,
208 6 stack ,elbuf_str,nlay ,drape_sh4n , nft ,
209 7 nel ,indx_drape,sedrape,numel_drape)
210 ENDIF
211 IF (idril>0) THEN
212 CALL cbaderirz(jft ,jlt ,area ,x13 ,x24 ,
213 2 y13 ,y24 ,bm0rz,bmkrz,bmerz,
214 3 vcore,nplat,iplat,ismstr)
215 ELSE
216 CALL cbabec3(jft ,jlt ,x13 ,x24 ,y13 ,y24 ,bm, nplat, iplat)
217 END IF
218C-----------------------------------------------
219C BOUCLE SUR POINTS D'INTEGRATION DE GAUSS
220C-----------------------------------------------
221 lenf = nel*gbuf%G_FORPG/npg
222 lenm = nel*gbuf%G_MOMPG/npg
223 nptr = elbuf_str%NPTR
224 npts = elbuf_str%NPTS
225 DO is = 1,npts
226 DO ir = 1,nptr
227 ng = nptr*(is-1) + ir
228 ptf = (ng-1)*lenf+1
229 ptm = (ng-1)*lenm+1
230 DO i=jft,jlt
231 cdet(i) = jac(i,ng)
232 vol0(i) = thk0(i)*cdet(i)
233 ENDDO
234C-----------------------------------------------
235C MATRICE [B]
236C-----------------------------------------------
237 CALL cbabe3(jft,jlt,ng,vcore,area,cdet,vqn,vqg,vjfi,
238 1 vnrm,vastn,hx,hy,veta,vksi,
239 2 bm,bmf,bf,bc,tc,bzz,nplat,iplat,
240 3 idril,brz )
241C-----------------------------------------------
242C IF [KT]
243C-----------------------------------------------
244 CALL cmatip3(jft ,jlt ,pm ,mat ,pid ,
245 1 mtn ,npt ,hm ,hf ,iorth ,
246 2 hmor ,hfor ,hmfor ,ng )
247C----------------------------------------------------------------------------
248C SUB-MATRICES [KE] LOCAL
249C----------------------------
250 CALL cbalke3(jft,jlt,cdet,thk0,thk2,hm,hf,hc,hz,
251 1 bm,bmf,bf,bc,tc,bzz,nplat,iplat,vol0,
252 2 ikgeo,gbuf%FORPG(ptf),gbuf%MOMPG(ptm),
253 3 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
254 4 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
255 5 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
256 6 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
257 7 iorth,hmor,hfor,idril,hmfor,
258 8 x13 ,x24 ,y13 ,y24,nel)
259 IF (idril>0) THEN
260 CALL cbaber3(jft ,jlt ,bm0rz,bmkrz,bmerz ,
261 2 bmrz ,brz ,bm ,nplat ,iplat,
262 3 ng )
263 CALL cbalkerz(jft ,jlt ,vol0 ,thk0 ,
264 2 hm ,hz ,bm ,
265 6 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
266 7 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
267 8 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
268 9 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
269 a iorth,hmor,hfor ,iplat,nplat,
270 b bmrz,brz ,gbuf%HOURG,ikgeo,ng ,hmfor,bf ,
271 c bmf ,nel)
272 END IF
273 ENDDO
274 ENDDO
275C---------FIN DE BOUCLE DE 4 POINTS DE GAUSS------------
276C-------------membrane shear traitement--------------------
277 IF (idril==0) THEN
278 CALL cbalkec3(jft,jlt,volg ,x13 ,x24 ,y13 ,y24, hm,
279 1 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
280 2 nplat,iplat,ikgeo,gbuf%FOR,m11,m22,m33,m44,
281 3 iorth,nel)
282 END IF
283C----------------------------
284C TRANSFORME [KE] LOCALE AU GLOBALE et ASSEMBLAGE----
285C----------------------------
286 CALL cbasumg3(
287 1 jft ,jlt ,vqn ,vq ,nplat ,iplat ,
288 2 k11,k12,k13,k14,k22,k23,k24,k33,k34,k44,
289 3 m11,m12,m13,m14,m22,m23,m24,m33,m34,m44,
290 4 mf11,mf12,mf13,mf14,mf22,mf23,mf24,mf33,
291 5 mf34,mf44,fm12,fm13,fm14,fm23,fm24,fm34,
292 6 ke11,ke22,ke33,ke44,ke12,ke13,ke14,ke23,
293 7 ke24,ke34,vcore,idril,iorth)
294
295 IF (neig>0) CALL c4eoff(
296 1 jft, jlt, ixc, etag, off)
297C
298 CALL assem_c4(
299 1 ixc ,nel ,iddl ,ndof ,k_diag ,
300 2 k_lt ,iadk ,jdik ,ke11 ,ke12 ,
301 3 ke13 ,ke14 ,ke22 ,ke23 ,ke24 ,
302 5 ke33 ,ke34 ,ke44 ,off )
303C
304 RETURN
305 END
306!||====================================================================
307!|| cbaini3 ../engine/source/elements/shell/coqueba/cbake3.F
308!||--- called by ------------------------------------------------------
309!|| cbake3 ../engine/source/elements/shell/coqueba/cbake3.F
310!||====================================================================
311 SUBROUTINE cbaini3(JFT,JLT,VKSI,VETA,
312 1 K11,K12,K13,K14,K22,K23,K24,K33,K34,K44,
313 2 M11,M12,M13,M14,M22,M23,M24,M33,M34,M44,
314 3 MF11,MF12,MF13,MF14,MF22,MF23,MF24,MF33,
315 4 MF34,MF44,FM12,FM13,FM14,FM23,FM24,FM34)
316C-----------------------------------------------
317C I M P L I C I T T Y P E S
318C-----------------------------------------------
319#include "implicit_f.inc"
320c-----------------------------------------------
321c g l o b a l p a r a m e t e r s
322c-----------------------------------------------
323C-----------------------------------------------
324C D U M M Y A R G U M E N T S
325C-----------------------------------------------
326 INTEGER JFT,JLT
327 MY_REAL
328 . VETA(4,4),VKSI(4,4)
329 my_real
330 . K11(9,*),K12(9,*),K13(9,*),K14(9,*),
331 . K22(9,*),K23(9,*),K24(9,*),K33(9,*),
332 . M11(9,*),M12(9,*),M13(9,*),M14(9,*),
333 . M22(9,*),M23(9,*),M24(9,*),M33(9,*),
334 . MF11(9,*),MF12(9,*),MF13(9,*),MF14(9,*),
335 . MF22(9,*),MF23(9,*),MF24(9,*),MF33(9,*),
336 . FM12(9,*),FM13(9,*),FM14(9,*),
337 . FM23(9,*),FM24(9,*),FM34(9,*),
338 . k34(9,*),k44(9,*),m34(9,*),m44(9,*),
339 . mf34(9,*),mf44(9,*)
340C-----------------------------------------------
341C L O C A L V A R I A B L E S
342C-----------------------------------------------
343 INTEGER I,J
344 my_real
345 . PG
346 DATA
347 . PG/.577350269189626/
348C-----------------------------------------------
349 DO J=1,9
350 do i=jft,jlt
351 k11(j,i)=zero
352 k12(j,i)=zero
353 k13(j,i)=zero
354 k14(j,i)=zero
355 k22(j,i)=zero
356 k23(j,i)=zero
357 k24(j,i)=zero
358 k33(j,i)=zero
359 k34(j,i)=zero
360 k44(j,i)=zero
361 m11(j,i)=zero
362 m12(j,i)=zero
363 m13(j,i)=zero
364 m14(j,i)=zero
365 m22(j,i)=zero
366 m23(j,i)=zero
367 m24(j,i)=zero
368 m33(j,i)=zero
369 m34(j,i)=zero
370 m44(j,i)=zero
371 ENDDO
372 ENDDO
373 DO j=1,9
374 DO i=jft,jlt
375 mf11(j,i)=zero
376 mf12(j,i)=zero
377 mf13(j,i)=zero
378 mf14(j,i)=zero
379 mf22(j,i)=zero
380 mf23(j,i)=zero
381 mf24(j,i)=zero
382 mf33(j,i)=zero
383 mf34(j,i)=zero
384 mf44(j,i)=zero
385 fm12(j,i)=zero
386 fm13(j,i)=zero
387 fm14(j,i)=zero
388 fm23(j,i)=zero
389 fm24(j,i)=zero
390 fm34(j,i)=zero
391 ENDDO
392 ENDDO
393C
394 vksi(1,1)=-fourth*(one+pg)
395 vksi(2,1)=-vksi(1,1)
396 vksi(3,1)= fourth*(one-pg)
397 vksi(4,1)=-vksi(3,1)
398 veta(1,1)=-fourth*(one+pg)
399 veta(2,1)=-fourth*(one-pg)
400 veta(3,1)=-veta(2,1)
401 veta(4,1)=-veta(1,1)
402 vksi(1,2)= vksi(1,1)
403 vksi(2,2)=-vksi(1,2)
404 vksi(3,2)= vksi(3,1)
405 vksi(4,2)=-vksi(3,2)
406 veta(1,2)= veta(2,1)
407 veta(2,2)= veta(1,1)
408 veta(3,2)=-veta(2,2)
409 veta(4,2)=-veta(1,2)
410 vksi(1,3)=-vksi(3,1)
411 vksi(2,3)=-vksi(1,3)
412 vksi(3,3)=-vksi(1,1)
413 vksi(4,3)=-vksi(3,3)
414 veta(1,3)= veta(1,2)
415 veta(2,3)= veta(2,2)
416 veta(3,3)=-veta(2,3)
417 veta(4,3)=-veta(1,3)
418 vksi(1,4)= vksi(1,3)
419 vksi(2,4)=-vksi(1,4)
420 vksi(3,4)= vksi(3,3)
421 vksi(4,4)=-vksi(3,4)
422 veta(1,4)= veta(1,1)
423 veta(2,4)= veta(2,1)
424 veta(3,4)=-veta(2,4)
425 veta(4,4)=-veta(1,4)
426C
427 RETURN
428 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
cbabec3
subroutine cbabec3(jft, jlt, x13, x24, y13, y24, bm, nplat, iplat)
Definition cbabe3.F:642
cbaber3
subroutine cbaber3(jft, jlt, bm0rz, bmkrz, bmerz, bmrz, brz, bm, nplat, iplat, ng)
Definition cbabe3.F:685
cbabe3
subroutine cbabe3(jft, jlt, ng, vcore, area, cdet, vqn, vq, vjfi, vnrm, vastn, hx, hy, veta, vksi, bm, bmf, bf, bc, tc, bzz, nplat, iplat, isrot, brz)
Definition cbabe3.F:38
cbacoork
subroutine cbacoork(jft, jlt, x, ixc, pm, offg, geo, area, vcore, jac, hx, hy, vqn, vqg, vq, vjfi, vnrm, vastn, nplat, iplat, x13_t, x24_t, y13_t, y24_t, elbuf_str, nlay, smstr, irep, npt, ismstr, dir_a, dir_b, pid, mat, ngl, off, isrot, nel)
Definition cbacoork.F:40
cbaderirz
subroutine cbaderirz(jft, jlt, area, x13, x24, y13, y24, bm0rz, bmkrz, bmerz, vcore, nplat, iplat, ismstr)
Definition cbadef.F:802
cbake3
subroutine cbake3(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 cbake3.F:57
cbaini3
subroutine cbaini3(jft, jlt, vksi, veta, 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 cbake3.F:316
cbalkec3
subroutine cbalkec3(jft, jlt, vol, x13, x24, y13, y24, hm, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, nplat, iplat, ikgeo, for, m11, m22, m33, m44, iorth, nel)
Definition cbalke3.F:1638
cbalke3
subroutine cbalke3(jft, jlt, cdet, thk0, thk2, hm, hf, hc, hz, bm, bmf, bf, bc, tc, bzz, nplat, iplat, vol, ikgeo, 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, iorth, hmor, hfor, idril, hmfor, x13, x24, y13, y24, nel)
Definition cbalke3.F:42
cbalkerz
subroutine cbalkerz(jft, jlt, vol, thk0, hm, hz, bm, 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, nplat, pmrz, brz, frz, ikgeo, ng, hmfor, bf, bmf, nel)
Definition cbalke3.F:1766
cbasumg3
subroutine cbasumg3(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, vcore, idril, iorth)
Definition cbasumg3.F:43
cmatip3
subroutine cmatip3(jft, jlt, pm, mat, pid, mtn, npt, hm, hf, iorth, hmor, hfor, hmfor, ipg)
Definition cmatc3.F:582
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
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
max
#define max(a, b)
Definition macros.h:21
drape_mod
Definition drape_mod.F:54
stack_mod
Definition stack_mod.F:72