SOLVE /POST1 ESEL,,LIVE PLVECT,U ALLSEL FINISH SAVE 5.4.1
!注意：两个常力产生的最大位移比一个常力产生的两倍还要大。
!当单元数目为50，100，200时计算得到的跨中最大位移是一样的，所以单元数目对位移响应没有多大的影响。
/config,nres,2000
/PREP7 !建模 ET,1,BEAM3 !梁单元 MP,EX,1,207E9 !弹性模量 MP,DENS,1,2000 !密度 MP,NUXY,1,0.3
R,1,0.1,0.0001/12,0.1 !横截面积、惯性距、截面高 NELEM=100
LBEAM=32 !梁长度 deltL=LBEAM/NELEM V=120 f=1000
w=10 !移动荷载速度 NNODE=NELEM+1 deltt=deltL/V*3.6 *do,i,1,NNODE,1
N,i,(i-1)*deltL !节点点矩为0.3 *enddo *do,i,1,NELEM E,i,i+1 *enddo finish /solu d,1,all d,NNODE,all antype,trans *do,ii,1,101+8,1 outres,all,all
!输出所有加载子步的所有结果 !每个荷载步的最终求解时间
time,ii*deltt *if,ii,lt,8,then
fdele,all,all
!删除所有的集中力荷载
nsubst,5 !每个荷载步包含五个荷载子步
nsel,s,loc,x,(ii-1)*deltl f,all,FY,-f*cos(w*ii*deltt) solve
*elseif,ii,lt,101 fdele,all,all
nsel,s,loc,x,(ii-1)*deltl nsel,a,loc,x,(ii-7-1)*deltl f,all,fy,-f*cos(w*ii*deltt) allsel solve *else
fdele,all,all
nsel,s,loc,x,(ii-7-1)*deltl f,all,fy,-f*cos(w*ii*deltt) allsel solve *endif *enddo finish /post26
!选择集中力作用节点
!施加集中荷载
allsel !选择所有
NSOL,2,nelem/2+1,U,Y,UY_2 !定义第二个变量为UY_2,其值为101号节点的Y方向的位移 xvar,1 !定义时间变量为坐标横轴 plvar,2 !输出位移图
deriv,3,2,1,,,,,1 !将位移对时间求导，得到速度，存为变量3 plvar,3 !输出速度图 fini 5.4.2
!注意：两个常力产生的最大位移比一个常力产生的两倍还要大。
!当单元数目为50，100，200时计算得到的跨中最大位移是一样的，所以单元数目对位移响应没有多大的影响。 /config,nres,2000
/PREP7 !建模 ET,1,BEAM3 !梁单元 MP,EX,1,207E9 !弹性模量 MP,DENS,1,2000 !密度 MP,NUXY,1,0.3
R,1,0.1,0.0001/12,0.1 !横截面积、惯性距、截面高 NELEM=100
LBEAM=32 !梁长度 deltL=LBEAM/NELEM f=1000
V=120 !移动荷载速度 NNODE=NELEM+1
deltt=deltL/V*3.6 *do,i,1,NNODE,1
N,i,(i-1)*deltL !节点点矩为0.3 *enddo *do,i,1,NELEM E,i,i+1 *enddo finish /solu d,1,all d,NNODE,all antype,trans *do,ii,1,nnode+8,1 outres,all,all
!输出所有加载子步的所有结果 !每个荷载步的最终求解时间
time,ii*deltt *if,ii,lt,8,then
fdele,all,all f,all,FY,-f solve
*elseif,ii,lt,nnode fdele,all,all
nsel,s,loc,x,(ii-1)*deltl nsel,a,loc,x,(ii-7-1)*deltl f,all,fy,-f allsel solve *else
fdele,all,all
nsel,s,loc,x,(ii-7-1)*deltl f,all,fy,-f allsel solve *endif *enddo finish /post26
NSOL,2,nelem/2+1,U,Y,UY_2 !定义第二个变量为UY_2,其值为101号节点的Y方向的位移 xvar,1 !定义时间变量为坐标横轴 plvar,2 !输出位移图
deriv,3,2,1,,,,,1 !将位移对时间求导，得到速度，存为变量3 plvar,3 !输出速度图
!删除所有的集中力荷载
!选择集中力作用节点
nsel,s,loc,x,(ii-1)*deltl
nsubst,5 !每个荷载步包含五个荷载子步
!施加集中荷载
allsel !选择所有
fini 5.5 连续梁 /prep7
/VIEW, 1 ,1,1,1 et,1,shell63 et,2,link8
r,1,0.25 !顶层与底层的厚度 r,2,0.6 !腹板的厚度 r,3,1 !V形支撑厚度
r,4,0.02,-0.005 !预应力钢筋的面积初始与初始应变 mp,ex,1,3.25e10 mp,prxy,1,0.2 mp,dens,1,2700 mp,ex,2,207e9 mp,prxy,2,0.3 mp,dens,2,7800 !桥面关键点 k,1 k,2,4 k,3,10 k,4,16 k,5,20 k,6,24 k,7,30 k,8,36 k,9,40 SAVE
!箱形梁底面关键点 k,10,4,-1.8 k,11,10,-1.8 k,12,16,-1.8 k,13,24,-1.8 k,14,30,-1.8 k,15,36,-1.8 SAVE SAVE
kgen,2,1,15,1,,,-16 !建立梁高为1.8m的关键点 kgen,2,1,15,1,,,-85 kgen,2,1,15,1,,,-154 kgen,2,1,15,1,,,-170 SAVE
k,101,,,-36 !建立梁高为1.8m的关键点 k,102,4,,-36 k,103,10,,-36
