两个荷载结构法ansys命令流实例

实例一:finish/clear/COM,Structural/prep7*AFUN,deg*set,Py,275461.0591*set,px1,71390*set,px2,116230.464*set,CQHD,0.5*set,CQDYCD,0.6*set,CQETXML,31e9*set,CQUBSB,0.2*set,WYTXKL,200e6*set,WYTXKLZGQ,850e6*set,CQMD,2300et,1,beam3r,1,0.5,0.5*0.5*0.5/12,0.5mp,ex,1,CQETXMLmp,dens,1,CQMDmp,prxy,1,CQUBSBk,1,311.419293,-97.854439k,2,308.756764,-98.346462k,3,306.445920,-99.757543k,4,304.948061,-101.812347k,5,304.418078,-104.299297k,6,304.452492,-105.061153k,7,304.555456,-105.816802k,8,304.958136,-106.652780k,9,305.731611,-107.165369k,10,308.539368,-107.845632k,11,311.419293,-108.074297k,12,314.299218,-107.845632k,13,317.106975,-107.165369k,14,317.880450,-106.652780k,15,318.283130,-105.816802k,16,318.386093,-105.061153k,17,318.420508,-104.299297k,18,317.890524,-101.812347k,19,316.392666,-99.757543k,20,314.081822,-98.346462!二衬轴线larc,1,3,2larc,3,5,4larc,5,7,6larc,7,9,8larc,9,11,10larc,11,13,12larc,13,15,14larc,15,17,16larc,17,19,18larc,19,1,20allsellsel,s,loc,x,311.419293,302.247076lcomb,allnummrg,alllesize,all,,,50allsellsel,s,loc,x,311.419293,320.122358lcomb,allnummrg,alllesize,all,,,50allslmesh,alleplotnummrg,all,,,,low!所有号数从1开始numcmp,node,eord!压缩节点号并排序@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@!施加水平弹簧支撑单元*do,i,16,50a1=abs(ny(i-1)-ny(i+1))*WYTXKLkx=a1/2PSPRNG,i,TRAN,kx,-0.3,,,,*enddo*do,i,52,61a1=abs(ny(i-1)-ny(i+1))*WYTXKLkx=a1/2PSPRNG,i,TRAN,kx,0.3,,,,*enddo*do,i,62,79a1=abs(ny(i-1)-ny(i+1))*WYTXKLZGQkx=a1/2PSPRNG,i,TRAN,kx,0.3,,,,*enddo*do,i,80,86a1=abs(ny(i-1)-ny(i+1))*WYTXKLkx=a1/2PSPRNG,i,TRAN,kx,0.3,,,,*enddo!施加竖直弹簧支撑单元*do,i,29,50b1=abs(nx(i-1)-nx(i+1))*WYTXKLky=b1/2PSPRNG,i,TRAN,ky,,-0.3,,,*enddob1=abs(nx(52)-nx(50))*WYTXKLky=b1/2PSPRNG,51,TRAN,ky,,-0.3,,,*do,i,52,61b1=abs(nx(i-1)-nx(i+1))*WYTXKLky=b1/2PSPRNG,i,TRAN,ky,,-0.3,,,*enddo*do,i,62,73b1=abs(nx(i-1)-nx(i+1))*WYTXKLZGQky=b1/2PSPRNG,i,TRAN,ky,,-0.3,,,*enddo!施加边界条件与围岩压力/solunsel,s,,,51!选择约束的节点,仰拱中间节点d,all,ux!施加水平方向的约束allselfcum,add,,!竖向荷载*do,i,2,28a=-Py*0.5*abs(nx(i-1)-nx(i+1))f,i,fy,a*enddo*do,i,80,99a=-Py*0.5*abs(nx(i-1)-nx(i+1))f,i,fy,a*enddoa=-Py*0.5*abs(nx(1)-nx(99))f,100,fy,aa=-Py*0.5*abs(nx(2)-nx(100))f,1,fy,a!水平向荷载（X方向压力，全部节点施加水平荷载，换算成节点荷载，在节点上施加）*do,i,52,61a=-(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5f,i,fX,a*enddo*do,i,80,99a=-(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5f,i,fX,a*enddoa=-(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(100)))*abs(ny(1)-ny(99))*0.5f,100,fX,a*do,i,2,50a=(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5f,i,fX,a*enddoacel,9.8,solvefinish!进入后处理，定义荷载工况并组合，输出图片和文本文件/post1ESEL,s,TYPE,,1!选择单元1（隧道轴力中心线单元）NSLE,s,ALL!选择所有节点etable,mi,smisc,6etable,mj,smisc,12plls,mi,mj,-1etable,fi,smisc,1etable,fj,smisc,7plls,fi,fj,1pldisp,1etable,reflpretab,mi,mj,fi,fj实例二:/prep7*set,shuxiang,276170!竖向荷载*set,zuoce,125500!水平荷载*set,youce,-125500*set,jizhongli,-100000!地震集中力*set,uxishu,136.966e6*0.5!水平弹性抗力系数*set,vxishu,87.4348e6*0.5!竖向弹性抗力系数!定义单元类型、实常数及材料性质et,1,beam3r,1,0.3,0.00225,0.3mp,ex,1,3.1e10!弹性模量mp,dens,1,2500!密度mp,prxy,1,0.167!泊松比!创建几何模型画点k,1,0,2.55k,2,-2.04,1.53k,3,-2.6,-0.15k,4,-2.319721,-2.023237k,5,-1.648852,-2.71152k,6,0,-2.98k,7,1.648852,-2.71152k,8,2.319721,-2.023237k,9,2.6,-0.15k,10,2.04,1.53k,1000,0,0!中心点!画隧道轮廓线larc,1,2,1000,2.55larc,2,3,1000,2.8larc,3,4,1000,6.4larc,4,5,1000,1.04969larc,5,6,1000,5.2larc,6,7,1000,5.2larc,7,8,1000,1.04969larc,8,9,1000,6.4larc,9,10,1000,2.8larc,10,1,1000,2.55!选择左边的所有线allsel!选择所有的实体lsel,s,loc,x,0,-5!选择左边的所有线:x从0到－5lcomb,all!合并所选的线lesize,all,,,30!把上面合并的线等分30段!选择右边的所有线allsellsel,s,loc,x,0,5lcomb,alllesize,all,,,30allsel!全选!划分单元lmesh,allnplot!显示节点!加弹簧!y方向*do,i,23,40PSPRNG,i,TRAN,vxishu,,-0.3,,,!tran-直线的*enddoPSPRNG,2,TRAN,vxishu,,-0.3,,,!固定结构用的!x方向*do,i,15,26PSPRNG,i,TRAN,uxishu,-0.3,,,,*enddo*do,i,37,48PSPRNG,i,TRAN,uxishu,0.3,,,,*enddoallselfinish!进入求解层，施加荷载，定义荷载步等/solu!施加约束nsel,s,,,2!选择约束的节点nsel,s,loc,x,0dnsel,r,loc,y,-2.98d,all,ux!施加水平方向的约束allselfcum,add,,!一定要，使荷载能叠加!竖向荷载*do,i,3,15a=-shuxiang*0.5*abs(nx(i)-nx(i+1))f,i+1,fy,af,i,fy,a*enddo*do,i,47,59a=-shuxiang*0.5*abs(nx(i)-nx(i+1))f,i+1,fy,af,i,fy,a*enddo!对未循环的节点施加竖向荷载a=-shuxiang*0.5*abs(nx(1)-nx(3))f,3,fy,af,1,fy,2*af,60,fy,a!水平荷载*do,i,3,30a=zuoce*0.5*abs(ny(i)-ny(i+1))f,i+1,fx,af,i,fx,a*enddo*do,i,32,59a=youce*0.5*abs(ny(i)-ny(i+1))f,i+1,fx,af,i,fx,a*enddof,1,fx,jizhongli!加集中力!对未循环的节点施加水平荷载a=youce*0.5*abs(ny(31)-ny(2))f,31,fx,-af,32,fx,a!施加重力acel,,9.8solve!求解finish!进入后处理，定义荷载工况并组合，输出图片和文本文件/post1!显示弯矩图etable,mi,smisc,6etable,mj,smisc,12plls,mi,mj,-1!显示轴力图etable,fi,smisc,1etable,fj,smisc,7plls,fi,fj,1!显示变形图pldisp,1