影响线的求法

fini
/clear
/title,yingxiangxian
/prep7
et,1,beam4

mp,ex,1,3e7
mp,dens,1,2.6
mp,prxy,1,0.3

r,1,1,1,1,1

n,1
n,127,12.6
fill,1,127

*do,i,1,126,1
e,i,i+1
*enddo

d,1,uy,,,,,uz,ux
d,127,uy,,,,,uz
fini

/solu
antype,trans
tim=0
delt=1
timint,on
*do,i,1,62,1.000
tim=tim+1
time,tim
fdele,all,all
allsel
nsel,s,loc,x,(i+64)*0.1

f,all,fy,-250

allsel
nsel,s,loc,x,(i+52)*0.1
f,all,fy,-250

allsel
nsel,s,loc,x,(i+12)*0.1
f,all,fy,-250

allsel
nsel,s,loc,x,i*0.1
f,all,fy,-250

allsel
nsubst,1
outres,all,all

solve
*enddo
fini

《Aogor、Ansys在桥梁工程中应用方法与实例》这本书是本非常经典的书，网上也流传了很多其命令流，可惜这个命令流还是不全，没有后处理部分！

上面的例子只是求得最不利荷载位置和该位置的值的大小，本人把书上地的例子修改了下，能够求左支点反力影响线
fini
/clear
/title,yingxiangxian
/prep7
et,1,beam4

mp,ex,1,3e7
mp,dens,1,2.6
mp,prxy,1,0.3

r,1,1,1,1,1

n,1
n,127,12.6
fill,1,127

*do,i,1,126,1
e,i,i+1
*enddo

d,1,uy,,,,,uz,ux
d,127,uy,,,,,uz
fini

/solu
antype,trans
tim=0
delt=1
timint,on
*do,i,1,62,1.000
tim=tim+1
time,tim
fdele,all,all
allsel
nsel,s,loc,x,(i+64)*0.1

f,all,fy,-250

allsel
nsel,s,loc,x,(i+52)*0.1
f,all,fy,-250

allsel
nsel,s,loc,x,(i+12)*0.1
f,all,fy,-250

allsel
nsel,s,loc,x,i*0.1
f,all,fy,-250

allsel
nsubst,1
outres,all,all

solve
*enddo
fini

/post26
/axlab,y,fy
esol,2,1,1,f,y,fy32
plvar,2,

[ 本帖最后由 水幽寒 于 2006-10-24 14:03 编辑 ]

　　对esol,2,1,1,f,y,fy32 命令的请教！请解释
一下好吗

ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 将结果存入变量
        NVAR: 变量号，2以上
ELEM: 单元号
NODE: 该单元的节点号，决定存储该单元的哪个量，如果空，则给出平均值
ITEM:
COMP:
NAME: 8字符的变量名， 缺省为ITEM加COMP

刚刚做的一个墩顶弯矩影响线

梁-----墩顶断面影响线

弯矩影响线就是书上的影响线的命令流只改一下后面的命令流吧，改一下变量　　　比如:  
/post26
/axlab,y,mz
nsol,2,96,u,y,wanju96
plvar,2

看了看出不来啊，怎么弄的啊你　

那本书上的命令流有些问题:lol
nsol,2,96,u,y,wanju96
提取的是位移,想画弯矩的影响线,应该改为
/post26
/axlab,y,mz32
esol,2,181,62,m,z,mz32
plvar,2
再尝试下 :D

esol,2,181,62,m,z,mz32
是提取的１８１单元的６２节点，应该不正确吧，是不是应该是
esol,2,64,64,m,z,mz64吧
另一个问题是好象你的时间都改变了
能不能把你的完整的命令流贴出来看一看啊

我这只是举个例子而已,具体的节点\单元编号,还要自己取,再有上面的图形并不是我上面的命令流做出来的,那是我最近刚刚做的一座连续梁桥的影响线.

好东西，这些东西很实用