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