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

Go to the source code of this file.

Functions/Subroutines

subroutine h3d_oned_vector (elbuf_tab, ifunc, iparg, geo, ixt, ixp, ixr, pm, anim, oned_vector, id_elem, ity_elem, info1, info2, is_written_oned, ipartt, ipartp, ipartr, h3d_part, keyword, x, d, tors)

Function/Subroutine Documentation

◆ h3d_oned_vector()

subroutine h3d_oned_vector ( type (elbuf_struct_), dimension(ngroup), target elbuf_tab,
integer ifunc,
integer, dimension(nparg,*) iparg,
geo,
integer, dimension(nixt,*) ixt,
integer, dimension(nixp,*) ixp,
integer, dimension(nixr,*) ixr,
pm,
anim,
oned_vector,
integer, dimension(*) id_elem,
integer, dimension(*) ity_elem,
integer info1,
integer info2,
integer, dimension(*) is_written_oned,
integer, dimension(*) ipartt,
integer, dimension(*) ipartp,
integer, dimension(*) ipartr,
integer, dimension(*) h3d_part,
character(len=ncharline100) keyword,
x,
d,
tors )

Definition at line 34 of file h3d_oned_vector.F.

40C-----------------------------------------------
41C M o d u l e s
42C-----------------------------------------------
43 USE elbufdef_mod
44 USE names_and_titles_mod, ONLY: ncharline100
45 use element_mod , only : nixt,nixr,nixp
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C C o m m o n B l o c k s
52C-----------------------------------------------
53#include "mvsiz_p.inc"
54#include "com01_c.inc"
55#include "com04_c.inc"
56#include "param_c.inc"
57C-----------------------------------------------
58C D u m m y A r g u m e n t s
59C-----------------------------------------------
60C REAL
61 my_real
62 . pm(npropm,*), geo(npropg,*),
63 . anim(*),oned_vector(3,*),x(3,*),d(3,*),tors(15,*)
64 INTEGER IPARG(NPARG,*),
65 . IXT(NIXT,*),IXP(NIXP,*),IXR(NIXR,*),IFUNC,
66 . NANIM1D_L,
67 . IS_WRITTEN_ONED(*),ID_ELEM(*),ITY_ELEM(*),
68 . IPARTT(*) ,IPARTP(*),IPARTR(*),H3D_PART(*)
69 INTEGER INFO1, INFO2
70C
71 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
72 CHARACTER(LEN=NCHARLINE100) :: KEYWORD
73C-----------------------------------------------
74C L o c a l V a r i a b l e s
75C-----------------------------------------------
76C REAL
77 my_real
78 . evar(mvsiz), mass(mvsiz) ,
79 . off, p, vonm2, vonm, s1, s2, s12, s3, value(3),
80 . a1,b1,b2,b3,yeq,f1,m1,m2,m3, xm,
81 . for, area, feq, eplas, rho0, a0, xx1, yy1, zz1, al0
82 INTEGER I, NG, NEL, NFT, ITY,
83 . IADD, N, J, LLT, MLW, NB1, NB2, NB3, NB4, NB5,
84 . NB6, NB7, NB8, NB9, NB10, NB11, NB12, NB13, NB14, NB15,
85 . NB16, LLL,NUVAR,IGTYP,
86 . ISTRAIN,NN, K1, K2,JTURB,MT,JALE, IMID, IALEL,IPID,
87 . NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9,NN10,NF,
88 . OFFSET,K,INC,KK,IHBE,ISROT,ILAYER,IR,IS,JJ(6),IOK_PART(MVSIZ),
89 . IS_WRITTEN_VALUE(MVSIZ),N1, N2
90
91C
92 TYPE(G_BUFEL_) ,POINTER :: GBUF
93
94C-----------------------------------------------
95C
96 nn1 = 1
97 nn3 = 1
98 nn4 = nn3
99 nn5 = nn4
100 nn6 = nn5
101 nn7 = nn6 + numelt
102 nn8 = nn7 + numelp
103 nn9 = nn8 + numelr
104 nn10= nn9
105C
106 DO ng=1,ngroup
107 mlw =iparg(1,ng)
108 nel =iparg(2,ng)
109 ity =iparg(5,ng)
110 igtyp =iparg(38,ng)
111C---
112 gbuf => elbuf_tab(ng)%GBUF
113C---
114 nft =iparg(3,ng)
115!
116 DO i=1,6
117 jj(i) = nel*(i-1)
118 ENDDO
119
120 value(1:3) = zero
121 DO i=1,nel
122 is_written_value(i) = 0
123 ENDDO
124c
125 IF (ity == 4) offset = 0
126 IF (ity == 5) offset = numelt
127 IF (ity == 6) offset = numelt+numelp
128c
129 DO i=1,nel
130 IF (ity == 4) THEN
131 id_elem(offset+nft+i) = ixt(nixt,nft+i)
132 ity_elem(offset+nft+i) = 4
133 IF( h3d_part(ipartt(nft+i)) == 1) iok_part(i) = 1
134 ELSEIF (ity == 5) THEN
135 id_elem(offset+nft+i) = ixp(nixp,nft+i)
136 ity_elem(offset+nft+i) = 5
137 IF( h3d_part(ipartp(nft+i)) == 1) iok_part(i) = 1
138 ELSEIF (ity == 6) THEN
139 id_elem(offset+nft+i) = ixr(nixr,nft+i)
140 ity_elem(offset+nft+i) = 6
141 IF( h3d_part(ipartr(nft+i)) == 1) iok_part(i) = 1
142 ENDIF
143 ENDDO
144
145 IF(ity==4 .OR. ity==5 .OR. ity==6)THEN
146 DO i=1,nel
147 oned_vector(1:3,offset+nft+i) = zero ! Default = zero in all cases !
148 ENDDO
149 ENDIF
150C-----------------------------------------------
151C TRUSS
152C-----------------------------------------------
153 IF(ity==4)THEN
154C--------------------------------------------------
155 IF (keyword == 'FINT') THEN
156C--------------------------------------------------
157 DO i=1,nel
158 value(1) = tors(1,offset+nft+i)
159 value(2) = tors(2,offset+nft+i)
160 value(3) = tors(3,offset+nft+i)
161 CALL h3d_write_vector(iok_part,is_written_oned,oned_vector,i,offset,nft,
162 . VALUE)
163 ENDDO
164C--------------------------------------------------
165c ELSEIF (KEYWORD == '') THEN
166C--------------------------------------------------
167c DO I=1,NEL
168c VALUE(1) =
169c VALUE(2) =
170c VALUE(3) =
171c CALL H3D_WRITE_VECTOR(IOK_PART,IS_WRITTEN_ONED,NODAL_VECTOR,I,0,0,
172c . VALUE)
173c ENDDO
174 ENDIF
175C-----------------------------------------------
176C POUTRES
177C-----------------------------------------------
178 ELSEIF(ity==5)THEN
179C--------------------------------------------------
180 IF (keyword == 'FINT') THEN
181C--------------------------------------------------
182 DO i=1,nel
183 value(1) = tors(1,offset+nft+i)
184 value(2) = tors(2,offset+nft+i)
185 value(3) = tors(3,offset+nft+i)
186 CALL h3d_write_vector(iok_part,is_written_oned,oned_vector,i,offset,nft,
187 . VALUE)
188 ENDDO
189C--------------------------------------------------
190c ELSEIF (KEYWORD == '') THEN
191C--------------------------------------------------
192c DO I=1,NEL
193c VALUE(1) =
194c VALUE(2) =
195c VALUE(3) =
196c CALL H3D_WRITE_VECTOR(IOK_PART,IS_WRITTEN_ONED,NODAL_VECTOR,I,0,0,
197c . VALUE)
198c ENDDO
199 ENDIF
200C-----------------------------------------------
201C RESSORTS
202C-----------------------------------------------
203 ELSEIF(ity==6)THEN
204C--------------------------------------------------
205 IF (keyword == 'FINT') THEN
206C--------------------------------------------------
207 DO i=1,nel
208 value(1) = tors(1,offset+nft+i)
209 value(2) = tors(2,offset+nft+i)
210 value(3) = tors(3,offset+nft+i)
211 CALL h3d_write_vector(iok_part,is_written_oned,oned_vector,i,offset,nft,
212 . VALUE)
213 ENDDO
214C--------------------------------------------------
215c ELSEIF (KEYWORD == '') THEN
216C--------------------------------------------------
217c DO I=1,NEL
218c VALUE(1) =
219c VALUE(2) =
220c VALUE(3) =
221c CALL H3D_WRITE_VECTOR(IOK_PART,IS_WRITTEN_ONED,NODAL_VECTOR,I,0,0,
222c . VALUE)
223c ENDDO
224 ENDIF
225 ENDIF
226 ENDDO
227
228 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
h3d_write_vector
subroutine h3d_write_vector(iok_part, is_written, vector, i, offset, nft, value)
Definition h3d_write_vector.F:36
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
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