!横梁参数
area_=1.26e3    !面积 mm^2
Izz=8.021e5     !mm^4
Iyy=1.705e6     !mm^4
TKZ=50          !mm
TKY=50          !mm

B=1400          !梁长度 mm

!划分间隔
num_=41        !等分数量

!荷载
shuiping_li = 1.236    !水平荷载 N/mm
shuxiang_li = 906.752    !竖向荷载 N



/prep7

Et,1,beam4
mp,ex,1,0.7e5   !铝合金
mp,nuxy,1,0.3
mp,dens,1,2800e-9
R,1,area_,Izz,Iyy,TKZ,TKY

k,1,0,0,0
k,2,B,0,0
l,1,2
k,3,0,B/5,0
k,4,B,B/5,0
l,3,4
k,5,0,B*2/5,0
k,6,B,B*2/5,0
l,5,6
k,7,0,B*3/5,0
k,8,B,B*3/5,0
l,7,8

allsel
latt,1,1,1         !选择面定义材料编号是1，实常数是编号1，单元类型是1
lesize,all,,,num_      !选择所有的面进行分割，40等分
lmesh,all          !划分



/solu

dk,1,ux,,,,uy,uz
dk,2,ux,,,,uy,uz,rotx
dk,3,ux,,,,uy,uz
dk,4,ux,,,,uy,uz,rotx
dk,5,ux,,,,uy,uz
dk,6,ux,,,,uy,uz,rotx
dk,7,ux,,,,uy,uz
dk,8,ux,,,,uy,uz,rotx


hezai_li = 1.236    !水平荷载 N/mm
singlelong = 40         !单根梁的元素总量,数量为双时，采用这个函数
start_e = 1              !梁的开始元素
end_e = start_e+singlelong-1               !梁的结束元素
!这个是添加（三角形荷载）
*do,i,start_e,(start_e+end_e-1)/2,1
SFBEAM, i, 1, PRES, ((i-start_e))*hezai_li/(singlelong/2), ((i-start_e)+1)*hezai_li/(singlelong/2) !定义各个单元上的线性分布荷载
*enddo
*do,i,end_e,(start_e+end_e-1)/2+1,-1
SFBEAM, i, 1, PRES, (end_e-i+1)*hezai_li/(singlelong/2), (end_e-i)*hezai_li/(singlelong/2)
*enddo


hezai_li = 1.236    !水平荷载 N/mm
singlelong = 41         !单根梁的元素总量,数量为单时，采用这个函数
start_e = 83              !梁的开始元素
end_e = start_e+singlelong-1               !梁的结束元素
!这个是添加（三角形荷载）
*do,i,start_e,(start_e+end_e-2)/2,1
SFBEAM, i, 1, PRES, ((i-start_e))*hezai_li/((singlelong-1)/2), ((i-start_e)+1)*hezai_li/((singlelong-1)/2) !定义各个单元上的线性分布荷载
*enddo
*do,i,end_e,(start_e+end_e)/2,-1
SFBEAM, i, 1, PRES, (end_e-i+1)*hezai_li/((singlelong+1)/2), (end_e-i)*hezai_li/((singlelong+1)/2)
*enddo



hezai_li = 1.236    !水平荷载 N/mm
singlelong = 41         !单根梁的元素总量
n_sj = 10                   !分配的三角型区域的元素
start_e = 1              !梁的开始元素
end_e = start_e+singlelong-1               !梁的结束元素
!这个是添加（梯形荷载）
*do,i,start_e,start_e+n_sj-1,1
SFBEAM, i, 1, PRES, ((i-start_e))*hezai_li/(n_sj), ((i-start_e)+1)*hezai_li/(n_sj) !定义各个单元上的线性分布荷载
*enddo
*do,i,start_e+n_sj,end_e-n_sj,1
SFBEAM, i, 1, PRES, hezai_li !定义各个单元上的线性分布荷载
*enddo
*do,i,end_e,end_e-n_sj+1,-1
SFBEAM, i, 1, PRES, (end_e-i+1)*hezai_li/(n_sj), (end_e-i)*hezai_li/(n_sj)
*enddo




allsel
SOLVE