[转载]ansys例子
deltim,0.1,0.05,0.2
autots,on !使用自动时间步
pred,on !打开时间步长预测器
lnsrch,on !打开线性搜索
nlgeom,on !打开大位移效果
nropt,full !设定牛顿-拉普森选项
cnvtol,f,,0.02,023.滞回分析
!定义单元类型，实常数，材料特性
ET,1,SHELL143
R,1,12, , , , ,
MP,EX,1,196784
MP,NUXY,1,0.3
!双线性随动强化模型
TB,BKIN,1,1,2,1
TBDATA,,310,600,,,,
!定义关键点、线、面
K,1,54,0,0
K,2,-54,0,0
K,3,54,0,1000
K,4,-54,0,1000
!定义边界荷强迫位移,划分网格
AESIZE,ALL,27,
MSHAPE,0,2D
CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CMSEL,S,_Y
OUTPR,BASIC,ALL,
OUTRES,ALL,ALL,
!第1荷载步
D,46,ux,60
NSUBST,10, , ,1
LSWRITE,01,
!第2荷载步
D,46,ux,-60
NSUBST,20, , ,1
LSWRITE,02,
!第3荷载步
D,46,ux,60
NSUBST,20, , ,1
LSWRITE,03,
!第4荷载步
D,46,ux,-60
NSUBST,20, , ,1
LSWRITE,04,
LSSOLVE,1,4,1,
!画出荷载位移曲线
NSOL,2,46,U,X,
RFORCE,3,46,F,X,
PLVAR,3, , , , , , , , , ,
1．ANSYS SOLID65环向布置钢筋的例子
! 一个管道，环向配筋率为1%，纵向配筋率为0.5%，径向配筋率为0.1%
! 单元属性
ET,1,SOLID65
KEYOPT,1,1,0
KEYOPT,1,5,0
KEYOPT,1,6,0
KEYOPT,1,7,1
!实参数1：不同方向配筋
R,1,2,.001, , ,2, .01,
RMORE, 90, ,2,.005 ,90 ,90 ,
!Adds real constants to a set.
!混凝土基本材料属性
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,30E2
MPDATA,PRXY,1,,.2
TB,MISO,1,1,5,
TBPT,,0.0005,15
TBPT,,0.001,21
TBPT,,0.0015,24
TBPT,,0.002,27
TBPT,,0.003,24
TB,CONC,1,1,9,
TBDATA,,.5,.9,3,30,,
!抗拉和抗压是10 倍的关系
TBDATA,,,,1,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
!钢材基本属性
MPDATA,EX,2,,200E3
MPDATA,PRXY,2,,.27
TB,BISO,2,1,2,
TBDATA,,310,2E3,,,,
! 管道内径和外径
CYL4,0,0,3000, , , ,10000
CYL4,0,0,2000, , , ,10000
VSBV, 1, 2
! 定义局部柱坐标
CSWPLA,11,1,1,1,
! Defines a local coordinate system at the origin of the working
KWPAVE, 11
!move to keypoint11
WPRO,,-90.000000,
WPCSYS,-1,0
!Defines the working plane location based on a coordinate
WPRO,,,-90.000000
ESIZE,500,0,
! 注意：设定单元局部坐标
VATT, 1, 1, 1, 11
VSWEEP,ALL
/DEVICE,VECTOR,1
/ESHAPE,1.0
!Displays elements with shapes determined from the real constants
or section
! 注意：红色代表最大配筋方向，绿色代表其次，蓝色表示最小配筋方向
11.四层弹簧－质点模型的地震分析
包括频率分析，谱分析和时程分析
作者：陆新征，清华大学土木工程系
DAMPRATIO=0.02 !振型阻尼系数都为0.02
ET,1,BEAM188
KEYOPT,1,7,1
KEYOPT,1,8,1
ET,2,MASS21,,,4
MP,EX,1,210E9
MP,NUXY,1,0.3
MP,DENS,1,7800
MPDATA,DAMP,1,,DAMPRATIO
SECTYPE,3,BEAM,I
SECDATA,0.2,0.2,0.5,0.01,0.01,0.008
R,1,1.6E2,1.6E2
R,2,1.2E2,1.2E2
K,5,0,12,0
K,6,0,0,10000
LSEL,,,,1,4
LATT,1,0,1,,6,,3
LESIZE,ALL, , ,3 , , , , ,1
NUMMRG,ALL
NUMCMP,ALL
*SET,NT,1001
*SET,DT,0.02
*DIM,AC,,NT
*VREAD,AC(1),RECORD,TXT
/SOLU !模态分析
MODOPT,SUBSP,8
MXPAND,8, , ,1
!得到自振频率1
*GET,FREQ1,MODE,1,FREQ
/ESHAPE,1.0
ANTYPE,SPECTR !谱分析
SPOPT,SPRS
MDAMP,1,DAMPRATIO, , , , , ,
SVTYPE,2,1
FREQ,0.167,0.25,0.333,0.,0.667,1,1.25,1.667
SV,,0.........
FREQ,2.222,10,100000
SV,,1...3528
!按7度多遇地震，地震影响系数0.08 第一组III场地TG＝0.45
SRSS,0.01,DISP
/INP,,MCOM
NSORT,U,SUM
*GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形
!开始地震计算
/CONFIG,NRES,20000
ANTYPE,TRANS
TRNOPT,FULL
ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926
BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)
*DO,I,1,1001
ACEL,AC(I),0,0
TIME,I*0.02
OUTRES,ALL,ALL
12.悬臂梁地震分析
!定义单元类型
ET,1,BEAM3
R,1,0.09,0..3, !0.3*0.3m截面
MP,EX,1,2.06e11
MP,NUXY,1,0.3
MP,DAMP,1,0.02 !阻尼比0.02
MP,dens,1,7.85e3
!--------------------------------------------------
!定义节点单元和位移约束
*do,i,1,11
n,i,(i-1)*0.5,0,0 !5m长,分10段
*do,i,1,10
!--------------------------------------------------
!定义和读入时程曲线
NT=1000 !时程曲线有NT个点
DT=0.01 !时间间隔
*dim,ac,,NT
!读入数据,这个数据文件可以用excel等软件来写
/input,tianjin,txt !天津波东西向共10秒
!--------------------------------------------------
NSUBST,1, , ,1 !1个子步
OUTRES,ALL,1 !输出每个子步的结果
ANTYPE,TRANS !时程分析
*do,i,1,NT
ACEL,0,ac ,0
!--------------------------------------------------
NSOL,2,11,U,Y,
PLVAR,2,3, , ,
下面是文件tianjin.txt，这个文件可以用excel等软件制作或编写很短的程序就能写出来先
贴前500个点
ac( 1)= -0.
ac( 2)= -0.
16.简支梁实体与预应力钢筋分析实例
/COM, Structural
egjx=2e5 !Ey
agjx=140 !单根钢绞线面积
ehnt=4e4 !Eh
xzxs=1.0e-5 !线胀系数
yjl=200000 !定义预加力
et,1,link8 !定义link8单元
et,2,solid95 !定义solid95单元
r,1,agjx !定义link8单元的面积
r,2 !定义第2种实常数
mp,ex,1,egjx !定义link8单元的弹性模量
mp,prxy,1,0.3 !定义link8单元的泊松系数
mp,alpx,1,1.0e-5 !定义线膨胀系数
mp,ex,2,ehnt !定义solid95单元的弹性模量
mp,prxy,2,0.3 !定义solid95单元的泊松系数
blc4, , ,100,200,3000 !定义梁体
/view,1,1,1,1 !定义ISO查看
vplot !绘制梁体
kwpave,6 !工作平面移动到关键点6
wpoff,-30 !工作平面移动-30mm(X)
wprot,0,0,90 !工作平面旋转
vsbw,1 !分割梁体
wpoff,0,0,-40 !工作平面移动-40mm(Z)
vsbw,2 !分割梁体
wpoff,0,40 !工作平面移动40mm(Y)
wprot,0,90 !工作平面旋转
vsbw,all !分割梁体
wpstyl !关闭工作平面显示
nummrg,all,,,,low !整理
numcmp,all !压缩编号
esize,30 !定义网分时边长控制
lsel,s,,,28,38,10 !定义line28和38为新的选择集
latt,1,1,1 !定义选择集的属性
lmesh,all !对线划分单元
allsel,all !新的选择集为所有的实体
gplot !绘制所有的实体
vsel,s,,,all !定义所有体为选择集
vatt,2,2,2 !定义选择集的属性
mshape,0,3d !将体划分单元的形状定位HEX
mshkey,1 !采用MAPPED划分器
vmesh,all !对体进行划分单元
dl,3,,all !对线line7施加约束(UX,UY,UZ)
dl,16,,all !对线line31施加约束(UX,UY,UZ)
dl,23,,all !对线line23施加约束(UX,UY,UZ)
dl,2,,uy !对线line4施加约束(UY)
dl,15,,uy !对线line30施加约束(UY)
dl,22,,uy !对线line23施加约束(UY)
dk,2,,,,,ux,uy !对关键点2约束(UX,UY)
bfl,28,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度
bfl,38,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度
solve !求解
plnsol,s,z,0,1 !绘制Z方向的应力
etable,sigi,ls,1 !定义钢筋单元数据表
plls,sigi,sigi,1 !绘制上述应力
!/exit,nosav
20.钢筋混凝土管的动力响应特性分析实例
在管内壁激振后，采用瞬态模块，模拟半埋设钢筋混凝土管在内部震源作用下的动力特性。
请多提宝贵意见，如果给分就太好了！
/COM,ANSYS MEDIA REL. 7.1 (03-13-2003) REF. VERIF. MANUAL: REL.
/CONFIG,NRES,2000
/VERIFY,ACOUSTIC_WELL_RADIUS
/TITLE,ACOUSTIC, ACOUSTIC_WELL_RADIUS
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!! 按平面问题分析
!! 钢板＋钢筋混凝土＋砂浆/泥浆＋原位岩/土（瞬态）
!! 6种介质成分，存放在介质属性数组中。
!! 计算边界外延一定宽度，以抑制边界干扰*****************
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
JH_HSGD=10& !计算高度
JH_GBHD=0.04& !钢板厚度
JH_HNTHD=0.7&& !混凝土厚度
JH_JYHD=0.6& !浆液厚度，高度方向上砂浆和泥浆各占一半，砂浆在下，泥浆在上
JH_YTTHD=0.6& !岩土体厚度
JH_RADWELL=4
ZHD=JH_RADWELL+JH_GBHD+JH_HNTHD+JH_JYHD+JH_YTTHD
BJCR1=ZHD+0.5
BJCR2=ZHD+1.0
BJCR3=ZHD+1.5
BJCR4=ZHD+2.0
SCBCJG=250
ENDTIME=0.001 !总计算时间
DST=1E-06 !记录的时间步长
ZYPL=5E4 !震源频率
PI=2*ASIN(1)
TCCS=6 !介质分类数目
!钢&& 钢筋混凝土 泥浆 砂浆 原位岩/土
*DIM,TC,ARRAY,TCCS,5
!该数组为介质参数，第一列为介质密度&&&
第二列为泊松比，
!&& 第三列为弹模
! 第四列为纵波波速，用于计算时间步长
! 第五列为横波波速，如为液体，横波波速为0
!土层自下而上排列，物理参数暂时按各土层相同考虑
TC(1,1)=00,00
TC(1,2)=0.285,0.32,0.48,0.35,0.35,0.35
!液体泊松比为0.5,计算出现错误，按提示改为0.499
TC(1,3)=210E9,36E9,2E6,6E9,20E9,20E9
!计算各种介质的P、S波速度,以计算时间步长
*DO,I,1,TCCS
LAME1=TC(I,3)*TC(I,2)/(1+TC(I,2))/(1-2*TC(I,2))
& LAME2=TC(I,3)/2/(1+TC(I,2))
TC(I,4)=((LAME1+2*LAME2)/TC(I,1))**.5
TC(I,5)=(LAME2/TC(I,1))**.5 !《土中应力波及其量测》P26
!#########生成各介质区的边界控制点############
JZFQ=4&&&&&&
!介质分区数目
R1=JH_RADWELL
R2=JH_RADWELL+JH_GBHD
R3=JH_RADWELL+JH_GBHD+JH_HNTHD
R4=JH_RADWELL+JH_GBHD+JH_HNTHD+JH_JYHD
R5=JH_RADWELL+JH_GBHD+JH_HNTHD+JH_JYHD+JH_YTTHD
*DIM,JHRB,ARRAY,JZFQ+BJCSM+1,1
JHRB(1,1)=R1,R2,R3,R4,R5,R6,R7,R8,R9!
*DIM,JHRA,ARRAY,JZFQ+BJCSM+1,4&&&&
!定义数组，用于存放
*DIM,AMAT,ARRAY,JZFQ+BJCSM+1,1&
*DIM,DYCC,ARRAY,JZFQ+BJCSM+1,1&&&
JHRA(1,1)=R1,R2,R3,R3,R4,R5,R6,R7,R8!
JHRA(1,2)=R2,R3,R4,R4,R5,R6,R7,R8,R9!
JHRA(1,3)=0,0,0,180,0,0,0,0,0
JHRA(1,4)=360,360,180,360,360,360,360,360,360
AMAT(1,1)=1,2,3,4,5,5,5,5,5&&&
!介质编号&&&&&
DYCC(1,1)=0.2,0.2,0.2,0.2,0.2,0.3,0.45,0.6,1!剖分尺寸
/REPLOT,RESIZE
*DO,I,1,JZFQ+BJCSM+1
CYL4,0,0,JHRA(I,1),JHRA(I,3),JHRA(I,2),JHRA(I,4)
NUMMRG,ALL
NUMCMP,ALL !压缩各定义项的编号,使之连续
!############计算时间步长##################
VMAX=TC(1,4) !钢板的纵波速度控制时间步长
DT=DYCC(1)/VMAX/2
*IF,DT,GE,DST,THEN !若时间步长大于设定的记录时间步长，采用记录时间步长记录
ZBC=ENDTIME/DT !时间总步数
!############震源处理##########
T0=1/ZYPL/2
!震源信号持续时间，采用正弦信号激震，也可采用RICKER信号，但要从数组文件中读取，可考虑设定一FLAG控制
XHCXSJ=NINT(1/ZYPL/DT)+1 !信号持续的时间步长
T0I=NINT(XHCXSJ/2)
!经MATLAB实际输出后获得的信号实际频率
!######信号记录步长处理#########
*IF,DST/DT,LT,NINT(DST/DT),THEN
!计算信号记录的时间步长间隔，实际记录的时间步长与设定值可能不同
& RECORDSTEP=NINT(DST/DT)-1
& RECORDSTEP=NINT(DST/DT)
!#################时间步计算结束#########################
!##################单元剖分################
ET,1,PLANE183
*DO,I,1,JZFQ+BJCSM+1
& LOCAL, 11, 1, 0, 0&
!定义局部坐标系，柱体
CSYS,11&&&
LSEL,S,LOC,X,JHRB(I),&&&
& LESIZE,ALL,,R1ANGLE
& ALLSEL,ALL
AREFINE,1,9,,1!精细剖分
AREFINE,1,4,,1!精细剖分
!#################赋予材料属性#################
*DO,I,1,TCCS
& MP,DENS,I,TC(I,1)
& MP,NUXY,I,TC(I,2)
& MP,EX,I,TC(I,3)
LOCAL, 11, 1, 0, 0& !定义局部坐标系，柱体
*DO,I,1,JZFQ+BJCSM+1
& NSEL,S,LOC,X,JHRA(I,1),JHRA(I,2)
& NSEL,R,LOC,Y,JHRA(I,3),JHRA(I,4)
ESLN,S&&&&&&&&
!边界节点单元
& MAT,AMAT(I)
EMODIFY,ALL&&&&&&&
!确认材料参数修改
ALLSEL,ALL&&&&&&&&
!重新选择ALL，准备下一步操作
/PNUM,MAT,1
LOCAL, 11, 1, 0, 0& !定义局部坐标系，柱体
CSYS,11&&&
LSEL,S,LOC,X,R1
*GET, R1NUM, NODE, 0, COUNT
*GET, R1NMIN, NODE, 0,NUM, MIN& !获取最小节点编号
*DIM,R1NODE,ARRAY,R1NUM,1& !定义内圈节点数组
*DIM,R1NODEXY,ARRAY,R1NUM,2& !定义内圈节点坐标数组，存储极坐标
R1NODE(1,1)=R1NMIN&&
!赋内圈节点编号的最小值
*GET,RNX, NODE, R1NODE(1,1), LOC, X&
!获取内圈的节点的坐标，
*GET,RNY, NODE, R1NODE(1,1), LOC, Y& !在此为极坐标
R1NODEXY(1,1)=RNX
R1NODEXY(1,2)=RNY
*DO,I,2,R1NUM
R1NODE(I,1)=NDNEXT(R1NODE(I-1))&
!对内圈节点赋值，NDNEXT的输出为选定的节点中编号大于且最接近该节点的节点编号
& *GET,RNX, NODE, R1NODE(I,1), LOC,
X& !获取内圈的节点的坐标，
& *GET,RNY, NODE, R1NODE(I,1), LOC,
Y& !在此为极坐标
& R1NODEXY(I,1)=RNX
& R1NODEXY(I,2)=RNY&
!###############按极坐标中角度，由小到大排列节点编号和对应坐标
*DO,I,1,R1NUM-1
*DO,J,I+1,R1NUM
& *IF,R1NODEXY(I,2),GT,R1NODEXY(J,2),THEN
RY=R1NODEXY(I,2)
R1NODEXY(I,2)=R1NODEXY(J,2)
R1NODEXY(J,2)=RY
RX=R1NODEXY(I,1)
R1NODEXY(I,1)=R1NODEXY(J,1)
R1NODEXY(J,1)=RX
RNODE=R1NODE(I)
R1NODE(I)=R1NODE(J)
R1NODE(J)=RNODE
*ENDIF&&&&&&&&&
!####################排序OVER#####################
*CFOPEN,R1NODEXY,DAT!输出井壁内表面节点编号和坐标
*DO,K,1,R1NUM
& NUM_NODE=R1NODE(K)
& RX=R1NODEXY(K,1)
& RY=R1NODEXY(K,2)
& *VWRITE,NUM_NODE,RX,RY
(F10.0,F13.2,F13.2)
ALLSEL,ALL
!#################求解过程##########################
!##与竖直方向采用直角坐标的不同：
!##1、震源位置变换方式
!##2、记录点位置
!##3、边界条件设定，不作专门的吸收边界，而以空间扩大来处理
*CFOPEN,USUM,DAT
!采用REIKER子波，在初始时段施加
NSUBST,1, , ,1 !1个子步
ANTYPE,TRANS !时程分析
OUTPR,BASIC,SCBCJG, !每隔一定时间步长输出基本值
NSEL,S,LOC,X,R9
D,ALL,UX,0
D,ALL,UY,0
ALLSEL,ALL
*DO,I,1,ZBC
*IF,I,GT,1,AND,I,LE,XHCXSJ,THEN&
!在信号持续时间段延续脉冲子波
*DO,J,1,R1NUM
*GET,RNX, NODE,R1NODE(J,1), LOC, X
*GET,RNY, NODE,R1NODE(J,1), LOC, Y
NN=R1NODE(J,1)
PIFT=PI*ZYPL*(I-T0I)*DT
PIFT=PIFT**2
AC=-(1-2*PIFT)*EXP(-PIFT)& !RECKER子波，何兵寿，P48
F,NN,FX,AC*RNX/R1&&
!施加X方向作用力
F,NN,FY,AC*RNY/R1& !施加Y方向作用力
&&& *ENDDO
& *IF,I,EQ,XHCXSJ+1,THEN&
!在XHCXSJ+1时刻，删除作用在节点上的位移&
FDELE,ALL& !在信号持续时间结束后，删除所有作用力
& TIME,I*DT
*DO,K,1,R1NUM
& NUM_NODE=R1NODE(K)
& RX=R1NODEXY(K,1)
& RY=R1NODEXY(K,2)
& *GET, USUMR1, NODE, NUM_NODE,U,SUM
!输出节点编号，极坐标R，角度，USUM
*VWRITE,I,K,NUM_NODE,USUMR1
(F10.0,F10.0,F10.0,E13.3)
BNODEN=NODE(0,R1,0)
BNODES=NODE(0,-R1,0)
FILE,'file','rst','.'
NSOL,2,BNODES,U,X,uxs
NSOL,4,BNODES,U,Y,uys
NSOL,3,BNODEN,U,X,uxn
NSOL,5,BNODEN,U,Y,uyn
PLVAR,uxs,uxn
/IMAGE,SAVE,UXS_N,JPG
PLVAR,uys,uyn, , , , , , ,
/IMAGE,SAVE,UYS_N,JPG
LOCAL, 11, 1, 0, 0& !定义局部坐标系，柱体
CSYS,11&&&
NSEL,S,LOC,X,R1,R5
PLNS,U,SUM
ANDATA,0.5, ,2,1,zbc,10,0,1
22.螺栓连接的模拟实现问题
!ANSYS6.1计算test3(对称模型，只建一半模型,试验前分析）
!/input,menust,tmp&&&&
,'',,,,,,,,,,,,,,,,1&
!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!参数定义!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!&
hc=400 ! 柱截面高度
bc=200 ! 柱截面宽度
tcf=10 ! 柱翼缘厚度
tcw=8& ! 柱腹板厚度
lc=1100 ! 柱构件伸出长度
tep=20 !端板厚度
bep=bc+20& !端板宽度
HEP1= hc+200&& !端板高度
!端板外伸部分加劲肋厚度
!端板外伸部分加劲肋高度
bst=bc/2-5&& !端板外伸部分加劲肋高度
lbt=2*tep&&&&
!螺栓杆长度
dbt=20&&&&&
!螺栓杆直径或有效直径
dbth=31.4&&&&
!螺栓头和螺母直径
lbth=12.5&&&
!螺栓头厚度
preten=155000&&
!螺栓施工预拉力&
miu=.4&& !端板间抗滑移系数
hb=400&&& !
梁截面高度&&&
bb=200&&& !
梁截面宽度&&&
tbf=8&&&&&
! 梁翼缘厚度&
! 梁腹板厚度&
lb1=870&&& !
梁构件伸出长度&
lb2=200&&&&
! 梁构件伸出的水平加载断长度
hb1=362&&&&&&
! 梁最左端高度
dh0=dbt+2&&&&
!螺栓孔直径
randa=0.05&&&&&
aa=50&&&&&
!螺栓中心到梁翼缘边缘（非受力方向）的距离
aa1=50&&&&&
!螺栓中心到梁翼缘边缘（受力方向的距离
ab=120&&&&&
!一二排螺栓间距
displa=-50&&&
!施加的位移荷载大小&&&
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
!!!!!!!!!!!建模开始!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!梁模型!!!!!!!!!!
block,o,tep, -hep1/2,hep1/2,0,bep/2
block,-tep,0,-hep1/2-100+4*tcf,hep1/2,0,bep/2
k,20,tep,hb/2
k,21,tep,hb/2-tbf
k,22,tep+lb1,hb/2-tbf+lb1*randa
k,23,tep+lb1,hb/2+lb1*randa
k,24,tep,hb/2,bb/2
k,25,tep,hb/2-tbf,bb/2
k,26,tep+lb1,hb/2-tbf+lb1*randa,bb/2
k,27,tep+lb1,hb/2+lb1*randa,bb/2
v,20,21,22,23,24,25,26,27
k,28,tep,hb/2-tbf
k,29,tep,-hb/2+tbf
k,30,tep+lb1,-hb/2+tbf+lb1*randa
k,31,tep+lb1,hb/2-tbf+lb1*randa
k,32,tep,hb/2-tbf,tbw/2&&&&&&&&&&&&&&&&&
k,33,tep,-hb/2+tbf
,tbw/2&&&&&&&&&&&&&&&
k,34,tep+lb1,-hb/2+tbf+lb1*randa ,tbw/2&
k,35,tep+lb1,hb/2-tbf+lb1*randa,tbw/2&&&
v,28,29,30,31,32,33,34,35
block,tep+lb1,tep+lb1+lb2,hb/2-tbf+lb1*randa,hb/2+lb1*randa,0,bb/2
block,tep+lb1,tep+lb1+lb2,-hb/2+tbf+lb1*randa,hb/2-tbf+lb1*randa,0,tbw/2
vgen,2,3,5,2,,-hb+tbf,, ,0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!颜色区分volumn!!!!!!!
/PNUM,KP,0&
/PNUM,LINE,0
/PNUM,AREA,0
/PNUM,VOLU,1
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,1&
/PNUM,ELEM,0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!柱模型!!!!!!!!!!
!!!!!!!端板对接斜坡!!!!!!!!!!!!!!!!!!!!
!!!!!!端板对接斜坡!!!!!!!!!!!!!!!!!!!!
k,80,-tep,-hep1/2-100+4*tcf
k,81,-tep,-hep1/2-100
k,82,-tep+tcf, -hep1/2-100
k,83,0,-hep1/2-100+4*tcf&
k,84,-tep,-hep1/2-100+4*tcf,bep/2&&&&&
k,85,-tep,-hep1/2-100,bep/2&&&&&&&&&&&&&
k,86,-tep+tcf,
-hep1/2-100,bep/2&&&&&&&
k,87,0,-hep1/2-100+4*tcf,bep/2&
v,80,81,82,83,84,85,86,87
!!!!!!柱翼缘!!!!!!!!!!!!!!!!!!!!
block,-tep,-tep+tcf,-hep1/2-100,-hb/2-lc,0,bc/2
block,-tep-hc+tcf,-tep-hc+2*tcf,-hb/2-lc,hb/2-(hc-2*tcf)*randa,0,bc/2&&&&&&&
!!!!!!柱腹板!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!
k,89,-tep-hc+2*tcf,-hb/2-lc
k,90,-tep,-hb/2-lc
k,92,-tep-hc+2*tcf,hb/2-(hc-2*tcf)*randa-tbf&
k,91,-tep,hb/2-tbf
k,93,-tep-hc+2*tcf,-hb/2-lc,tcw/2&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
k,94,-tep,-hb/2-lc,tcw/2&&&&&&&&&&&&&&&&&&&&&&&&&&&
k,96,-tep-hc+2*tcf,hb/2-(hc-2*tcf)*randa-tbf,tcw/2&
k,95,-tep,hb/2-tbf,tcw/2
v,89,90,91,92,93,94,95,96&&&&&&&&&&&&&&&&&&&&&&&
k,98, -tep,hb/2
k,97, -tep-hc+2*tcf,hb/2-(hc-2*tcf)*randa
-tep,hb/2,bb/2&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
-tep-hc+2*tcf,hb/2-(hc-2*tcf)*randa,bb/2&
v,91,92,97,98,95,96,99,100
block,-tep-hc+2*tcf,-tep,-hb/2,-hb/2+tst,tcw/2,tcw/2+bst&
!!!!!!端板外伸加劲肋!!!!!!!!!!!!!!!!!!!!
k,109,tep,hb/2
k,110,tep+hst,hb/2+hst*randa
k,111,tep,hb/2+hst
k,112,tep,hb/2,tst/2&&&&&&&&&&&&&&&
k,113,tep+hst,hb/2+hst*randa,tst/2&
k,114,tep,hb/2+hst,tst/2&&&&&&&&&&&
V, 109,110,111,112,113,114
k,115,tep,-hb/2&&&&&&&&&&&&&&&&&&&&&
k,116,tep+hst,-hb/2+hst*randa&&&&&&&
k,117,tep,-hb/2-hst&&&&&&&&&&&&&&&&&
k,118,tep,-hb/2,tst/2&&&&&&&&&&&&&&&&&&&&&
k,119,tep+hst,-hb/2+hst*randa
,tst/2&&&&&
k,120,tep,-hb/2-hst,tst/2&&&&&&&&&&&&&&&&&
V, 115,116,117,118,119,120
k,121,-tep,hb/2&&&&&&&&&&&&&&&&&&&&&&&&&
k,122,-tep-hst,hb/2-hst*randa&&&&&&&&&&&&&&&&&&&&&
k,123,-tep,hb/2+hst&&&&&&&&&&&&&&&&&&&&&
k,124,-tep,hb/2,tst/2&&&&&&&&&&&&&&&&&
k,125,-tep-hst,hb/2-hst*randa
k,126,-tep,hb/2+hst
,tst/2&&&&&&&&&&&&&
121,122,123,124,125,126&&&&&&&&&&&
!!!!!!梁加劲肋!!!!!!!!!!!!!!!!!!!!
wpoff,tep+lb1+lb2/2,-hb/2+tbf+lb1*randa,0
block,-tst/2,tst/2,0,hb-2*tbf,tbw/2,tbw/2+bst
block,75-tst/2,tst/2+75,0,hb-2*tbf,tbw/2,tbw/2+bst
block,-75-tst/2,-75+tst/2,0,hb-2*tbf,tbw/2,tbw/2+bst
wplane,,0,hb/2+100-aa,bep/2-aa1,0,hb/2+aa+10,bep/2-aa1,0,hb/2-aa+100+10,bep/2-aa1+10
!!!!!!端板打孔!!!!!!!!!!!!!!!!!!!
cylind,dh0/2,,-tep-10,tep+10,0,360
vsel,S,LOC,Y,hb/2+100-aa&
vgen,2,all, , ,,-ab , , ,0
vgen,2,all, , ,,-(hep1-2*aa-ab) , , ,0
cm,bolt,volu
!!!!布尔运算!!!!!!!
vsel,s,loc,x,-tep+5,tep&
VSEL,U, , ,bolt
cm,endplate,volu
VSBV,endplate,bolt,keep
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!建螺栓!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
cylind,dbt/2,,-tep,tep,0,360
cylind,dbth/2,,-tep-lbth,-tep,0,360
cylind,dbth/2,,tep,tep+lbth,0,360
wpoff,-ab,0
cylind,dbt/2,,-tep,tep,0,360
cylind,dbth/2,,-tep-lbth,-tep,0,360
cylind,dbth/2,,tep,tep+lbth,0,360
wpoff,-(hep1-2*aa-ab),0
cylind,dbt/2,,-tep,tep,0,360
cylind,dbth/2,,-tep-lbth,-tep,0,360
cylind,dbth/2,,tep,tep+lbth,0,360
cylind,dbt/2,,-tep,tep,0,360
cylind,dbth/2,,-tep-lbth,-tep,0,360
cylind,dbth/2,,tep,tep+lbth,0,360
allsel,all
NUMMRG,ALL, , , ,LOW
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!粘贴柱腹板翼缘端板和相应的加劲肋!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,x,-tep-hc,0&
vsel,u,loc,x,0,tep
VGLUE,all&
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!粘贴梁腹板翼缘端板和相应的加劲肋!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,x,0,lb1+lb2+tep
vsel,u,loc,x,-tep,0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!粘贴螺栓!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,z,bep/2-aa1
cm,bolt,volu
ALLSEL,ALL&
NUMMRG,ALL, , , ,LOW
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!定义材料和单元属性!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!下面定义材料1，Q345钢材用于梁柱端板!!!!!!!!!
!!!!!!!!!!!!!!!!
et,1,solid92
mp,ex,1,206e3
mp,prxy,1,0.3
tb,biso,1,1,2
tbdata,,345,0,,,,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!下面定义材料2，用于高强螺栓!!!!!!!!!
!!!!!!!!!!!!!!!!
mp,ex,2,206e3
mp,prxy,2,0.3
tb,biso,2,1,2
tbdata,,940,0,,,,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!给单元附材料号和属性!!!!!!!!!
!!!!!!!!!!!!!!!!
!!!!高强螺栓附属性!!!!!!!!
cmsel,s,bolt
VATT, 2,,1,0&
!!!!!!!!!!!!!!!!
!!!!梁柱端板附属性!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL,ALL
cmsel,u,bolt
VATT, 1,,1,0&
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!网格划分!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh螺栓杆!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!
cmsel,s,bolt
vsel,r,loc,x,0
LSEL,s,radius,,dbt/2
LESIZE,all, , ,4, , , , ,1&
LSEL,s,length,,2*tep
LESIZE,all, , ,10, , , , ,1&
MSHAPE,1,3D
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh螺栓头和螺母!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!
cmsel,s,bolt
vsel,u,loc,x,0
LSEL,s,radius,,dbth/2
LESIZE,all, , ,4, , , ,
LSEL,s,radius,,dh0/2
LESIZE,all, , ,8, , , , ,1
LSEL,s,length,,lbth
LESIZE,all, , ,2, , , , ,1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh端板!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!先建endplate的component!!!!!!!!!!!!!!!!!!!
allsel,all
vsel,s,loc,x,-tep+tcf/2,tep
vsel,u,loc,x,0
cm,endplate,volu
ESIZE,bep/10
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh柱!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,x,-tep-hc,0
vsel,r,loc,z,0,tcf
vsel,u,loc,y,hb/2,hb/2+hst
cm,col_web,volu
eSIZE,hc/8
vsel,s,loc,x,-tep-hc,0
cmsel,u,bolt
cmsel,u,endplate
cmsel,u,col_web
eSIZE,bc/6
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!mesh梁!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,x,tep,tep+lb1+lb2
cmsel,u,bolt
lsel,s,loc,x,tep+lb1+lb2
LSEL,r,length,,bb/2
LESIZE,all,bb/6
lsel,s,loc,x,tep+lb1+lb2
LSEL,r,length,,hb-2*tbf
LESIZE,all,hb/8
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!接触建立!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
allsel,all
cmsel,s,endplate
asel,r,loc,x
asel,u,loc,y,-hep/2,-hep/2-10
cm,target1,area
allsel,all
cmsel,s,endplate
asel,r,loc,x
asel,r,loc,y,-hep/2,-hep/2-10
cm,contact1,area
cmsel,s,target1
cmsel,a,contact1
MP,MU,1,miu
KEYOPT,4,9,0
ASEL,S,,,target1
ESURF,ALL&
ASEL,S,,,contact1
ESURF,ALL&
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!对高强螺栓psmesh!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
cmsel,s,bolt
vsel,r,loc,x,
cm,boltshank1,volu
PSMESH,1,bolt1, ,v,1, 0,x,
PSMESH,2,bolt2, ,V,22, 0,x,
PSMESH,3,bolt3, ,V,27, 0,x,
PSMESH,4,bolt4, ,V,30, 0,x,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!下面施加荷载和位移条件进行求解
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!下面施加loadstep 1,边界约束条件、柱轴向压力和螺栓预拉力
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL,ALL&
NSUBST,10,0,0&
NSUBST,10,0,0&
OUTRES,alL,1
NEQIT,200&
CNVTOL,F, ,0.01,2, ,
ASEL,S,LOC,Z,&
DA,all,SYMM
nsel,s,loc,y,-(hb/2+lc+1)
D,all, ALL
ALLSEL,ALL
SLOAD,ALL,9,LOCK,FORC,preten, 1,2
LSWRITE,1,
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!下面施加loadstep 2,粱端荷载
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
ALLSEL,ALL
ANTYPE,,REST,,,0&
NSUBST,20,5000,1
OUTRES,alL,1
ALLSEL,ALL&
num=node(tep+lb1+lb2/2,hb/2+(lb1+lb2/2)*randa,0)&&&
D,num, uy,displa&
LSWRITE,2,&
ALLSEL,ALL
27.在任意面施加任意方向任意变化的压力
在某些特殊的应用场合，可能需要在结构件的某个面上施加某个坐标方向的随坐标位置变化的压力载荷，当然，这在一定程度上可以通过ANSYS表面效应单元实现。如果利用ANSYS的参数化设计语言，也可以非常完美地实现此功能，下面通过一个小例子描述此方法。
!!!在执行如下加载命令之前,请务必用选择命令asel将需要加载的几何面选择出来
et,500,shell63
press=100e6
! 如果载荷的反向是一个特殊坐标系的方向,可在此建立局部坐标系,并将
! 所有节点坐标系旋转到局部坐标系下.
*get,enmax,elem,,num,max !把当前选中的对象的最大单元号赋予enmax
dofsel,s,fx,fy,fz
fcum,add !!!将力的施加方式设置为"累加",而不是缺省的"替代"
*do,i,1,enmax
*if,esel(i),eq,1,then !i若在所选择单元中，esel（i）＝1，否则为0。
*get,ae,elem,i,area !此命令用单元真实面积，如用投影面积，请用下几条命令
! *get,ae,elem,i,aproj,x !此命令用单元X投影面积，
! *get,ae,elem,i,aproj,y&&
!此命令用单元Y投影面积
! *get,ae,elem,i,aproj,z&&
!此命令用单元Z投影面积
xe=centrx(i)&&&
!单元i中心X坐标(用于求解压力值)
ye=centry(i)&&&
!单元i中心Y坐标(用于求解压力值)
ze=centrz(i)&&&
!单元i中心Z坐标(用于求解压力值)
! 下面输入压力随坐标变化的公式,本例的压力随X和Y坐标线性变化.
p_e=(xe-10)*press+(ye-5)*press
f_tot=p_e*ae
esel,s,elem,,i
nsle,s,corner
*get,nn,node,,count
f_n=f_tot/nn
*do,j,1,nn
f,nelem(i,j),fx,f_n !压力的作用方向为X方向
! f,nelem(i,j),fy,f_n !压力的作用方向为Y方向
! f,nelem(i,j),fz,f_n !压力的作用方向为Z方向
aclear,all
fcum,repl !!!将力的施加方式还原为缺省的"替代"
dofsel,all
33.等截面杆单元生死应用实例
等截面杆两端固定，承受均匀的温度载荷时将其中间1/3段移去。过程是将其应变自由化并移去均匀温度。分析其热应力和应变情况。
& /TITLE, ELEMENT BIRTH/DEATH IN A FIXED BAR WITH
THERMAL LOADING
& ET,1,LINK1 !二维杆单元
& MP,EX,1,30E6 !材料特性
& MP,ALPX,1,.00005
& MP,EX,2,30E6
& MP,ALPX,2,.00005 !重新‘出生’单元的特性
& MP,REFT,2,100 !单元出生的参考温度
& EGEN,3,1,-1 !生成三个单元
& ANTYPE,STATIC
& D,1,ALL,,,4,3 !固定杆的两端
& TREF,0 !参考温度0度
& TUNIF,100 !均匀温度载荷
& NROPT,FULL
& OUTPR,BASIC,ALL
& EKILL,2 !‘杀死’中间的单元
& EALIVE,2 !重新激活中间单元
& MPCHG,2,2 !将材料特性改为2以生成自由应变单元
& TUNIF,0 !删除温度载荷
36.简单的例子－如何对结构的振动控制分析
第一步：首先建立模型，如图，一个很简单的角形结构梁
其中上面的号码是节点号，连接1，10的就是调频器，节点2与1均有质量。
/clear,start
et,1,beam44
R,1,78,680,,
RMODIF,1,21,,
RMODIF,1,25,,0,0,7.,
UIMP,1,EX, , ,210000,
UIMP,1,DENS, , ,1.1261e-5,
et,2,mass21
第二步：进行模态分析
antype,modal
modopt,subsp,10
自己进入后处理器，察看结果
&&& TIME 1
&&&&&&&&&&&&&&&&&
再进入扩展模态，找出对应想抑制震动的主自由度
在上面就以二阶和三阶为例吧，他们的振型是在y，z平面，这也好理解，这个方向刚度很低，所以低频率就出现在这个振型上。所以上面加的tmd的方向也很讲究，也就是节点10的位置是垂直于y，z平面，不然，在其他方向就没有好的效果了。
et,3,combin14
r,3,100,10
r,4,100,100,100
n,10,0,0,-1000
然后进入阻尼模态分析，由于关心的频率范围是前两阶，因此，设置了频率范围在0.1到0.2，
在此强调一下，在damp中频率范围一定要设置，不然感兴趣的结果有很多遗漏。
antype,modal
modopt,damp,10
MXPAND,10, , ,0
MODOPT,DAMP,10,0.1,0.2, ,OFF
然后进入后处理器
有兴趣者可以对tmd命令流中的r,3,100,10中的阻尼进行变化
就可发现阻尼的作用了
对于低频作用不大，但对高频有很好的抑制效果
弹簧刚度的作用相反。
这个简单的例子是对其方法的一个简单描述，可能有不少错误和不详细的地方，希望大家共同探讨指正，谢谢
37.模态分析结果的输出实例
在gui里看到的既不能修改，也不能copy
现用命令流导出
考虑斑竹以积分鼓励大家，以一分为限
例子如下：
/clear,start
mass=83500
et,1,beam44
et,2,mass21
R,1,78,680,,
RMODIF,1,21,,
RMODIF,1,25,,0,0,7.,
R,2,mass-j,mass-j,mass-j,0,0,0,
UIMP,1,EX, , ,210000,
UIMP,1,DENS, , ,1.1261e-5,
UIMP,2,EX, , ,210000,
UIMP,2,DENS, , ,0.78e-5,
latt,1,1,1
lesize,1,,,5,1
katt,2,2,2
d,1,,,,,,ux,uy,uz,rotx,roty,rotz
antype,modal
modopt,subsp,5
FILE !确定数据源
FILE,'file','rst','.'
/UI,COLL,1 !打开对话框
NUMVAR,200 !post26中最大允许存100个
SOLU,191,NCMIT !确定每时间解的存储
STORE,MERGE ! 存数据
FILLDATA,191,,,,1,1 !以斜坡函数来填充数据
REALVAR,191,191 !只使用实部
! Save time history variables to file b.csv
*DEL,_P26_EXPORT !删除参数
*DIM,_P26_EXPORT,TABLE,5,1 ! 定义数组大小
VGET,_P26_EXPORT(1,0),1 !移动数据
*if,i,eq,1,then
/OUTPUT,'a1','csv','.' ! 生成文件
*elseif,i,eq,2,then
/OUTPUT,'a2','csv','.'
*elseif,i,eq,3,then
/OUTPUT,'a3','csv','.'
*elseif,i,eq,4,then
/OUTPUT,'a4','csv','.'
*elseif,i,eq,5,then
/OUTPUT,'a5','csv','.'
*VWRITE,'TIME' ! 按一定次序写入数据
*VWRITE,_P26_EXPORT(1,0)
/OUTPUT,TERM
! End of time history save
38.火车过桥动态加载实例（部分）
! 用对话框提示用户输入基本数据
!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
*if,ttype,eq,0, then
&& q1 = 147000 !中-活载
&& q2 = 92000
&& q3 = 80000
&& L1 = 7.5
&& L2 = 30
*elseif,ttype,eq,1,then
&& *msg,ui
&& error value of ttype(1 or
2),please enter again
&& *go,: reread
! 识别桥面过列车的单元，并将单元按节点坐标从小到大排序
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 ! 存放桥面过车的节点
*do,i,1,elmax
& *if,esel ,gt,0, then
&&& eldk(j) =
*enddo ! 记录桥面单元的标号，保存于数组eldk中
eldknum = j ! 总桥面单元数
! 对eldk数组排序，按节点坐标从小到大排列
*do,i,1,eldknum
&& ii = eldk
&& ndi = nelem(ii,1)
&& ndj = nelem(ii,2)
&& i1 = nx(ndi)
&& i2 = nx(ndj)
&& ndx = i1
&& *if,i2,lt,i1,then
*do,i,1,eldknum-1
& *do,j,i+1,eldknum
*if,ndx(j),lt,ndx ,then
temp = ndx
ndx = ndx(j)
ndx(j) = temp
temp = eldk
eldk = eldk(j)
eldk(j)= temp
&&& *endif
nsle,s ! 选中桥面单元的节点
*do,i,1,ndmax
& *if,nsel ,gt,0,then
nodes(j)=i
ndnum = j ! 总桥面节点数
xmin = nx(nelem(eldk(1),1))
temp = nx(nelem(eldk(1),2))
*if,temp,lt,xmin,then
& xmin = temp
xmax = nx(nelem(eldk(eldknum),1))
temp = nx(nelem(eldk(eldknum),2))
*if,temp,gt,xmax,then
& xmax = temp
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
*do,kk,0,lth,5,
time = time+1
!/uis,msgpop,3 ! 不显示数组被重新定义的警告
!/uis,msgpop,2
/solu ! 进入求解器
antype,static
sfedele,all,1,pres, !先删除上一次加上的所有单元荷载
! 计算列车各部分在桥上的单元号
tpos1 = kk ! 火车头所在的位置
*if,tpos1,lt,(qth+L1+L2), then !火车进桥情况
&& xlth1 = 0
&& *do,i,1,eldknum
ndi = nelem(ii,1)
ndj = nelem(ii,2)
elth1 = abs(nx(ndi)-nx(ndj))
xlth1 = xlth1 + elth1
*if,xlth1,ge,tpos1,then
&& tele1 = i ! 火车头所在的单元
&& *if,tpos1,gt,L1,then
tpos2 =tpos1 -L1 !煤水车头所在的位置
*do,i,1,eldknum
ndi = nelem(ii,1)
ndj = nelem(ii,2)
elth2 = abs(nx(ndi)-nx(ndj))
xlth2 = xlth2 + elth2
*if,xlth2,ge,tpos2,then
tele2 = i !旅客列车头所在的单元
*if,tpos1,gt,(L1+L2),then
tpos3 =tpos1 -(L1+L2) !旅客车头所在的位置
*do,i,1,eldknum
ndi = nelem(ii,1)
ndj = nelem(ii,2)
elth3 = abs(nx(ndi)-nx(ndj))
xlth3 = xlth3 + elth3
*if,xlth3,ge,tpos3,then
&&&&&&&&&&&&
tele3 = i !旅客列车头所在的单元
*if,tpos1,gt,(qth+L1+L2+182.5),then ! 火车过桥情况(182.5为任意取的一个常数)
& *if,tpos1,lt,lth1, then
&&& tpos4 =
tpos1-(qth+L1+L2+182.5) !旅客列车尾部所在的位置
&&& xlth4 =
*do,i,1,eldknum
ndi = nelem(ii,1)
ndj = nelem(ii,2)
elth4 = abs(nx(ndi)-nx(ndj))
xlth4 = xlth4 + elth4
*if,xlth4,ge,tpos4,then
&&& *enddo
&&& tele4 = i
!旅客列车尾部在的单元
! 加荷载于相应的单元
! SFBEAM, ELEM, LKEY, Lab, VALI, VALJ, VAL2I, VAL2J, IOFFST,
*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,,,,,
*do,j,1,tele1-1
sfbeam,eldk(j),1,pres,q1,q1,,,,,
qq = q1*c2/elth1
sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
*elseif,tpos1,le,(L1+L2),then ! 煤水车进桥
& c1 = xlth1 - tpos1
& *if,tele1,eq,tele2,then !火车头在同一个单元上
&&& c2 = xlth2 -
&&& c3 = elth2 -
(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
&&& 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,,,,,
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,,,,,
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,,,,,
*elseif,tpos1,le,qth,then ! 旅客列车进桥
& c1 = xlth1 - tpos1
& *if,tele1,eq,tele2, then
&&& c2 = xlth2 -
&&& c3 = elth2 -
(q1*L1+q2*c3)/elth1
sfbeam,eldk(tele1),1,pres,qq,qq,,,,,
&&& c4 = xlth3 -
&&& c5 = elth3 -
&&& 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,,,,,
*do,j,tele3+1,tele2-1
sfbeam,eldk(j),1,pres,q2,q2,,,,,
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,,,,,
&&& 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,,,,,
*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,,,,,
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,,,,,
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,,,,,
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,,,,,
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,,,,,
*do,j,tele3+1,tele2-1
&&&&&&&&&&&&
sfbeam,eldk(j),1,pres,q2,q2,,,,,
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,,,,,
&&& *endif
*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,,,,,
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,,,,,
&& *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,,,,,
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,,,,,
&& *elseif,dtele1,gt,1,then
*do,j,tele2+1,eldknum
sfbeam,eldk(j),1,pres,q1,q1,,,,,
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,,,,,
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,,,,,
*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,,,,,
*do,j,tele3+1,eldknum
sfbeam,eldk(j),1,pres,q2,q2,,,,,
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,,,,,
*elseif,tpos1,le,(qth+L1+L2+182.5),then ! 全桥均布列车荷载(182.5为任意取的
&& *do,j,1,eldknum
sfbeam,eldk(j),1,pres,q3,q3,,,,,
*elseif,tpos1,le,lth1,then ! 均布列车荷载出桥
&& c1 = xlth4 - tpos4
&& c2 = elth4 - c1
&& *do,j,tele4+1,eldknum
sfbeam,eldk(j),1,pres,q3,q3,,,,,
&& qq = q3*c1/elth4
sfbeam,eldk(tele4),1,pres,qq,qq,,,,,
/psf,pres,norm,1
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))
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
save,,,,all
! 结果处理
*do,i,1,ndnum
& nd = nodes
& disymax = abs(disy(1,nd))
& *do,j,2,dim
*if,abs(disy(j,nd)),gt,disymax ,then
disymax = abs(disy(j,nd))
&&& *endif
*do,i,1,ndnum-1
&& *do,j,i+1,ndnum
*if,disymax(j),gt,disymax ,then
temp = disymax
disymax = disymax(j)
disymax(j) = temp
temp = nodes
nodes = nodes(j)
nodes(j) = temp
*do,i,1,ndnum ! 前5个位移最大的节点
&& numnd = nodes
&& disy1(0,i) = numnd
&& *do,j,1,dim
kk = (j-1)*5
disy1(j,0) = kk
disy1(j,i) = disy(kk,numnd)
*do,i,1,elmax
&& fmin = ffx(i,1)
&& fmax = ffx(i,1)
&& *do,j,2,dim
*if,ffx(i,j),gt,fmax ,then
fmax = ffx(i,j)
*if,ffx(i,j),lt,fmin ,then
fmin = ffx(i,j)
detaf = fmax - fmin ! 每个单元的内力变化值
detaf1 = detaf
*do,i,1,elmax
&& elnum = i
*do,i,1,elmax-1
&& *do,j,i+1,elmax
*if,detaf1(j),gt,detaf1 ,then
temp = detaf1
detaf1 = detaf1(j)
detaf1(j) = temp
temp = elnum
elnum = elnum(j)
elnum(j)=temp
*do,i,1,20
& detaf2 = detaf1 !前20个内力变化最大的单元
& elnum2 = elnum
*do,i,1,20 ! 前20个内力变化最大的单元
&& numel = elnum2
&& 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) ! 相应的弯矩
! 前20个内力最大的单元
*do,i,1,elmax
&& fabs =abs(fmax )
&& *if,fabs ,lt,abs(fmin
fabs =abs(fmin )
*do,i,1,elmax
&& elnum3 =i
*do,i,1,elmax-1
&& *do,j,i+1,elmax
*if,fabs(j),gt,fabs ,then
temp = fabs
fabs = fabs(j)
fabs(j) =temp
temp = elnum3
elnum3 = elnum3(j)
elnum3(j) = temp
&&& *enddo
*do,i,1,20
& elnum4 = elnum3
& nmax = fabs
*do,i,1,20 ! 前20个内力最大的单元
&& numel = elnum4
&& 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) ! 相应的弯矩
save,,,,all
!/title,the 91 node displament vs train position ! 显示图形的名称
!*vplot,diy(1,0),diy(1,1)
!/ui,copy,save,bmp,graph
*do,i,1,30
cc = nodes
/title, y displacement of node ?% VS. train position
*cfopen,result,txt,,,
('**********************南京桥计算结果*****************************')
(' ****节点最大位移**** ')
ndnum1 = 5*nint(ndnum/5)
*do,i,1,ndnum1,5
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(
(f5.0,,5(3x,f10.5))
(' ****前20个轴力变化最大的单元和相应的轴力值**** ')
*vwrite,elnum2(1),detaf2(1)
(f15.0,' ', f15.1)
*do,i,1,18,3
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))
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))
(' ***前20个轴力最大的单元和相应的弯矩值*** ')
*vwrite,elnum4(1),detaf2(1)
(f15.0,' ', f15.1)
*do,i,1,18,3
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))
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))
39.悬索结构的找形和计算的例题
取自&悬索结构设计&沈世钊，中国建筑工业出版设，1997.8
p149例题3-1
进行了找形计算和荷载作用下的静力计算
找形时取较小的弹性模量，静力计算时采用弹性模量的真值
这个例题说明，ANSYS还是可以进行悬索结构的分析的
ET,1,LINK10
!找形时采用很小的弹性模量
MP,EX,1,0.60551e9
MP,NUXY,1,0.3
*do,j,1,2*i-1
k=(i-1)*(i-1)+j
n,k,(j-i)*9.15,-36.6+(i-1)*9.15,0
*do,j,1,2*i-1
k=(i-1)*(i-1)+j+25
n,k,(j-i)*9.15,36.6-(i-1)*9.15,0
*do,j,1,2*i
*do,j,1,2*i
k=i*i+j+25
*do,j,1,2*i-1
k=(i-1)*(i-1)+j
*do,j,1,2*i-1
k=(i-1)*(i-1)+j+25
*do,i,18,24
!施加位移约束
d,1,uz,3.66
d,2,uz,1.83
d,4,uz,1.83
d,10,uz,-1.83
d,16,uz,-1.83
d,17,uz,-3.66
d,25,uz,-3.66
d,35,uz,-1.83
d,41,uz,-1.83
d,27,uz,1.83
d,29,uz,1.83
d,26,uz,3.66
!设定时间步
NSUBST,20, , ,1
!输出每个子步的结果
OUTPR,BASIC,ALL,
OUTRES,ALL,ALL,
!设定大变形，应力刚化
/STATUS,SOLU
PRNSOL,DOF,
PRESOL,SMISC,1
PLESOL,SMISC,1
PLNSOL,U,X,0,1
!更新几何形状
!a.rst为计算结果文件名，最后一个为目录
!这两个参数应根据你的计算情况定
UPGEOM,1,LAST,LAST,a,rst,E:JZD1
!弹性模量恢复为真值
R,1,0..0027248,
MP,EX,1,2.0e11
!重新施加位移约束
/STATUS,SOLU
!观察施加荷载之前结构的位移
!可见位移基本为零,预张力基本不变
PLNSOL,U,Z,0,1
PLESOL,SMISC,1
!写荷载工况文件01
LSWRITE,01,
!施加节点荷载
*do,i,1,41
f,i,fz,-167445
!设置第2荷载步
NSUBST,20, , ,1
LSWRITE,02,
LSSOLVE,1,2,1,
!可以看到，该结果和书中结果一致
NSOL,2,21,U,Z,
PLVAR,2, , , , , , , , , ,
ESOL,3,33,1,F,Y,
PLVAR,3, , , , , , , , , ,
44.一个例子（含地震影响，求振兴与频率）
!地处Ⅱ类场地,抗震设防烈度为八度,近震
!反应谱最大地震影响系数αｍａｘ=0.16,
!特征周期Ｔｇ=0.30.结构阻尼比取为0.05.
BLOCK,0,0.15,0,2.85,0,-0.15,
BLOCK,0.15,3,0,2.85,0,-0.15,
BLOCK,3,6,0,2.85,0,-0.15
BLOCK,6,9,0,2.85,0,-0.15
BLOCK,0,0.15,0,2.85,-0.15,-4.5,
BLOCK,0,0.15,2.85,3,0,-0.15,
BLOCK,0.15,3,2.85,3,0,-0.15
BLOCK,3,6,2.85,3,0,-0.15
BLOCK,6,9,2.85,3,0,-0.15
BLOCK,0,0.15,2.85,3,-0.15,-4.5,
BLOCK,0.15,3,2.85,3,-0.15,-4.5
Block,3,6,2.85,3,-0.15,-4.5
Block,6,9,2.85,3,-0.15,-4.5
BLOCK,0,0.15,3,6,0,-0.15,
BLOCK,0.15,3,3,6,0,-0.15,
BLOCK,3,6,3,6,0,-0.15,
BLOCK,6,9,3,6,0,-0.15,
BLOCK,0,0.15,3,6,-0.15,-4.5,
vsel,s,,,all
/pnum,area,on
vsel,s,,,all
cm,f_sub,volum !前面部分
cmsel,s,f_sub
vsymm,z,all,,,,0,0
vsel,s,,,all
vsel,u,,,f_sub
vgen,,all,,,,,-9,,,1
vsel,s,,,all
et,1,solid45
mp,ex,1,2.5e8
mp,nuxy,1,0.2
mp,dens,1,2500
lsel,s,,,all
lsel,s,length,,,0.15
lesize,all,,3
lsel,s,,,all
lsel,u,length,,0.15
lesize,all,0.15
mshape,0,3D
vsel,s,,,all
asel,s,loc,y,0
cm,base_f,area !基础面
asel,s,,,all
asel,s,loc,x,9
cm,l_wall,area !x对称面
antype,modal !求振型和频率
modopt,subsp,4 !用子空间迭代法求前3阶振型
asel,s,,,all
cmsel,s,base_f
asel,s,,,all
asel,s,,,l_wall
dsym,symm,x
set,list,2
/view,1,1,1,1,
anmode,10,0.5e-1,20
/anfile,save,'1','avi',''
/view,1,1,1,1,
anmode,10,0.5e-1,20
/anfile,save,'2','avi',''
/view,1,1,1,1,
anmode,10,0.5e-1,20
/anfile,save,'3','avi',''
antype,spectr !求振型和频率
spopt,sprs,3,yes !单点谱分析
svtype,2 !加速度谱
sed,1,0,0 !X方向地震作用
!定义反应谱(抗震规范GBJ11-89)：Tg=0.4s,amax=0.16(8度) !F=a*G=a*m*g
a值乘以9.82
freq,0.........,
sv,,0.........,
freq,0.......5,10,1000
sv,,0.......2,0.70704
antype,modal !求振型和频率
mxpand,3,,,yes !扩展前3阶振型，计算单元解和反力
antype,spectr
srss,1,disp !应用SRSS方法计算单元应力,该命令将写*.mcom文件
/inp,,mcom
pletab,avg
prnsol,dof
presol,elem
48.一个圆形水池的静力分析
!-------------------------------------------
ET,1,SHELL63 !壳单元
MP,EX,1,2.1e5 !C20混凝土
MP,NUXY,1,0.2
R,1,200, , , , , , !壁厚200
R,2,150, , , , , , !肋厚150
!-------------------------------------------
k,2,325,0,0
k,3,325,520.4250,0
k,10,0,5975,0
AROTAT,1,2,3,4,5,6,1,10,360,4,
AROTAT,7,8, , , , ,1,10,360,4,
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
!-------------------------------------------
ASEL,S,,,1,24,1,,
AATT,1,1,1,0,
ASEL,S,,,25,32,1,,
AATT,1,2,1,0,
ALLSEL,ALL
AESIZE,ALL,400,
!-------------------------------------------
LSEL,S,,,27,29,2,,
LSEL,A,,,33,35,2,,
LSEL,A,,,39,41,2,,
LSEL,A,,,45,47,2,,
DL,ALL,,ALL,
ALLSEL,ALL
!-------------------------------------------
!施加静水压力
ASEL,S,,,1,5,1,,
ASEL,A,,,7,11,1,,
ASEL,A,,,13,17,1,,
ASEL,A,,,19,23,1,,
SFGRAD,PRES,0,Y,e-5
SF,ALL,PRES,0
ALLSEL,ALL
!-------------------------------------------
!-------------------------------------------
PLNSOL,U,SUM,0,
49.ANSYS中混凝土模式预应力模拟的算例
a.初应力方式&
b.温度方式（建议使用）&
c.等效荷载法（细部不便模拟）&
下面有一小例，使用solid95和link8,或许对你有参考！&
!简支梁实体与预应力钢筋分析
/COM, Structural
egjx=2e5 !Ey
agjx=140 !单根钢绞线面积
ehnt=4e4 !Eh
xzxs=1.0e-5 !线胀系数
yjl=200000 !定义预加力
et,1,link8 !定义link8单元
et,2,solid95 !定义solid95单元
r,1,agjx !定义link8单元的面积
r,2 !定义第2种实常数
mp,ex,1,egjx !定义link8单元的弹性模量
mp,prxy,1,0.3 !定义link8单元的泊松系数
mp,alpx,1,1.0e-5 !定义线膨胀系数
mp,ex,2,ehnt !定义solid95单元的弹性模量
mp,prxy,2,0.3 !定义solid95单元的泊松系数
blc4, , ,100,200,3000 !定义梁体
/view,1,1,1,1 !定义ISO查看
vplot !绘制梁体
kwpave,6 !工作平面移动到关键点6
wpoff,-30 !工作平面移动-30mm(X)
wprot,0,0,90 !工作平面旋转
vsbw,1 !分割梁体
wpoff,0,0,-40 !工作平面移动-40mm(Z)
vsbw,2 !分割梁体
wpoff,0,40 !工作平面移动40mm(Y)
wprot,0,90 !工作平面旋转
vsbw,all !分割梁体
wpstyl !关闭工作平面显示
nummrg,all,,,,low !整理
numcmp,all !压缩编号
esize,30 !定义网分时边长控制
lsel,s,,,28,38,10 !定义line28和38为新的选择集
latt,1,1,1 !定义选择集的属性
lmesh,all !对线划分单元
allsel,all !新的选择集为所有的实体
gplot !绘制所有的实体
vsel,s,,,all !定义所有体为选择集
vatt,2,2,2 !定义选择集的属性
mshape,0,3d !将体划分单元的形状定位HEX
mshkey,1 !采用MAPPED划分器
vmesh,all !对体进行划分单元
dl,3,,all !对线line7施加约束(UX,UY,UZ)
dl,16,,all !对线line31施加约束(UX,UY,UZ)
dl,23,,all !对线line23施加约束(UX,UY,UZ)
dl,2,,uy !对线line4施加约束(UY)
dl,15,,uy !对线line30施加约束(UY)
dl,22,,uy !对线line23施加约束(UY)
dk,2,,,,,ux,uy !对关键点2约束(UX,UY)
bfl,28,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度
bfl,38,temp,-yjl/(xzxs*egjx*agjx) !对钢绞线施加温度
solve !求解
plnsol,s,z,0,1 !绘制Z方向的应力
etable,sigi,ls,1 !定义钢筋单元数据表
plls,sigi,sigi,1 !绘制上述应力
!/exit,nosav
50.悬臂梁受重力作用发生大变形求其固有频率
基本过程：
2、静力分析
&& NLGEOM,ON
&& PSTRES,ON
3、求静力解
4、开始新的求解：modal
&& PSTRES,ON
&& UPCOORD,1,ON 修正坐标
&& PSOLVE...
5、扩展模态解
6、察看结果
ET,1,BEAM189 !使用beam189梁单元
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,210e9
MPDATA,PRXY,1,,0.3
MPDATA,DENS,1,,7850
SECTYPE, 1, BEAM, RECT, secA, 0 !定义梁截面secA
SECOFFSET, CENT
SECDATA,0.005,0.01,0,0,0,0,0,0,0,0
K, ,,,, !建模与分网
LSTR, 1, 2
LATT,1, ,1, , 3, ,1
LESIZE,1, , ,20, , , , ,1
/SOL !静力大变形求解
PSTRES,ON !计及预应力效果
DK,1, , , ,0,ALL, , , , , ,
ACEL,0,9.8,0, !只考虑重力作用
NSUBST,20, , ,1
ANTYPE,2 !进行模态求解
MODOPT,LANB,10
MXPAND,10, , ,0 !取前十阶模态
PSTRES,1 !打开预应力效应
MODOPT,LANB,10,0,0, ,OFF
UPCOORD,1,ON !修正坐标以得到正确的应力
PSOLVE,TRIANG !三角化矩阵
PSOLVE,EIGLANB !提取特征值和特征向量
EXPASS,1 !扩展模态解
PSOLVE,EIGEXP
SET,LIST !观察结果
51.循环对称结构模态分析
这是ANSYS HELP里的例子，但那个命令流似乎有些问题，下面是整理过的命令流。
基本过程：
2、define cyclic
3、定义约束
4、求模态解
5、展开并察看结果
r1=5 !建模
alpha_deg=360/nsect
alpha_rad=2*acos(-1)/nsect
k,2,0,0,d1
k,3,r1,0,0
k,4,r1,0,d1
arotat,1,,,,,,1,2,alpha_deg/2
k,7,r2,0,0
k,8,r2,0,d1
arotat,5,,,,,,1,2,alpha_deg/2
arotat,2,,,,,,1,2,alpha_deg/2
arotat,6,,,,,,1,2,alpha_deg/2
a,5,6,10,9
mp,ex,1,10e6
mp,prxy,1,0.3
mp,dens,1,1e-4
asel,,,,1,4
CYCLIC, , , ,'CYCLIC' !定义基本扇区
AMESH,all !分网
/SOL !求模态解
MODOPT,LANB,10
MXPAND,10, , ,0
MODOPT,LANB,10,0,0, ,OFF
DK,5, , , ,0,UZ, , , , , ,
SET, , ,1, , , ,4,
/CYCEXPAND, ,ON
PLNS,U,SUM !观察扩展模态
52.三角平台受谐波载荷作用的结构响应
谐波载荷为作用在平台上面一点的集中力，幅值为10N，频率范围5Hz~70Hz
基本过程：
2. 求模态解
3、用模态叠加法作谐分析
4. 观察结果
/PREP7 !建模
ET,1,BEAM189
ET,2,SHELL93
R,1,0.01, , , , , ,
MP,EX,1,210e9
MP,PRXY,1,0.3
MP,DENS,1,7850
SECTYPE, 1, BEAM, RECT, secA, 0
SECOFFSET, CENT
SECDATA,0.005,0.008,0,0,0,0,0,0,0,0
K, ,-0.5,,,
K, ,0.5,,,
K, ,-0.5,1,,
K, ,0.5,1,,
K, ,,,0.5,
LSTR, 1, 5
LSTR, 3, 4
LSTR, 2, 6
LSEL,S,LINE,,4,6
LATT,1,1,1, , 7, ,1
LSEL, , , ,ALL
LESIZE,ALL, , ,10, , , , ,1
MSHAPE,0,2D
/SOL !为了使用模态叠加法谐分析
MODOPT,LANB,20 !先取结构前20阶模态
MXPAND,20, , ,0
MODOPT,LANB,20,0,0, ,OFF
KSEL,S,KP,,1,3
DK,ALL, , , ,0,ALL, , , , , ,
!/POST1 !最好事先看一下模态结果
!SET,LIST !以了解模态频率范围
/SOL !模态叠加法谐响应分析
HROPT,MSUP
HROPT,MSUP,20, ,
HROUT,ON,OFF,0
F,177,FY,-10, !施加10N的力在节点177的负Y方向
HARFRQ,5,70, !载荷的强制频率范围
NSUBST,200,
/POST26 !察看位移响应
FILE,,rfrq
NSOL,2,139,U,Y,
53.三角平台受一地震谱激励的应力分布和支反力
基本过程：
2、求模态解
4、扩展模态
5、模态合并
6、观察结果
/PREP7 !建模
ET,1,BEAM189
ET,2,SHELL93
R,1,0.01, , , , , ,
MP,EX,1,210e9
MP,PRXY,1,0.3
MP,DENS,1,7850
SECTYPE, 1, BEAM, RECT, secA, 0
SECOFFSET, CENT
SECDATA,0.005,0.008,0,0,0,0,0,0,0,0
K, ,-0.5,,,
K, ,0.5,,,
K, ,-0.5,1,,
K, ,0.5,1,,
K, ,,,0.5,
LSTR, 1, 5
LSTR, 3, 4
LSTR, 2, 6
LSEL,S,LINE,,4,6
LATT,1,1,1, , 7, ,1
LSEL, , , ,ALL
LESIZE,ALL, , ,10, , , , ,1
MSHAPE,0,2D
/SOL !取前十阶模态
MODOPT,LANB,10
KSEL,S,KP,,1,3
DK,ALL, , , ,0,ALL, , , , , ,
/SOL !谱分析
SPOPT,SPRS,10,1
SED,0,1,0, !给出激励方向
FREQ,0.25,1.34,6.73,14.6,28.9,0,0,0,0 !激励谱
SV,0,0.16,0.34,0.00052,
/SOL !扩展模态
MXPAND,10,0,0,1,0.001,
/SOL !模态合并
SRSS,0.001,DISP
/POST1 !观察结果
/INPUT,,mcom
54.三角平台受时程载荷作用的应力分布和变形过程
基本过程：
2. 施加随时间变化载荷，定义载荷步
4. 在POST1和POST26中观察结果
/PREP7 !建模
ET,1,BEAM189
ET,2,SHELL93
R,1,0.01, , , , , ,
MP,EX,1,210e9
MP,PRXY,1,0.3
MP,DENS,1,7850
SECTYPE, 1, BEAM, RECT, secA, 0
SECOFFSET, CENT
SECDATA,0.005,0.008,0,0,0,0,0,0,0,0
K, ,-0.5,,,
K, ,0.5,,,
K, ,-0.5,1,,
K, ,0.5,1,,
K, ,,,0.5,
LSTR, 1, 5
LSTR, 3, 4
LSTR, 2, 6
LSEL,S,LINE,,4,6
LATT,1,1,1, , 7, ,1
LSEL, , , ,ALL
LESIZE,ALL, , ,10, , , , ,1
MSHAPE,0,2D
TRNOPT,FULL !完全法瞬态分析
KSEL,S,KP,,1,3
DK,ALL, , , ,0,ALL, , , , , ,
OUTRES,ALL,1
ALPHAD,5, !α阻尼
TIME,2 !定义载荷曲线
NSUBST,50, , ,1
SFA,ALL,1,PRES,500
LSWRITE,1,
LSWRITE,2,
SFA,ALL,1,PRES,150
LSWRITE,3,
SFA,ALL,1,PRES,
LSWRITE,4,
LSSOLVE,1,4,1, !求解
/POST26 !观察变形随时间的变化
NSOL,2,177,U,Y,uyy
55.经典层合板理论
下面是我用APDL写的程序。程序基于经典层合板理论。写的比较冗长，希望牛人能够给点意见，有没有哪位大侠用ANSYS做过层合板缺口件强度。希望交流交流程序写的很辛苦。
拿出来共享，希望能加分。
*AFUN,DEG !参数定义为角度
OVER=1 !定义控制结束的参数,OVER=0,循环结束
N=8 !定义总层数
H=1 !定义总厚度
TK=H/N !定义每层的厚度
PX=100 !定义x方向外载荷
PY=0 !定义y方向外载荷
PXY=0 !定义xy方向剪切载荷
*dim,ARG,array,N !定义铺层的角度
ARG(3)=-45
Strainx=0 !定义中面x方向应变
strainy=0 !定义中面y方向应变
strainxy=0 !定义中面xy方向应变
*dim,strainSTEP,array,2*N !定义各个阶段的应变
*dim,pxstep,array,2*N !定义各个阶段的应力
*dim,strain1,array,N !定义各层x方向主应变
*DIM,Strain2,array,N !定义各层y方向主应变
*dim,Strain12,array,N !定义各层xy方向主应变
*dim,s1,array,N !定义各层X方向主应力
*dim,s2,array,N !定义各层Y方向主应力
*dim,s12,array,N !定义各层XY方向的剪应力
!****************************************************************
!*Define the strength parameter of lamia
!****************************************************************
XT=2250 !定义纤维方向拉伸强度
XC=1600 !定义纤维方向压缩强度
YT=34.5 !定义基体方向拉伸强度
YC=200 !定义基体方向压缩强度
s=110 !定义剪切强度
F1=1/XT-1/XC
F2=1/YT-1/YC
F11=1/(XT*XC)
F22=1/(YT*YC)
F66=1/(S**2)
F12=(-1/2)*sqrt(F11*F22)
*dim,a,array,N !定义系数a
*Dim,b,array,N !定义系数b
*dim,R,array,N !定义各层的强度比
Rmin=1000000 !定义最小强度比
RminN=0 !定义最小强度比层号&
!****************************************************************
!*Define the parameter of lamia
!****************************************************************
V12=E2*V21/E1
VV=1/(1.0-V21*V12)
*dim,Q11,array,N !定义退化刚度系数
*dim,Q22,array,N
*dim,Q12,array,N
*dim,Q66,array,N
*dim,HQ11,array,N !定义偏轴刚度系数
*dim,HQ12,array,n
*dim,HQ16,array,n
*dim,HQ22,array,n
*dim,HQ26,array,n
*dim,HQ66,array,n
*dim,z,array,N+1
*do,i,1,nint(N/2+1-0.5)
z =-H/2+(i-1)*Tk
z(N+2-i)=-z
*do,i,1,n !给退化刚度系数赋初值
Q11 =VV*E1
Q22 =VV*E2
Q12 =V21*VV*E2
!*****************************************************************
!***********************进入循环计算层合板强度********************
!*****************************************************************
*do,J,1,2*n+1
*if,over,NE,0,then !判断层合板是否失效
RM=COS(ARG )
RN=SIN(ARG )
RM4=RM2*RM2
RN4=RN2*RN2
RMN2=RMN*RMN
HQ11 =Q11 *RM4+2.0*(Q12 +2.0*Q66 )*RMN2+Q22 *RN4
HQ12 =(Q11 +Q22 -4.0*Q66 )*RMN2+Q12 *(RM4+RN4)
HQ16 =-RMN*RN2*Q22 +RM2*RMN*Q11 -RMN*(RM2-RN2)*(Q12 +2.0*Q66
HQ22 =Q11 *RN4+2.0*(Q12 +2.0*Q66 )*RMN2+Q22 *RM4
HQ26 =-RMN*RM2*Q22 +RMN*RN2*Q11 +RMN*(RM2-RN2)*(Q12 +2.0*Q66
HQ66 =(Q11 +Q22 -2*Q12 )*RMN2+Q66 *(RM2-RN2)*(RM2-RN2)
!*****************************************************************
!*Dedine Z Coordinate of Each Lamia
!*****************************************************************
!*****************************************************************
!* Define The Stiffness of Laminate
!*****************************************************************
*Do,i,1,N !积分求解层合板的刚度系数
A11=A11+HQ11 *(z(i+1)-z )
A12=A12+HQ12 *(z(i+1)-z )
A22=A22+HQ22 *(z(i+1)-z )
A26=A26+HQ26 *(z(i+1)-z )
A16=A16+HQ16 *(z(i+1)-z )
A66=A66+HQ66 *(z(i+1)-z )
B11=B11+0.5*HQ11 *(Z(i+1)**2-z **2)
B12=B12+0.5*HQ12 *(Z(i+1)**2-z **2)
B22=B22+0.5*HQ22 *(Z(i+1)**2-z **2)
B26=B26+0.5*HQ26 *(Z(i+1)**2-z **2)
B16=B16+0.5*HQ16 *(Z(i+1)**2-z **2)
B66=B66+0.5*HQ66 *(Z(i+1)**2-z **2)
D11=D11+1/3*HQ11 *(z(i+1)**3-z **3)
D12=D12+1/3*HQ12 *(z(i+1)**3-z **3)
D22=D22+1/3*HQ22 *(z(i+1)**3-z **3)
D16=D16+1/3*HQ16 *(z(i+1)**3-z **3)
D26=D26+1/3*HQ11 *(z(i+1)**3-z **3)
D66=D11+1/3*HQ11 *(z(i+1)**3-z **3)
!*****************************************************************
!* Define The roughness of Laminate
!*****************************************************************
DA=(A11*A22*A66+2*A12*A26*A16-A16**2*A22) !求解层合板的柔度系数
DA=DA-A11*A26**2-A12**2*A66
HA11=(A22*A66-A26**2)/DA
HA12=(A16*A26-A12*A66)/DA
HA22=(A11*A66-A16**2)/DA
HA66=(A11*A22-A12**2)/DA
HA16=(A12*A26-A12*A16)/DA
HA26=(A12*A16-A11*A26)/DA
!*****************************************************************
!* Caculate The Strain and stress of each Lamina
!*****************************************************************
STRAINX=HA11*PX+HA12*PY+HA16*PXY !计算层合板的应变
STRAINY=HA12*PX+HA22*PY+HA26*PXY
STRAINXY=HA16*PX+HA26*PY+HA66*PXY
strainstep(j)=STRAINX !记录下各个阶段的应变
pxstep(j)=PX !记录下各个阶段的应力
*Do,i,1,N !计算各层的主应变
RM=COS(ARG )
RN=SIN(ARG )
strain1 =RM2*STRAINX+RN2*STRAINY+(RMN)*STRAINXY
Strain2 =RN2*STRAINX+RM2*STRAINY+(-RMN)*STRAINXY
STrain12 =-2*RMN*STRAINX+2*RMN*STRAINY+(RM2-RN2)*STRAINXY
S1 =strain1 *Q11 +strain2 *Q12 !计算各层的主应力
S2 =STrain1 *Q12 +STRAIN2 *Q22
s12 =strain12 *Q66
a =F11*(s1 **2)+F22*(S2 **2)
a =a +2*F12*s1 *s2 +F66*(S12 **2)
b =F1*S1 +F2*S2
R =-b /(2*a )+sqrt(b **2+4*a )/(2*a )
*Vscfun,Rmin,min,R(1) !求出最小强度比
*Vscfun,RminN,lmin,R(1) !求出最小强度所在的层号
&&&&&&&&&&&&&&&&&&&&&&
*if,s1(RminN)*(Rmin),gt,0,then !判断破坏层的失效模式
*if,s1(RminN)*(Rmin),lt,xt,then !并进行刚度退化
Q22(RminN)=Q22(RminN)*(10e-14)
Q12(RminN)=Q12(RminN)*(10e-14)
Q66(RminN)=Q66(RminN)*(10e-14)
Q22(RminN)=Q22(RminN)*(10e-14)
Q12(RminN)=Q12(RminN)*(10e-14)
Q66(RminN)=Q66(RminN)*(10e-14)
Q11(RminN)=Q11(RminN)*(10e-14)
*ELSEIF,s1(RminN)*(Rmin),le,0,then
*if,(-1)*s1(RminN)*(Rmin),lt,xc,then
Q22(RminN)=Q22(RminN)*(10e-14)
Q12(RminN)=Q12(RminN)*(10e-14)
Q66(RminN)=Q66(RminN)*(10e-14)
Q22(RminN)=Q22(RminN)*(10e-14)
Q12(RminN)=Q12(RminN)*(10e-14)
Q66(RminN)=Q66(RminN)*(10e-14)
Q11(RminN)=Q11(RminN)*(10e-14)
*do,i,1,n !检查是否还有其他层破坏
*if,R ,eq,Rmin,then
Q11 =Q11(RminN)
Q22 =Q22(RminN)
Q12 =Q12(RminN)
Q66 =Q66(RminN)
PX=PX*(RMIN) !外载荷按比例增加继续加载
PY=PY*(RMIN)
PXY=PXY*(RMIN)
*IF,Rmin,lt,1,THEN !判断层合板是否能继续承载
*else !全部单层完全破坏，退去循环
58．钢筋混凝土整体式模型例子
! 单元属性
ET,1,SOLID65
KEYOPT,1,1,0
KEYOPT,1,5,0
KEYOPT,1,6,0
KEYOPT,1,7,1
!实参数1：未配筋
R,1,2, , , ,2, ,
RMORE, , ,2, , , ,
!实参数2：X方向配筋
R,2,2,.05, , ,2, ,
RMORE, , ,2, , , ,
!混凝土基本材料属性
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,30e9
MPDATA,PRXY,1,,.2
TB,MISO,1,1,5,
TBPT,,0.001,21e6
TBPT,,0.002,27e6
TBPT,,0.003,24e6
TB,CONC,1,1,9,
TBDATA,,.5,.9,3e6,30e6,,
TBDATA,,,,1,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
!钢材基本属性
MPDATA,EX,2,,200e9
MPDATA,PRXY,2,,.27
TB,BISO,2,1,2,
TBDATA,,310e6,2e9,,,,
KGEN,2,1,2 ,1 , ,.05, , ,0
KGEN,2,1,2 ,1 , ,0.5, , ,0
VEXT,1,2 ,1 ,0,0,.2,,,,
VSEL, , , , 1
VATT, 1, 2, 1, 0
VSEL, , , , 2
VATT, 1, 1, 1, 0
LESIZE,all,.2, , , , , , ,1
!力二范数收敛，2%
CNVTOL,F, ,0.02,2, ,
D,3, , , , , ,UY, , , , ,
D,35, , , , , ,UY, , , , ,
D,16, , , , , ,UY, , , , ,
D,48, , , , , ,UY, , , , ,
D,35, , , , , ,ALL, , , , ,
D,48, , , , , ,UZ, , , , ,
!加竖向位移荷载
D,75, ,-0.01, , , ,UY, , , , ,
D,123, ,-0.01, , , ,UY, , , , ,
!加载最少分为20步
NSUBST,20,
!输出每步结果
OUTRES,ALL,1
!每步至少迭代25次
CUTCONTROL,NOITERPRED,1
59．在荷载步之间改变材料属性例子
! 材料泊松比随荷载增加而逐步增大
FORCE=1. !初始荷载
FC=30. !极限荷载
NSTEP=30 !加载步数
EMU0=0.2 !初始泊松比为0.2
EMUU=0.499 !最终泊松比为0.499
SVM=0. !VON MISES应力
ET,1,SOLID45
MP,EX,1,30E3
MP,NUXY,1,EMU0
BLC4,0,0,100,100,100
ESIZE,100,0,
!输出RESTART文件
RESCONTRL,DEFINE,ALL,-1,1
*DO,I,1,NSTEP
!使用重启动功能
*IF,I,GT,1,THEN
ANTYPE,,REST,
PARRES, CHANGE , PARAM, TXT,
! 如果荷载超过强度的50%，则线性提高泊松比
*IF,SVM,GE,FC*0.5,THEN
MP,EX,1,30E3
MP,NUXY,1,EMU0+(EMUU-EMU0)*(SVM/FC-0.5)/0.5
!得到下一步荷载
FORCE=FORCE+1
SFE,ALL,4,PRES, , FORCE, , ,
!得到VON MISES应力
*GET,SVM,ELEM,1,NMISC, 4
PARSAV, ALL, PARAM, TXT,
62．柱子稳定分析算（预应力，特征值屈曲，初始缺陷）
! 要点：预应力，特征值屈曲，添加初始缺陷，几何非线性分析
OFFSET=0.1 !初始缺陷为0.1
! 建立模型
ET,1,BEAM4
ET,2,LINK10
R,1,0.1*0.12,0.12*0.1**3/12,0.1*0.12**3/12,0.12,0.1, ,
R,2,0.002*0.002,2e-3, !预应力
MPTEMP,1,0
MPDATA,EX,1,,200e9
MPDATA,PRXY,1,,0.27
MPDATA,DENS,1,,7800
k,3,0,0,-5
k,11,0.2,0,0
k,12,-0.2,0,0
k,13,0,0.2
k,15,0,-0.2
lsel,,,,1,6
latt,1,1,1
ALLSEL,ALL
lsel,,,,7,14
latt,1,2,2
ALLSEL,ALL
lsel,,,,1,6
LESIZE,all,0.3, , , , , , ,1
lsel,,,,7,14
LESIZE,all, , ,1 , , , , ,1
ALLSEL,ALL
/ESHAPE,1.0
! 求解特征值屈曲荷载
*DO,I,1,100
DK,3, , , ,0,UX,UY,UZ, , ,ROTZ ,
DK,2, , , ,0,UX,UY, , , , ,
FK,2,FZ,-FORCE
!设定时间步
NSUBST,1, , ,1
ANTYPE,BUCKLE ! Buckling analysis
BUCOPT,LANB,1 ! Use Block Lanczos solution method, extract 1
MXPAND,1 ! Expand 1 mode shape
PSTRES,ON ! INCLUDE PRESTRESS EFFECTS
*GET,FREQ1,MODE,1,FREQ
*IF,FREQ1,LT,0.01,THEN
FREQ1=0.01
*IF,ABS(FREQ1-1),LT,0.01,THEN !如果频率误差小于1%
FORCE=FORCE*FREQ1
! 添加初始缺陷
NSORT,U,SUM
*GET,DMAX,SORT, ,MAX
upgeom,OFFSET/DMAX,1,1,beammodel,rst
! 非线性求解
FORCE=FORCE*3.
DK,3, , , ,0,UX,UY,UZ, , ,ROTZ ,
DK,2, , , ,0,UX,UY, , , , ,
FK,2,FZ,-FORCE
NSUBST,20,0,0
OUTRES,ALL,1
ARCLEN,1,0,0
AUTOTS,-1.0
! 绘制荷载位移曲线
FILE,'beammodel','rst','.'
/UI,COLL,1
NUMVAR,200
FILLDATA,191,,,,1,1
REALVAR,191,191
NUMVAR,200
FILLDATA,191,,,,1,1
REALVAR,191,191
NSOL,2,2,U,Z, UZ_2
STORE,MERGE
RFORCE,3,19,F,Z, FZ_3
STORE,MERGE
63． module MConcrete !混凝土模板
implicit none
& type::typ_Concrete
!混凝土抗拉强度,抗压强度,初始弹性模量,初始泊松比
&& !最终泊松比,初始剪切模量
Ft,Fc,E0,MU0,MUU,G0
&& !抗拉下降段参数,裂面剪力折减系数
&& real*8 A1,A2
!Crack=1,一条裂缝,=2两条裂缝,AddLoad=1加载,=0,卸载
&& integer*4 Crack, AddLoad
&& !裂缝角度
&& real*8 ANGLE
&& !t时的应力,主应力,应力增量,t+dt时的应力
SIG(3),SIGP(3),dSIG(3),Stress(3)
EPS(3),EPSP(3),dEPS(3),Strain(3)
&& !非线性指标, 破坏面,最大非线性指标
&& real*8 Beta,J2f,BetaMax
&& !弹性本构矩阵,割线本构矩阵,本构矩阵
&& real*8 De(3,3), Ds(3,3),
&& !坐标转换矩阵
&& real*8 N(3,3)
&& integer(4) INC, NCycle
& end type typ_Concrete
& subroutine Con_Initial(C) !初始化混凝土参数
67．导出刚度矩阵
其原理很简单，即使用ansys的超单元即可解决问题。定义超单元，然后列出超单元的刚度矩阵即可。
面是一个小例题，自可明白。
et,1,beam3
mp,ex,1,2e5
mp,prxy,1,0.3
lesize,all,,,10
!----以上正常建立模型，不必施加约束和荷载
antype,7 !substructuring分析类型
seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)
nsel,all !选择所有节点
m,all,all !定义所有节点自由度为主自由度
solve !求解
selist,matname,3 !列出整体刚度矩阵
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。