曲柄滑块机构基于ANSYS的运动特性分析Word下载.docx
《曲柄滑块机构基于ANSYS的运动特性分析Word下载.docx》由会员分享，可在线阅读，更多相关《曲柄滑块机构基于ANSYS的运动特性分析Word下载.docx（12页珍藏版）》请在冰豆网上搜索。
ThispaperusingfiniteelementmethodinvestigatetheCrank-slidermechanismanalysisofmovementcharacteristics.Firstly,byusinganalyticalmethodtosolveproblems，butsolvingprocessisrelativelycomplexandlargeamountofcalculation.ThereisacertaintypeofANSYSwhichisatypeofstructuredynamicsanalysis,inaccordancewiththerequiredanalysiswhichcanbeused.ThetopicchoicesthreedimensionalhingeunitCOMBIN7ontheanalysissolution.Intheend,theanalyticalsolutionandANSYSfiniteelementsolutioniscompared,wecangettheoptimalsolutionfromthefiniteelement.
KeyWords：
Finiteelementmethod;
Movementcharacteristics；
ANSYS;
Transientdynamic
第一章引言
1.1本文的主要研究内容
ANSYS是大型的通用有限元软件，其功能强大，可靠性好，具有强大的结构分析能力和优化设计模块，因而被国内外大多数机械行业所采用。
本文将基于ANSYS建立曲柄滑块机构的有限元模型，对曲柄滑块机构进行动力学分析。
首先，对ANSYS进行了简要的介绍，为曲柄滑块机构的有限元分析做好准备工作；
然后，利用ANSYS建立了曲柄摇杆机构实体单元模型，然后利用ANSYS对结构的动态特性进行研究，并给出分析步骤。
最后将利用解析法计算出的结果与ANSYS运算出的结果进行比较，从而对曲柄滑块机构的位移、速度、加速度进行分析。
第二章曲柄滑块机构的求解与分析
2.1问题描述
分析如图所示曲柄滑块机构的位移、速度、加速度随时间变化的情况。
已知：
=120mm,
=300mm,e=50mm,曲柄AB转速ω=100rad/s.
2.2有限元解
2.2.1定义参量
分析前要将分析过程中所需的参量进行定义，需要定义PI的值以及铰链A、B、C、D的坐标，OMGA1（曲柄转速），T（曲柄转动一周所需时间，单位为s），具体定义值如下。
参数
参数值或表达式
PI
3.1415926
OMGA1
100/（2*PI）
T
1/OMGA1
0.12
0.3
FI0
ASIN（E/（R+L））
E
0.05
（AX,AY）
（0,0）
（BX,BY）
（R*COS（FI0）,-R*SIN（FI0））
（CX,CY）
（（R+L）*COS（FI0）,-E）
2.2.2创建单元类型
本次课题中的俯仰结构是三维弹性材料，既承受拉压力也有弯矩存在，因此可以选用BEAM4单元和COMBIN7单元。
COMBIN7是个三维铰接连接单元,是三维销钉（或旋转）铰链单元，可用于在公共点上连接模型的两个或多个部分，适用于运动学静力分析和运动学动力分析。
BEAM4是一种可用于承受拉、压、弯、扭的单轴受力单元。
这种单元在每个节点上有六个自由度：
x、y、z三个方向的线位移和绕x,y,z三个轴的角位移。
可用于计算应力硬化及大变形的问题。
通过一个相容切线刚度矩阵的选项用来考虑大变形（有限旋转）的分析。
在Mainmenu中选择单元类型命令，选择Combination中的Revolutejoint7（COMBIN7），选择StructuralBeam中的3Delastic4（BEAM4），这样就完成了BEAM4单元和COMBIN7单元的创建。
2.2.3定义材料特性
在系统的材料属性设置对话框中依次选取Structural（结构）、Linear（线性）、Elastic（弹性）、Isotropic（各向同性），弹性模量（EX）为2E11和泊松比（PRXY）为0.3；
再依次选取Structural（结构）、Density（密度），不考虑各杆的惯性力，密度近似为0，将密度（DENS）定义为1e-14
，完成材料特性的定义。
2.2.4定义实常数
由于本课题使用了两种单元类型，需将两个单元分别加以定义。
Type1COMBIN7单元的实常量定义，需要定义X-Y平动刚度、Z方向刚度、X-Y旋转刚度、扭转刚度4个实常数，分别为1E9、1E3、1E3、0。
Type2BEAM4单元的实常量定义，需定义横截面积、Z轴惯性矩、y轴惯性矩、沿Z轴单元厚度、沿Y轴单元厚度5个实常数，分别为4E-4、1.3333E-8、1.3333E-8、0.02、0.02。
2.2.5创建节点
创建1~5号节点，1~5号节点的坐标值如表3.1所示。
表3.1关键点的坐标值
节点号
活动坐标系的坐标值
1
（AX，AY，0）
2
（BX，BY，0）
3
4
（CX，CY，0）
5
（BX,BY,-1）
2.2.6指定单元属性
拾取菜单MainMenuElem→Preprocessor→modeling→create→elements→elemAttributes选项，确定“COMBIN7”单元的单元属性，表示选定使用的“COMBIN7”单元，其材料属性编号为1，实常数编号为1。
2.2.7创建铰链单元
拾取节点2,3,5，创建了一个铰链单元，其中心在节点2,3处（节点2、3重合于B点）。
2.2.8指定单元属性
确定“BEAM4”单元的单元属性，表示选定使用的“BEAM4”单元，其材料属性编号为1，实常数编号为2。
2.2.9创建梁单元
分别拾取节点1,2，点击apply，再在拾取窗口的文本框中输入3,4，创建2个梁单元。
2.2.10指定分析类型
打开Solution中的NewAnalysis命令。
指定分析类型为“Transient”，即指定分析类型为瞬态动力学分析。
2.2.11打开大变形选项
本课题属于非线性稳定问题，因此需打开大变形选项。
打开Solution中的AnalysisOptions选项。
将“NLGEOM”打开，打开大变形选项。
2.2.12确定第一个载荷步时间和时间步长
打开Solution中的Time-TimeStep选项，弹出对话框，将第一个载荷步时间定义为T即曲柄转动一周所需的时间；
将时间步长定义为T/25，选择“KBC”为“Ramped”，指定载荷为坡度载荷。
2.2.13确定数据库和结果文件中包含的内容
在对话框中，选择下拉表框“Item”为“AllItems”，选中“Everysubstep”。
表明在结果中记录所有项目的各个分步。
2.2.14设定非线性分析的收敛值
在对话框中，分别选择“Structural”和“ForceF”，将受力的收敛值的参考值定义为1，将收敛参考值定义为0.1，再分别选择“Structural”和“MomentM”，将力矩的收敛值定义为1，将力矩收敛参考值定义为0.1，完成非线性分析收敛值的设定。
2.2.15施加约束
在建好的模型中选择所有节点，在对话框中，约束“UZ”“ROTX”、“ROTY”方向的自由度；
拾取节点1，在“ROTZ”方向上施加值为2*PI的约束；
再次拾取节点1，在“UX”、“UY”方向上施加值为0的约束；
再次拾取节点4，在“UY”方向上施加值为0的约束，完成约束定义。
2.2.16求解
选择Solve下的CurrentLS命令，进行问题的求解，出现“Solutionisdone！
”提示时，求解结束，进行结果的查看。
2.2.17定义变量
在对话框中，将“TypeofVariable”定义为“NodalDOFresult”，即变量类型为节点的自由度。
拾取节点4，在对话框中，在右侧列表中选择“TranslationUX”，完成变量2的定义，表示滑块的位移
。
2.2.18对变量进行数学操作
把变量2对时间t微分，得到滑块的速度
；
速度
对时间t微分，得到滑块的角加速度
在对话框中，在“IR”文本框中输入3，在“IY”文本框中输入2，在“IX”文本框中输入1，点击apply，弹出对话框，在“IR”文本框中输入4，在“IY”文本框中输入3，在“IX”文本框中输入1，完成对变量的数学操作。
经过以上操作，得到两个新的变量3和4.其中，变量3是变量2对变量1的微分，而变量2是位移
，变量1是时间T（系统设定），所以，变量3就是角速度
同样可知，变量4就是角加速度
2.2.19用曲线图显示角位移、角速度和角加速度
打开TimeHistPostpro下GraphVariables显示曲线图，在“NVAR1”文本框中输入2，点击OK，结果如图2.2所示。
再重复执行两次以上命令，在对话框的“NVAR1”文本框中分别输入3和4，结果如图2.3、2.4所示。
图3.2位移曲线
图3.3速度曲线
图3.4加速度曲线
3.3.20列表显示角位移、角速度和角加速度
打开TimeHistPostpro下的ListVariables命令。
在对话框中的“NVAR1”和”“NVAR2”文本框中分别输入2、3，在得到的列表中可以看出变量2即位移
的最大值为0.244047s，此即滑块的行程，该值对应的时间为1.8096s，即空回行程经历的时间，速度为0.411117E-02rad/s。
*****ANSYSPOST26VARIABLELISTING*****
TIME4UX3DERI
UX3
0.15080-0.536064E-02-0.825746E-01
0.30159-0.211211E-01-0.134514
0.45239-0.459290E-01-0.186453
0.60319-0.773539E-01-0.219561
0.75398-0.112147-0.230054
0.90478-0.146736-0.218103
1.0556-0.177925-0.188346
1.2064-0.203540-0.148371
1.3572-0.222673-0.105665
1.5080-0.235408-0.651425E-01
1.6588-0.242319-0.286455E-01
1.8096-0.2440470.411117E-02
1.9604-0.2410790.343714E-01
2.1112-0.2336810.635222E-01
2.2619-0.2219220.925966E-01
2.4127-0.2057550.121884
2.5635-0.1851620.150528
2.7143-0.1603570.176190
2.8651-0.1320250.195063
3.0159-0.1015270.202529
3.1667-0.709432E-010.194426
3.3175-0.428894E-010.168422
3.4683-0.201483E-010.124856
3.6191-0.523362E-020.668064E-01
3.76990.233469E-120.875674E-02
图3.6
第四章总结和展望
4.1本文所做工作总结
基于ANSYS的曲柄滑块机构运动分析使我在机械原理以及有限元软件分析上受到了很好的训练。
对于所学知识，有了一次综合运用的机会。
学习ANSYS的时间较短，加之英文操作界面，使得在ANSYS的运用上遇到不少困难。
首先对课题进行研究，认识到这属于瞬态动力学分析。
在对问题认识清楚后，采用了梁单元BEAM4和三维铰链单元COMBIN7建立了曲柄摇杆机构模型，施加约束后对其进行了运动特性分析。
通过这次毕业设计，我熟悉了ANSYS操作步骤，从定义单元类型到建模到结果分析，有了明确的思路和整体的认识。
了解了一般问题进行ANSYS求解的步骤方法。
而且，我还认识到不同的工程问题，ANSYS在定义时是根据实际情况定义的，没有固定的方式、方法。
4.2工作的展望
本次毕设是一个工程问题的实际应用型课题，在我国，机械产业还远远没有达到国际领先水平，并且与之相关的检测、质保、维护产业也没有很好的发展。
近年来，国家对于机械制造业有了大力的发展和推动，同时对于其相关产业也有大幅度的推进。
ANSYS作为工程问题解决的重要手段之一，也被得到重视、推广。
ANSYS可以模拟结构力学、热力学、流体力学、电磁学、声学等领域的各种复杂物理现象。
因此，ANSYS软件在我国的机械设计与制造、航空航天、交通、化工和能源、电子等行业得到了巨大的应用，将来会为各领域的科学研究和产品开出做出巨大贡献
参考文献
[1]程社联等.有限元法在机械设计中的应用[A].杨凌职业技术学院报，2009
[2]龚曙光，黄云清.有限元分析与ANSYSAPDL编程及高级应用[M]，北京，机械工业出版社，2009
[3]王富耻等.ANSYS10.0有限元分析理论与工程应用[M].北京:
电子工业出版社,2006
[4]杨咸启等.现代有限元理论技术与工程应用[M]。
北京：
北京航空航天大学出版社，2007
[5]黄云清.有限元分析与ANSYSAPDL编程及高级应用[M]。
机械工业出版社，2009
[6]盛和太等.ANSYS有限元原理与工程应用实例大全[M].北京:
清华大学出版社,2006
[7]杜平安等.有限元法—原理、建模及应用，北京:
国防工业出版社,2004
[8]邢静忠.ANSYS应用实例与分析[M].北京:
科技出版社,2006
[9]王富耻等.ANSYS10.0有限元分析理论与工程应用[M].北京:
电子工业出版社,2006
升级会员 