两个荷载结构法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,2300

et,1,beam3 r,1,0.5,0.5*0.5*0.5/12,0.5 mp,ex,1, CQETXML mp,dens,1, CQMD mp,prxy,1, CQUBSB

k,1,311.419293,-97.854439 k,2,308.756764,-98.346462 k,3,306.445920,-99.757543 k,4,304.948061,-101.812347 k,5,304.418078,-104.299297 k,6,304.452492,-105.061153 k,7,304.555456,-105.816802 k,8,304.958136,-106.652780 k,9,305.731611,-107.165369 k,10,308.539368,-107.845632 k,11,311.419293,-108.074297 k,12,314.299218,-107.845632 k,13,317.106975,-107.165369 k,14,317.880450,-106.652780 k,15,318.283130,-105.816802

k,16,318.386093,-105.061153 k,17,318.420508,-104.299297 k,18,317.890524,-101.812347 k,19,316.392666,-99.757543 k,20,314.081822,-98.346462

!二衬轴线 larc,1,3,2 larc,3,5,4 larc,5,7,6 larc,7,9,8 larc,9,11,10 larc,11,13,12 larc,13,15,14 larc,15,17,16 larc,17,19,18 larc,19,1,20

allsel

lsel,s,loc,x,311.419293,302.247076 lcomb,all nummrg,all lesize,all,,,50

allsel

lsel,s,loc,x,311.419293,320.122358 lcomb,all nummrg,all lesize,all,,,50 alls

lmesh,all eplot

nummrg,all,,,,low !所有号数从1开始 numcmp,node,eord !压缩节点号并排@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

!施加水平弹簧支撑单元 *do,i,16,50

序

a1=abs(ny(i-1)-ny(i+1))*WYTXKL kx=a1/2

PSPRNG,i,TRAN,kx,-0.3,,,, *enddo

*do,i,52,61

a1=abs(ny(i-1)-ny(i+1))*WYTXKL kx=a1/2

PSPRNG,i,TRAN,kx,0.3,,,, *enddo

*do,i,62,79

a1=abs(ny(i-1)-ny(i+1))*WYTXKLZGQ kx=a1/2

PSPRNG,i,TRAN,kx,0.3,,,, *enddo

*do,i,80,86

a1=abs(ny(i-1)-ny(i+1))*WYTXKL kx=a1/2

PSPRNG,i,TRAN,kx,0.3,,,, *enddo

!施加竖直弹簧支撑单元

*do,i,29,50

b1=abs(nx(i-1)-nx(i+1))*WYTXKL ky=b1/2

PSPRNG,i,TRAN,ky,,-0.3,,, *enddo

b1=abs(nx(52)-nx(50))*WYTXKL ky=b1/2

PSPRNG,51,TRAN,ky,,-0.3,,,

*do,i,52,61

b1=abs(nx(i-1)-nx(i+1))*WYTXKL ky=b1/2

PSPRNG,i,TRAN,ky,,-0.3,,, *enddo

*do,i,62,73

b1=abs(nx(i-1)-nx(i+1))*WYTXKLZGQ ky=b1/2

PSPRNG,i,TRAN,ky,,-0.3,,, *enddo

!施加边界条件与围岩压力

/solu

nsel,s,,,51 !选择约束的节点,仰拱中间节点 d,all,ux !施加水平方向的约束

allsel

fcum,add,,

!竖向荷载

*do,i,2,28 a=-Py*0.5*abs(nx(i-1)-nx(i+1)) f,i,fy,a *enddo

*do,i,80,99 a=-Py*0.5*abs(nx(i-1)-nx(i+1)) f,i,fy,a *enddo

a=-Py*0.5*abs(nx(1)-nx(99)) f,100,fy,a

a=-Py*0.5*abs(nx(2)-nx(100)) f,1,fy,a

!水平向荷载（X方向压力，全部节点施加水平荷载，换算成节点荷载，在节点上施加）

*do,i,52,61

