OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
c3init3.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!|| c3init3 ../starter/source/elements/sh3n/coque3n/c3init3.F
25!||--- called by ------------------------------------------------------
26!|| initia ../starter/source/elements/initia/initia.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| c1buf3 ../starter/source/elements/shell/coque/c1buf3.F
30!|| c3coori ../starter/source/elements/sh3n/coque3n/c3coori.F
31!|| c3derii ../starter/source/elements/sh3n/coque3n/c3derii.F
32!|| c3eps11_ini ../starter/source/elements/sh3n/coque3n/c3init3.F
33!|| c3epschk ../starter/source/elements/sh3n/coque3n/c3init3.F
34!|| c3epsini ../starter/source/elements/sh3n/coque3n/c3epsini.F
35!|| c3evec3 ../starter/source/elements/sh3n/coque3n/c3evec3.F
36!|| c3inmas ../starter/source/elements/sh3n/coque3n/c3inmas.F
37!|| c3veok3 ../starter/source/elements/sh3n/coque3n/c3veok3.F
38!|| cbufxfe ../starter/source/elements/xfem/cbufxfe.F
39!|| cdkevec3 ../starter/source/elements/sh3n/coquedk/cdkevec3.F
40!|| cfailini ../starter/source/elements/shell/coque/cfailini.F
41!|| cm27in3 ../starter/source/materials/mat/mat027/cm27in3.F
42!|| cm35in3 ../starter/source/materials/mat/mat035/cm35in3.F
43!|| cm58in3 ../starter/source/materials/mat/mat058/cm58in3.f
44!|| cmatini ../starter/source/materials/mat_share/cmatini.F
45!|| corth3 ../starter/source/elements/shell/coque/corth3.F
46!|| corthdir ../starter/source/elements/shell/coque/corthdir.F
47!|| corthini ../starter/source/elements/shell/coque/corthini.F
48!|| csigini ../starter/source/elements/shell/coque/csigini.F
49!|| cuserini ../starter/source/elements/shell/coque/cuserini.F
50!|| fail_windshield_init ../starter/source/materials/fail/windshield_alter/fail_windshield_init.F
51!|| fractal_dmg_init ../starter/source/materials/fail/fractal/fractal_dmg_init.F90
52!|| fretitl2 ../starter/source/starter/freform.F
53!|| initemp_shell ../starter/source/materials/therm/initemp_shell.F90
54!|| law158_init ../starter/source/materials/mat/mat158/law158_init.F
55!|| law42c_ini ../starter/source/materials/mat/mat042/law42c_ini.F90
56!|| layini1 ../starter/source/elements/shell/coqueba/layini1.F
57!|| thickini ../starter/source/elements/shell/coqueba/thickini.F
58!||--- uses -----------------------------------------------------
59!|| drape_mod ../starter/share/modules1/drape_mod.F
60!|| fractal_dmg_init_mod ../starter/source/materials/fail/fractal/fractal_dmg_init.F90
61!|| initemp_shell_mod ../starter/source/materials/therm/initemp_shell.F90
62!|| law42c_ini_mod ../starter/source/materials/mat/mat042/law42c_ini.F90
63!|| message_mod ../starter/share/message_module/message_mod.F
64!|| stack_mod ../starter/share/modules1/stack_mod.F
65!||====================================================================
66 SUBROUTINE c3init3(ELBUF_STR,
67 1 IXTG ,PM ,X ,GEO ,IGRSH4N ,
68 2 XMAS ,IN ,NVC ,DTELEM ,IGRSH3N ,
69 3 XREFTG ,OFFSET ,NEL ,ITHK ,THK ,
70 4 ISIGSH ,SIGSH ,STIFN ,STIFR ,PARTSAV ,
71 5 V ,IPART ,MSTG ,INTG ,PTG ,
72 8 SKEW ,IPARG ,NSIGSH ,IGEO ,IUSER ,
73 9 ETNOD ,NSHNOD ,STTG ,PTSH3N ,IPM ,
74 A BUFMAT ,SH3TREE,MCP ,MCPTG ,TEMP ,
75 B CPT_ELTENS,PART_AREA,ITAGE,ITAGN,IXFEM ,
76 C NPF , TF ,SH3TRIM ,XFEM_STR,ISUBSTACK,
77 D STACK ,RNOISE ,DRAPE ,SH3ANG ,IDDLEVEL,
78 E GEO_STACK,IGEO_STACK,STRTG,PERTURB ,ISH3N,
79 F IYLDINI ,ELE_AREA,NLOC_DMG,NG,GROUP_PARAM,
80 G IDRAPE , DRAPEG,MAT_PARAM,FAIL_FRACTAL,FAIL_BROKMANN,GLOB_THERM)
81C-----------------------------------------------
82C M o d u l e s
83C-----------------------------------------------
84 USE my_alloc_mod
85 USE mat_elem_mod
86 USE message_mod
87 USE stack_mod
88 USE groupdef_mod
89 USE nlocal_reg_mod
90 USE group_param_mod
91 USE drape_mod
92 USE random_walk_def_mod
93 USE fractal_dmg_init_mod
94 USE names_and_titles_mod , ONLY : nchartitle
95 use brokmann_random_def_mod
96 use glob_therm_mod
97 use initemp_shell_mod
98 use element_mod , only : nixtg
99 use law42c_ini_mod
100C-----------------------------------------------
101C I m p l i c i t T y p e s
102C-----------------------------------------------
103#include "implicit_f.inc"
104C-----------------------------------------------
105C G l o b a l P a r a m e t e r s
106C-----------------------------------------------
107#include "mvsiz_p.inc"
108C-----------------------------------------------
109C C o m m o n B l o c k s
110C-----------------------------------------------
111#include "com01_c.inc"
112#include "com04_c.inc"
113#include "com_xfem1.inc"
114#include "param_c.inc"
115#include "scr03_c.inc"
116#include "scr17_c.inc"
117#include "vect01_c.inc"
118#include "scry_c.inc"
119C-----------------------------------------------
120C V a r i a b l e s
121C-----------------------------------------------
122 INTEGER OFFSET,NEL,ITHK,ISIGSH,NSIGSH,IUSER,
123 . CPT_ELTENS,ISUBSTACK,IYLDINI,ISH3N,NG,IDRAPE
124 INTEGER IXTG(NIXTG,*),IPART(*),IGEO(NPROPGI,*),IPM(NPROPMI,*),
125 . IPARG(*),NSHNOD(*), PTSH3N(*),NPF(*),
126 . SH3TREE(*),ITAGE(*),ITAGN(*),IXFEM,SH3TRIM(*),
127 . IGEO_STACK(*),PERTURB(NPERTURB)
128 INTEGER ,INTENT(IN) :: IDDLEVEL
129C REAL
130 my_real
131 . PM(*),X(3,*),GEO(NPROPG,*),XMAS(*),XREFTG(3,3,*),
132 . IN(*),DTELEM(*),THK(*),SIGSH(NSIGSH,*),
133 . STIFN(*),STIFR(*),PARTSAV(20,*), V(*), SKEW(LSKEW,*),
134 . MSTG(*),INTG(*),PTG(3,*),
135 . ETNOD(*), STTG(*),BUFMAT(*),MCP(*),MCPTG(*),TEMP(*),
136 . part_area(*),tf(*),rnoise(*),sh3ang(*),
137 . geo_stack(*),strtg(*),ele_area(*)
138 TYPE(elbuf_struct_), TARGET :: ELBUF_STR
139 TYPE(ELBUF_STRUCT_), TARGET ,DIMENSION(NGROUP,*):: XFEM_STR
140 ! when XFEM is ON, XFEM_STR's dimension = NGROUP,NXEL
141 TYPE (STACK_PLY) :: STACK
142 TYPE (NLOCAL_STR_) :: NLOC_DMG
143 TYPE (GROUP_PARAM_) :: GROUP_PARAM
144 TYPE (DRAPE_) :: DRAPE(NUMELC_DRAPE + NUMELTG_DRAPE)
145 TYPE (DRAPEG_) :: DRAPEG
146 TYPE (MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(INOUT) :: MAT_PARAM
147 TYPE (FAIL_FRACTAL_) ,INTENT(IN) :: FAIL_FRACTAL
148 TYPE (FAIL_BROKMANN_) ,INTENT(IN) :: FAIL_BROKMANN
149 TYPE (glob_therm_) ,intent(in) :: glob_therm
150C-----------------------------------------------
151 TYPE (GROUP_) , DIMENSION(NGRSHEL) :: IGRSH4N
152 TYPE (GROUP_) , DIMENSION(NGRSH3N) :: IGRSH3N
153C-----------------------------------------------
154C L o c a l V a r i a b l e s
155C-----------------------------------------------
156 INTEGER I,J,NDEPAR,IGTYP,NVC,NUVAR,NLAY,IR,IS,IL,IFAIL,NUPARAM,
157 . nptr,npts,nptt,ixel,ii,it,ilaw,imat,iprop,irep,itg,igmat,
158 . ifram_old,npt_all,mpt,laynpt_max,lay_max
159 INTEGER IORTHLOC(MVSIZ),MAT(MVSIZ),PID(MVSIZ),NGL(MVSIZ),JJ(6),
160 . ix1(mvsiz),ix2(mvsiz),ix3(mvsiz)
161 my_real
162 . vx(mvsiz),vy(mvsiz),vz(mvsiz),aldt(mvsiz),area(mvsiz)
163 my_real,
164 . DIMENSION(MVSIZ) :: px1g,py1g,py2g,x2s,x3s,y3s,dt
165 my_real x1(mvsiz), x2(mvsiz), x3(mvsiz) ,x4(mvsiz),
166 . y1(mvsiz), y2(mvsiz), y3(mvsiz),y4(mvsiz),
167 . z1(mvsiz), z2(mvsiz), z3(mvsiz),z4(mvsiz),
168 . e1x(mvsiz), e1y(mvsiz), e1z(mvsiz),
169 . e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
170 . e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),
171 . x31(mvsiz), y31(mvsiz), z31(mvsiz),
172 . x2l(mvsiz), x3l(mvsiz), y3l(mvsiz)
173 INTEGER ID
174 CHARACTER(LEN=NCHARTITLE)::TITR
175 my_real,
176 . ALLOCATABLE, DIMENSION(:) :: dir_a,dir_b,phi1,phi2,
177 . coor1,coor2,coor3,coor4
178 INTEGER, ALLOCATABLE, DIMENSION(:) :: INDX
179 parameter(laynpt_max = 10)
180 parameter(lay_max = 100)
181 INTEGER, DIMENSION(:),ALLOCATABLE::MATLY
182 my_real, DIMENSION(:,:),ALLOCATABLE :: POSLY
183C-----------------------------------------------
184 my_real,
185 . DIMENSION(:) ,POINTER :: UVAR,DIR1,DIR2
186 TYPE(G_BUFEL_) ,POINTER :: GBUF
187 TYPE(l_bufel_) ,POINTER :: LBUF
188 TYPE(buf_lay_) ,POINTER :: BUFLY
189C=======================================================================
190 ! Allocate buffers
191 CALL my_alloc(matly,mvsiz*lay_max)
192 CALL my_alloc(posly,mvsiz,lay_max*laynpt_max)
193
194 gbuf => elbuf_str%GBUF
195 bufly => elbuf_str%BUFLY(1)
196c
197 imat = ixtg(1,1+nft) ! mat N
198 iprop = ixtg(nixtg-1,1+nft) ! property N
199 igtyp = nint(geo(12,iprop))
200 id = igeo(1,iprop)
201 igmat = igeo(98,iprop)
202 irep = iparg(35)
203 ifail = iparg(43)
204c
205 IF (ish3n==3.AND.ish3nfram==0) THEN
206 ifram_old =0
207 ELSE
208 ifram_old =1
209 END IF
210C
211 CALL fretitl2(titr,igeo(npropgi-ltitr+1,iprop),ltitr)
212 nuvar = ipm(8,ixtg(1,1+nft))
213 vx = zero
214 vy = zero
215 vz = zero
216 iorthloc = 0
217 itg = 1+numelc
218C---
219 ir = 1
220 is = 1
221 nlay = elbuf_str%NLAY
222 nxel = elbuf_str%NXEL
223 nptt = elbuf_str%NPTT
224C
225 npt_all = 0
226 DO il=1,nlay
227 npt_all = npt_all + elbuf_str%BUFLY(il)%NPTT
228 ENDDO
229 mpt = max(1,npt_all)
230 IF(npt_all == 0) npt_all = nlay
231 IF (iparg(6) == 0.OR.npt==0) mpt=0
232 IF((igtyp == 51 .OR. igtyp == 52) .AND. idrape > 0) THEN
233 ALLOCATE(dir_a(npt_all*nel*2))
234 ALLOCATE(dir_b(npt_all*nel*2))
235 dir_a = zero
236 dir_b = zero
237 ALLOCATE(phi1(mvsiz*npt_all))
238 ALLOCATE(phi2(nvsiz*npt_all))
239 phi1 = zero
240 phi2 = zero
241 ALLOCATE(coor1(npt_all*mvsiz))
242 ALLOCATE(coor2(npt_all*mvsiz))
243 ALLOCATE(coor3(npt_all*mvsiz))
244 ALLOCATE(coor4(npt_all*mvsiz))
245 coor1 = zero
246 coor2 = zero
247 coor3 = zero
248 coor4 = zero
249 ELSE
250 ALLOCATE(dir_a(nlay*nel*2))
251 ALLOCATE(dir_b(nlay*nel*2))
252 dir_a = zero
253 dir_b = zero
254 ALLOCATE(phi1(nlay*mvsiz))
255 ALLOCATE(phi2(nlay*mvsiz))
256 phi1 = zero
257 phi2 = zero
258 ALLOCATE(coor1(nlay*mvsiz))
259 ALLOCATE(coor2(nlay*mvsiz))
260 ALLOCATE(coor3(nlay*mvsiz))
261 ALLOCATE(coor4(nlay*mvsiz))
262 coor1 = zero
263 coor2 = zero
264 coor3 = zero
265 coor4 = zero
266 npt_all = nlay
267 ENDIF
268 !
269C
270 DO j=1,6
271 jj(j) = nel*(j-1)
272 ENDDO
273C
274 DO i=lft,llt
275 mat(i) = imat
276 pid(i) = iprop
277 ENDDO
278C
279 IF (ixfem > 0) THEN
280 DO i=lft,llt
281 itagn(ixtg(2,i+nft)) =1
282 itagn(ixtg(3,i+nft)) =1
283 itagn(ixtg(4,i+nft)) =1
284 itage(i+nft) = 1
285 ENDDO
286 ENDIF
287C
288 CALL c3coori(x,xreftg(1,1,nft+1),ixtg(1,nft+1),ngl,
289 . x1 ,x2 ,x3 ,y1 ,y2 ,y3 ,
290 . z1 ,z2 ,z3 ,ix1 ,ix2 ,ix3 )
291 CALL c3veok3(nvc ,ix1 ,ix2 ,ix3 )
292 CALL c3evec3(lft ,llt ,area,
293 . x1 ,x2 ,x3 ,y1 ,y2 ,y3 ,
294 . z1 ,z2 ,z3 ,e1x ,e2x ,e3x ,
295 . e1y ,e2y ,e3y ,e1z ,e2z ,e3z ,
296 . x31, y31, z31 ,x2l ,x3l ,y3l )
297C------------
298C Tags total area of the part (needed in /ADMAS for shells)
299C------------
300 IF ((imasadd > 0).OR.(nloc_dmg%IMOD > 0)) THEN
301 DO i=lft,llt
302 j = ipart(i+nft)
303C PART_AREA(J) = PART_AREA(J) + AREA(I)
304 ele_area(numelc+i+nft) = area(i)
305 IF (gbuf%G_AREA > 0) gbuf%AREA(i) = area(i)
306 ENDDO
307 ENDIF
308C-----------------------------------------------------------------------
309! Initialize element temperature from /initemp
310!
311 IF (jthe == 0 .and. glob_therm%NINTEMP > 0) THEN
312 CALL initemp_shell(elbuf_str,temp,nel,numnod,numeltg,3,nixtg,ixtg)
313 END IF
314C------------
315 IF(idrape > 0 ) THEN
316 ALLOCATE(indx(numeltg))
317 indx = 0
318 indx(1:numeltg) = drapeg%INDX(numelc + 1 : numelc + numeltg)
319 ELSE
320 ALLOCATE(indx(0))
321 ENDIF
322 CALL c3inmas(x,xreftg(1,1,nft+1),ixtg,geo,pm,xmas,in,thk,
323 . partsav,v,ipart(nft+1),mstg(nft+1),intg(nft+1),
324 . ptg(1,nft+1),igeo ,imat ,iprop ,area ,
325 . etnod,nshnod,sttg(nft+1) ,sh3tree ,mcp ,
326 . mcptg(nft+1),temp ,sh3trim,isubstack,nlay ,
327 . elbuf_str ,stack ,gbuf%THK_I,rnoise,drape,
328 . perturb,ix1 ,ix2 ,ix3 ,glob_therm%NINTEMP,
329 . x2l ,x3l ,y3l ,idrape,indx)
330C-----------------------------------------------
331 CALL c3derii(lft,llt,pm,geo,px1g,py1g,py2g,
332 . stifn ,stifr ,ixtg(1,nft+1),
333 . thk,sh3tree,aldt ,bufmat , ipm ,igeo,
334 . stack%PM,isubstack,strtg(nft+1),imat,iprop,
335 . area ,dt ,x31 ,y31 ,z31 ,
336 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,
337 . e1z ,e2z ,e3z ,x2l ,x3l ,y3l ,
338 . group_param)
339C
340 CALL c1buf3(geo,gbuf%THK,gbuf%OFF,thk,ksh3tree,sh3tree)
341C
342 IF (ixfem > 0) THEN
343 DO ixel=1,nxel
344 DO i=lft,llt
345 xfem_str(ng,ixel)%GBUF%THK(i) = thk(i)
346 xfem_str(ng,ixel)%GBUF%OFF(i) = -one
347 END DO
348 ENDDO
349 ENDIF
350C-------new local system for ortho or aniso init
351 IF (ifram_old ==0 )
352 . CALL cdkevec3(lft ,llt ,area,
353 . x1 ,x2 ,x3 ,y1 ,y2 ,y3 ,
354 . z1 ,z2 ,z3 ,e1x ,e2x ,e3x ,
355 . e1y ,e2y ,e3y ,e1z ,e2z ,e3z )
356C----
357C PHI, COOR used only with dimension(NLAY,MVSIZ)
358c
359 CALL corthini(
360 . lft ,llt ,nft ,nlay ,numeltg ,
361 . nsigsh ,nixtg ,ixtg(1,nft+1),igeo ,geo ,
362 . skew ,sigsh ,ptsh3n ,phi1 ,phi2 ,
363 . vx ,vy ,vz ,coor1 ,coor2 ,
364 . coor3 ,coor4 ,iorthloc ,isubstack ,stack ,
365 . irep ,elbuf_str ,drape ,sh3ang(nft+1),x ,
366 . geo_stack ,e3x ,e3y ,e3z ,
367 . gbuf%BETAORTH,x1 ,x2 ,y1 ,y2 ,
368 . z1 ,z2 ,nel ,gbuf%G_ADD_NODE,gbuf%ADD_NODE,
369 . npt_all ,idrape ,indx)
370c---
371 IF(igtyp == 51 .OR. igtyp == 52 .OR. igmat > 0) THEN
372c
373 CALL corthdir(elbuf_str,
374 . igeo ,geo ,vx ,vy ,vz ,
375 . phi1 ,phi2 ,coor1 ,coor2 ,coor3 ,
376 . coor4 ,iorthloc ,nlay ,irep ,isubstack,
377 . stack ,geo_stack ,igeo_stack ,ir ,is ,
378 . nel ,imat ,iprop ,
379 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
380 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
381 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z ,
382 . npt_all ,idrape)
383C
384 ELSEIF (mtn == 27) THEN
385 CALL cm27in3(elbuf_str,
386 . geo ,igeo ,pm ,ipm ,ixtg(1,1+nft) ,nixtg,
387 . nlay,ir ,is ,imat )
388 ELSEIF (mtn == 35) THEN
389 nptr = elbuf_str%NPTR
390 npts = elbuf_str%NPTS
391 nptt = elbuf_str%NPTT
392 CALL cm35in3(elbuf_str,thk,area,nel,nlay,
393 . nptr,npts,nptt,igtyp)
394 ELSEIF (mtn == 15 .or. mtn == 19 .or. mtn == 25 .or. mtn >= 28)THEN
395 IF (mtn == 19 .AND. igtyp /= 9) THEN
396 CALL ancmsg(msgid=5,
397 . anmode=aninfo,
398 . msgtype=msgerror,
399 . i1=igeo(1,ixtg(nixtg-1,nft+1)))
400 ENDIF
401c
402 CALL corthdir(elbuf_str,
403 . igeo ,geo ,vx ,vy ,vz ,
404 . phi1 ,phi2 ,coor1 ,coor2 ,coor3 ,
405 . coor4 ,iorthloc ,nlay ,irep ,isubstack,
406 . stack ,geo_stack ,igeo_stack ,ir ,is ,
407 . nel ,imat ,iprop ,
408 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
409 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
410 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z ,
411 . npt_all ,idrape)
412 ENDIF
413c
414 IF ((mtn == 58 .or. mtn == 158) .AND.
415 . igtyp /= 16 .AND. igtyp /= 51 .AND. igtyp /= 52) THEN
416 CALL ancmsg(msgid=658,
417 . msgtype=msgerror,
418 . anmode=aninfo_blind_1,
419 . i1=id,
420 . c1=titr,
421 . i2=mtn,
422 . i3=igtyp)
423 ELSEIF (mtn == 58 .or. mtn == 158 .OR. igtyp == 51 .OR. igtyp == 52) THEN
424 IF (idrape == 0 ) THEN
425 DO il = 1,nlay
426 nptt = elbuf_str%BUFLY(il)%NPTT
427 imat = elbuf_str%BUFLY(il)%IMAT
428 ilaw = elbuf_str%BUFLY(il)%ILAW
429 nuvar = elbuf_str%BUFLY(il)%NVAR_MAT
430 dir1 => elbuf_str%BUFLY(il)%DIRA
431 dir2 => elbuf_str%BUFLY(il)%DIRB
432 nuparam = mat_param(imat)%NUPARAM
433!
434 IF (ilaw == 58) THEN
435 DO it=1,nptt
436 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
437 uvar => elbuf_str%BUFLY(il)%MAT(ir,is,it)%VAR
438 CALL cm58in3(
439 . irep ,dir1 ,dir2 ,mat_param(imat)%UPARAM,
440 . uvar ,aldt ,nel ,nuvar ,lbuf%ANG ,
441 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
442 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
443 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
444 ENDDO
445 ELSE IF (ilaw == 158) THEN
446 DO it=1,nptt
447 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
448 uvar => elbuf_str%BUFLY(il)%MAT(ir,is,it)%VAR
449 CALL law158_init(dir1 ,dir2 ,
450 . uvar ,aldt ,nel ,nuvar ,lbuf%ANG ,
451 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
452 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
453 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
454 ENDDO
455 ENDIF
456!
457 ENDDO ! DO IL = 1,NLAY
458 ELSE
459 DO il = 1,nlay
460 nptt = elbuf_str%BUFLY(il)%NPTT
461 imat = elbuf_str%BUFLY(il)%IMAT
462 ilaw = elbuf_str%BUFLY(il)%ILAW
463 nuvar = elbuf_str%BUFLY(il)%NVAR_MAT
464 nuparam = mat_param(imat)%NUPARAM
465!
466 IF (ilaw == 58) THEN
467 DO it=1,nptt
468 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
469 uvar => elbuf_str%BUFLY(il)%MAT(ir,is,it)%VAR
470 dir1 => elbuf_str%BUFLY(il)%LBUF_DIR(it)%DIRA
471 dir2 => elbuf_str%BUFLY(il)%LBUF_DIR(it)%DIRB
472 CALL cm58in3(
473 . irep ,dir1 ,dir2 ,mat_param(imat)%UPARAM,
474 . uvar ,aldt ,nel ,nuvar ,lbuf%ANG ,
475 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
476 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
477 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
478 ENDDO
479 ELSE IF (ilaw == 158) THEN
480 DO it=1,nptt
481 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
482 uvar => elbuf_str%BUFLY(il)%MAT(ir,is,it)%VAR
483 dir1 => elbuf_str%BUFLY(il)%LBUF_DIR(it)%DIRA
484 dir2 => elbuf_str%BUFLY(il)%LBUF_DIR(it)%DIRB
485 CALL law158_init(dir1 ,dir2 ,
486 . uvar ,aldt ,nel ,nuvar ,lbuf%ANG ,
487 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
488 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
489 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
490 ENDDO
491 ENDIF
492!
493 ENDDO ! DO IL = 1,NLAY
494 ENDIF
495 ENDIF ! IF (MTN == 58) THEN
496 IF (mtn == 42 .OR. mtn == 69 .OR. igtyp == 51 .OR. igtyp == 52) THEN
497 DO il = 1,nlay
498 nptt = elbuf_str%BUFLY(il)%NPTT
499 ilaw = elbuf_str%BUFLY(il)%ILAW
500 nuvar = elbuf_str%BUFLY(il)%NVAR_MAT
501!
502 IF (ilaw == 42 .OR. mtn == 69) THEN
503 DO it=1,nptt
504 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
505 uvar => elbuf_str%BUFLY(il)%MAT(ir,is,it)%VAR
506 CALL law42c_ini(nuvar,uvar,llt)
507 ENDDO
508 ENDIF ! IF (ILAW == 42) THEN
509 ENDDO ! DO IL = 1,NLAY
510 ENDIF
511C-----------------------------------------------------------------------
512C CALCULATION OF INITIAL STRESSES
513C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
514 IF (isigsh /= 0 .OR. ithkshel == 2) THEN
515C to be checked for IGTYP = 51
516 IF (mpt>0) THEN
517 CALL layini1(
518 . elbuf_str ,lft ,llt ,geo ,igeo ,
519 . mat ,pid ,matly ,posly ,igtyp ,
520 . nlay ,mpt ,isubstack ,stack ,drape ,
521 . nft ,gbuf%THK ,nel ,idrape , stdrape ,
522 . indx)
523 CALL corth3(elbuf_str,dir_a ,dir_b ,lft ,llt ,
524 . nlay ,irep ,nel ,
525 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
526 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
527 . e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z ,
528 . idrape , igtyp)
529 END IF
530 CALL csigini(elbuf_str,
531 1 lft ,llt ,nft ,mpt ,istrain ,
532 2 gbuf%THK,gbuf%EINT,gbuf%STRA,gbuf%HOURG,gbuf%PLA ,
533 3 gbuf%FOR,gbuf%MOM ,sigsh ,nlay ,gbuf%G_HOURG,
534 4 numeltg ,ixtg ,nixtg ,nsigsh ,numsh3n ,
535 5 ptsh3n ,igeo ,thk ,nel ,e1x ,
536 6 e2x ,e3x ,e1y ,e2y ,e3y ,
537 7 e1z ,e2z ,e3z ,isigsh ,dir_a ,
538 8 dir_b ,posly ,igtyp )
539 IF (mpt == 0) THEN
540 DO i=lft,llt
541 gbuf%FOR_G(i+jj(1:5))=gbuf%FOR(i+jj(1:5))
542 END DO
543 END IF
544 ELSEIF ( ithkshel == 1 ) THEN
545 CALL thickini(lft ,llt ,nft ,ptsh3n,numeltg,
546 2 gbuf%THK,thk ,ixtg ,nixtg ,nsigsh ,
547 3 sigsh )
548 ENDIF
549c-----------------------------------------------------------------------
550c Failure model initialisation
551c-----------------------------------------------------------------------
552c set element damage initialized in /FAIL/FRACTAL_DMG
553 IF (fail_fractal%NFAIL > 0) THEN
554 CALL fractal_dmg_init(elbuf_str,mat_param,fail_fractal,
555 . nummat ,numeltg ,nel ,nft ,ngl ,ity )
556 ENDIF
557c tag edge elements in local UVAR for /FAIL/ALTER (XFEM)
558 IF (ifail > 0 .and. iddlevel == 1) THEN
559 CALL fail_windshield_init(elbuf_str,mat_param,fail_brokmann,
560 . nel ,nft ,ity ,igrsh4n ,igrsh3n ,
561 . aldt ,thk ,ngl )
562 ENDIF
563c
564C to be checked for IGTYP = 51
565c IF(NPERTURB /= 0 .OR. NVSHELL1 /= 0) THEN
566 CALL cfailini(elbuf_str,mat_param,
567 . nptt ,nlay ,sigsh ,nsigsh ,ptsh3n ,
568 . rnoise ,perturb ,aldt ,thk )
569c ENDIF
570C-----------------------------------------------------------------------
571C CALCUL DES DEFORMATIONS INITIALES (MEMBRANE)
572C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
573 IF (istrain == 1 .AND. nxref > 0) THEN
574 uvar => elbuf_str%BUFLY(1)%MAT(1,1,1)%VAR
575 imat = elbuf_str%BUFLY(1)%IMAT
576C to be checked for IGTYP = 51
577 CALL c3epsini(elbuf_str,mat_param(imat),
578 . lft ,llt ,ismstr ,mtn ,ithk ,
579 . pm ,geo ,ixtg(1,nft+1),x ,xreftg(1,1,nft+1),
580 . gbuf%FOR,gbuf%THK ,gbuf%EINT ,gbuf%STRA,nlay ,
581 . px1g ,py1g ,py2g ,x2s ,x3s ,
582 . y3s ,uvar ,ipm ,
583 . nel ,dir_a ,dir_b ,gbuf%SIGI,npf ,
584 . tf ,irep ,ifram_old ,imat )
585c
586 CALL c3epschk(lft, llt,nft, pm, geo,ixtg(1,nft+1), gbuf%STRA,thk,
587 . nel,cpt_eltens)
588c
589 IF (ismstr == 1 .AND. mtn==19) iparg(9) = 11
590c
591 ELSEIF (ismstr == 11 .OR. (ismstr==1 .AND. mtn==19)) THEN
592C to be checked for IGTYP = 51
593 CALL c3eps11_ini(
594 . lft ,llt ,ixtg(1,nft+1),x ,x2s ,
595 . x3s ,y3s )
596 ENDIF
597c
598 IF (ismstr == 10) THEN
599 DO i=lft,llt
600 ii = nft + i
601 elbuf_str%GBUF%SMSTR(jj(1)+i) = x(1,ixtg(3,ii))-x(1,ixtg(2,ii))
602 elbuf_str%GBUF%SMSTR(jj(2)+i) = x(2,ixtg(3,ii))-x(2,ixtg(2,ii))
603 elbuf_str%GBUF%SMSTR(jj(3)+i) = x(3,ixtg(3,ii))-x(3,ixtg(2,ii))
604 elbuf_str%GBUF%SMSTR(jj(4)+i) = x(1,ixtg(4,ii))-x(1,ixtg(2,ii))
605 elbuf_str%GBUF%SMSTR(jj(5)+i) = x(2,ixtg(4,ii))-x(2,ixtg(2,ii))
606 elbuf_str%GBUF%SMSTR(jj(6)+i) = x(3,ixtg(4,ii))-x(3,ixtg(2,ii))
607 ENDDO
608 ELSEIF (ismstr == 11 .OR.(ismstr==1 .AND. mtn==19)) THEN
609 DO i=lft,llt
610 elbuf_str%GBUF%SMSTR(jj(1)+i) = x2s(i)
611 elbuf_str%GBUF%SMSTR(jj(2)+i) = x3s(i)
612 elbuf_str%GBUF%SMSTR(jj(3)+i) = y3s(i)
613 ENDDO
614 ENDIF
615C
616 IF (iuser == 1 .and. mtn > 28) THEN
617C to be checked for IGTYP = 51
618 CALL cuserini(elbuf_str,
619 1 lft ,llt ,nft ,nel ,npt ,
620 2 istrain,sigsh ,numeltg ,ixtg ,nixtg ,
621 3 nsigsh ,numsh3n,ptsh3n ,ir ,is ,
622 4 nlay )
623 ENDIF
624C-------------------------------------------
625 IF (iyldini == 1 .AND. (mtn== 36.OR. mtn==87))THEN
626 CALL cmatini(elbuf_str,
627 1 lft ,llt ,nft ,nel ,npt ,
628 2 istrain,sigsh ,numeltg ,ixtg ,nixtg ,
629 3 nsigsh ,numsh3n,ptsh3n ,ir ,is ,
630 4 nlay )
631 ENDIF
632C-------------------------------------------
633C CALCULATION OF ELEMENTARY TIMESTEPS
634C-------------------------------------------
635c IGTYP=GEO(12,IXTG(5,I+NFT))
636 IF (igtyp /= 0 .AND. igtyp /= 1 .AND.
637 . igtyp /= 9 .AND. igtyp /= 10 .AND.
638 . igtyp /= 11 .AND. igtyp /= 16 .AND.
639 . igtyp /= 17 .AND. igtyp /= 51 .AND.
640 . igtyp /= 52 ) THEN
641 CALL ancmsg(msgid=25,
642 . anmode=aninfo,
643 . msgtype=msgerror,
644 . i1=id,
645 . c1=titr,
646 . i2=iprop)
647 ENDIF
648 ndepar=numels+numelc+numelt+numelp+numelr+nft
649 DO i=lft,llt
650 dtelem(ndepar+i) = dt(i)
651 END DO
652C---
653 IF (ixfem > 0) THEN
654 CALL cbufxfe(elbuf_str,xfem_str,isubstack,stack ,
655 . igeo ,geo ,lft ,llt ,mat,
656 . pid ,npt ,nptt ,nlay,ir ,
657 . is ,ixfem,mtn ,ng)
658 ENDIF
659C------------
660 ! Compute the initial volume
661 DO i=lft,llt
662 IF (gbuf%G_VOL > 0) gbuf%VOL(i) = area(i)*gbuf%THK(i)
663 ENDDO
664 IF (ixfem > 0) THEN
665 DO ixel=1,nxel
666 DO i=lft,llt
667 IF (xfem_str(ng,ixel)%GBUF%G_VOL > 0)
668 . xfem_str(ng,ixel)%GBUF%VOL(i) = area(i)*gbuf%THK(i)
669 END DO
670 ENDDO
671 ENDIF
672C---
673 IF (ALLOCATED(dir_b)) DEALLOCATE(dir_b)
674 IF (ALLOCATED(dir_a)) DEALLOCATE(dir_a)
675 IF (ALLOCATED(indx)) DEALLOCATE(indx)
676
677 DEALLOCATE(matly)
678 DEALLOCATE(posly)
679C------------
680 RETURN
681 END
682C
683!||====================================================================
684!|| c3epschk ../starter/source/elements/sh3n/coque3n/c3init3.F
685!||--- called by ------------------------------------------------------
686!|| c3init3 ../starter/source/elements/sh3n/coque3n/c3init3.F
687!|| cdkinit3 ../starter/source/elements/sh3n/coquedk/cdkinit3.F
688!||--- calls -----------------------------------------------------
689!|| ancmsg ../starter/source/output/message/message.F
690!||--- uses -----------------------------------------------------
691!|| message_mod ../starter/share/message_module/message_mod.F
692!||====================================================================
693 SUBROUTINE c3epschk(JFT,JLT,NFT,PM,GEO,IXTG,GSTR,THK,
694 . NEL,CPT_ELTENS)
695 USE message_mod
696 use element_mod , only : nixtg
697C-----------------------------------------------
698C I m p l i c i t T y p e s
699C-----------------------------------------------
700#include "implicit_f.inc"
701C-----------------------------------------------
702C G l o b a l P a r a m e t e r s
703C-----------------------------------------------
704#include "mvsiz_p.inc"
705C-----------------------------------------------
706C C o m m o n B l o c k s
707C-----------------------------------------------
708#include "param_c.inc"
709#include "scr03_c.inc"
710C-----------------------------------------------
711C D u m m y A r g u m e n t s
712C-----------------------------------------------
713 INTEGER JFT, JLT,NFT ,NEL ,IXTG(NIXTG,*),CPT_ELTENS
714 my_real
715 . PM(NPROPM,*), GSTR(NEL,8),THK(*), GEO(NPROPG,*)
716C------------------------------------------------------
717C L o c a l V a r i a b l e s
718C------------------------------------------------------
719 INTEGER I
720 my_real
721 . DELT(MVSIZ),X1(MVSIZ),X2(MVSIZ)
722C-------------------------------------------------------------
723 DO I=jft,jlt
724 delt(i)=((gstr(i,1)+gstr(i,2))*(gstr(i,1)+gstr(i,2)))-four*
725 . (gstr(i,1)*gstr(i,2)-fourth*gstr(i,3)*gstr(i,3))
726 IF(delt(i)>=zero)THEN
727 x1(i) = (gstr(i,1)+gstr(i,2)-sqrt(delt(i)))/two
728 x2(i) = (gstr(i,1)+gstr(i,2)+sqrt(delt(i)))/two
729 IF((x1(i)>zero.AND.(x1(i)>em10)).OR.(x2(i)>zero
730 . .AND.(x2(i)>em10)))THEN
731 IF(ipri == 5) THEN
732 CALL ancmsg(msgid=607,
733 . msgtype=msgwarning,
734 . anmode=aninfo_blind_1,
735 . r1=min(x1(i),x2(i)),
736 . r2=max(x1(i),x2(i)),
737 . i1=ixtg(6,i))
738 ELSE
739 cpt_eltens = cpt_eltens + 1
740 ENDIF
741 ENDIF
742 ENDIF
743 ENDDO
744c-----------
745 RETURN
746 END
747!||====================================================================
748!|| c3eps11_ini ../starter/source/elements/sh3n/coque3n/c3init3.f
749!||--- called by ------------------------------------------------------
750!|| c3init3 ../starter/source/elements/sh3n/coque3n/c3init3.F
751!||--- uses -----------------------------------------------------
752!||====================================================================
753 SUBROUTINE c3eps11_ini(
754 1 JFT ,JLT ,IXTG ,X ,X2S ,
755 2 X3S ,Y3S )
756 use element_mod , only : nixtg
757C-----------------------------------------------
758C I m p l i c i t T y p e s
759C-----------------------------------------------
760#include "implicit_f.inc"
761C-----------------------------------------------
762C G l o b a l P a r a m e t e r s
763C-----------------------------------------------
764#include "mvsiz_p.inc"
765C-----------------------------------------------
766C D u m m y A r g u m e n t s
767C-----------------------------------------------
768 INTEGER JFT, JLT
769 INTEGER IXTG(NIXTG,*)
770 my_real x(3,*),x2s(*), x3s(*), y3s(*)
771C------------------------------------------------------
772C L o c a l V a r i a b l e s
773C-----------------------------------------------
774 INTEGER I
775 my_real
776 . X1(MVSIZ) , X2(MVSIZ) , X3(MVSIZ) ,
777 . Y1(MVSIZ) , Y2(MVSIZ) , Y3(MVSIZ) ,
778 . Z1(MVSIZ) , Z2(MVSIZ) , Z3(MVSIZ)
779 my_real SUM
780 my_real
781 . E1X(MVSIZ), E1Y(MVSIZ), E1Z(MVSIZ),
782 . E2X(MVSIZ), E2Y(MVSIZ), E2Z(MVSIZ),
783 . E3X(MVSIZ), E3Y(MVSIZ), E3Z(MVSIZ),
784 . X31(MVSIZ), Y31(MVSIZ), Z31(MVSIZ),
785 . x32(mvsiz), y32(mvsiz), z32(mvsiz),
786 . x21(mvsiz), y21(mvsiz), z21(mvsiz)
787C---------------------------------------------------------------
788 DO i=jft,jlt
789 x1(i)=zero
790 y1(i)=zero
791 z1(i)=zero
792 x2(i)=x(1,ixtg(3,i))-x(1,ixtg(2,i))
793 y2(i)=x(2,ixtg(3,i))-x(2,ixtg(2,i))
794 z2(i)=x(3,ixtg(3,i))-x(3,ixtg(2,i))
795 x3(i)=x(1,ixtg(4,i))-x(1,ixtg(2,i))
796 y3(i)=x(2,ixtg(4,i))-x(2,ixtg(2,i))
797 z3(i)=x(3,ixtg(4,i))-x(3,ixtg(2,i))
798 ENDDO
799C
800 DO i=jft,jlt
801 x21(i)=x2(i)
802 y21(i)=y2(i)
803 z21(i)=z2(i)
804 x31(i)=x3(i)
805 y31(i)=y3(i)
806 z31(i)=z3(i)
807 x32(i)=x3(i)-x2(i)
808 y32(i)=y3(i)-y2(i)
809 z32(i)=z3(i)-z2(i)
810 ENDDO
811 DO i=jft,jlt
812 e1x(i)= x21(i)
813 e1y(i)= y21(i)
814 e1z(i)= z21(i)
815 sum = one/sqrt(e1x(i)*e1x(i)+e1y(i)*e1y(i)+e1z(i)*e1z(i))
816 e1x(i)=e1x(i)* sum
817 e1y(i)=e1y(i)* sum
818 e1z(i)=e1z(i)* sum
819 ENDDO
820C
821 DO i=jft,jlt
822 e3x(i)=y31(i)*z32(i)-z31(i)*y32(i)
823 e3y(i)=z31(i)*x32(i)-x31(i)*z32(i)
824 e3z(i)=x31(i)*y32(i)-y31(i)*x32(i)
825 sum = one/sqrt(e3x(i)*e3x(i)+e3y(i)*e3y(i)+e3z(i)*e3z(i))
826 e3x(i)=e3x(i)* sum
827 e3y(i)=e3y(i)* sum
828 e3z(i)=e3z(i)* sum
829 ENDDO
830C
831 DO i=jft,jlt
832 e2x(i) = e3y(i) * e1z(i) - e3z(i) * e1y(i)
833 e2y(i) = e3z(i) * e1x(i) - e3x(i) * e1z(i)
834 e2z(i) = e3x(i) * e1y(i) - e3y(i) * e1x(i)
835 sum = one/sqrt(e2x(i)*e2x(i)+e2y(i)*e2y(i)+e2z(i)*e2z(i))
836 e2x(i)=e2x(i)*sum
837 e2y(i)=e2y(i)*sum
838 e2z(i)=e2z(i)*sum
839 ENDDO
840C
841 DO i=jft,jlt
842 y3s(i)=e2x(i)*x31(i)+e2y(i)*y31(i)+e2z(i)*z31(i)
843 x3s(i)=e1x(i)*x31(i)+e1y(i)*y31(i)+e1z(i)*z31(i)
844 x2s(i)=e1x(i)*x21(i)+e1y(i)*y21(i)+e1z(i)*z21(i)
845 ENDDO
846C-----------
847 RETURN
848 END
c1buf3
subroutine c1buf3(geo, thk, off, thke, kshtree, shtree)
Definition c1buf3.F:32
c3coori
subroutine c3coori(x, xreftg, ixp, ngl, x1, x2, x3, y1, y2, y3, z1, z2, z3, ix1, ix2, ix3)
Definition c3coori.F:39
c3derii
subroutine c3derii(jft, jlt, pm, geo, px1, py1, py2, stifn, stifr, ixtg, thk, sh3tree, aldt, uparam, ipm, igeo, pm_stack, isubstack, strtg, imat, iprop, area, dt, x31g, y31g, z31g, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, x2, x3, y3, group_param)
Definition c3derii.F:38
c3epsini
subroutine c3epsini(elbuf_str, mat_param, jft, jlt, ismstr, mtn, ithk, pm, geo, ixtg, x, xreftg, for, thk, eint, gstr, nlay, px1g, py1g, py2g, x2s, x3s, y3s, uvar, ipm, nel, dir_a, dir_b, sigi, npf, tf, irep, ish3nfr, imat)
Definition c3epsini.F:45
c3epschk
subroutine c3epschk(jft, jlt, nft, pm, geo, ixtg, gstr, thk, nel, cpt_eltens)
Definition c3init3.F:695
c3init3
subroutine c3init3(elbuf_str, ixtg, pm, x, geo, igrsh4n, xmas, in, nvc, dtelem, igrsh3n, xreftg, offset, nel, ithk, thk, isigsh, sigsh, stifn, stifr, partsav, v, ipart, mstg, intg, ptg, skew, iparg, nsigsh, igeo, iuser, etnod, nshnod, sttg, ptsh3n, ipm, bufmat, sh3tree, mcp, mcptg, temp, cpt_eltens, part_area, itage, itagn, ixfem, npf, tf, sh3trim, xfem_str, isubstack, stack, rnoise, drape, sh3ang, iddlevel, geo_stack, igeo_stack, strtg, perturb, ish3n, iyldini, ele_area, nloc_dmg, ng, group_param, idrape, drapeg, mat_param, fail_fractal, fail_brokmann, glob_therm)
Definition c3init3.F:81
c3eps11_ini
subroutine c3eps11_ini(jft, jlt, ixtg, x, x2s, x3s, y3s)
Definition c3init3.F:756
c3inmas
subroutine c3inmas(x, xreftg, ixtg, geo, pm, ms, tiner, thke, partsav, v, ipart, mstg, intg, ptg, igeo, imat, iprop, area, etnod, nshnod, sttg, sh3tree, mcp, mcptg, temp, sh3trim, isubstack, nlay, elbuf_str, stack, thki, rnoise, drape, perturb, ix1, ix2, ix3, nintemp, x2, x3, y3, idrape, indx)
Definition c3inmas.F:46
c3veok3
subroutine c3veok3(nvc, ix1, ix2, ix3)
Definition c3veok3.F:36
cbufxfe
subroutine cbufxfe(elbuf_str, xfem_str, isubstack, stack, igeo, geo, lft, llt, mat, pid, npt, nptt, nlay, ir, is, ixfem, mtn, ng)
Definition cbufxfe.F:39
cdkevec3
subroutine cdkevec3(jft, jlt, area, x1, x2, x3, y1, y2, y3, z1, z2, z3, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition cdkevec3.F:38
cfailini
subroutine cfailini(elbuf_str, mat_param, nptt, nlay, sigsh, nsigsh, ptsh, rnoise, perturb, aldt, thk)
Definition cfailini.F:36
cm27in3
subroutine cm27in3(elbuf_str, geo, igeo, pm, ipm, ix, nix, nlay, ir, is, imat)
Definition cm27in3.F:38
cm35in3
subroutine cm35in3(elbuf_str, thk, area, nel, nlay, nptr, npts, nptt, igtyp)
Definition cm35in3.F:34
cm58in3
subroutine cm58in3(irep, dir1, dir2, uparam, uvar, aldt, nel, nuvar, tan_phi, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition cm58in3.F:35
cmatini
subroutine cmatini(elbuf_str, jft, jlt, nft, nel, npt, istrain, sigsh, numel, ix, nix, nsigsh, numsh, ptsh, ir, is, nlay)
Definition cmatini.F:36
corth3
subroutine corth3(elbuf_str, dir_a, dir_b, jft, jlt, nlay, irep, nel, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, idrape, igtyp)
Definition corth3.F:42
corthdir
subroutine corthdir(elbuf_str, igeo, geo, vx, vy, vz, phi1, phi2, coor1, coor2, coor3, coor4, iorthloc, nlay, irep, isubstack, stack, geo_stack, igeo_stack, ir, is, nel, imat, iprop, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, npt_all, idrape)
Definition corthdir.F:42
corthini
subroutine corthini(jft, jlt, nft, nlay, numel, nsigsh, nix, ix, igeo, geo, skew, sigsh, ptsh, phi1, phi2, vx, vy, vz, coor1, coor2, coor3, coor4, iorthloc, isubstack, stack, irep, elbuf_str, drape, angle, x, geo_stack, e3x, e3y, e3z, betaorth, x1, x2, y1, y2, z1, z2, nel, g_add_node, add_node, npt_all, idrape, indx)
Definition corthini.F:48
my_real
#define my_real
Definition cppsort.cpp:32
csigini
subroutine csigini(elbuf_str, jft, jlt, nft, npt, istrain, thk, eint, gstr, hh, plas, for, mom, sigsh, nlay, g_hourg, numel, ix, nix, nsigsh, numsh, ptsh, igeo, thke, nel, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, isigsh, dir_a, dir_b, posly, igtyp)
Definition csigini.F:46
cuserini
subroutine cuserini(elbuf_str, jft, jlt, nft, nel, npt, istrain, sigsh, numel, ix, nix, nsigsh, numsh, ptsh, ir, is, nlay)
Definition cuserini.F:36
fail_windshield_init
subroutine fail_windshield_init(elbuf_str, mat_param, fail_brokmann, nel, nft, ity, igrsh4n, igrsh3n, aldt, thk, ngl)
Definition fail_windshield_init.F:40
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
law158_init
subroutine law158_init(dir1, dir2, uvar, aldt, nel, nuvar, tan_phi, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition law158_init.F:35
layini1
subroutine layini1(elbuf_str, jft, jlt, geo, igeo, mat, pid, matly, posly, igtyp, nlay, npt, isubstack, stack, drape, nft, thk, nel, idrape, numel_drape, indx)
Definition layini1.F:40
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
drape_mod
Definition drape_mod.F:54
drape_mod::stdrape
integer stdrape
Definition drape_mod.F:92
groupdef_mod
Definition groupdef_mod.F:668
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008
nlocal_reg_mod
Definition nlocal_reg_mod.F:119
stack_mod
Definition stack_mod.F:75
c3evec3
subroutine c3evec3(jft, jlt, area, x1, x2, x3, y1, y2, y3, z1, z2, z3, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, x31, y31, z31, x2l, x3l, y3l)
Definition c3evec3.F:39
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895
fretitl2
subroutine fretitl2(titr, iasc, l)
Definition freform.F:799
starter
program starter
Definition starter.F:39
thickini
subroutine thickini(jft, jlt, nft, ptsh, numel, thk, thke, ix, nix, nsigsh, sigsh)
Definition thickini.F:36