OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
h3d_shell_vector.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!|| h3d_shell_vector ../engine/source/output/h3d/h3d_results/h3d_shell_vector.F
25!||--- called by ------------------------------------------------------
26!|| genh3d ../engine/source/output/h3d/h3d_results/genh3d.F
27!||--- calls -----------------------------------------------------
28!|| h3d_shell_vector_1 ../engine/source/output/h3d/h3d_results/h3d_shell_vector_1.F
29!||--- uses -----------------------------------------------------
30!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
31!|| initbuf_mod ../engine/share/resol/initbuf.F
32!|| multi_fvm_mod ../common_source/modules/ale/multi_fvm_mod.F90
33!|| names_and_titles_mod ../common_source/modules/names_and_titles_mod.F
34!|| schlieren_mod ../engine/share/modules/schlieren_mod.F
35!|| stack_mod ../engine/share/modules/stack_mod.F
36!||====================================================================
37 SUBROUTINE h3d_shell_vector(
38 . ELBUF_TAB ,SHELL_VECTOR,IFUNC ,IPARG ,GEO ,
39 . IXQ ,IXC ,IXTG ,PM ,
40 . EL2FA ,NBF ,IADP ,
41 . NBF_L ,EHOUR ,ANIM ,NBPART ,IADG ,
42 . IPM ,IGEO ,THKE ,ERR_THK_SH4 ,ERR_THK_SH3,
43 . INVERT ,X ,V ,W ,
44 . NV46 ,NERCVOIS ,NESDVOIS ,LERCVOIS ,LESDVOIS,
45 . STACK ,ID_ELEM ,ITY_ELEM ,INFO1 ,INFO2 ,
46 . IS_WRITTEN_SHELL,IPARTC,IPARTTG ,LAYER_INPUT ,IPT_INPUT ,
47 . PLY_INPUT ,GAUSS_INPUT,IUVAR_INPUT,H3D_PART ,KEYWORD ,
48 . D ,MULTI_FVM)
49C-----------------------------------------------
50C M o d u l e s
51C-----------------------------------------------
52 USE initbuf_mod
53 USE elbufdef_mod
55 USE stack_mod
56 USE multi_fvm_mod
58C-----------------------------------------------
59C I m p l i c i t T y p e s
60C-----------------------------------------------
61#include "implicit_f.inc"
62C-----------------------------------------------
63C C o m m o n B l o c k s
64C-----------------------------------------------
65#include "mvsiz_p.inc"
66#include "com01_c.inc"
67#include "com04_c.inc"
68#include "param_c.inc"
69C-----------------------------------------------
70C D u m m y A r g u m e n t s
71C-----------------------------------------------
73 . shell_vector(3,*),x(3,numnod),v(3,numnod),w(3,numnod),d(3,numnod),thke(*),ehour(*),geo(npropg,numgeo),
74 . anim(*),pm(npropm,nummat),err_thk_sh4(*), err_thk_sh3(*)
75 INTEGER IPARG(NPARG,NGROUP),IXC(NIXC,NUMELC),IXTG(NIXTG,NUMELTG),EL2FA(*),
76 . IXQ(NIXQ,NUMELQ),IFUNC,NBF,
77 . IADP(*),NBF_L, NBPART,IADG(NSPMD,*),IPM(NPROPMI,NUMMAT),
78 . IGEO(NPROPGI,NUMGEO),INVERT(*), NV46,ID_ELEM(*),ITY_ELEM(*),
79 . INFO1,INFO2,IS_WRITTEN_SHELL(*),IPARTC(*),IPARTTG(*),H3D_PART(*),
80 . LAYER_INPUT ,IPT_INPUT,GAUSS_INPUT,PLY_INPUT,IUVAR_INPUT
81 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
82 TYPE (STACK_PLY) :: STACK
83 CHARACTER(LEN=NCHARLINE100)::KEYWORD
84 TYPE (MULTI_FVM_STRUCT), INTENT(IN) :: MULTI_FVM
85C-----------------------------------------------
86C L o c a l V a r i a b l e s
87C-----------------------------------------------
89 . evar(mvsiz),dam1(mvsiz),dam2(mvsiz),
90 . wpla(mvsiz),dmax(mvsiz),wpmax(mvsiz),fail(mvsiz),
91 . epst1(mvsiz),epst2(mvsiz),epsf1(mvsiz),epsf2(mvsiz),
92 . value(mvsiz),vg(5),vly(5),ve(5),mass(mvsiz)
94 . off, p,vonm2,s1,s2,s12,s3,dmgmx,fac,
95 . dir1_1,dir1_2,dir2_1,dir2_2,aa,bb,v1,v2,v3,x21,x32,x34,
96 . x41,y21,y32,y34,y41,z21,z32,z34,z41,suma,vr,vs,x31,y31,
97 . z31,e11,e12,e13,e21,e22,e23,sum,area,x2l,var,
98 . e1x,e1y,e1z,e2x,e2y,e2z,e3x,e3y,e3z,rx,ry,rz,sx,sy,sz,
99 . rho0,thk0,xx1,xx2,xx3,yy1,yy2,yy3,zz1,zz2,zz3,a0
100 INTEGER I,I1,II,J,NG,NEL,NPTR,NPTS,NPTT,NLAY,L,IFAIL,ILAY,
101 . IR,IS,IT,IL,MLW, NUVAR,IUS,LENF,PTF,PTM,PTS,NFAIL,
102 . N,NN,K,K1,K2,JTURB,MT,IMID,IALEL,IPID,ISH3N,NNI,
103 . NN1,NN2,NN3,NN4,NN5,NN6,NN9,NF,BUF,NVARF,
104 . offset,ihbe,nptm,npg, mpt,ipt,iadd,iadr,ipmat,ifailt,
105 . iigeo,iadi,isubstack,ithk,nercvois(*),nesdvois(*),
106 . lercvois(*),lesdvois(*),id_ply,nb_plyoff,iok,n1,n2,n3,n4
107 INTEGER PID(MVSIZ),MAT(MVSIZ),MATLY(MVSIZ*100),FAILG(100,MVSIZ),
108 . PTE(4),PTP(4),PTMAT(4),PTVAR(4),NPT_ALL,IPLY,
109 . ID_ELEM_TMP(MVSIZ),NIX,IOK_PART(MVSIZ),JJ(5),NPGT,IUVAR,
110 . IS_WRITTEN_VALUE(MVSIZ)
111 REAL R4
112 TYPE(G_BUFEL_) ,POINTER :: GBUF
113 TYPE(L_BUFEL_) ,POINTER :: LBUF
114 TYPE(BUF_LAY_) ,POINTER :: BUFLY
115 TYPE(BUF_FAIL_) ,POINTER :: FBUF
116 my_real,
117 . DIMENSION(:), POINTER :: UVAR
118 TYPE(l_bufel_) ,POINTER :: LBUF1,LBUF2,LBUF3,LBUF4
119C-----------------------------------------------
120c print *,'LAYER_INPUT',LAYER_INPUT
121c print *,'IPT_INPUT',IPT_INPUT
122c print *,'GAUSS_INPUT',GAUSS_INPUT
123c print *,'PLY_INPUT',PLY_INPUT
124c print *,'IUVAR_INPUT',IUVAR_INPUT
125c print *,'KEYWORD',KEYWORD,'KEYWORD'
126C
127 DO i=1,numelc+numeltg
128 is_written_shell(i) = 0
129 ENDDO
130c a corriger
131 nn3 = 0
132C
133 DO ng=1,ngroup
135 . elbuf_tab ,shell_vector,ifunc ,iparg ,geo ,
136 . ixq ,ixc ,ixtg ,pm ,
137 . el2fa ,nbf ,iadp ,
138 . nbf_l ,ehour ,anim ,nbpart ,iadg ,
139 . ipm ,igeo ,thke ,err_thk_sh4 ,err_thk_sh3,
140 . invert ,x ,v ,w ,
141 . nv46 ,nercvois ,nesdvois ,lercvois ,lesdvois,
142 . stack ,id_elem ,ity_elem ,info1 ,info2 ,
143 . is_written_shell,ipartc,iparttg ,layer_input ,ipt_input ,
144 . ply_input ,gauss_input,iuvar_input,h3d_part ,keyword ,
145 . d ,ng ,multi_fvm)
146 ENDDO
147C-----------------------------------------------
148 RETURN
149 END
150
#define my_real
Definition cppsort.cpp:32
subroutine h3d_shell_vector(elbuf_tab, shell_vector, ifunc, iparg, geo, ixq, ixc, ixtg, pm, el2fa, nbf, iadp, nbf_l, ehour, anim, nbpart, iadg, ipm, igeo, thke, err_thk_sh4, err_thk_sh3, invert, x, v, w, nv46, nercvois, nesdvois, lercvois, lesdvois, stack, id_elem, ity_elem, info1, info2, is_written_shell, ipartc, iparttg, layer_input, ipt_input, ply_input, gauss_input, iuvar_input, h3d_part, keyword, d, multi_fvm)
subroutine h3d_shell_vector_1(elbuf_tab, shell_vector, ifunc, iparg, geo, ixq, ixc, ixtg, pm, el2fa, nbf, iadp, nbf_l, ehour, anim, nbpart, iadg, ipm, igeo, thke, err_thk_sh4, err_thk_sh3, invert, x, v, w, nv46, nercvois, nesdvois, lercvois, lesdvois, stack, id_elem, ity_elem, info1, info2, is_written_shell, ipartc, iparttg, layer_input, ipt_input, ply_input, gauss_input, iuvar_input, h3d_part, keyword, d, ng, multi_fvm)
subroutine area(d1, x, x2, y, y2, eint, stif0)
integer, parameter ncharline100