[求助]车、桥在风的作用下的响应问题

绿荷雨丝

　　

我在Ansys中国上看到好像是长沙铁道学院做了一个用流固耦合的方法分析车、桥在风的作用下的响应问题，不知哪位高手能不吝赐教。

对于车桥共振问题我想能不能采用函数加载的方法来判断结构的响应，对于谐响应分析在这上面的应用请有教于我。不胜感谢
不知道用单点半幅的正弦波怎么样?

hbchen

车桥耦合振动问题本身就是时变问题如果考虑车惯性ansys可能不能模拟吧，只能简化为移动荷载，你去仿真论坛看看吧，那里有一些例子。

linda

　　看看下面的例子：
ANSYS一个火车过桥动态加载的例子

! 用对话框提示用户输入基本数据
!:reread
!multipro,'start',3
! *cset,1,3,y_bot,'enter the value of y_bot',0
! *cset,4,6,y_top,'enter the value of y_bot',0
! *cset,7,9,ttype,'enter the value of train type',0 !输入火车活载的类型

!multipro,'end'
y_bot=13.9
y_top=14.1
ttype = 0
*if,ttype,eq,0, then
q1 = 147000 !中-活载
q2 = 92000
q3 = 80000
L1 = 7.5
L2 = 30
*elseif,ttype,eq,1,then
q1 = 0
q2 = 0
q3 = 0
L1 = 0
L2 = 0
*else
*msg,ui
error value of ttype(1 or 2),please enter again
*go,: reread
*endif
! 识别桥面过列车的单元，并将单元按节点坐标从小到大排序
lsel,s,loc,y,y_bot,y_top
esll,s ! 选中桥面单元
elmax = elmiqr(0,14) ! 最大单元号
ndmax = ndinqr(0,14) ! 最大节点号
*dim,eldk,,elmax ! 存放桥面过车单元
*dim,ndx,,elmax ! 存放过车单元的节点x坐标(两节点坐标中的较小值)
*dim,nodes,,ndmax ! 存放桥面过车的节点
j=0
*do,i,1,elmax
*if,esel(i),gt,0, then
j = j+1
eldk(j) = i
*endif
*enddo ! 记录桥面单元的标号，保存于数组eldk中
eldknum = j ! 总桥面单元数
! 对eldk数组排序，按节点坐标从小到大排列
*do,i,1,eldknum
ii = eldk(i)
ndi = nelem(ii,1)
ndj = nelem(ii,2)
i1 = nx(ndi)
i2 = nx(ndj)
ndx(i) = i1
*if,i2,lt,i1,then
ndx(i) = i2
*endif
*enddo
*do,i,1,eldknum-1
*do,j,i+1,eldknum
*if,ndx(j),lt,ndx(i),then
temp = ndx(i)
ndx(i) = ndx(j)
ndx(j) = temp
temp = eldk(i)
eldk(i)= eldk(j)
eldk(j)= temp
*endif
*enddo
*enddo
nsle,s ! 选中桥面单元的节点
j=0
*do,i,1,ndmax
*if,nsel(i),gt,0,then
j=j+1
nodes(j)=i
*endif
*enddo
ndnum = j ! 总桥面节点数
xmin = nx(nelem(eldk(1),1))
temp = nx(nelem(eldk(1),2))
*if,temp,lt,xmin,then
xmin = temp
*endif
xmax = nx(nelem(eldk(eldknum),1))
temp = nx(nelem(eldk(eldknum),2))
*if,temp,gt,xmax,then
xmax = temp
*endif
qth = xmax -xmin ! 桥面过车长度
lth1 = 2*qth + L1+L2 +182.5
lth = (nint(lth1/5)+8)*5 ! 总长度
dim = lth/5 +1
*dim,disy,table,dim,ndnum ! 桥面节点的y方向的位移值
*dim,disy1,table,dim,ndnum ! 桥面节点的y方向的位移值(从大到小)
*dim,disx,table,dim,ndnum ! 桥面节点的x方向的位移值
*dim,disa,table,dim,ndnum ! 桥面节点的转角
*dim,disymax,,ndnum
*dim,diy,table,dim
*dim,ffx,table,elmax,dim ! 单元轴力
*dim,fmz,table,elmax,dim ! 单元弯矩
*dim,detaf,,elmax ! 每个单元轴力变化值
*dim,detaf1,,elmax ! 每个单元轴力变化值(从小到大)
*dim,elnum,,elmax ! 轴力变化最大的单元号
*dim,detaf2,,20 ! 前20个轴力变化最大值
*dim,elnum2,,20 ! 前20个轴力变化最大值的单元
*dim,fmax,,elmax ! 每个单元的轴力最大值
*dim,fmin,,elmax ! 每个单元的轴力最小值
*dim,fabs,,elmax ! 每个单元轴力最大值(绝对值)
*dim,elnum3,,elmax ! 轴力最大的单元号
*dim,elnum4,,20 ! 前20个轴力最大的单元号
*dim,nmax,,20
*dim,dfmaxel,table,dim,20 ! 前20各轴力变化最大单元的轴力
*dim,dmmaxel,table,dim,20 ! 前20各轴力变化最大单元的弯矩
*dim,fmaxel,table,dim,20 ! 前20各轴力最大单元的轴力
*dim,mmaxel,table,dim,20 ! 前20各轴力最大单元的弯矩
eplot ! plot element
allsel
time = 0
*do,kk,0,lth,5,
time = time+1
!/uis,msgpop,3 ! 不显示数组被重新定义的警告
!/uis,msgpop,2
/solu ! 进入求解器
antype,static
esel,all
sfedele,all,1,pres, !先删除上一次加上的所有单元荷载
! 计算列车各部分在桥上的单元号
tpos1 = kk ! 火车头所在的位置
*if,tpos1,lt,(qth+L1+L2), then !火车进桥情况
xlth1 = 0
*do,i,1,eldknum
ii = eldk(i)
ndi = nelem(ii,1)
ndj = nelem(ii,2)
elth1 = abs(nx(ndi)-nx(ndj))
xlth1 = xlth1 + elth1
*if,xlth1,ge,tpos1,then
*exit
*endif
*enddo
tele1 = i ! 火车头所在的单元
*if,tpos1,gt,L1,then
tpos2 =tpos1 -L1 !煤水车头所在的位置
xlth2 = 0
*do,i,1,eldknum
ii = eldk(i)
ndi = nelem(ii,1)
ndj = nelem(ii,2)
elth2 = abs(nx(ndi)-nx(ndj))
xlth2 = xlth2 + elth2
*if,xlth2,ge,tpos2,then
*exit
*endif
*enddo
tele2 = i !旅客列车头所在的单元
*if,tpos1,gt,(L1+L2),then
tpos3 =tpos1 -(L1+L2) !旅客车头所在的位置
xlth3 = 0
*do,i,1,eldknum
ii = eldk(i)
ndi = nelem(ii,1)
ndj = nelem(ii,2)
elth3 = abs(nx(ndi)-nx(ndj))
xlth3 = xlth3 + elth3
*if,xlth3,ge,tpos3,then
*exit
*endif
*enddo
tele3 = i !旅客列车头所在的单元
*endif
*endif
*endif
*if,tpos1,gt,(qth+L1+L2+182.5),then ! 火车过桥情况(182.5为任意取的一个常数)


　　

