OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
lec_inistate.F File Reference
#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"
#include "scry_c.inc"
#include "scr16_c.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine lec_inistate (ixs, ixq, ixc, ixt, ixp, ixr, geo, pm, kxsp, ixtg, index, itri, nsigsh, igeo, ipm, nsigs, nsigsph, ksysusr, ptshel, ptsh3n, ptsol, ptquad, ptsph, numel, nsigrs, unitab, isolnodd00, lsubmodel, rtrans, idrape, nsigi, ptspri, nsigbeam, ptbeam, nsigtruss, pttruss, sigi, sigsh, sigsp, sigsph, sigrs, sigbeam, sigtruss, strsglob, straglob, orthoglob, isigsh, iyldini, ksigsh3, fail_ini, iusolyld, iuser, igrbric, map_tables, iparg, stack, iworksh, mat_param, numsph, nisp)

Function/Subroutine Documentation

◆ lec_inistate()

subroutine lec_inistate ( integer, dimension(nixs,*) ixs,
integer, dimension(nixq,*) ixq,
integer, dimension(nixc,*) ixc,
integer, dimension(nixt,*) ixt,
integer, dimension(nixp,*) ixp,
integer, dimension(nixr,*) ixr,
geo,
pm,
integer, dimension(*) kxsp,
integer, dimension(nixtg,*) ixtg,
integer, dimension(*) index,
integer, dimension(*) itri,
integer nsigsh,
integer, dimension(npropgi,*) igeo,
integer, dimension(npropmi,*) ipm,
integer nsigs,
integer nsigsph,
integer, dimension(*) ksysusr,
integer, dimension(*) ptshel,
integer, dimension(*) ptsh3n,
integer, dimension(*) ptsol,
integer, dimension(*) ptquad,
integer, dimension(*) ptsph,
integer numel,
integer nsigrs,
type (unit_type_), intent(in) unitab,
integer, dimension(*) isolnodd00,
type(submodel_data), dimension(*) lsubmodel,
rtrans,
integer, dimension(nplymax,*) idrape,
integer nsigi,
integer, dimension(*) ptspri,
integer nsigbeam,
integer, dimension(*) ptbeam,
integer nsigtruss,
integer, dimension(*) pttruss,
sigi,
sigsh,
sigsp,
sigsph,
sigrs,
sigbeam,
sigtruss,
integer, dimension(*) strsglob,
integer, dimension(*) straglob,
integer, dimension(*) orthoglob,
integer isigsh,
integer iyldini,
integer ksigsh3,
integer, dimension(5) fail_ini,
integer iusolyld,
integer iuser,
type (group_), dimension(ngrbric) igrbric,
type(mapping_struct_) map_tables,
integer, dimension(nparg,ngroup), intent(in) iparg,
type (stack_ply) stack,
integer, dimension(3,numelc + numeltg), intent(in) iworksh,
type (matparam_struct_), dimension(nummat), intent(inout) mat_param,
integer, intent(in) numsph,
integer, intent(in) nisp )

Definition at line 36 of file lec_inistate.F.

