ANSYS在结构中的仿真模拟与分析Word下载.docx
- 文档编号:22457447
- 上传时间:2023-02-04
- 格式:DOCX
- 页数:45
- 大小:919.53KB
ANSYS在结构中的仿真模拟与分析Word下载.docx
《ANSYS在结构中的仿真模拟与分析Word下载.docx》由会员分享,可在线阅读,更多相关《ANSYS在结构中的仿真模拟与分析Word下载.docx(45页珍藏版)》请在冰豆网上搜索。
LESIZE,ALL,2
LMESH,ALL
LSEL,U,,,ALL
建立外环梁模型
L,61,65
L,65,89
L,89,92
L,92,68
L,68,72
L,72,120
L,120,109
L,109,61
LGEN,10,161,168,1,,,3,60
LATT,1,,1,,4000,,2
建立内框架梁模型
L,63,111
L,85,89
L,89,113
L,70,118
L,92,96
L,92,116
LGEN,10,241,246,1,,,3,60
LATT,1,,1,,3000,,3
!
建立次梁模型
L,62,110
L,64,112
L,69,117
L,71,119
L,97,108
L,73,77
L,80,84
L,90,102
L,91,103
LGEN,10,301,309,1,,,3,60
LATT1,,1,,3000,,4
ALLS
建立楼层板模型
*DO,I,61,64
A,I,I+1,I+13,I+12
AGEN,4,1,4,1,,4,,12
*DO,I,68,71,1
AGEN,4,17,20,1,,4,,12
*DO,I,89,91,1
AGEN,10,1,35,1,,,3,60
A,653,656,644,641
AATT,2,2,2
AESIZE,ALL,2
AMESH,ALL
ASEL,U,,,ALL
建立外墙模型并划分网格
*DO,I,61,64,1
A,I,I+1,I+61,I+60
*ENDDO
*DO,I,109,112,1
*DO,I,116,119,1
*DO,I,61,97,12
A,I,I+12,I+72,I+60
*DO,I,65,77,12
*DO,I,68,80,12
*DO,I,72,108,12
AGEN,9,352,382,,,,3,60
AATT,3,2,2
建立筒体模型
A,53,56,656,653
A,41,44,644,641
AATT,1,3,2
ESIZE,1,0
MSHAPE,0,2D
MSHKEY,1
ASEL,U,,,1
A,41,53,653,641
A,44,56,656,644
NUMMRG,ALL
NUMCMP,ALL
ALLSEL
SAVE
FINISH
至此,模型已经全建立完毕
进入求解模块
/SOLU
ANtYPE,STATIC
DELTIM,0.1,0.05,0.2
AUTOTS,ON
NLGEOM,ON
PRED,ON
LNSRCH,ON
NSEL,S,LOC,Z,0
D,ALL,ALL
ACEL,,,10
开始施工过程的模拟
ESEL,ALL
EKILL,ALL
NSEL,ALL
D,ALL,ALL
首先修建柱子和梁,同事再修建楼板
将柱、梁和楼板所属的单元激活,并删除其全部约束
ESEL,S,TYPE,,1
NSLE,S
NSEL,R,LOC,Z,0,3
ESLN,R,1
EALIVE,ALL
DDELE,ALL,ALL
ESEL,S,MAT,,2
NSEL,R,LOC,Z,0
SOLVE
激活第二至第十层的柱、梁和楼板单元
*DO,I,2,10,1
NSEL,R,LOC,Z,3*(I-1),3*I
至此,梁、柱和楼板施工完毕
以下进行外墙及隔墙的施工
ESEL,S,MAT,,3
ESEL,S,MAT,,1
NSEL,S,LOC,Z,0
激活第二至第十层的墙面单元
求解过程结束
进入通用后处理器查看结果
/POST1
PLNSOL,U,X
PLNSOL,U,y
PLNSOL,U,z
PLNSOL,S,X
PLNSOL,S,Y
PLNSOL,S,Z
查看柱子和梁的轴力,剪力和弯矩
ETABLE,IF,SMISC,1!
定义轴力单元列表1
ETABLE,JF,SMISC,7!
定义轴力单元列表2
ETABLE,II,SMISC,2!
定义剪力单元列表1
ETABLE,JI,SMISC,8!
定义剪力单元列表2
ETABLE,IM,SMISC,6!
定义弯矩单元列表1
ETABLE,JM,SMISC,12!
定义弯矩单元列表2
PLETAB,IF,NOAV!
显示轴力单元列表1
PLETAB,JF,NOAV!
显示轴力单元列表2
PLETAB,II,NOAV!
显示剪力单元列表1
PLETAB,JI,NOAV!
显示剪力单元列表2
PLETAB,IM,NOAV!
显示弯矩单元列表1
PLETAB,JM,NOAV!
显示弯矩单元列表2
2.斜拉桥三维仿真分析
/filname,cable-stayedbridge,1
keyw,pr_struc,1
/prep7
定义单元类型
et,1,beam4
et,2,link10
定义材料属性
mp,ex,1,3.5e10
mp,prxy,1,0.17
mp,dens,1,2500
mp,ex,2,10e15
mp,prxy,2,0
mp,dens,2,0
mp,ex,3,1.9e10
mp,prxy,3,0.25
mp,dens,3,1200
mp,damp,3,0.5
定义实常数
r,1,25.6,5.46,546.133,16,1.6
r,2,16,29.417,15.394,3.4,4.7
r,3,54,364.5,162,6,9
r,4,40,213.3,83.3,5,8
r,5,1,1/12,1/12,1,1
r,6,0.012,0.012!
索的
创建节点和单元
建立主梁节点
/view,1,1,1,1
/angle,1,270,xm,0
/replot
*do,i,1,59!
此循环用于建立主梁的半跨节点
x=-174*2+(i-1)*6!
最左端x=174*2,x=0左边的节点x坐标值,间距为6
y1=-14!
桥面宽28米,故左边节点为-14
y2=14!
桥面宽28米,故右边节点为-14
n,3*(i-1)+1,x!
建立主梁节点3*(i-1)+1为节点号
n,3*(i-1)+2,x,y1!
以下两行建立桥面两边节点
n,3*i,x,y2!
能想出这种建模命令的绝对是编程高手,哈哈
*enddo!
完全可以先建立端部的三个节点,然后用这三个节点在x方向上复制59份,间距为6
建立主梁单元
type,1
real,1
mat,1
*do,i,1,58,1!
以下循环建立建立桥面中线主梁单元
j=3*(i-1)+1
e,j,j+3
*enddo
建立鱼刺刚横梁
real,5
mat,2
*do,i,1,59,1!
以下循环用于建立桥面鱼刺横梁的节点
j1=3*(i-1)+2
j2=3*i
e,j,j1
e,j,j2
建立半跨主塔
i=59*3!
变量用于记录桥面的节点数,即至此已经建立了59*3个节点了,用于指导以后设定节点的编号
n,i+1,-174,-10,-30!
以下两行记录塔脚节点
n,i+2,-174,10,-30
n,i+3,-174,-15!
以下两行用于建立与桥面齐高的主塔节点
n,i+4,-174,15
*do,j,1,5,1!
以下循环用于建立索塔在桥面以上的节点
k=i+4+j
n,k,-174,0,60+(j-1)*18
建立下索塔单元
real,4
e,i+1,i+3!
以下用于建立主塔在桥面以下的两根塔柱单元
e,i+2,i+4
建立中索塔单元
real,3
e,i+3,i+5!
以下用于建立倒Y分叉点到桥面间的两根塔柱单元
e,i+4,i+5
建立上索塔单元
real,2
*do,j,1,4,1!
以下用于建立倒Y分叉点以上的塔柱单元
e,k,k+1
建立与塔的倒Y分叉点链接的索单元
type,2
real,6
mat,3
e,i+5,89
e,i+5,90
建立主塔倒Y分叉点以上第一个张拉点连接的索单元
*do,j,1,8,1
此循环用于建立主塔倒Y分叉点以上第一个张拉点连接的所有索单元,共32个
e,i+6,89+3*j
e,i+6,89-3*j
e,i+6,90+3*j
e,i+6,90-3*j
建立与主塔的其他三个张拉点连接的单元
*do,k,1,3,1
*do,j,1,7,1
e,i+6+k,113+(k-1)*21+3*j!
一共有28个索单元连接在每个张拉点上
e,i+6+k,65-(k-1)*21-3*j
e,i+6+k,114+(k-1)*21+3*j
e,i+6+k,66-(k-1)*21-3*j
生成全桥模型节点
i=i+9!
记录半跨的所有节点数
nsym,x,i,all!
用映射法直接建立另半跨节点
esym,,i,all!
用映射法直接建立另半跨单元
nummrg,all!
合并所有节点和单元
建立索塔连接横梁单元
j=i
i=i-9
n,1000,-174
e,1000,i+3
e,1000,i+4
n,2000,174
e,2000,i+3+j
e,2000,i+4+j
施加主塔的四个脚上的全约束
nsel,s,loc,z,-30
d,all,all
allsel
在左桥端施加y,z约束
nsel,s,loc,x,-348!
仅给左端主梁施加约束
nsel,r,loc,y,0
d,all,uy
d,all,uz
在右桥端施加y约束
nsel,s,loc,x,348!
仅给右端主梁施加约束
numcmp,all
施加重力场
acel,,,9.8
耦合节点,耦合跨中由于对称而重复的单元节点以及两主塔上塔横梁和主梁的重合节点,
cpintf,uy
cpintf,uz
cpintf,rotx
cpintf,rotz
成桥状态的确定
静力的初步计算
直接进行静力计算
/solu
solve
finish
得到最大位移为1.288m,由于偏差太大需要重新计算(与事实不符)
修改实常数后重新计算:
令r,6,0.012,0.012,即给索以预应变0.012
/solu!
为了将计算应力用于下面的动力分析,这里打开预应力和集中质量设置开关
lumpm,on
pstres,on
求得最大位移为0.0329m,说明已经达到成桥状态要求,(与事实不符)
模态分析
分析设置
antype,2
MODOPT,LANB,20
EQSLV,SPAR
MXPAND,10,,,0
LUMPM,1
PSTRES,1
MODOPT,LANB,20,0,100,,OFF
UPCOORD,1,ON!
更新模型计算坐标,目的是为了将预应力效应准确的应用能够到模态分析上来
分析设置完毕
激励耦合分析
恢复成桥静力计算结果
RESUME,cable-stayedbridge,db,!
路径:
utilitymenu>
file>
resumejobname.db
生成地震激励向量,将记事本格式的地震波数据调入到工作目录下,并执行以下命令
*DIM,aay1,ARRAY,2,50,1
*DIM,aaz1,ARRAY,2,50,1
*CREATE,ansuitmp
*VREAD,aay1(1,1),'
tjx'
'
txt'
'
50
(e9.3,e11.3)
*END
/INPUT,ansuitmp
*VREAD,aaz1(1,1),'
tjy'
/INPUT,ansuitmp
输入车辆激励波,Z=1000*cos(10*t).车子跨越一个单元的时间为0.1秒(速度216公里/小时),全桥共696米,历时11.6秒,纵向共116个单元,117个节点
n=117!
定义向量维数
*dim,fcar,array,n!
定义车载荷向量
*do,i,1,n,1!
以下循环为生成车载荷数据,即向“farc”矩阵中写入数据
tt=(i-1)*0.1
fcar(i)=1000*cos(10*tt)
输入风载荷激励p=50*sin(1.5*t)。
由于梁长6米,高3米,故每个节点的受力可表示为:
*dim,fwind,array,n,n!
定义风荷载向量,***注意要与车辆荷载步即时间相等***
*do,i,1,n,1
tt=(i-1)*0.1!
按不同时刻,每一次都是给列赋值,即各节点赋予荷载值,第a行第b列表示第b时刻第a个节点的荷载值
p=50*sin(1.5*tt)
fwind(1,i)=3*3*p!
两端的每根梁的面积为6*3,两个节点上所受的力应为一半面积与压力的乘积:
3*3*p
*do,j,1,n-2
fwind(j+1,i)=6*3*p!
中间的每根梁的面积为6*3,两个节点上所受的力应为两个相邻梁各一半面积与压力的乘积:
6*3*p
fwind(n,i)=3*3*p
加载步耦合
由于瞬态分析是以荷载步为阶段来求解的,所以应该将117个载荷步分别写入荷载文件中
finish!
完成上一阶段的设置
/config,nres,2000!
设置最大子步数为2000
antype,trans!
定义分析类型为瞬态
trnopt,full!
求解方式为完全法
m=n!
m变量表示向量的需要存储的时间点数m=n=117
*do,i,1,m,1!
总循环开始,开始激励的耦合循环
time,i*0.1!
给每个载荷步文件设置结束点时间
auto,off!
关闭自动时间步选项
deltim,0.05!
设置子步求解的时间间隔
kbc,0!
选择ramped选项
outres,basic!
选择基本结果输出
cnvtol,f,,0.01!
选择力收敛准则,容差为0.01
施加地震激励
*if,i,le,50,then!
如果i小于等于50(le为littleequalto)
acel,,-aay1(2,i),10-aaz1(2,i)!
其中最后一项为重力加速度与地震竖向加速度的差值
*else
acel,0,0,10
*endif
施加移动载荷激励
*if,i,eq,1,then!
如果力作用在最左端的主梁上时,单独定义施加的载荷,目的为了以下删除载荷的语句可以使用
nsel,s,loc,x,-174*2
nsel,r,loc,z,0
f,all,fz,-fcar(i)
*else!
如果力不是作用在最左端时
nsel,s,loc,x,-174*2+i*6!
选择受力节点施加载荷
nsel,s,loc,x,-174*2+(i-1)*6!
删除前一节点上的载荷
fdele,all,fz
*endif
施加风载荷,(子循环)
*do,j,1,m,1
nsel,s,loc,y,-14!
选择有风作用一侧的节点
nsel,r,loc,x,-174*2+(j-1)*6
f,all,fy,fwind(j,i)!
给所选节点施加载荷
allsel
*enddo
lswrite,i!
写数据到载荷步文件i
总循环结束
lssolve,1,m,1
3.用ANSYS对钢筋混凝土梁进行计算模拟
分离式模型
FINISH
/CLEAR
/PREP7
1.定义单元与材料属性
ET,1,SOLID65,,,,,,,1
ET,2,LINK8
MP,EX,1,13585!
混凝土材料的初始弹模以及泊松比
MP,PRXY,1,0.2
FC=14.3!
混凝土单轴抗压强度和单轴抗拉强度
FT=1.43
TB,CONCR,1
TBDATA,,0.5,0.95,FT,-1!
定义混凝土材料及相关参数,关闭压碎
TB,KINH,1,,11!
定义混凝土应力应变曲线,用KINH模型
TBPT,,0.0002,FC*0.19
TBPT,,0.0004,FC*0.36
TBPT,,0.0006,FC*0.51
TBPT,,0.0008,FC*0.64
TBPT,,0.0010,FC*0.75
TBPT,,0.0012,FC*0.84
TBPT,,0.0014,FC*0.91
TBPT,,0.0016,FC*0.96
TBPT,,0.0018,FC*0.99
TBPT,,0.002,FC
TBPT,,0.0033,FC*0.85
MP,EX,2,2.0E5!
钢筋材料的初始弹模以及泊松比
MP,PRXY,2,0.3
TB,BISO,2
TBDATA,,300,0!
钢筋的应力应变关系,用BISO模型
PI=ACOS(-1)
R,1
R,2,0.25*PI*22*22
R,3,0.25*PI*10*10
TBPLOT,MISO,1!
混凝土材料的数据表绘图
TBPLOT,BISO,2!
钢筋材料的数据表绘图
2.创建几何模型
BLC4,,,150,300,2000
*
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 结构 中的 仿真 模拟 分析