最新ansys屈曲分析练习模型.docx
《最新ansys屈曲分析练习模型.docx》由会员分享，可在线阅读，更多相关《最新ansys屈曲分析练习模型.docx（28页珍藏版）》请在冰点文库上搜索。
最新ansys屈曲分析练习模型
ansys屈曲分析练习模型：
边界条件：
底端固定
几何：
长为100mm，截面：
10mm×10mm
惯性矩：
Izz=833.333
材料性质：
E=2.0e5MPa,v=0.3
分析压力的临界值
分析过程：
特征值屈曲分析方法：
1、建立关键点1（000）,2（01000）
2、在关键点1、2之间建立直线
3、定义单元类型（Beam3）
4、定义单元常数
5、定义材料属性
6、定义网格大小，指定单元边长为10
7、划分网格
（首先此处应该做一次模态分析，有模态数据文件，后出来才可以看屈曲模态。
）
8、定义分析类型（static）
9、激活预应力效应。
要进行屈曲分析，必须激活预应力效应。
10、施加位移约束（关键点1固定）
11、施加集中荷载，Fy=-1N
12、求解
13、结束求解，
14、重新定义分析类型（EigenBuckling）
15、设置屈曲分析选项，提取1阶模态（菜单路径：
Solution-->AnalysisType-->Analysisoptions
16、求解，结束后退出
17、解的展开
1）设置expansionpass“on”
2）设置展开模态为1（LoadStepOptions>ExpansionsPass>SingleExpand>ExpandModes
3）重新求解
18、查看结果（临界载荷和屈曲模态等）
二、非线性分析方法
前8步与上述过程相同
9、设置分析控制（主要黄色高亮部分区域需要修改）
10、施加位移约束（关键点1固定）
11、施加集中荷载，Fy=-50000N，Fx=-250N
12、求解
13、查看变形和位移
14、定义时间-历史变量
1）进入时间历程后处理器（TimeHistPostproc）
2）在弹出的对话框中选择左上角的+号，添加一个监控变量（节点2的Y方向位移）
15、查看位移-载荷曲线
屈曲分析是一种用于确定结构开始变得不稳定时的临介荷载和屈曲结构发生屈曲响应时的模态形状的技术。
ANSYS提供两种结构屈曲荷载和屈曲模态分析方法：
非线性屈曲分析和特征值屈曲分析。
非线性屈曲分析是在大变形效应开关打开的情况下的一种非线性静力学分析，该分析过程一直进行到结构的极限荷载或最大荷载。
非线性屈曲分析的方法是，逐步地施加一个恒定的荷载增量，直到解开始发散为止。
尤其重要的是，要一个足够小的荷载增量，来使荷载达到预期的临界屈曲荷载。
若荷载增量太大，则屈曲分析所得到的屈曲荷载就可能不准确，在这种情况下打开自动时间步长功能，有助于避免这类问题，打开自动时间步长功能，ANSYS程序将自动寻找屈曲荷载。
特征值屈曲分析步骤为：
1.建模
2.获得静力解：
与一般静力学分析过程一致，但必须激活预应力影响，通常只施加一个单位荷载就行了
3.获得特征屈曲解:
A.进入求解
B.定义分析类型
C.定义分析选项
D.定义荷载步选项
E.求解
4.扩展解
之后就可以察看结果了
示例1：
!
ansys7.0有限元分析实用教程
!
3．命令流求解
!
ANSYS命令流：
!
EigenvalueBuckling
FINISH !
这两行命令清除当前数据
/CLEAR
/TITLE,EigenvalueBucklingAnalysis
/PREP7 !
进入前处理器
ET,1,BEAM3 !
选择单元
R,1,100,833.333,10 !
定义实常数
MP,EX,1,200000 !
弹性模量
MP,PRXY,1,0.3 !
泊松比
K,1,0,0!
创建梁实体模型
K,2,0,100
L,1,2!
创建直线
ESIZE,10 !
单元边长为1mm
LMESH,ALL,ALL !
划分网格
FINISH !
退出前处理
！
屈曲特征值部分
/SOLU!
进入求解
ANTYPE,STATIC !
在进行屈服分析之前，ANSYS需要从静态分析提取数据
PSTRES,ON !
屈服分析中采用预应力
DK,1,ALL !
定义约束
FK,2,FY,-1 !
顶部施加载荷
SOLVE !
求解
FINISH !
退出求解
/SOLU !
重新进入求解模型进行屈服分析
ANTYPE,BUCKLE !
屈服分析类型
BUCOPT,LANB,1 !
1阶模态，子空间法
SOLVE !
求解
FINISH !
退出求解
/SOLU !
重新进入求解展开模态
EXPASS,ON !
模态展开打开
MXPAND,1 !
定义需要展开的阶数
SOLVE !
求解
FINISH !
退出求解
/POST1 !
进入通用后处理
SET,LIST !
列出特征值求解结果
SET,LAST !
入感兴趣阶数模态结果
PLDISP !
显示变形后图形
！
NonLinearBuckling!
非线性分析部分
FINISH!
这两行命令清除当前数据
/CLEAR
/TITLE,NonlinearBucklingAnalysis
/PREP7 !
进入前处理
ET,1,BEAM3 !
选择单元
MP,EX,1,200000 !
弹性模量
MP,PRXY,1,0.3 !
泊松比
R,1,100,833.333,10 !
定义实常数
K,1,0,0,0 !
底端节点
K,2,0,100,0 !
顶点
L,1,2 !
连成线
ESIZE,1 !
网格尺寸参数设定
LMESH,ALL !
划分网格
FINISH !
退出前处理
/SOLU !
进入求解
ANTYPE,STATIC !
静态分析类型（非屈服分析）
NLGEOM,ON !
打开非线性大变形设置
OUTRES,ALL,ALL !
选择输出数据
NSUBST,20 !
5个子步加载
NEQIT,1000 !
20步迭代
AUTOTS,ON !
自动时间步长
LNSRCH,ON !
激活线搜索选项
/ESHAPE,1 !
显示二维状态下变形图
DK,1,ALL,0 !
约束底部节点
FK,2,FY,-50000 !
顶部载荷稍微比特征值分析结果大
FK,2,FX,-250 !
施加水平扰动载荷
SOLVE !
求解
FINISH !
退出求解
/POST26 !
进入时间-历程后处理器
RFORCE,2,1,F,Y !
2#变量表示力
NSOL,3,2,U,Y !
3#变量表示y方向位移
XVAR,2 !
将x轴显示2＃变量
PLVAR,3 !
y轴显示3＃变量数据
/AXLAB,Y,DEFLECTION !
修改y轴标签
/AXLAB,X,LOAD !
修改x轴标签
/REPLOT !
重新显示图形
示例2：
!
悬臂梁受端部轴向压力作用的屈曲分析
!
先进行静力分析，在进行特征值屈曲分析，最后进行非线性分析
!
静力分析
fini
/cle
/filname,beam-flexure
/tittle,beam-flexure
/prep7
*set,f1,-1e6 !
设置轴向压力荷载参数
et,1,beam189
mp,dens,1,7.85e3 !
设置材料参数
mp,ex,1,2.06e11
mp,nuxy,1,0.2
sectype,1,beam,I,,2 !
设置截面参数
secoffset,cent
secdata,0.15,0.15,0.25,0.015,0.015,0.015,0,0,0,0
k,1,0
k,2,2.5,0
k,3,1.25,1
lstr,1,2
latt,1,,1,,3,,1
lesize,1,,,10
lmesh,1
/view,1,1,1,1
/eshape,1.0
dk,1,,,,0,all,
fk,2,fx,f1 !
施加关键点压力
finish
/solu
antype,0
eqslv,spar !
求解器设置 稀疏矩阵直接法
pstres,on !
打开预应力开关
solve
finish
!
特征值屈曲分析
/solu
antype,1
bucopt,lanb,6,0 !
取前六阶模态分析
mxpand,6,0,0,1,0.001
solve
finish
/post1
*do,i,1,6
set,i
pldisp,i
*enddo
*get,freq1,mode,1,freq
finish
!
非线性屈曲分析
/config,nres,200 !
只记录两百步的结果
/prep7
tb,biso,1,1,2 !
定义材料非线性
tbtemp,0
tbdata,,2.0e8,0
upgeom,0.01,1,1,'beam-flexure','rst'
!
对有限元模型进行一阶模态的位移结果0.01倍的修改
save,beam-flexure,db
finish
resu,beam-flexure,db
/solu
antype,0
nlgeom,1 !
打开大变形
outres,all,all
arclen,1,0 !
弧长法设置
arctrm,l !
弧长法终止准则达到第一个峰值时终止计算
nsubst,200,,,1
fk,2,fx,f1*freq1
！
fk,2,fx,f1*freq1*1.2 !
将轴向压力值放大，放大系数为第一阶模态的主频
solve
finish
问题描述
一根直的细长悬臂梁，一端固定一端自由。
在自由端施加载荷。
本模型做特征值屈曲分析，并进行非线性载荷和变形研究。
研究目标为确定梁发生分支点失稳（标志为侧向的大位移）的临界载荷。
问题特性参数
本例使用如下材料特性：
杨氏模量=1.0X10e4psi
泊松比=0.0
本例使用如下的几何特性：
L=100in
H=5in
B=2in
本例的载荷为：
P=1lb
问题示意图
特征值屈曲分析是线性化的计算过程，通常用于弹性结构。
屈曲一般发生在小于特征值屈曲分析得到的临界载荷时。
这种分析比完全的非线性屈曲分析需要的求解时间要少。
用户还可以做非线性载荷和位移研究，这时用弧长法确定临界载荷。
对于更通用的分析，一般要进行崩溃分析。
在模型中有缺陷时一定要做非线性崩溃分析，因为此时模型不会表现出屈曲。
可以通过使用特征值分析求解的特征向量来添加缺陷。
特征向量是最接近于实际屈曲模态在预测值。
添加的缺陷应该比梁的标准厚度要小。
缺陷删除了载荷-位移曲线的突变部分。
通常情况下，缺陷最大不小于10%的梁厚度。
UPGEOM命令在前一步分析的基础上添加位移并更新变形的几何特征。
第一步：
设置分析名称和图形选项
1．选择菜单UtilityMenu>File>ChangeTitle。
2．输入“LateralTorsionalBucklingAnalysis”并单击OK。
3．确认PowerGraphics正在运行。
选择菜单UtilityMenu>lotCtrls>Style>Hidden-LineOptions。
确认PowerGraphics选项打开并单击OK。
4．将GraphicalSolutionTracking打开。
选择菜单MainMenu>Solution>-LoadStepOpts-OutputCtrls>GrphSoluTrack并确认对话框中radio按钮设置为ON。
单击OK。
5．生成屈曲分析图的输出文件。
选择菜单UtilityMenu>lotCtrls>RedirectPlots>ToGRPHFile。
将文件名改为buckle.grph并单击OK。
第二步：
定义几何模型
1．进入前处理器并生成梁的关键点。
选择菜单MainMenu>reprocessor>-Modeling-Create>Keypoints>InActiveCS，然后输入下列关键点号和坐标值：
关键点号：
1坐标值：
0，0，0
关键点号：
2坐标值：
100，0，0
关键点号：
3坐标值：
50，5，0
2．在关键点1和2之间生成一条直线。
选择菜单MainMenu>reprocessor>-Modeling-Create>-Lines-Lines>StraightLine。
将弹出生成直线对话框。
在图形窗口选择关键点1和2并单击OK。
3．存储模型。
选择菜单UtilityMenu>File>SaveAs。
在“SaveDatabaseto”对话框中输入buckle.db作为文件名并单击OK。
第三步：
定义单元类型和横截面信息
1．选择菜单MainMenu>reprocessor>ElementType>Add/Edit/Delete。
将弹出单元类型对话框。
2．单击Add。
将出现单元类型库对话框。
3．在左列选择“StructuralBeam”。
4．在右列选择“3Dfinitestrain,3node189”以选中BEAM189。
5．单击OK，然后的单元类型对话框中单击Close。
6．定义梁的矩形截面。
选择菜单MainMenu>Preprocessor>Sections>-Beam-CommonSects。
将出现梁工具对话框。
缺省时ANSYS将截面号设置为1，将子类型设置为RECT（在子类型处图示一个矩形）。
因为要生成一个矩形横截面，在子类型处不作修改。
7．在梁工具对话框的底部，可以看到横截面形状和尺寸的图示。
在B标志的部分输入0.2作为横截面的宽度；在H标志的部分输入5.0作为横截面的高度。
单击OK确定设置。
8．列出当前截面特性。
选择菜单MainMenu>Preprocessor>Sections>ListSections。
ANSYS缺省选择特性号1。
单击OK显示横截面信息。
在浏览过以后，在SLIST窗口单击Close。
2006-12-112:
55#1
冰彬鱼
工程师
精华4
积分61
帖子82
水位149
技术分16
状态离线第四步：
定义材料特性并定位结点
1．选择菜单MainMenu>reprocessor>MaterialProps>-Constant-Isotropic。
2．单击OK确认材料号为1。
将出现各向同性材料特性对话框。
3．在杨氏模量框输入1E4。
4．在泊松比（minor）处输入0.0，并单击OK。
5．选择菜单UtilityMenu>Select>Entities来选择线。
选择下列选项：
Lines，ByNum/Pick，FromFull并单击OK。
6．出现选择线对话框。
在图形窗口单击线实体。
在对话框中单击OK。
7．作为线的属性定义结点定位。
选择菜单MainMenu>reprocessor>-Attributes-Define>AllLines。
单击PickOrientationKeypointradio按钮旁边的radio按钮将其改变为Yes并单击OK。
ANSYS将材料特性号指向1，将单元类型号指向1并将截面特性号指向1。
8．出现线属性对话框。
在图形窗口选择关键点3并在对话框中单击OK。
9．存储模型。
选择菜单UtilityMenu>File>SaveAs。
选择OK，当ANSYS询问是否覆盖时，单击OK。
第五步：
对线划分网格并确认梁的定位
1．定义网格大小和分段数。
选择菜单MainMenu>reprocessor>-Meshing-SizeCntrls>-Lines-AllLines。
在No.OfElementDivisions框中输入10并单击OK。
2．对线划分网格。
选择菜单MainMenu>reprocessor>-Meshing-Mesh>Lines。
确认在“MeshLines”对话框中Pick和Single选定，然后在图形窗口选择线。
在对话框中单击OK对线划分网格。
3．旋转划分好网格的线。
选择菜单UtilityMenu>lotCtrls>Pan,
Zoom,Rotate。
弹出Pan,Zoom,Rotate对话框。
选择ISO并单击Close。
图形窗口中梁将旋转。
4．确认梁的定位。
选择菜单UtilityMenu>PlotCtrls>Style>
Size&Shape。
选择/ESHAPE旁边的radio按钮并单击OK。
5．显示横截面形状。
选择菜单MainMenu>Preprocessor>Sections>
PlotSection并单击OK。
6．重新显示网格。
选择菜单UtilityMenu>Plot>Elements。
第六步：
定义边界条件
1．定义固定端的边界条件。
选择菜单MainMenu>Solution>-Loads-Apply>-Structural-Displacement>OnKeypoints。
将弹出ApplyU,ROTonKPs对话框。
2．定义关键点1为固定端。
在ANSYS输入窗口，输入1并回车，然后单击OK。
3．在对话框中选择“AllDOF”，然后单击OK。
在ANSYS图形窗口将显示边界条件。
4．在自由端施加集中力。
选择菜单MainMenu>Solution>-Loads-Apply>-Structural-Force/Moment>OnKeypoints。
将出现ApplyF/MonKPs对话框。
5．定义关键点2为自由端。
在ANSYS输入窗口，输入2并回车，然后单击OK。
6．在Directionofforce/mom框中选择FY。
7．在数值处输入1并单击OK。
在ANSYS图形窗口将出现集中力标志。
8．存储模型。
选择菜单UtilityMenu>File>SaveAs。
选择OK，当ANSYS询问是否覆盖时，单击Yes。
9．选择菜单MainMenu>F
升级会员 