ansys子结构分析实例Word格式文档下载.docx
- 文档编号:151264
- 上传时间:2023-04-28
- 格式:DOCX
- 页数:16
- 大小:1.02MB
ansys子结构分析实例Word格式文档下载.docx
《ansys子结构分析实例Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《ansys子结构分析实例Word格式文档下载.docx(16页珍藏版)》请在冰点文库上搜索。
然后定义单元类型shell63;
定义实常数1为:
2(板厚度)。
材料性能:
弹性模量E=201000;
波松比μ=0.3;
密度ρ=7.8e-9;
单位为mm-s-N-MPa。
采用边长1划分单元;
一端设置位移约束all,另一端所有(21个)节点各承受Z向力5。
计算模型如下图:
静力分析的计算结果如下:
为了后面比较的方便,分别给出两个area上的结果:
超单元部分,按照上述步骤操作如下:
(1)生成超单元
选择后半段作为超单元,前半段作为非超单元(主单元)。
按照ANSYS使用超单元的要求,超单元与非超单元部分的界面节点必须一致(重合),且最好分别的节点编号也相同,否则需要分别对各节点对建立耦合方程,操作比较麻烦。
实际上,利用ANSYS中提供的mesh200单元,对超单元和非超单元的界面实体,按照同一顺序,先于所有其它实体划分单元,很容易满足界面节点编号相同的要求。
对于多级超单元的情况,则还要结合其它操作(如偏移节点号等)以满足这一要求。
对于本例,采用另一办法,即先建立整个模型,然后再划分超单元和非超单元。
即:
将上述模型分别存为se_1.db(超单元部分)和se_main.db(非超单元部分)两个文件,然后分别处理。
对于se_1.db模型,按照超单元方式进行处理。
由于模型及边界条件已建立,只需删除前半段上的划分,结果就是超单元所需的模型。
然后直接进入创建超单元矩阵的操作,首先说明一下创建超单元矩阵的一般步骤:
A进入求解模块:
命令:
/Solu
GUI:
Mainmenu->
Solution
●
B设置分析类型为“子结构或部件模态综合“
ANTYPE
GUI:
Solution->
AnalysisType->
NewAnalysis
选择Substructuring/CMS(子结构或部件模态综合)
C设置子结构选项
SEOPT
Mainmenu->
AnalysisOptions
设置内容有:
方程求解方法(波前法–默认,或稀疏矩阵法–推荐);
要求生成的矩阵,可以是:
刚度矩阵;
刚度和质量矩阵;
刚度、质量和阻尼矩阵(稀疏矩阵法不能生成阻尼矩阵);
要求输出到output窗口的矩阵(矩阵和载荷矢量或仅输出载荷矢量,默认不输出);
扩展超单元的方法(后续操作时,对本超单元扩展时使用的方法),默认是回代法(backsubstitutionmethod),需要存储主自由度的三角矩阵(波前法为sebname.tri,稀疏矩阵法为sename.LNxx);
完全求解方法(Fullresolvemethod不存储任何三角矩阵)。
随着问题规模加大,三角矩阵可能很大;
完全求解方法不需要三角矩阵,但花费时间较长,且不能重启动子结构分析(对于不同的载荷矢量)。
质量矩阵公式:
是否使用集中质量矩阵,推荐用一致质量矩阵。
D选择主自由度
M
MasterDOFs->
UserSelected->
Define
主自由度是超单元和其它单元之间的界面,必须确保将超单元与其它单元的界面上的所有节点自由度定义为主自由度。
假如模型中除了超单元外没有其它单元,也必须定义主自由度。
在将超单元用于动力分析时,主自由度的动力特性应能代表超单元的动力特性。
对于在使用超单元时将被定义为约束或施加载荷的节点,应该定义为主自由度(如果它们也位于当前的超单元中)。
对于大挠度问题或使用SETRAN命令的情况,主节点的6个自由度都应作为主自由度。
本例中选择超单元外端面上所有节点的所有自由度作为主自由度。
一般情况均推荐选择相应主节点的All自由度。
E在超单元上施加载荷
作为超单元的结构可以施加任意载荷,但要记住:
程序会生成一个载荷矢量,其中包含了施加的所有载荷。
对于每个载荷步,在超单元文件中写一个等效的载荷矢量,最多可以写31个载荷矢量。
在超单元结构中可以使用非零位移约束,并作为载荷矢量的一部分。
求解后对超单元进行结果扩展时,对于有非零位移的载荷步,其数据库必须匹配,否则需再次施加非零位移约束后再进行结果扩展。
可以在使用超单元时再定义载荷和约束,此时应注意将相应节点选为主自由度节点。
同样也可以在使用超单元时再定义惯性力(加速度或角速度),但在超单元部分要生成质量矩阵。
如果在使用超单元时需要对其进行旋转,推荐这样做,因为超单元的载荷矢量是随着超单元一起转动的。
注意:
如果生成质量矩阵,推荐在使用超单元时再施加约束,但要注意将约束节点选为主自由度。
如果要施加加速度载荷,建议再生成超单元时施加,而不要施加到简缩的质量矩阵上。
对于有大转动的问题,推荐在使用超单元时再施加约束。
F设置载荷步选项
在生成超单元的阶段,唯一有效的载荷步选项是设置阻尼。
以上内容设置完毕,存储模型并求解。
对于本例,首先选择新的求解类型:
设置分析选项:
定义主自由度:
然后将模型存盘并求解。
求解后将生成超单元文件se_1.sub。
ANSYS中的超单元-静力分析
(2)
(2)使用超单元
使用超单元的过程是:
A建立新模型
对于本例,将se_1.db文件存盘,点击File->
Clear&
StartNew...菜单项,修改Jobname为se_main,,Resume模型。
此时读入的模型仍是个完整的模型,在前处理的Mesh下,清除后半段area上的网格,只保留前半段area上的网格(以及已经施加的节点Z向载荷向载荷)。
图10非超单元部分的模型
为了使用超单元,除了其它单元类型外,还需要定义单元类型Matrix50(超单元)。
图11定义超单元类型
在对非超单元部分建模时,必须确保其界面节点与超单元模型的界面节点精确匹配且最好节点编号也完全一致,否则需要偏移节点编号或在超单元节点和非超单元节点之间建立耦合,或由老的超单元模型创建一个新的超单元模型(SETRAN命令),其界面节点与非超单元模型相匹配。
对于本例,由于非超单元模型与超单元模型都是由同一个整体模型修改而来,故其界面节点的位置和编号都是完全相同的。
B读入超单元矩阵
首先设置单元属性为Superelement:
Type
Preprocessor->
Modeling->
Create->
Elements->
ElemAttributes->
Type->
Superelements
图12定义单元属性为超单元
然后读入超单元矩阵se_1.sub:
SE
Superelements->
From.SUBFile
读入se_1.sub文件,自动与非超单元部分组装到一起。
可以通过Utility菜单的List->
Other->
SuperelemData列出超单元se_1的数据:
现在应该施加超单元载荷,然后将模型存盘并求解(一般的静力求解)。
求解后存盘,再行后处理:
(3)将结果扩展到超单元部分
在生成超单元阶段产生的.EMAT、、.ESAV、.SUB、、.TRI、.DB、以及.SEID文件,和使用该超单元的模型求解时生成的.DSUB文件都是可用的,则可以将结果扩展到超单元部分。
清除当前模型,将Jobname改为超单元模型文件名(se_1),读入超单元模型,转换为Solution模块,激活ExpansionPass:
EXPASS
MainMenu->
ExpansionPass
在EXPASS表单中,设置:
ExpansionPass->
ON
然后读取超单元矩阵.SUB、以及非超单元模型求解时,为使用的超单元模型生成的数据文件.DSUB:
命令命令:
SEEXP
LoadStepOpts->
SingleExpand->
ExpandSuperelements
然后在SEEXP表单中分别输入超单元文件的.sub文件名和使用该超单元的主结构在求解时生成的.dsub文件名:
然后进行求解->
solve,求解完毕,可以到后处理查看计算结果。
注意此时并不生成.rst文件,但可以直接查看计算结果。
对于超单元部分,相应的位移和VonMises应力云图分别见图18和图19。
如果存在完整的结构模型(例如se_all.db),则可以显示整个结构的结果,具体操作方法是:
在将Jobname改为se_1后,不是读取其本身,而是Resumese_all.db,然后只选择超单元部分的单元和节点,并按照前面相同的办法对结果进行扩展。
结果扩展完成后,显示整个结构的单元和节点,然后进行后处理。
进入Post1后,首先执行:
Data&
FileOpts->
ReadSingleResultFile->
选择se_1.rst;
再点击ReadResult->
First,从se_1.rst中读取结果。
然后执行:
选择se_2.rst;
First,从se_2.rst中读取结果。
读入这一结果不会删除从se_1.rst中读入的结果。
然后可以显示或列表整个模型的结果:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ansys 结构 分析 实例