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

Go to the source code of this file.

Functions/Subroutines

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)
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)

Function/Subroutine Documentation

◆ cbaini3()

subroutine cbaini3 ( integer jft,
integer 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 at line 313 of file cbake3.F.

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
#define my_real
Definition cppsort.cpp:32

◆ cbake3()

subroutine cbake3 ( 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 49 of file cbake3.F.

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
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
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
subroutine c4eoff(jft, jlt, ixc, etag, off)
Definition c4eoff.F:32
subroutine cbabec3(jft, jlt, x13, x24, y13, y24, bm, nplat, iplat)
Definition cbabe3.F:642
subroutine cbaber3(jft, jlt, bm0rz, bmkrz, bmerz, bmrz, brz, bm, nplat, iplat, ng)
Definition cbabe3.F:685
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
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
subroutine cbaderirz(jft, jlt, area, x13, x24, y13, y24, bm0rz, bmkrz, bmerz, vcore, nplat, iplat, ismstr)
Definition cbadef.F:779
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
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
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
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
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
subroutine cmatip3(jft, jlt, pm, mat, pid, mtn, npt, hm, hf, iorth, hmor, hfor, hmfor, ipg)
Definition cmatc3.F:582
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
subroutine area(d1, x, x2, y, y2, eint, stif0)
#define max(a, b)
Definition macros.h:21