OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cmaini3.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "vect01_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine cmaini3 (elbuf_str, pm, geo, nel, nlay, skew, igeo, ix, nix, numel, nsigsh, sigsh, ptsh, igtyp, iorthloc, ipm, propid, aldt, mat_param, ir, is, isubstack, stack, irep, drape, shang, geo_stack, igeo_stack, igmat, imat, iprop, nummat, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, x, npt_all, idrape, numel_drape, indx)

Function/Subroutine Documentation

◆ cmaini3()

subroutine cmaini3 ( type(elbuf_struct_), target elbuf_str,
pm,
geo,
integer nel,
integer nlay,
skew,
integer, dimension(npropgi,*) igeo,
integer, dimension(nix,*) ix,
integer nix,
integer numel,
integer nsigsh,
sigsh,
integer, dimension(*) ptsh,
integer igtyp,
integer, dimension(mvsiz) iorthloc,
integer, dimension(npropmi,*) ipm,
integer propid,
aldt,
type (matparam_struct_), dimension(nummat), intent(in) mat_param,
integer ir,
integer is,
integer isubstack,
type (stack_ply) stack,
integer irep,
type(drape_), dimension(numelc_drape + numeltg_drape) drape,
shang,
geo_stack,
integer, dimension(*) igeo_stack,
integer igmat,
integer imat,
integer iprop,
integer, intent(in) nummat,
intent(in) x1,
intent(in) x2,
intent(in) x3,
intent(in) x4,
intent(in) y1,
intent(in) y2,
intent(in) y3,
intent(in) y4,
intent(in) z1,
intent(in) z2,
intent(in) z3,
intent(in) z4,
intent(in) e1x,
intent(in) e2x,
intent(in) e3x,
intent(in) e1y,
intent(in) e2y,
intent(in) e3y,
intent(in) e1z,
intent(in) e2z,
intent(in) e3z,
x,
integer npt_all,
integer idrape,
integer numel_drape,
integer, dimension(numel_drape) indx )

Definition at line 40 of file cmaini3.F.