49C-----------------------------------------------
50C M o d u l e s
51C-----------------------------------------------
52 USE unitab_mod
53 USE groupdef_mod
54 USE submodel_mod
55 USE mapping_option_mod
56 USE stack_mod
57 USE matparam_def_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 "com01_c.inc"
66#include "com04_c.inc"
67#include "param_c.inc"
68#include "scry_c.inc"
69#include "scr16_c.inc"
70#include "units_c.inc"
71C-----------------------------------------------
72C D u m m y A r g u m e n t s
73C-----------------------------------------------
74 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
75 INTEGER IXS(NIXS,*), IXQ(NIXQ,*),IXC(NIXC,*),
76 . IGEO(NPROPGI,*), IXT(NIXT,*),IXP(NIXP,*), IXR(NIXR,*),
77 . IXTG(NIXTG,*),INDEX(*),ITRI(*),KXSP(*),IPM(NPROPMI,*),
78 . KSYSUSR(*),PTSHEL(*),PTSH3N(*),PTSOL(*),PTQUAD(*),PTSPH(*),
79 . IDRAPE(NPLYMAX,*),PTSPRI(*),PTBEAM(*),PTTRUSS(*)
80 INTEGER NSIGI,NSIGSH,NSIGS, NSIGSPH, NSIGRS,
81 . NUMEL,ISOLNODD00(*),NSIGBEAM,NSIGTRUSS,STRSGLOB(*),
82 . STRAGLOB(*),ORTHOGLOB(*),ISIGSH,IYLDINI,KSIGSH3,FAIL_INI(5),
83 . IUSOLYLD,IUSER,VARMAX
84 INTEGER, DIMENSION(NPARG,NGROUP) ,INTENT(IN):: IPARG
85 my_real
86 . geo(*),pm(npropm,*),rtrans(ntransf,*),
87 . sigi(nsigs,*),sigsh(max(1,nsigsh),*),sigtruss(nsigtruss,*),
88 . sigsp(nsigi,*),sigsph(nsigsph,*),sigrs(nsigrs,*),sigbeam(nsigbeam,*)
89C
90 TYPE(SUBMODEL_DATA) LSUBMODEL(*)
91 TYPE (GROUP_) , DIMENSION(NGRBRIC) :: IGRBRIC
92C
93 TYPE(MAPPING_STRUCT_) :: MAP_TABLES
94 TYPE (STACK_PLY) :: STACK
95 INTEGER, INTENT(IN) :: IWORKSH(3,NUMELC + NUMELTG)
96 TYPE (MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(INOUT) :: MAT_PARAM
97 INTEGER, INTENT(IN) :: NUMSPH
98 INTEGER, INTENT(IN) :: NISP
99C-----------------------------------------------
100C L o c a l V a r i a b l e s
101C-----------------------------------------------
102 INTEGER I, J, N, stat
103 INTEGER NIBRICK, NIQUAD, NISHELL, NISH3N, NISPRING, NIBEAM, NITRUSS, NISPHCEL
104 INTEGER, DIMENSION(:), ALLOCATABLE :: ID_SIGSH
105 INTEGER, DIMENSION(:), ALLOCATABLE :: ID_SOLID_SIGI
106 INTEGER, DIMENSION(:), ALLOCATABLE :: ID_QUAD_SIGI
107 INTEGER, DIMENSION(:), ALLOCATABLE :: ID_SIGSPRI
108 INTEGER, DIMENSION(:), ALLOCATABLE :: ID_SIGBEAM
109 INTEGER, DIMENSION(:), ALLOCATABLE :: ID_SIGTRUSS
110 INTEGER, DIMENSION(:), ALLOCATABLE :: ID_SIGSPH
111 INTEGER, DIMENSION(:), ALLOCATABLE :: WORK
112 LOGICAL IS_STATE
113C=======================================================================
114C
115C -- LECTURE OF INITIAL STATE DATA - EXTRACTED FROM INITIA.F
116C
117C=======================================================================
118 iuser = 0
119 iyldini = 0
120 iusolyld = 0
121 fail_ini(1:5) = 0
122 isigsh =0
123 ksigsh3 = 0
124 iortshel = 0
125 inispri = 0
126 varmax = max(nsigsh,nsigi,nsigi,nsigtruss,nsigbeam,nsigrs)
127 is_state = .false.
128C
129 ALLOCATE (id_sigsh(numshel+numsh3n) ,stat=stat)
130 ALLOCATE (id_solid_sigi(numsol) ,stat=stat)
131 ALLOCATE (id_quad_sigi(numquad) ,stat=stat)
132 ALLOCATE (id_sigspri(numspri) ,stat=stat)
133 ALLOCATE (id_sigbeam(numbeam) ,stat=stat)
134 ALLOCATE (id_sigtruss(numtrus) ,stat=stat)
135 ALLOCATE (id_sigsph(numsph) ,stat=stat)
136 ALLOCATE (work(70000) ,stat=stat)
137C
138 IF(numshel+numsh3n > 0) id_sigsh(1:numshel+numsh3n) = 0
139 IF(numsol > 0 )id_solid_sigi(1:numsol) = 0
140 IF(numquad > 0 )id_quad_sigi(1:numquad) = 0
141 IF(numspri > 0 )id_sigspri(1:numspri) = 0
142 IF(numbeam > 0 )id_sigbeam(1:numbeam) = 0
143 IF(numtrus > 0 )id_sigtruss(1:numtrus) = 0
144 work(1:70000) = 0
145C
146 IF (abs(isigi) == 3.OR.abs(isigi) == 4.OR.abs(isigi) == 5) THEN
147 DO i=1,numshel+numsh3n
148 DO j=1,nsigsh
149 sigsh(j,i)=zero
150 ENDDO
151 ENDDO
152 ENDIF
153
154C--------------------------------------------------------
155C CONTRAINTES INITIALES + ENERGIES DENSITES EPS-PLAST
156C SOLIDE-QUAD-SPRING READ ON FILE
157C--------------------------------------------------------
158
159 IF (isigi == 1.OR.isigi == 2) THEN
160C
161C FICHIER S00 (Obsolete)
162C
163 IF (ioutp_fmt == 2) THEN
164 DO i=1,numels+numelq
165 READ(iin4,'(I8,3F16.0/8X,3F16.0)') n,(sigi(j,i),j=1,6)
166 sigi(7,i) = n
167 ENDDO
168 ELSE
169 DO i=1,numels+numelq
170 READ(iin4,'(I10,3F20.0/8X,3F20.0)') n,(sigi(j,i),j=1,6)
171 sigi(7,i) = n
172 ENDDO
173 ENDIF
174
175 ELSEIF (isigi == 3.OR.isigi == 4.OR.isigi == 5) THEN
176C
177C FICHIER Y000
178C
179 CALL lec_inistate_yfile(
180 1 nsigsh ,nsigs ,nsigsph ,nsigrs ,nsigi ,
181 2 sigsh ,sigi ,sigsph ,sigrs ,sigsp ,
182 3 isigsh ,iuser ,
183 4 id_sigsh , id_solid_sigi, id_quad_sigi )
184
185 ENDIF
186
187C-----------------------------------------
188C CONTRAINTES INITIALES FICHIER D00
189C-----------------------------------------
190 IF (isigi == -3.OR.isigi == -4.OR.isigi == -5) is_state = .true.
191
192 IF (isigi == -3.OR.isigi == -4.OR.isigi == -5) THEN
193 isigi = -isigi
194! CALL LEC_INISTATE_D00 (
195! 1 IXS ,IXQ ,IXC ,IXT ,IXP ,
196! 2 IXR ,GEO ,PM ,IXTG ,INDEX ,
197! 3 ITRI ,NSIGSH ,IGEO ,
198! 4 IPM ,NSIGS ,NSIGSPH ,KSYSUSR ,NSIGRS ,
199! 5 UNITAB ,ISOLNODD00 ,LSUBMODEL ,RTRANS ,IDRAPE ,
200! 6 NSIGI ,NSIGBEAM ,NSIGTRUSS ,
201! 7 SIGI ,SIGSH ,SIGSP ,SIGSPH ,SIGRS ,
202! 8 SIGBEAM ,SIGTRUSS ,STRSGLOB ,STRAGLOB ,ORTHOGLOB ,
203! 9 ISIGSH ,IYLDINI ,FAIL_INI ,IUSOLYLD ,IUSER ,
204! A ID_SIGSH ,ID_SOLID_SIGI,ID_QUAD_SIGI ,ID_SIGSPRI ,ID_SIGBEAM,
205! b id_sigtruss,work ,igrbric )
206 CALL hm_read_inistate_d00 (
207 1 ixs ,ixq ,ixc ,ixt ,ixp ,
208 2 ixr ,geo ,pm ,ixtg ,index ,
209 3 itri ,nsigsh ,igeo ,
210 4 ipm ,nsigs ,nsigsph ,ksysusr ,nsigrs ,
211 5 unitab ,isolnodd00 ,lsubmodel ,rtrans ,idrape ,
212 6 nsigi ,nsigbeam ,nsigtruss ,
213 7 sigi ,sigsh ,sigsp ,sigsph ,sigrs ,
214 8 sigbeam ,sigtruss ,strsglob ,straglob ,orthoglob ,
215 9 isigsh ,iyldini ,fail_ini ,iusolyld ,iuser ,
216 a id_sigsh ,id_solid_sigi,id_quad_sigi ,id_sigspri ,id_sigbeam,
217 b id_sigtruss,work ,igrbric ,nibrick ,niquad ,
218 c nishell ,nish3n ,nispring ,nibeam ,nitruss ,
219 d map_tables ,varmax ,iparg ,ptshel ,ptsh3n ,
220 e stack ,iworksh ,iout ,mat_param ,nisphcel ,
221 f numsph ,nisp ,kxsp ,id_sigsph)
222 ENDIF
223C------------------------------------------------------------------------------------------
224 CALL lec_inistate_tri(
225 1 ixs ,ixq ,ixc ,ixt ,ixp ,
226 2 ixr ,kxsp ,ixtg ,index ,itri ,
227 3 nsigsh ,nsigs ,nsigsph ,ksysusr ,ksigsh3 ,
228 4 nsigrs ,nsigi ,nsigbeam ,nsigtruss ,
229 5 ptshel ,ptsh3n ,ptsol ,ptquad ,ptsph ,
230 6 ptspri ,ptbeam ,pttruss ,sigi ,sigsh ,
231 7 sigsp ,sigsph ,sigrs ,sigbeam ,sigtruss ,
232 8 id_sigsh ,id_solid_sigi,id_quad_sigi ,id_sigspri ,id_sigbeam ,
233 9 id_sigtruss,work ,id_sigsph ,is_state)
234
235C
236 IF(numsol > 0) DEALLOCATE (id_solid_sigi)
237 IF(numquad > 0 )DEALLOCATE (id_quad_sigi)
238 IF(numshel+numsh3n > 0 )DEALLOCATE (id_sigsh)
239 IF(numspri > 0 )DEALLOCATE (id_sigspri)
240 IF(numbeam > 0 )DEALLOCATE (id_sigbeam)
241 IF(numtrus > 0 )DEALLOCATE (id_sigtruss)
242 IF(numsph > 0 )DEALLOCATE (id_sigsph)
243 DEALLOCATE (work)
244C
245 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
hm_read_inistate_d00
subroutine hm_read_inistate_d00(ixs, ixq, ixc, ixt, ixp, ixr, geo, pm, ixtg, index, itri, nsigsh, igeo, ipm, nsigs, nsigsph, ksysusr, nsigrs, unitab, isolnodd00, lsubmodel, rtrans, idrape, nsigi, nsigbeam, nsigtruss, sigi, sigsh, sigsp, sigsph, sigrs, sigbeam, sigtruss, strsglob, straglob, orthoglob, isigsh, iyldini, fail_ini, iusolyld, iuser, id_sigsh, id_solid_sigi, id_quad_sigi, id_sigspri, id_sigbeam, id_sigtruss, work, igrbric, nibrick, niquad, nishell, nish3n, nispring, nibeam, nitruss, map_tables, varmax, iparg, ptshel, ptsh3n, stack, iworksh, iout, mat_param, nisphcel, numsph, nisp, kxsp, id_sigsph)
Definition hm_read_inistate_d00.F:65
lec_inistate_tri
subroutine lec_inistate_tri(ixs, ixq, ixc, ixt, ixp, ixr, kxsp, ixtg, index, itri, nsigsh, nsigs, nsigsph, ksysusr, ksigsh3, nsigrs, nsigi, nsigbeam, nsigtruss, ptshel, ptsh3n, ptsol, ptquad, ptsph, ptspri, ptbeam, pttruss, sigi, sigsh, sigsp, sigsph, sigrs, sigbeam, sigtruss, id_sigsh, id_solid_sigi, id_quad_sigi, id_sigspri, id_sigbeam, id_sigtruss, work, id_sigsphcel, is_state)
Definition lec_inistate_tri.F:43
lec_inistate_yfile
subroutine lec_inistate_yfile(nsigsh, nsigs, nsigsph, nsigrs, nsigi, sigsh, sigi, sigsph, sigrs, sigsp, isigsh, iuser, id_sigsh, id_solid_sigi, id_quad_sigi)
Definition lec_inistate_yfile.F:36
max
#define max(a, b)
Definition macros.h:21
groupdef_mod
Definition groupdef_mod.F:662
mapping_option_mod
Definition dichotomy_mod.F:36
stack_mod
Definition stack_mod.F:72
submodel_mod
Definition submodel_mod.F:712
unitab_mod
Definition unitab_mod.F:513