offz = -10
/com,
/com, Define a beam channel section. Unequal lengths and thicknesses are
/com, used to demonstrate an unsymmetric section with noncoincident
/com, centroid and shear center.
/com,
/com, An arbitrary nodal offset is used to demonstrate coupling between
/com, axial and bending stiffness, as well as coupling between shear
/com, and torsional stiffness.
/com,
sectype,1,beam,chan
secdata,1,1.5,2, 0.1, 0.15, 0.2
secoffset,user,offy,offz
/com, For the channel section, generate elements along the x-axis
secnum,1
n,1
do,i,1,nelem
loc = ileng/nelem
n,i+1,loc
e,i,i+1
*enddo
/com,
/com, Retrieve the essential calculated and geometric properties from
/com, the channel section
/com,
get,xarea,secp,1,prop,area
get,iyy,secp,1,prop,iyy
get,izz,secp,1,prop,izz
get,iyz,secp,1,prop,iyz
get,J,secp,1,prop,tors
get,iw,secp,1,prop,warp
get,scyy,secp,1,prop,scyy
get,sczz,secp,1,prop,sczz
get,cgy,secp,1,prop,cgy
get,cgz,secp,1,prop,cgz
get,offy,secp,1,prop,offy
get,offz,secp,1,prop,offz
get,shcy,secp,1,prop,shcy
get,shcz,secp,1,prop,shcz
yc = offy - cgy! Distance between nodal offset and centroid
zc = offz - cgz
ys = offy - shcy! Distance between nodal offset and shear center
zs = offz - shcz
yssq = ysys
zssq = zszs
/com,
/com, Adjust the preintegrated torsional stiffness based on nodal offset
k44 = GxyJ
k44 = k44 + ysysGxyxareasczz
k44 = k44 + zszsGxyxarea*scyy
/com,
/com, Use the material and channel section properties to define the 7x7
/com, upper triangle of the preintegrated section stiffness matrix
/com,
sectype,2,comb,matrix
cbmx,1, Exxxarea, -zcExxxarea, ycExxxarea
cbmx,2, Exxiyy+zczcExxxarea, -ExxIyz-yczcExxxarea
cbmx,3, Exxizz+ycycExxxarea
cbmx,4, k44, -ysGxyxareasczz, zsGxyxareascyy
cbmx,5, Gxyxareasczz
cbmx,6, Gxyxareascyy
cbmx,7, Exxiw
slist
/com,
Release 15.0 - © SAS IP, Inc. All rights reserved. - Contains proprietary and confidential information
Using Preintegrated Composite Beam Sections