*if,tpos1,lt,lth1, then
tpos4 = tpos1-(qth+L1+L2+182.5) !旅客列车尾部所在的位置
xlth4 = 0
*do,i,1,eldknum
ii = eldk(i)
ndi = nelem(ii,1)
ndj = nelem(ii,2)
elth4 = abs(nx(ndi)-nx(ndj))
xlth4 = xlth4 + elth4
*if,xlth4,ge,tpos4,then
*exit
*endif
*enddo
tele4 = i !旅客列车尾部在的单元
*endif
*endif
! 加荷载于相应的单元
! SFBEAM, ELEM, LKEY, Lab, VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST
*if,tpos1,le,L1, then ! 火车头进桥
c1 = xlth1 -tpos1
c2 = elth1 - c1
*if,tele1,eq,1,then
qq = q1*c2/elth1
sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
*else
*do,j,1,tele1-1
sfbeam,eldk(j),1,pres,q1,q1,,,,,
*enddo
qq = q1*c2/elth1
sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
*endif

linda

　　*elseif,tpos1,le,(L1+L2),then ! 煤水车进桥
c1 = xlth1 - tpos1
*if,tele1,eq,tele2,then !火车头在同一个单元上
c2 = xlth2 - tpos2
c3 = elth2 - c2
qq = (q1*L1+q2*c3)/elth1
sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
*do,j,1,tele1-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
*else
dtele = tele1 - tele2
*if,dtele,eq,1, then
c2 = xlth2 - tpos2
c3 = elth2 - c2
qq0 = q1*(elth1-c1)/elth1
qq = (q1*c2+q2*c3)/elth2
sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
*do,j,1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
*else
c2 = xlth2 - tpos2
c3 = elth2 - c2
qq0 = q1*(elth1-c1)/elth1
sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
*do,j,tele2+1,tele1-1
sfbeam,eldk(j),1,pres,q1,q1,,,,,
*enddo
qq = (q1*c2+q2*c3)/elth2
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
*do,j,1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
*endif
*endif
*elseif,tpos1,le,qth,then ! 旅客列车进桥
c1 = xlth1 - tpos1
*if,tele1,eq,tele2, then
c2 = xlth2 - tpos2
c3 = elth2 - c2
qq = (q1*L1+q2*c3)/elth1
sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
c4 = xlth3 - tpos3
c5 = elth3 - c4
dtele1 = tele2 - tele3
*if,dtele1,eq,1,then
qq1 = (q2*c4+q3*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*else
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q3*c5)/(elth3)
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*else
dtele1 = tele1 -tele2
*if,dtele1,eq,1,then
c2 = xlth2 - tpos2
c3 = elth2 -c2
qq0 = q1*(elth1-c1)/elth1
sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
c4 = xlth3- tpos3
c5 = elth3 - c4
dtele2 = tele2 - tele3
*if,dtele2,eq,0,then
qq0 = (q1*c2+q2*L2+q3*c3)/elth2
sfbeam,eldk(tele2),1,pres,qq0,qq0,,,,,
*do,j,1,tele2-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*elseif,dtele2,eq,1,then
qq0 = (q1*c2+q2*c3)/elth2
sfbeam,eldk(tele2),1,pres,qq0,qq0,,,,,
qq1 = (q2*c4+c3*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*else
qq0 = (q1*c2+q2*c3)/elth2
sfbeam,eldk(tele2),1,pres,qq0,qq0,,,,,
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q3*c5)/(elth3)
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*else
c2 = xlth2 - tele2
c3 = elth2 - c2
qq0 = q1*(elth1-c1)/elth1
sfbeam,eldk(tele1),1,pres,qq0,qq0,,,,,
*do,j,tele2+1,tele1-1,
sfbeam,eldk(j),1,pres,q1,q1,,,,,
*enddo
qq = (q1*c2+q2*c3)/(elth2)
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
dtele2 = tele2 - tele3
*if,dtele2,eq,1,then
qq1 = (q2*c4+q3*c5)/(elth3)
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*else
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q3*c5)/(elth3)
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*endif
*endif

