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