a=-(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5 f,i,fX,a *enddo

*do,i,80,99

a=-(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5 f,i,fX,a *enddo

a=-(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(100)))*abs(ny(1)-ny(99))*0.5 f,100,fX,a

*do,i,2,50

a=(px1+(px2-px1)/(ny(51)-ny(1))*(ny(51)-ny(i)))*abs(nY(i-1)-nY(i+1))*0.5 f,i,fX,a *enddo

acel,9.8, solve finish

!进入后处理，定义荷载工况并组合，输出图片和文本文件 /post1

ESEL,s,TYPE,,1 !选择单元1（隧道轴力中心线单元） NSLE,s,ALL !选择所有节点

etable,mi,smisc,6

etable,mj,smisc,12 plls,mi,mj,-1

etable,fi,smisc,1 etable,fj,smisc,7 plls,fi,fj,1

pldisp,1

etable,refl

pretab, 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,beam3

r,1,0.3,0.00225,0.3

mp,ex,1,3.1e10 !弹性模量 mp,dens,1,2500 !密度 mp,prxy,1,0.167 !泊松比

!创建几何模型 画点 k,1,0,2.55 k,2,-2.04,1.53 k,3,-2.6,-0.15

k,4,-2.319721,-2.023237 k,5,-1.648852,-2.71152 k,6,0,-2.98

k,7,1.648852,-2.71152

k,8,2.319721,-2.023237 k,9,2.6,-0.15 k,10,2.04,1.53

k,1000,0,0 !中心点

!画隧道轮廓线 larc,1,2,1000,2.55 larc,2,3,1000,2.8 larc,3,4,1000,6.4 larc,4,5,1000,1.04969 larc,5,6,1000,5.2 larc,6,7,1000,5.2 larc,7,8,1000,1.04969 larc,8,9,1000,6.4 larc,9,10,1000,2.8 larc,10,1,1000,2.55

!选择左边的所有线

allsel !选择所有的实体

lsel,s,loc,x,0,-5 !选择左边的所有线:x从0到－5 lcomb,all !合并所选的线

lesize,all,,,30 !把上面合并的线等分30段

!选择右边的所有线 allsel

lsel,s,loc,x,0,5 lcomb,all lesize,all,,,30

allsel !全选

!划分单元 lmesh,all

nplot !显示节点

!加弹簧 !y方向 *do,i,23,40

PSPRNG,i,TRAN,vxishu,,-0.3,, , !tran-直线的 *enddo

PSPRNG,2,TRAN,vxishu,,-0.3,, , !固定结构用的

!x方向 *do,i,15,26

PSPRNG,i,TRAN,uxishu,-0.3,,, , *enddo

*do,i,37,48

PSPRNG,i,TRAN,uxishu,0.3,,, , *enddo

allsel finish

!进入求解层，施加荷载，定义荷载步等 /solu

!施加约束

nsel,s,,,2 !选择约束的节点 nsel,s,loc,x,0d nsel,r,loc,y,-2.98

d,all,ux !施加水平方向的约束

allsel

fcum,add,, !一定要，使荷载能叠加

!竖向荷载 *do,i,3,15

a=-shuxiang*0.5*abs(nx(i)-nx(i+1)) f,i+1,fy,a f,i,fy,a *enddo *do,i,47,59

a=-shuxiang*0.5*abs(nx(i)-nx(i+1)) f,i+1,fy,a f,i,fy,a *enddo

!对未循环的节点施加竖向荷载 a=-shuxiang*0.5*abs(nx(1)-nx(3)) f,3,fy,a f,1,fy,2*a f,60,fy,a

!水平荷载 *do,i,3,30

a=zuoce*0.5*abs(ny(i)-ny(i+1)) f,i+1,fx,a f,i,fx,a *enddo *do,i,32,59

a=youce*0.5*abs(ny(i)-ny(i+1)) f,i+1,fx,a f,i,fx,a *enddo

f,1,fx,jizhongli !加集中力

!对未循环的节点施加水平荷载 a=youce*0.5*abs(ny(31)-ny(2)) f,31,fx,-a f,32,fx,a

!施加重力 acel,,9.8

solve !求解 finish

!进入后处理，定义荷载工况并组合，输出图片和文本文件 /post1

!显示弯矩图

etable,mi,smisc,6 etable,mj,smisc,12 plls,mi,mj,-1

!显示轴力图 etable,fi,smisc,1 etable,fj,smisc,7 plls,fi,fj,1

!显示变形图 pldisp,1

本文来源：https://www.bwwdw.com/article/nbqh.html