ansys模态分析及其详细过程.docx
《ansys模态分析及其详细过程.docx》由会员分享,可在线阅读,更多相关《ansys模态分析及其详细过程.docx(27页珍藏版)》请在冰豆网上搜索。
ansys模态分析及其详细过程
压电变换器的自振频率分析及详细过程
1.模态分析的定义及其应用
模态分析用于确定设计结构或机器部件的振动特性(固有频率和振型),即结构的固有频率和振型,它们是承受动态载荷结构设计中的重要参数。
同时,也可以作为其它动力学分析问题的起点,例如瞬态动力学分析、谐响应分析和谱分析,其中模态分析也是进行谱分析或模态叠加法谐响应分析或瞬态动力学分析所必需的前期分析过程。
ANSYS的模态分析可以对有预应力的结构进行模态分析和循坏对称结构模态分析。
前者有旋转的涡轮叶片等的模态分析,后者则允许在建立一部分循坏对称结构的模型来完成对整个结构的模态分析。
ANSYS提供的模态提取方法有:
子空间法(subspace)>分块法(blocklancets),缩减法(reduced/householder)动态提取法(powerdynamics)、非对称法(unsymmetric),阻尼法(damped),QR阻尼法(QRdamped)等,大多数分析都可使用子空间法、分块法、缩减法。
ANSYS的模态分析是线形分析,任何非线性特性,例如塑性、接触单元等,即使被定义了也将被忽略。
2.模态分析操作过程
一个典型的模态分析过程主要包括建模、模态求解、扩展模态以及观察结果四个步骤。
⑴•建模
模态分析的建模过程与其他分析类型的建模过程是类似的,主要包括定义单元类型、单元实常数、材料性质、建立几何模型以及划分有限元网格等基本步骤。
(2).施加载荷和求解
包括指定分析类型、指定分析选项、施加约束、设置载荷选项,并进行固有频率的求解等。
指定分析类型,MainMenu-Solution-AnalysisType-NewAnalysis,选择Modalo
指定分析选项>MainMenu-Solution-AnalysisType-AnalysisOptions.选择MODOPT(模态提取方法),设置模态提取数量MXPAND.
定义主自由度,仅缩减法使用。
施力fl约束,MainMenu-Solution-DefineLoads-Apply-Structural-Displacement。
求解,MainMenu-Solution-Solve-CurrentLS。
(3).扩展模态
如果要在POSTI中观察结果,必须先扩展模态,即将振型写入结果文件。
过程包括重新进入求解器、激话扩展处理及其选项、指定载荷步选项、扩展处理等。
激活扩展处理及其选项,MainMenu-Solution-LoadStepOpts-Expansionpass-SingleExpand-Expandmodeso
指定载荷步选项。
扩展处理,MainMenu-solution~So1ve-CurrentLS。
注意:
扩展模态可以如前述办法单独进行,也可以在施加载荷和求解阶段同时进行。
本例即采用了后面的方法
(4).查看结果
模态分析的结呆包括结构的频率、振型、相对应力和力等
所有图片人小可以自己调节
实例
1.在开始菜单中启动MechanicalAPDLProductlauncher,如图1所示其中WorkingDirectory
FileProfilesOptionsToolsLinksHelp
SimulationEnvironment
小ANSYSH
License:
Add-onModules1
ULS-DYNA(-DYN)UANSYSDesignXpiorer(-DVT)
wurapnysics1
File厂CustomizatK)rV/HighPerformance斤
ManagementPreferencesComputingSetup
WorkingDirectory[F:
\modei2
Browse.
JobName.|model2
Browse...
Run
ProductHelp
A12.1:
ANSYSMechanicalAPDLProductLauncher[Profile:
***LastANSYSRun***)HostnameTw...
代表你的工作路径,也就是所有ansys的生成文件都存在这个文件加下。
而jobName就代表你所做实例的名字,实例中我们将目录设为F盘,工作名设为model2.M后点击run运行进入ansys界面。
如图2所示。
图2
2.定义单元类型
拾取菜单MainM--Preprocessor-ElementType~Add/Edit/Delete,如图3所示。
然后单击,
图3
然后单击add按钮。
如图4所示,单击ok,选择solid5单元,这样就完成了对单元的选择。
然后再单击图3中的close按钮。
ALibraryof口ementTypes
LibraryofElementTypes
Elementtypereferencenumber
OK
3.
Apply|
图4
定义材料属性
r»=n,w
拾取菜单MainMenu-Preprocessor--MaterialProps~MaterialModels.弹出对话框,如图5所示。
定义相对介电常数。
在右侧列表中依次双击"electromagnetic","relativepermittivity^,"orthotropic".输入相关参数。
图5
在右侧再点击Density,如图6所示。
定义材料的密度。
密度为7500,单击ok。
/\DensityforMaterialNumber1
DensityforMaterialNumber1
TemperaturesDENS
7SOO
GraphI
AddTemperatureIDeleteTemperatureI
OK
GenoaI
图6
在右侧在依次点击“Piezoelectrics”“Piezoelectricmatrix"定义压电应力矩阵,如图7所示。
输入相关参数。
单击ok。
/\PiezoelectricMatrixforMaterialNumber1
PiezoelectricMatrixforMaterialNumber1
PiezoelectricMatrixOptionsIPiezoelectricstressmatrix[u]
Concel
图7
在右侧列表中依次双击"Structural","Linear","Elastic","Anisotropic"如图8所示,定义各向异性弹性矩阵设置,输入参数,单击ok。
AAnisotropicElastic:
tyforMaterialNumber1
AnisotropicElasticityforMaterialNumber1
AnisotropicElasticMatrixOptions
[stiffnessform
d
Temperature
Dll
012
013
D14
D15
016
022
D23
D24
D25
026
D33
D34
D35
D36
D44
D45
D36
044
D45
D46
D55
D56
T1
AddTemperatureDeleteTemperature]
AddRow
DeleteRow
Cancel
Graph|
Help
D66
点击Material-Exit完成对材料的定义。
4健模
在selection中输入L=10E-3然后点击accept,再输
选取菜单Parameters-ScalarParameters
入h=20e込点击accept,按照相同的步骤输入,w=10e-3,a3=1000.如图9所示。
liems
A3
H・2OOOOOOOOOE-D2
L=1.000000000E-02W=1.000000000E-02
Selection
Accept
Delete
Close
Help
图9
点击closeo
拾取菜单MainMenu一Preprocessor一Modeling一Create一Keypoints~InActiveCS弹出对话框,如图10所示,输入1,创建关键点1,坐标为0,0,0.点击applyo
再输入如图11所示
图11
点击applyo输入图12所示。
图12
图13
建立四个关键点。
复制关键点选取modeling-copy-keypoints弹出选框,如图14所示。
点击pickall
改变视角选取plotctrls-panzoomrotate-iso然后点击close
是原始坐标系不显示选取plotctrls-windowscontrols-windowsoptions如图16所示
选取noshown点击ok。
AWindowOptions
[/PLOPTS]WindowOptions
INFODisplayoflegend
LEG1Legendheader
LEG2Viewportionoflegend
LEGSContourlegend
FRAMEWindowframe
TITLETitle
MINMMin-Maxsymbols
FILEJobname
LOGOANSYSlogodisplay
WINSAutomaticwindowsizing・
・whenentirelegendturnedonoroffWPWPdrawnaspartofplot?
DATEDATE/TIMEdisplay
[/TRIAD]Locationoftriad
[/REPLOT]ReplotUponOK/Apply?
(MultiLegend▼
“On
FOn
&On
lOn
FOn
POn
[Graphicallogo
厂No
DateandTime
Notshown
IReplot▼|
Apply_JCancel|
图16
仓ij建直纟戈MainMenu一Preprocessor一Modeling~Create一line-1ine~straightline选取点1和5,点击ok。
划分直线网格,分为4份。
MainMenu-Preprocessor-meshing一sizecntrls~manualsize~line-alllines,如图17所示,点击oko
图17
通过关键点创建面,MainMenu-Preprocessor-Modeling~Create~
areas-arbitrary-throughkps弹出选框,选取点1,2,3,4<>
通过面拉彳eR体,MainMenu一Preprocessor-Modeling-operate~extrude~areas"alonglines弹出选框,选取刚才那个形成的面,单击ok,然后再选取前面我们建立的线,点击。
k,形成的体如图18所示。
P0IWT5
TYPEWUN
ANSYSl
JUM262011
00:
36:
53
图18
5.网格划分
显加线号plotctrls~numberinglinenumbers前打钩。
点击ok。
图19所不
图19
显示线plot-lines
选取meshing;-meshtool如图20所示。
选取line后边的set弹岀选框,选取线11.和13,单击。
心弹出如21所示,输入2,单击。
k。
划分为两部分。
MeshTool
£10.0心
ok花e
l^bbsl
-JSet1
11SmartSire
Jj
上
Fnc6
Coarse
SizeCofttro勺
Gbbol
Set|
Oeer|
Atobs
Set|
C3eer|
Lincs
屜1
Copy|
Oeer|
FlipI
Layer
Set|
C1"r|
KeypU
Sut|
Meeh:
1ValumoE
Shape
CTei
俺Hex
CHee
(9MafifXfd
r
|3or7sided
I06CJ
AE^tertSuesoiP^dli^
ND”No.ofdem:
divisions
ANGSEDwiicndK(degree)
(LKEqDtrnenisizesonpictedlinesSIZEEhnentedgelength
(NC(VisusedorlyJS12E6blarkorzero)
KYNCP/SEWcanbechenged
SPACESpacngratio
°KIAppW|
(i^eANGSEonlyifnuirberof/visions(N^andelementedgelength[SIZE)areblankorzero)
Cle-srattachedare^sandvclun^s
No
No
5el|
图20图21
选取hex,mapped,点击meshtool中的mesh弹出选框,选取体,单击ok。
如图22所示。
图22
6.施加约束
选取节点select-entities弹出如图23所示的对话框,点击ok。
/VSelectEntities
图23
扌合取菜单MainMenu一一Solution~DefineLoads一Apply一Structural"Displacement一symmetryb・c-onnodes弹出图24选框,单击ok。
图24
同理约束y轴,将图24中改为y轴,单击。
k。
7.分析类型
拾取菜单HainMenu-Solution-AnalysisType-NewAnalysis。
弹出对话框,选择"Typ已ofAnalysis"为"Modal",单击"OK"按钮。
8.指定分析选项
拾取菜单HainMenu-Solution-AnalysisType-AnalysisOptionso弹出对话框,在"No.ofmodestoextractzz文本框中输入10.扩展模态数为10如图25所示,单击ok。
AModalAnalysis
[MODOPT]Modeextractionmethod
BlockLanczos
PCGLanczos
Reduced
UnsymmetricDampedQRDampedSupernode
No.ofmodestoextract
(mustbespecifiedforallmethodsexcepttheReducedmethod)[MXPAND]
Expandmodeshapes
NMODENo.ofmodestoexpand
▽Yei
10
ElcalcCalculateelemresults?
[LUMPM]Uselumpedmassapprox?
[PSTRES]Inclprestresseffects?
OK
图25频率范围50000到150000,图26所示,单击oko
图26
约束z二0和z二h的面的节点电压为0.
选取节点select-entities弹出如图27所示的对话框,点击ok。
/\电瘁&
|Nodes▼]
|Byl-ocwtion
CXcoorsCYcoordinates皿NcoordinatesMln,Max
1°1
••FromFull
「AlsoSeloct
U»-»select
Selezxii
Invert
3e?
le?
None?
1c?
Det1c>
OK
zxppty
&
PlotZplM
Canuof|XEp|
•t
•
[Nodes
1ByLocation▼]
厂XcoordhmtoK
Ycoordinato8
<•Zcoordinates
Min,Max
h|
(»FromFull
「Rcselect
「、/XlsoSelect
「Unselect
ScleAll
Invert
SqIqNono
SoloBalo
X
ApplyI
■
Plot
Replot|
Cancel
Help
AApplyVOLTonnodes
Constantvalue
拾取菜单MainMenu—Solution-DefineLoads~Apply-e1ectric~boundary~voltage-onnode如图28所示,单击ok。
[D]ApplyVOLTonnodesasa
IfConstantvaluethen:
VALUELoadVOLTvalue
OKI
Apply
CancelIHelp
图28
同理定义z=h,电压为0.如图29.所示
施加约束时,先选取该面上的节点,然后直接施加就是该面的约束。
在每一步施加约束时,我们都先选择该面上的节点后再施加。
图23,27,29都是选择相应面上的节点。
然后是假载荷,不然施加的是全部节点。
9.求解
拾取菜单MainMenu-Solution~Solve~CurrentLS。
单击"SolveCurrentLoadStep"对话框的"OK”按钮。
出现“Solutionisdone!
"提示时,求解结束,即可査看结果了。
10列表固有频率
拾取菜单MainMe-GeneralPostproc-ResultsSummaryo弹出窗1丨,列表中显刀£了模型的前10率,可以看出结果虽然存在一定的误差,但与解析解是基本符合的。
查看完毕后,关闭该窗口。
11.拾取菜单MainMenu-GeneralPostproc-ReadResults-FirstSet
12.用动画观察模型的一阶模态
拾取菜单Ut订ityMenu-PlotCtrls-Animate-ModeShape。
弹出对话框,单击"0K”按钮。
观察完毕,单击“AnimationControllerv对话框的"close"按钮
13.观察其余各阶模态
拾取菜单MainMenu-GeneralPostproc-ReadResults-NextSet.依次将其余各阶模态的结果读入,然后重复步骤1
观察完模型的各阶模态后,请读者自行分析频率结果产生误差的原因,并改进以上分析过程。
KXXXX
INDEXOFDATASETSONRESULTSFILExxxxx
SET
TIME/FREQ
64944・
LOADSTEP
SUBSTEPCUMULATIUE
65772.
68952・
73930.
579656・1
685755.1
792109.1
893052・1
90.11862E+06
10Q.119Q7E+GG
1
2
3
U
5
6
7
8
9
10
1
2
3
4
5
6
7
8
9
1Q
图30
图30为各阶自振频率
DISPLACEMDH
STEF=15UB=1FREQ64944DMX-14.45B
•Foncuord/0Dclw/Drd
CForwardOnly
St6M
1DEI
Nqxi
Previous|
Close
Help|
图31
图31为1阶振型
总结,再施加约束时,仅仅是几个面,如果忽略图23,27,29这几个过程,那结果可想而知。
如果出错,仔细检查,你添加的边界约束是否正确,重复修改,结果就是这个。
我只是
把别人的例子详细的做了一遍。
比较适合初学者。