ANSYS动力学分析Word下载.docx
《ANSYS动力学分析Word下载.docx》由会员分享,可在线阅读,更多相关《ANSYS动力学分析Word下载.docx(21页珍藏版)》请在冰豆网上搜索。
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,,,,1,
K,,,1,1,
K,,-0.5,1,,
K,,0.5,1,,
K,,,,0。
5,
A,4,5,6
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
LMESH,4,6
TYPE,2
MSHAPE,0,2D
MSHKEY,1
AMESH,1
FINISH
/SOL!
取前十阶模态
ANTYPE,2
MODOPT,LANB,10
KSEL,S,KP,,1,3
DK,ALL,,,,0,ALL,,,,,,
SOLVE
谱分析
ANTYPE,8
SPOPT,SPRS,10,1…………0是no,1是yes包括应力响应
SVTYP,3
SED,0,1,0,!
给出激励方向
FREQ,0.25,1。
34,6。
73,14.6,28.9,0,0,0,0!
激励谱
SV,0,0。
00073,0。
00016,0。
00034,0.00034,0.00052,
扩展模态
EXPASS,1
MXPAND,10,0,0,1,0。
001,
模态合并
SRSS,0.001,DISP
/POST1!
观察结果
SET,LIST
/INPUT,,mcom
PRRSOL,F
问题一:
悬臂梁受重力作用发生大变形,求其固有频率。
基本过程:
2、静力分析
 
 
NLGEOM,ON
&
nbsp;
nbsp;
PSTRES,ON
3、求静力解
4、开始新的求解:
modal
PSTRES,ON
UPCOORD,1,ON修正坐标
PSOLVE。
..
5、扩展模态解
6、察看结果
/PREP7
ET,1,BEAM189!
使用beam189梁单元
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,210e9
MPDATA,PRXY,1,,0。
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,,,,,!
建模与分网
K,,2,,,
K,,2,1,,
LSTR,1,2
LATT,1,,1,,3,,1
LESIZE,1,,,20,,,,,1
LMESH,1
/SOL!
静力大变形求解
ANTYPE,0
NLGEOM,1
PSTRES,ON!
计及预应力效果
DK,1,,,,0,ALL,,,,,,
ACEL,0,9.8,0,!
只考虑重力作用
TIME,1
AUTOTS,1
NSUBST,20,,,1
KBC,0
/SOLUTION
ANTYPE,2!
进行模态求解
MSAVE,0
MODOPT,LANB,10
MXPAND,10,,,0!
PSTRES,1!
打开预应力效应
MODOPT,LANB,10,0,0,,OFF
UPCOORD,1,ON!
修正坐标以得到正确的应力
PSOLVE,TRIANG!
三角化矩阵
PSOLVE,EIGLANB!
提取特征值和特征向量
/SOLU
EXPASS,1!
扩展模态解
PSOLVE,EIGEXP
/POST1
SET,LIST!
问题二:
循环对称结构模态分析
这是ANSYSHELP里的例子,但那个命令流似乎有些问题,下面是整理过的命令流.
2、definecyclic
3、定义约束
4、求模态解
5、展开并察看结果
r1=5!
r2=10
d1=2
nsect=24
alpha_deg=360/nsect
alpha_rad=2*acos(-1)/nsect
/prep7
csys,1
k,1,0,0,0
k,2,0,0,d1
k,3,r1,0,0
k,4,r1,0,d1
l,3,4
arotat,1,,,,,,1,2,alpha_deg/2
k,7,r2,0,0
k,8,r2,0,d1
l,7,8
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
mshkey,1
et,1,181
r,1,0。
20
r,2,0.1
mp,ex,1,10e6
mp,prxy,1,0.3
mp,dens,1,1e—4
esize,0.5
asel,,,,1,4
aatt,,1
asel,,,,5
aatt,,2
allsel
CYCLIC,,,,'
CYCLIC'
!
定义基本扇区
AMESH,all!
分网
求模态解
MODOPT,LANB,10
MXPAND,10,,,0
PSTRES,0
MODOPT,LANB,10,0,0,,OFF
DK,5,,,,0,UZ,,,,,,
SET,LIST
SET,,,1,,,,4,
/CYCEXPAND,,ON
PLNS,U,SUM!
观察扩展模态
问题三:
三角平台受谐波载荷作用,求结构响应
谐波载荷为作用在平台上面一点的集中力,幅值为10N,频率范围5Hz~70Hz
2.求模态解
3、用模态叠加法作谐分析
观察结果
/PREP7!
ET,1,BEAM189
ET,2,SHELL93
R,1,0。
01,,,,,,
MP,EX,1,210e9
MP,PRXY,1,0。
MP,DENS,1,7850
SECDATA,0.005,0.008,0,0,0,0,0,0,0,0
K,,-0。
5,,,
K,,0.5,,,
K,,,,1,
K,,,1,1,
K,,—0.5,1,,
K,,0.5,1,,
K,,,,0.5,
A,4,5,6
LSTR,2,6
LSEL,S,LINE,,4,6
LATT,1,1,1,,7,,1
LSEL,,,,ALL
LESIZE,ALL,,,10,,,,,1
LMESH,4,6
TYPE,2
MSHAPE,0,2D
AMESH,1
为了使用模态叠加法谐分析
ANTYPE,2
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!
以了解模态频率范围
模态叠加法谐响应分析
ANTYPE,3
HROPT,MSUP
HROUT,ON
HROPT,MSUP,20,,
HROUT,ON,OFF,0
F,177,FY,-10,!
施加10N的力在节点177的负Y方向
HARFRQ,5,70,!
载荷的强制频率范围
NSUBST,200,
KBC,1
ALPHAD,5,
/POST26!
察看位移响应
FILE,,rfrq
NUMVAR,20
NSOL,2,139,U,Y,
PLVAR,2
问题四:
/PREP7!
ET,1,BEAM189
ET,2,SHELL93
R,1,0。
01,,,,,,
SECTYPE,1,BEAM,RECT,secA,0
SECDATA,0.005,0。
008,0,0,0,0,0,0,0,0
K,,—0.5,,,
K,,0。
K,,,,1,
K,,,1,1,
K,,—0。
5,1,,
5,1,,
K,,,,0.5,
A,4,5,6
LSTR,1,5
LSEL,S,LINE,,4,6
LATT,1,1,1,,7,,1
LSEL,,,,ALL
LESIZE,ALL,,,10,,,,,1
LMESH,4,6
MSHAPE,0,2D
MSHKEY,1
KSEL,S,KP,,1,3
DK,ALL,,,,0,ALL,,,,,,
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.00073,0。
00016,0。
00034,0.00034,0。
00052,
MXPAND,10,0,0,1,0。
ANTYPE,8
/INPUT,,mcom
PRRSOL,F
问题五:
三角平台受时程载荷作用,求应力分布和变形过程
2.施加随时间变化载荷,定义载荷步
3、求解
4.在POST1和POST26中观察结果
R,1,0。
01,,,,,,
MP,EX,1,210e9
MP,PRXY,1,0.3
MP,DENS,1,7850
SECTYPE,1,BEAM,RECT,secA,0
SECDATA,0。
005,0。
008,0,0,0,0,0,0,0,0
5,,,
K,,,,1,
K,,,1,1,
K,,—0.5,1,,
5,1,,
K,,,,0.5,
A,4,5,6
LSTR,3,4
LSTR,2,6
LSEL,S,LINE,,4,6
LATT,1,1,1,,7,,1
LSEL,,,,ALL
LESIZE,ALL,,,10,,,,,1
/SOL
ANTYPE,4
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,
TIME,3
LSWRITE,2,
TIME,4
SFA,ALL,1,PRES,150
LSWRITE,3,
TIME,5
SFA,ALL,1,PRES,
LSWRITE,4,
LSSOLVE,1,4,1,!
求解
观察变形随时间的变化
NSOL,2,177,U,Y,uyy
PLVAR,2
素混凝土圆柱体轴压试验数值模拟LOG
ANSYS2009-10-1914:
27阅读12评论0字号:
大大中中小小/PREP7
ET,1,SOLID65,0,,,,1,3,!
材料1:
混凝土
KEYOPT,1,1,1!
不包含大位移
KEYOPT,1,3,0
KEYOPT,1,5,1
KEYOPT,1,6,3
KEYOPT,1,7,1!
裂后应力不释放
KEYOPT,1,8,0
MP,EX,1,2。
9397E4!
砼弹模
MP,NUXY,1,0。
2!
砼泊松比
MP,DENS,1,2.5E-6!
砼密度
TB,CONCR,1,,,,!
定义混凝土五参数破坏准则
TBDATA,1,0.5,0。
8,2。
51,32。
7,,,!
可以关闭压碎将32。
7改为—1
TBDATA,7,,,0。
7
TB,MKIN,1!
定义混凝土等效单轴压曲线
fc=32。
7!
峰值应力
Strnmax=0。
002!
峰值应变
NPOIN=4!
描述曲线上升段的点数,不包括原点(0,0)
*DIM,XE,ARRAY,NPOIN!
strain/峰值应变
*DIM,Strn,ARRAY,NPOIN!
strain
*DIM,Strs,ARRAY,NPOIN!
stress
XE
(1)=0。
3,0。
5,0。
8,1。
*VOPER,Strn
(1),XE
(1),MULT,Strnmax
TBMODIF,1,1,0!
曲线原点
TBMODIF,2,1,0
*DO,I,1,NPOIN
Strs(I)=fc*(2.2*XE(I)-1。
4*XE(I)**2+0.2*XE(I)**3)!
《规范》采用的曲线
TBMODIF,1,I+1,Strn(I)
TBMODIF,2,I+1,Strs(I)
*ENDDO
TBMODIF,1,NPOIN+2,0。
004!
曲线下降段
TBMODIF,2,NPOIN+2,fc*0.85
TBPLOT,MKIN,1
R,1,
et,2,mesh200!
辅助平面网格划分
KEYOPT,2,1,6!
4nodequadelement
ET,3,SOLID45
MP,EX,2,2。
0E5!
钢板弹模
MP,NUXY,2,0.2!
钢板泊松比
参数定义
RR=75!
圆柱半径mm
HH=300!
圆柱高mm
THICK=10!
钢垫板厚度mm
建模
K,1,0,0,0
K,2,RR,0,0
K,3,0,RR,0
L,1,2
L,1,3
CIRCLE,1,RR,,,90,,
NUMMRG,ALL,,,,LOW
NUMCMP,ALL
ALLSEL,ALL
AL,ALL
lsel,s,,,1,2
LESIZE,all,,,10,,,,,1!
分10份
lsel,s,,,3
LESIZE,all,,,20,,,,,1!
分20份
type,2
allsel,all
amesh,all
由面拉伸成圆柱体
TYPE,1
MAT,1
EXTOPT,ESIZE,20,0,!
20份
EXTOPT,ACLEAR,1
VEXT,all,,,0,0,HH,,,,!
拉伸长度HH
需要拉伸作垫板的面(A1,A2)
TYPE,3
MAT,2
EXTOPT,ESIZE,1,0,!
1份
EXTOPT,ACLEAR,1
asel,s,,,1
VEXT,all,,,0,0,—THICK,,,,!
下垫板厚THICK
EXTOPT,ESIZE,1,0,!
asel,s,,,2
VEXT,all,,,0,0,THICK,,,,!
上垫板厚THICK
施加对称面约束
asel,s,loc,x,0。
DA,all,SYMM
asel,s,loc,y,0.0
DA,all,SYMM
施加底面钢板约束(z=-THICK)
asel,s,loc,z,—THICK
aplot
da,all,all,0
将面约束转移到有限元模型上
DTRAN
求解过程
allsel,all
/solu
在顶面施加面力(A1)
ASEL,S,LOC,Z,HH+THICK
SFA,ALL,1,PRES,100!
均布荷载=100N
SFTRAN!
将面荷载转化到有限元模型上
非线性参数设置
antype,static
非线性参