!   beam_cs.mac
!
mod_assemb_prep
mod_coord_system
!
create_cs,0,'np',0,5,5,5,0,0,0
!
angweb=arg1
pload=100
wid=5
hght=5
L=10*wid
cswid=0.25
!
et,1,shell63
et,2,solid45
mat_select,'c410'
!
/title, Angled Web Beam Cross-Section / Web Angle = %angweb%
!
k,,-(1.2*(wid/2)),hght/2,0
k,,(1.2*(wid/2)),hght/2,0
k,,(1.2*(wid/2)),(hght/2)+cswid,0
k,,-(1.2*(wid/2)),(hght/2)+cswid,0
!
create_area
!
k,,-(1.2*(wid/2)),-(hght/2)-cswid,0
k,,(1.2*(wid/2)),-(hght/2)-cswid,0
k,,(1.2*(wid/2)),-hght/2,0
k,,-(1.2*(wid/2)),-hght/2,0
!
create_area
!
create_cs,0,'vweb',0,0,0,0,angweb,0,0
!
k,,-cswid/2,-(hght/2)/cos(angweb*rconv),0
k,,cswid/2,-(((hght/2)/cos(angweb*rconv))+(cswid*tan(angweb*rconv))),0
k,,cswid/2,(hght/2)/cos(angweb*rconv),0
k,,-cswid/2,((hght/2)/cos(angweb*rconv))+(cswid*tan(angweb*rconv)),0
!
create_area
!
aadd,all
!
init_entity_num
!
extrude_aset,L
!
numcmp,all
init_entity_num
!
csys,np
!
k,,-(1.2*(wid/2)),-(hght/2)-cswid,0
k,,(1.2*(wid/2)),-(hght/2)-cswid,0
k,,(1.2*(wid/2)),(hght/2)+cswid,0
k,,-(1.2*(wid/2)),(hght/2)+cswid,0
!
create_volume,,cswid
!
!k,,-(1.2*(wid/2)),-(hght/2)-cswid,L-cswid
!k,,(1.2*(wid/2)),-(hght/2)-cswid,L-cswid
!k,,(1.2*(wid/2)),(hght/2)+cswid,L-cswid
!k,,-(1.2*(wid/2)),(hght/2)+cswid,L-cswid
!
!create_volume,,cswid
!
numcmp,all
init_entity_num
!
vovlap,all
!
init_entity_num
!
lplot
!
!------------------------------------------------------------------------
!
!
!                (1)     (2)     (3)      (4)      (5)   (6)   (7)    (8)   (9)
! MSH_STANDRD, perfmsh, mstyp, mshdens, elmszovr, etpa, etpv, elmnl, elmxl, mat,
!
!              real, perflsz
!              (10)   (11)
!
csys,np
!
asel,s,loc,z,0
msh_standrd,3,0,0,1,1,2,0.2,0.2
!
!asel,s,loc,z,L-(cswid/3),L+(cswid/3)
!msh_standrd,3,0,0,1,1,2,0.2,0.2
!
!asel,s,area,,am-19
!msh_standrd,2,0,0,-2,1,2,0.2,0.2
!
asel,s,area,,am-12
msh_standrd,2,0,0,-2,1,2,0.2,0.2
!


!*go,:skip1

nummrg,all
init_entity_num
!
lplot
!
!------------------------------------------------------------------------
!
csys,np
!
nsel,s,loc,z,0
d,all,all
!
init_entity_num
!
/PSF,PRES,NORM,2,0
asel,s,loc,y,(hght/2)+cswid
sfa,all,,pres,pload
!
init_entity_num
!
lplot
!
!------------------------------------------------------------------------
!
/solu
eqslv,iter,5,3
solve
!
lplot
!

:skip1