均布荷载.txt1.均布荷载(1)y方向荷载施加*set,q,-350000*do,i,3,20a=q*abs(nx(i)-nx(i+1))/2注意:此时布置的均布荷载是实际均布荷载的一半f,i,fy,af,i+1,fy,a*enddo*do,i,43,60a=q*abs(nx(i)-nx(i+1))/2f,i,fy,af,i+1,fy,a*enddoa=q*abs(nx(1)-nx(3))/2f,1,fy,af,3,fy,aa=q*abs(nx(1)-nx(61))/2f,1,fy,af,61,fy,a(2)x方向荷载施加!左侧*set,e,45000*do,i,3,30a=e*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a*enddoa=e*abs(ny(2)-ny(31))/2,2,fx,af,31,fx,aa=e*abs(ny(1)-ny(3))/2f,1,fx,a第1页均布荷载.txtf,3,fx,a*do,i,63,70a=e*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a*enddoa=e*abs(ny(2)-ny(63))/2f,2,fx,af,63,fx,a*do,i,82,88a=e*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a*enddoa=e*abs(ny(62)-ny(82))/2f,62,fx,af,82,fx,a右侧*do,i,32,60a=-e*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a*enddoa=-e*abs(ny(61)-ny(1))/2f,61,fx,af,1,fx,a*do,i,73,80a=-e*abs(ny(i)-ny(i+1))/2*do,i,73,80a=-e*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a*enddoa=-e*abs(ny(73)-ny(32))/2f,73,fx,a第2页均布荷载.txtf,32,fx,aa=-e*abs(ny(72)-ny(81))/2f,72,fx,af,81,fx,a*do,i,89,95a=-e*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a1.均布荷载(1)y方向荷载施加*do,i,1,10a=q*abs(nx(i)-nx(i+1))/2f,i,fy,af,i+1,fy,a*enddo(2)x方向荷载施加*do,i,1,10a=q*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a*enddo2.梯形荷载*set,shuxiang(竖向),292608.9*set,shangce(上侧),78149.5*set,xiace(下侧),127057.6!竖向荷载*do,i,73,80a=-e*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a*enddoa=-e*abs(ny(73)-ny(32))/2f,73,fx,af,32,fx,aa=-e*abs(ny(72)-ny(81))/2f,72,fx,af,81,fx,a第3页均布荷载.txt*do,i,89,95a=-e*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a1.均布荷载(1)y方向荷载施加*do,i,1,10a=q*abs(nx(i)-nx(i+1))/2f,i,fy,af,i+1,fy,a*enddo(2)x方向荷载施加*do,i,1,10a=q*abs(ny(i)-ny(i+1))/2f,i,fx,af,i+1,fx,a*enddo2.梯形荷载*set,shuxiang(竖向),292608.9*set,shangce(上侧),78149.5*set,xiace(下侧),127057.6!竖向荷载*do,i,5,18a=-shuxiang*0.5*abs(nx(i)-nx(i+1))f,i+1,fy,af,i,fy,a*enddo!左水平荷载(11.82为梯形的高度)!3指梯形上底对应的节点号!12指梯形下底对应的节点号*do,i,3,11a=abs(ny(i)-ny(12))*(xiace-shangce)/11.82+shangceb=abs(ny(i+1)-ny(12))*(xiace-shangce)/11.82+shangcef,i+1,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))f,i,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))*enddo*do,i,22,26a=abs(ny(i)-ny(12))*(xiace-shangce)/11.82+shangce第4页均布荷载.txtb=abs(ny(i+1)-ny(12))*(xiace-shangce)/11.82+shangcef,i+1,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))f,i,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))*enddoa=abs(ny(1)-ny(12))*(xiace-shangce)/11.82+shangceb=abs(ny(3)-ny(12))*(xiace-shangce)/11.82+shangcef,3,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))f,1,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))a=abs(ny(1)-ny(12))*(xiace-shangce)/11.82+shangceb=abs(ny(22)-ny(12))*(xiace-shangce)/11.82+shangcef,22,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))f,1,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))!右水平荷载*do,i,12,20a=abs(ny(i)-ny(12))*(xiace-shangce)/11.82+shangceb=abs(ny(i+1)-ny(12))*(xiace-shangce)/11.82+shangcef,i+1,fx,-(a/6+b/3)*abs(ny(i)-ny(i+1))f,i,fx,-(a/3+b/6)*abs(ny(i)-ny(i+1))*enddoa=abs(ny(2)-ny(12))*(xiace-shangce)/11.82+shangceb=abs(ny(21)-ny(12))*(xiace-shangce)/11.82+shangcef,21,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))f,2,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))a=abs(ny(2)-ny(12))*(xiace-shangce)/11.82+shangceb=abs(ny(32)-ny(12))*(xiace-shangce)/11.82+shangcef,32,fx,(a/6+b/3)*abs(ny(i)-ny(i+1))f,2,fx,(a/3+b/6)*abs(ny(i)-ny(i+1))上述程序中,i为节点号,q为竖向围岩压力,a为等效节点力计算中间值,nx(i)为i节点的水平坐标值,竖直坐标用ny(i)表示。第5页