51C-----------------------------------------------
52C M o d u l e s
53C-----------------------------------------------
54 USE mat_elem_mod
55 USE message_mod
56 USE stack_mod
57 USE drape_mod
58 USE names_and_titles_mod , ONLY : nchartitle
59C-----------------------------------------------
60C I m p l i c i t T y p e s
61C-----------------------------------------------
62#include "implicit_f.inc"
63C-----------------------------------------------
64C G l o b a l P a r a m e t e r s
65C-----------------------------------------------
66#include "mvsiz_p.inc"
67C-----------------------------------------------
68C C o m m o n B l o c k s
69C-----------------------------------------------
70#include "param_c.inc"
71#include "vect01_c.inc"
72C-----------------------------------------------
73C D u m m y A r g u m e n t s
74C-----------------------------------------------
75 INTEGER ,INTENT(IN) :: NUMMAT
76 INTEGER NEL,NLAY,NIX,NSIGSH,NUMEL,IGTYP,PROPID,IR,IS,IREP,IGMAT,IMAT,IPROP,
77 . IDRAPE ,NPT_ALL,NUMEL_DRAPE
78 INTEGER IGEO(NPROPGI,*),IX(NIX,*),PTSH(*),IORTHLOC(MVSIZ),
79 . IPM(NPROPMI,*),ISUBSTACK,IGEO_STACK(*)
80 INTEGER , DIMENSION(NUMEL_DRAPE) :: INDX
81 my_real
82 . pm(*),geo(npropg,*),skew(lskew,*),aldt(*),
83 . sigsh(nsigsh,*),shang(*),geo_stack(*),x(3,*)
84 TYPE(ELBUF_STRUCT_), TARGET :: ELBUF_STR
85 TYPE (STACK_PLY) :: STACK
86 my_real, DIMENSION(MVSIZ), INTENT(IN) :: e1x,e2x,e3x,e1y,e2y,e3y,e1z,e2z,e3z,
87 . x1,x2,x3,x4,y1,y2,y3,y4,z1,z2,z3,z4
88 TYPE (MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(IN) :: MAT_PARAM
89C-----------------------------------------------
90C L o c a l V a r i a b l e s
91C-----------------------------------------------
92 INTEGER IT,IL,NUVAR,NPTT,ILAW,NUPARAM
93 my_real
94 . vx(mvsiz),vy(mvsiz),vz(mvsiz),phi1(npt_all,mvsiz),phi2(npt_all,mvsiz)
95 my_real
96 . coor1(npt_all,mvsiz),coor2(npt_all,mvsiz),coor3(npt_all,mvsiz),
97 . coor4(npt_all,mvsiz)
98 CHARACTER(LEN=NCHARTITLE) :: TITR,TITR1
99 my_real,
100 . DIMENSION(:),POINTER :: uvar,dir1,dir2
101 TYPE(G_BUFEL_) ,POINTER :: GBUF
102 TYPE(L_BUFEL_) ,POINTER :: LBUF
103 TYPE(BUF_LAY_) ,POINTER :: BUFLY
104 TYPE(DRAPE_) :: DRAPE(NUMELC_DRAPE + NUMELTG_DRAPE)
105C=======================================================================
106 coor1 = zero
107 coor2 = zero
108 coor3 = zero
109 coor4 = zero
110C
111 gbuf => elbuf_str%GBUF
112C----
113 CALL corthini(
114 . lft ,llt ,nft ,nlay ,numel ,
115 . nsigsh ,nix ,ix ,igeo ,geo ,
116 . skew ,sigsh ,ptsh ,phi1 ,phi2 ,
117 . vx ,vy ,vz ,coor1 ,coor2 ,
118 . coor3 ,coor4 ,iorthloc ,isubstack ,stack ,
119 . irep ,elbuf_str ,drape ,shang ,x ,
120 . geo_stack ,e3x ,e3y ,e3z ,
121 . gbuf%BETAORTH,x1 ,x2 ,y1 ,y2 ,
122 . z1 ,z2 ,nel ,gbuf%G_ADD_NODE,gbuf%ADD_NODE,
123 . npt_all ,idrape ,indx)
124c---
125C----
126 IF(igtyp == 51 .OR. igtyp == 52 .OR. igmat > 0) THEN
127 CALL corthdir(elbuf_str,
128 . igeo ,geo ,vx ,vy ,vz ,
129 . phi1 ,phi2 ,coor1 ,coor2 ,coor3 ,
130 . coor4 ,iorthloc ,nlay ,irep ,isubstack,
131 . stack ,geo_stack ,igeo_stack,ir ,is ,
132 . nel ,imat ,iprop ,
133 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
134 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
135 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z ,
136 . npt_all ,idrape)
137
138 ELSEIF (mtn == 27) THEN
139 CALL cm27in3(elbuf_str,
140 . geo ,igeo ,pm ,ipm ,ix ,nix,
141 . nlay,ir ,is ,imat )
142 ELSEIF (mtn==15 .or. mtn==19 .or. mtn==25 .or. mtn >= 28) THEN
143 IF (mtn == 19 .and. igtyp /= 9) THEN
144 CALL ancmsg(msgid=5,
145 . anmode=aninfo,
146 . msgtype=msgerror,
147 . i1=ix(1,1))
148 ENDIF
149c
150 CALL corthdir(elbuf_str,
151 . igeo ,geo ,vx ,vy ,vz ,
152 . phi1 ,phi2 ,coor1 ,coor2 ,coor3 ,
153 . coor4 ,iorthloc ,nlay ,irep ,isubstack,
154 . stack ,geo_stack ,igeo_stack,ir ,is ,
155 . nel ,imat ,iprop ,
156 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
157 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
158 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z ,
159 . npt_all ,idrape )
160 ENDIF
161C-----
162 IF ((mtn == 58 .or. mtn == 158) .AND.
163 . igtyp /= 16 .AND. igtyp /= 51 .AND. igtyp /= 52) THEN
164 CALL ancmsg(msgid=658,
165 . msgtype=msgerror,
166 . anmode=aninfo_blind_1,
167 . i1=propid,
168 . c1=titr,
169 . i2=mtn,
170 . i3=igtyp)
171 ELSEIF (mtn == 58 .or. mtn == 158 .OR. igtyp == 51 .OR. igtyp == 52) THEN
172C
173 IF (idrape == 0) THEN
174 DO il = 1,nlay
175 nptt = elbuf_str%BUFLY(il)%NPTT
176 imat = elbuf_str%BUFLY(il)%IMAT
177 ilaw = elbuf_str%BUFLY(il)%ILAW
178 nuvar = elbuf_str%BUFLY(il)%NVAR_MAT
179 dir1 => elbuf_str%BUFLY(il)%DIRA
180 dir2 => elbuf_str%BUFLY(il)%DIRB
181 nuparam = mat_param(imat)%NUPARAM
182!
183 IF (ilaw == 58) THEN
184 DO it=1,nptt
185 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
186 uvar => elbuf_str%BUFLY(il)%MAT(ir,is,it)%VAR
187 CALL cm58in3(irep ,dir1 ,dir2 ,mat_param(imat)%UPARAM,
188 . uvar ,aldt ,nel ,nuvar ,lbuf%ANG ,
189 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
190 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
191 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
192 ENDDO
193 ELSE IF (ilaw == 158) THEN
194 DO it=1,nptt
195 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
196 uvar => elbuf_str%BUFLY(il)%MAT(ir,is,it)%VAR
197 CALL law158_init(dir1 ,dir2 ,
198 . uvar ,aldt ,nel ,nuvar ,lbuf%ANG ,
199 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
200 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
201 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
202 ENDDO
203 ENDIF ! ILAW
204 ENDDO ! DO IL = 1,NLAY
205 ELSE
206 DO il = 1,nlay
207 nptt = elbuf_str%BUFLY(il)%NPTT
208 imat = elbuf_str%BUFLY(il)%IMAT
209 ilaw = elbuf_str%BUFLY(il)%ILAW
210 nuvar = elbuf_str%BUFLY(il)%NVAR_MAT
211 nuparam = mat_param(imat)%NUPARAM
212!
213 IF (ilaw == 58) THEN
214 DO it=1,nptt
215 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
216 uvar => elbuf_str%BUFLY(il)%MAT(ir,is,it)%VAR
217 dir1 => elbuf_str%BUFLY(il)%LBUF_DIR(it)%DIRA
218 dir2 => elbuf_str%BUFLY(il)%LBUF_DIR(it)%DIRB
219 CALL cm58in3(irep ,dir1 ,dir2 ,mat_param(imat)%UPARAM,
220 . uvar ,aldt ,nel ,nuvar ,lbuf%ANG ,
221 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
222 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
223 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
224 ENDDO
225 ELSE IF (ilaw == 158) THEN
226 DO it=1,nptt
227 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
228 uvar => elbuf_str%BUFLY(il)%MAT(ir,is,it)%VAR
229 dir1 => elbuf_str%BUFLY(il)%LBUF_DIR(it)%DIRA
230 dir2 => elbuf_str%BUFLY(il)%LBUF_DIR(it)%DIRB
231 CALL law158_init(dir1 ,dir2 ,
232 . uvar ,aldt ,nel ,nuvar ,lbuf%ANG ,
233 . x1 ,x2 ,x3 ,x4 ,y1 ,y2 ,
234 . y3 ,y4 ,z1 ,z2 ,z3 ,z4 ,
235 . e1x, e2x, e3x, e1y, e2y, e3y ,e1z, e2z, e3z )
236 ENDDO
237 ENDIF ! ILAW
238 ENDDO ! DO IL = 1,NLAY
239 ENDIF
240 ENDIF
241C-----------
242 RETURN
cm27in3
subroutine cm27in3(elbuf_str, geo, igeo, pm, ipm, ix, nix, nlay, ir, is, imat)
Definition cm27in3.F:38
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
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
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
drape_mod
Definition drape_mod.F:54
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
stack_mod
Definition stack_mod.F:72
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:889