linda

　　*elseif,tpos1,le,(qth+L1),then !火车头出桥
c1 = tpos1 - qth !火车头出桥长度
dtele1 = eldknum - tele2
*if,dtele1,eq,0,then
c2 = xlth2 - tpos2
c3 = elth2 - c2
qq = (q1*c2+q2*c3)/(elth2)
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
dtele2 = tele2 - tele3
c4 = xlth3 -tpos3
c5 = elth3 - c4
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q3*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*elseif,dtele1,eq,1,then
sfbeam,eldk(eldknum),1,pres,q1,q1,,,,,
qq = (q1*c2+q2*c3)/(elth2)
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
dtele2 = tele2- tele3
c4 = xlth3 - tpos3
c5 = elth3 - c4
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q2*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*elseif,dtele1,gt,1,then
*do,j,tele2+1,eldknum
sfbeam,eldk(j),1,pres,q1,q1,,,,,
*enddo
qq = (q1*c2+q2*c3)/(elth2)
sfbeam,eldk(tele2),1,pres,qq,qq,,,,,
c4 = xlth3 - tpos3
c5 = elth3 - c4
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq1 = (q2*c4+q2*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq1,qq1,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*elseif,tpos1,le,(qth+L1+L2),then !煤水车头出桥
c4 = xlth3 - tpos3
c5 = elth3 - c4
dtele1 = eldknum - tele3
*if,dtele1,eq,0,then
qq = (q2*c4+q3*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq,qq,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*else
*do,j,tele3+1,eldknum
sfbeam,eldk(j),1,pres,q2,q2,,,,,
*enddo
qq = (q2*c4+q3*c5)/elth3
sfbeam,eldk(tele3),1,pres,qq,qq,,,,,
*do,j,1,tele3-1
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*endif
*elseif,tpos1,le,(qth+L1+L2+182.5),then ! 全桥均布列车荷载(182.5为任意取的
常数)
*do,j,1,eldknum
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
*elseif,tpos1,le,lth1,then ! 均布列车荷载出桥
c1 = xlth4 - tpos4
c2 = elth4 - c1
*do,j,tele4+1,eldknum
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*enddo
qq = q3*c1/elth4
sfbeam,eldk(tele4),1,pres,qq,qq,,,,,
*endif
/psf,pres,norm,1
/replot
solve
finish

linda

　　/post1
diy(time,0) = tpos1
diy(time) = uy(91) ! 得到中间节点的y坐标
*vplot,diy(1,0),diy(1,1) ! 动态显示中间节点位移结果
disx(time,0) = tpos1
disy(time,0) = tpos1
disa(time,0) = tpos1
*do,j,1,ndnum
disx(0,j) = nodes(j)
disy(0,j) = nodes(j)
disa(0,j) = nodes(j)
disx(time,j) = ux(nodes(j))
disy(time,j) = uy(nodes(j))
disa(time,j) = rotz(nodes(j))
*enddo
etable,forcex,smisc,1 ! 定义单元轴力表
etable,momentz,smisc,6 ! 定义单元弯矩表
*vget,ffx(1,time),elem,1,etab,forcex,,,,
*vget,fmz(1,time),elem,1,etab,momentz,,,,
ffx(0,time) = tpos1
fmz(0,time) = tpos1
*if,time,eq,1,then
*do,j,1,elmax
ffx(j,0) = j
fmz(j,0) = j
*enddo
*endif
*enddo
save,,,,all
! 结果处理
*do,i,1,ndnum
nd = nodes(i)
disymax(i) = abs(disy(1,nd))
*do,j,2,dim
*if,abs(disy(j,nd)),gt,disymax(i),then
disymax(i) = abs(disy(j,nd))
*endif
*enddo
*enddo
*do,i,1,ndnum-1
*do,j,i+1,ndnum
*if,disymax(j),gt,disymax(i),then
temp = disymax(i)
disymax(i) = disymax(j)
disymax(j) = temp
temp = nodes(i)
nodes(i) = nodes(j)
nodes(j) = temp
*endif
*enddo
*enddo
*do,i,1,ndnum ! 前5个位移最大的节点
numnd = nodes(i)
disy1(0,i) = numnd
*do,j,1,dim
kk = (j-1)*5
disy1(j,0) = kk
disy1(j,i) = disy(kk,numnd)
*enddo
*enddo
*do,i,1,elmax
fmin(i) = ffx(i,1)
fmax(i) = ffx(i,1)
*do,j,2,dim
*if,ffx(i,j),gt,fmax(i),then
fmax(i) = ffx(i,j)
*endif
*if,ffx(i,j),lt,fmin(i),then
fmin(i) = ffx(i,j)
*endif
*enddo
detaf(i) = fmax(i) - fmin(i) ! 每个单元的内力变化值
detaf1(i) = detaf(i)
*enddo
*do,i,1,elmax
elnum(i) = i
*enddo
*do,i,1,elmax-1
*do,j,i+1,elmax
*if,detaf1(j),gt,detaf1(i),then
temp = detaf1(i)
detaf1(i) = detaf1(j)
detaf1(j) = temp
temp = elnum(i)
elnum(i) = elnum(j)
elnum(j)=temp
*endif
*enddo
*enddo
*do,i,1,20
detaf2(i) = detaf1(i) !前20个内力变化最大的单元
elnum2(i) = elnum(i)
*enddo
*do,i,1,20 ! 前20个内力变化最大的单元
numel = elnum2(i)
dfmaxel(0,i) = i
dmmaxel(0,i) = i
*do,j,1,dim
kk = (j-1)*5
dfmaxel(j,0) = kk
dfmaxel(j,i) = ffx(numel,kk) ! 轴力变化最大单元的轴力
dmmaxel(j,0) = kk
dmmaxel(j,i) = fmz(numel,kk) ! 相应的弯矩
*enddo
*enddo
! 前20个内力最大的单元
*do,i,1,elmax
fabs(i) =abs(fmax(i))
*if,fabs(i),lt,abs(fmin(i)),then
fabs(i) =abs(fmin(i))
*endif
*enddo
*do,i,1,elmax
elnum3(i) =i
*enddo
*do,i,1,elmax-1
*do,j,i+1,elmax
*if,fabs(j),gt,fabs(i),then
temp = fabs(i)
fabs(i) = fabs(j)
fabs(j) =temp
temp = elnum3(i)
elnum3(i) = elnum3(j)
elnum3(j) = temp
*endif
*enddo
*enddo
*do,i,1,20
elnum4(i) = elnum3(i)
nmax(i) = fabs(i)
*enddo
*do,i,1,20 ! 前20个内力最大的单元
numel = elnum4(i)
fmaxel(0,i) = i
mmaxel(0,i) = i
*do,j,1,dim
kk = (j-1)*5
fmaxel(j,0) = kk
fmaxel(j,i) = ffx(numel,kk) ! 轴力最大单元的轴力
mmaxel(j,0) = kk
mmaxel(j,i) = fmz(numel,kk) ! 相应的弯矩
*enddo
*enddo
save,,,,all
!/title,the 91 node displament vs train position ! 显示图形的名称
!*vplot,diy(1,0),diy(1,1)
!/ui,copy,save,bmp,graph
!/title
!/replot
*do,i,1,30
cc = nodes(i)
/title, y displacement of node %cc% VS. train position
*vplot,disy1(1,0),disy1(1,i)
/ui,copy,save,bmp,graph
*enddo
/title
/sys,rename njq*.bmp aaa*.bmp
*do,i,31,ndnum
cc = nodes(i)
/title, y displacement of node %cc% VS. train position
*vplot,disy1(1,0),disy1(1,i)
/ui,copy,save,bmp,graph
*enddo
/title
/sys,rename njq*.bmp aa*.bmp
!*do,i,1,20
!cc = elnum2(i)
!/title,axis of element %cc% VS. train position
! *vplot,dfmaxel(1,0),dfmaxel(1,i)
!/ui,copy,save,bmp,graph
!/title,moment of element %cc% VS. train position
! *vplot,dmmaxel(1,0),dmmaxel(1,i)
!/ui,copy,save,bmp,graph
!*enddo
!/title
*do,i,1,20
cc = elnum4(i)
/title,axis of element %cc% vs train position
*vplot,fmaxel(1,0),fmaxel(1,i)
/ui,copy,save,bmp,graph
/title,moment of element %cc% vs train position
*vplot,mmaxel(1,0),mmaxel(1,i)
/ui,copy,save,bmp,graph
*enddo
/title
*cfopen,result,txt,,,
*vwrite
('' '')
*vwrite

linda

　　(''**********************南京桥计算结果*****************************'')
*vwrite
('' '')
*vwrite
('' ****节点最大位移**** '')
ndnum1 = 5*nint(ndnum/5)
*do,i,1,ndnum1,5
a=nodes(i)
b=nodes(i+1)
c=nodes(i+2)
d=nodes(i+3)
e=nodes(i+4)
*vwrite,a,b,c,d,e
(''荷载位置'',5(f8.0,5x))
*vwrite,disy1(1,0),disy1(1,i),disy1(1,i+1),disy1(1,i+2),disy1(1,i+3),disy1(
1,i+4)
(f5.0,,5(3x,f10.5))
*enddo
*vwrite
('' ****前20个轴力变化最大的单元和相应的轴力值**** '')
*vwrite,elnum2(1),detaf2(1)
(f15.0,'' '', f15.1)
*vwrite
('' '')
*do,i,1,18,3
*vwrite
('' '')
a= elnum2(i)
b= elnum2(i+1)
c= elnum2(i+2)
*vwrite,a,b,c
(''荷载位置'',9x,''单元'',f4.0,14x,''单元'',f4.0,14x,''单元'',f4.0)
*vwrite,dfmaxel(1,0),dfmaxel(1,i),dmmaxel(1,i),dfmaxel(1,i+1),dmmaxel(1,i+1)
,dfmaxel(1,i+2),dmmaxel(1,i+2)
(f5.0,3(f11.1,'' '',f11.1))
*enddo
*vwrite
('' '')
a= elnum2(19)
b= elnum2(20)
c= elnum2(20)
*vwrite,a,b,c
(''荷载位置'',9x,''单元'',f4.0,14x,''单元'',f4.0,14x,''单元'',f4.0)
*vwrite,dfmaxel(1,0),dfmaxel(1,19),dmmaxel(1,19),dfmaxel(1,20),dmmaxel(1,20)
,dfmaxel(1,20),dmmaxel(1,20)
(f5.0,3(f11.1,'' '',f11.1))
*vwrite
('' '')
*vwrite
('' '')
*vwrite
('' ***前20个轴力最大的单元和相应的弯矩值*** '')
*vwrite,elnum4(1),detaf2(1)
(f15.0,'' '', f15.1)
*do,i,1,18,3
*vwrite
('' '')
a= elnum4(i)
b= elnum4(i+1)
c= elnum4(i+2)
*vwrite,a,b,c
(''荷载位置'',9x,''单元'',f4.0,14x,''单元'',f4.0,14x,''单元'',f4.0)
*vwrite,fmaxel(1,0),fmaxel(1,i),mmaxel(1,i),fmaxel(1,i+1),mmaxel(1,i+1),fmax
el(1,i+2),mmaxel(1,i+2)
(f5.0,3(f11.1,'' '',f11.1))
*vwrite
('' '')
*enddo
*vwrite
('' '')
a= elnum4(19)
b= elnum4(20)
c= elnum4(20)
*vwrite,a,b,c
(''荷载位置'',9x,''单元'',f4.0,14x,''单元'',f4.0,14x,''单元'',f4.0)
*vwrite,fmaxel(1,0),fmaxel(1,19),mmaxel(1,19),fmaxel(1,20),mmaxel(1,20),fmax
el(1,20),mmaxel(1,20)
(f5.0,3(f11.1,'' '',f11.1))
*cfclos
补充说明:此贴来源于水木,原帖没有完全版，只是提供了与用户互动的APDL怎么操作(如用户可以随意改变参数)，以及过程是什么样的,建模需自己完成。