基于ANSYS／LSDYNA的平行分度凸轮机构的应力分析

1、 题目 基于ls-dyna的平行分度凸轮机构的应力分析 专 业 机械设计制造及其自动化 年 级 学生姓名 学生学号 指导教师 目 录第一章 绪论 11.1平行分度凸轮机构的来源及其优缺点11.2 本论文的设计思路21.3 平行分度凸轮机构的应用及发展趋势21.4 本论文研究的内容2第二章 平行分度凸轮机构的载荷计算32.1 转盘的运动规律计算32.2 主、从动件的负载转矩计算5第三章 凸轮材料属性及各参数83.1 凸轮参数设置83.2 *vread命令读取数据8第四章 建模及划分有限元网格114.1 建立工作轮廓曲线114.2 凸轮安装及整个模型的建立124.3 定义材料属性124.4 划分有

2、限元网格13第五章 施加载荷155.1创建part155.2 定义接触155.3施加载荷175.4设置求解选项18第六章 结果分析 196.1 通用后处理器/post1中滚子的受力分析 196.2 通用后处理器/post1中凸轮的受力分析 216.3 时间历程处理器/post26中节点的受力分析23第七章 结论 257.1 小结 25第八章 apdl方式建模及处理268.1 宏程序shuju1.mac 268.2 宏程序yundong.mac 268.3 命令流文件 26致谢 34参考文献 35英文摘要36基于ls-dyna的平行分度凸轮机构的应力分析贺国军 重庆三峡学院应用技术学院系机械设计

3、制造及其自动化专业2005级 重庆万州 404000摘要 本文主要论述了平行分度凸轮机构的设计过程，并根据所选择的改进正弦加速度运动规律计算了凸轮机构的载荷。通过ansys/ls-dyna有限元分析软件，完成了ansys/ls-dyna与vb的数据接口，成功导入了平行分度凸轮机构轮廓曲线的坐标值，正确建立了平行分度凸轮机构的有限元模型，并通过采用ansys/ls-dyna软件中的参数模块，将从动件运动规律和载荷分别计算得出具体的数值并绘制成曲线图。本文采用通用有限元分析的方法，对平行分度凸轮机构进行了正确的有限元分析，并采用post1和post26两个后处理器对有限元分析结果进行了定点定位的分

4、析，得出了平行分度凸轮机构的最大应力状态和最大应力值。此外，本文通过两种方式建模和分析，重点介绍了gui方式的分析过程和方法，简要阐述了apdl方式建模和分析，并给出了相应的分析命令流文件。通过ansys/ls-dyna的分析，得出平行分度凸轮机构接触强度的薄弱环节在凸轮轮廓面，最大应力出现在凸轮轮廓面的棱边上；用ansys/ls-dyna对平行分度凸轮机构进行应力分析，为平行分度凸轮机构的强度设计提供了理论参考；采用apdl建模，实现了以更改凸轮基本参数而任意更改凸轮大小和以更换vb数据源而任意设计凸轮运动规律的模块化建模，大大提高了分析的效率。关键词：平行凸轮机构 apdl参数建模 应力分

5、析 数组 曲线载荷 函数加载2009届 机械设计制造及其自动化专业毕业设计（论文）第一章 绪论1.1平行分度凸轮机构的来源及其优缺点平行分度凸轮机构，如图1-1所示，最早诞生于美国，1961年美国第一次公布了外接平行分度凸轮机构的专利，随后由美国商业凸轮公司（commercial cam and machine co.）和菲固索公司（ferguson machine co.）试制。上世纪七十年代初，日本著名机构学家牧野洋教授研究并发表了关于外接平行分度凸轮机构的论文，并编入自动机械机构学一书。1979年前后，他又同前田隆发表了关于内接和直动式平行分度凸轮机构的研究论文。同时日本的东芝精机、大塚

6、、三共、椿本等凸轮专业公司纷纷投入外接、内接和直动平行分度凸轮机构的系列化生产。从1983年起，西北轻工业学院的彭国勋开始研究平行分度凸轮机构，目前该机构在我国的应用还比较少(特别在数控刀架上)，而日本、韩国和我国台湾的数控车床刀架，普遍采用液压马达驱动，用平行共轭分度凸轮机构进行分度和预定位。图1-1 平行分度凸轮机构平行分度凸轮机构是使从动件作预期规律运动的高副机构，其主要的优点是：从动件的运动规律可以任意拟定，只要设计了相应的凸轮轮廓曲线，就可以使从动件按拟定的规律运动，此外，与圆柱凸轮机构相比，它具有输入输出轴平行、刚性好、分度数多、分度精度高、传动转矩大、运动平稳、易于制造等优点，既

7、可以采用内啮合，也可以采用外啮合，与齿轮传动配合可产生多种形式的传动方案，满足多种应用场合的使用要求。其主要的缺点是：设计过程复杂，设计质量不易保证，在凸轮接触处难以保证良好的润滑，凸轮叶片和滚子容易磨损，寿命不长，不能有过大的传动力，高速凸轮机构的动力特性较为复杂，难以计算等等。鉴于此，为保证凸轮良好的传动性能，更加广泛的为各行各业服务，这就需要根据实际工况需求对凸轮工作轮廓曲线进行精确设计，以及准确的受力分析和寿命分析，来对凸轮机构的设计进行综合考虑。1.2 本论文的设计思路平行分度凸轮机构的工作轮廓曲线的设计，必须满足精确的转位停歇转位工程需要，且转盘在转位的过程中具有一定的运动规律，所

8、以凸轮工作轮廓曲线的精确设计是相当重要的。本文采用vb 软件计算的构成凸轮轮廓曲线的各动态点的坐标值，在ansys/ls-dyna的前处理器中分别采用gui鼠标操作和apdl参数化两种方式进行建模及受力分析。ls-dyna软件最早是由美国john hallquist发布的，直到ansys公司购买了ls-dyna的使用权，才形成了ansys/ls-dyna产品，ansys/ls-dyna继承了ansys中强大的建模功能，尤其近几年来，apdl的开发已相当成熟，在前处理的建模过程中，极为方便，大大提高了设计效率。1.3 平行分度凸轮机构的应用及发展趋势凸轮机构广泛应用于食品机械、填充机械、药品机械

9、、封罐封盖机械、灌装机械、输送机构、冲床自动送料机构、玻璃机械、陶瓷机械、烟草机械、化工机械、电子设备，数控机床加工中心，旋转定位及直线定位间歇机构中。近年来，计算机辅助设计为人们提供了极大的方便，为人们的设计工作节省了大量的时间，凸轮轮廓曲线的设计达到了相当精确的地步，凸轮的应用越来越广泛，尤其是平行分度凸轮的应用更为人们所追崇，这是因为平行分度凸轮机构，占据空间较小，具有多用性和灵活性，且能满足用户的任意运动规律，因此随着凸轮轮廓曲线的设计精度的提高和有限元分析软件对凸轮的受力分析、寿命分析的深入，平行分度凸轮机构将向高精度，高承载能力，高适用性能等方向发展。1.4 本论文研究的内容在进行

10、平行分度凸轮机构强度设计时，如果用赫兹理论校核滚子和凸轮轮廓面间的接触强度，存在下述的困难：其一凸轮轮廓曲线的曲率半径计算十分复杂，其二凸轮在哪个运动状态时接触应力最大难以确定，再就是当平行分度凸轮与超过两个滚子同时接触时，属超静定问题，计算滚子与平行分度凸轮在接触处的法向力非常困难。鉴于此，本文通过利用ansys/ls-dyna软件，建立了改进正弦加速运动规律的平行分度凸轮机构模型，并采用通用有限元的分析思路，对凸轮机构进行了划分网格、加载、求解、post1通用后处理器和post26时间历程后处理器对计算结果进行了分析，为其他凸轮机构模型的建立和有限元分析，提供了良好的思路和方法。通过本论文

11、的研究，将提高平行分度凸轮机构的设计效率和设计精度。通过对平行分度凸轮机构的有限元分析，确定平行分度凸轮机构在实际工况下的受力情况，达到定量计算和评价平行分度凸轮机构的寿命和接触强度，改进平行分度凸轮机构的目的。第二章 平行分度凸轮机构的载荷计算凸轮机构的工作条件分为三类：低速轻载、中速中载、高速重载或靠模凸轮，不同的工作条件需要选用不同的材料，而不同的材料所能承受的压力、冲击载荷是不一样的，所以在设计平行分度凸轮机构时，要先选好平行分度凸轮机构的运动规律。不同的运动规律根据实际情况而定，常用的运动规律有：等速、等加速、等减速、余弦加速度、正弦加速度、改进正弦加速度、改进梯形加速度等。本文以从

12、动件作改进正弦加速度规律运动为研究对象。2.1 转盘的运动规律计算改进正弦加速运动规律在工程上应用得较为广泛，为了设计的广泛应用性，在此选取改进正弦加速运动规律，凸轮转速为n=100转/分,合n=弧度/s，凸轮在工作的一个周期内保持匀速转动，转盘在前半周期内按照改进的正弦加速度规律运动，在后半周期内停歇，以便工作台的动作，如此循环往复以达到凸轮机构的实际工况要求。在一个周期内研究凸轮的运动规律，即分度期时间为 （2-1）停歇期时间为 （2-2）其中：凸轮角速度为，转盘分度期转位角，凸轮转角=·t, 为转盘角位移，为转盘角速度，为转盘角加速度，计算如下：转盘分度期角位移： t0，) （

13、2-3） t，) （2-4） t， （2-5）转盘角速度： t0，) （2-6） t，) （2-7） t， （2-8）转盘角加速度： t0，) （2-9） t，) （2-10） t， （2-11）ansys/ls-dyna前处理器模块中的函数编辑器提供了强大的计算功能，可实现各种复杂函数、矩阵的计算并输出数组、数据和线形图，利用函数编辑器对凸轮的角速度、角加速度进行计算并绘制曲线如图2-1、图2-2：图2-1 转盘角速度图2-2 转盘角加速度2.2 主、从动件的负载转矩计算 凸轮和转盘在转动过程中，分别受驱动转矩和阻力转矩，由于转盘的运动规律的复杂性和不规则性，使得驱动转矩随着时间的变化而不断

14、变化，这对凸轮的性能有较高的要求。根据凸轮和转盘所受驱动转矩和阻力矩公式计算：转盘阻力矩取 凸轮的驱动转矩为 （2-12）式（2-12）中，为转盘的转动惯量，为转盘的角加速度，为转盘角速度，为凸轮角速度，在以上已计算。采用函数编辑器来运算并绘制凸轮的驱动转矩，如图2-3：图2-3 凸轮驱动转矩根据以上计算结果，创建四个数组：时间数组(time）,主转速数组(zzhuansu)，主转矩数组(zzhuanju)，阻力转矩数组(czhuanju)，并将计算结果分别读入此四个数组，以便施加载荷。注意在定义数组时，时间数组要与其他数组的行数列数要保持一致。转盘的运动规律曲线是由三条分段曲线组成，为了提高

15、曲线插值精度，把每条曲线分成20个时间段，加上起始时间段和终止时间段，共62个时间段，在每个时间段通过函数编辑器计算输出数值，并分别存放在time.txt，zzhuansu.txt，zzhuanju.txt，czhuanju.txt,四个文本文档中，以便ansys读取。综上所述，在定义数组时，需要定义四个1列62行数组，按钮、，添加zzhuansu数组，设置内容如图2-4，apply后，用同样的方法定义其他三个数组。 图2-4 定义zzhuansu数组然后读取数据，在按钮、下，读取数据，弹出如图2-5示图，设置如图，用同样的方法读取其他三个数组数据。图2-5 读取主转速 数据注意，数组与数据源

16、的格式要保持一致，即均为1列62行，（1f12.9）是表示读取数据的格式，不能省略或填写错误，否则将不能读取数据或读取的数据将不正确。第三章 凸轮材料属性及各参数根据从动件的运动规律来确定满足工况需要的材料，在ansys/ls-dyna中，对不同的材料，是通过定义不同的弹性模量、泊松比、密度等参数来表征的。在选择材料时应注意：相同金属材料比不同金属材料的粘着倾向大，单相材料、塑性材料比多相材料、脆性材料的粘着倾向大。禁忌的材料匹配是：非淬硬钢-青铜、非淬硬钢、尼龙及积层热压树脂；淬硬钢-硬青铜；淬硬镍钢-淬硬镍钢。3.1 凸轮参数设置由于数据均从vb中获得，所以为了建模的方便，本文一律采用国际

17、单位制。考虑到凸轮的适用环境，凸轮选合金结构钢，弹性模量e取pa，泊松比v取0.28,密度dens取7850 ，从动件选碳素工具钢，弹性模量e取 pa，泊松比v取0.25,密度dens取7850。 具体的凸轮参数见表3-1：名称代号值单位 名称代号值单位 中心距c0.1m转盘厚度a10.015m凸轮转速n100转/分转盘节圆半径rp0.046m分度转位i4个转盘半径rz1rp+0.012m凸轮头数h2个转盘中心半径rz20.02m滚子总数zh*i个凸轮总长ht0.09m滚子半径r0.01m凸轮轴外半径hr10.035m滚子高度b2*rm凸轮轴半径hrrz2m凸轮分度转角cf180度单片凸轮厚度

18、hl0.015m最大压力角amax50.4786度凸轮叶片间距离hh0.035m安装相位角cp145.973度数据行数hs181行 表3-1 平行分度凸轮各参数本文采用gui方式和apdl方式两种方式来完成整个模型的建立和分析，用户可以根据任何需要而任意更改上表中的任何值，从而改变凸轮的尺寸大小和基本参数，这对今后类似的凸轮建模提供了极大的方便，节省了大量的重复建模时间，且便于同业人士之间的交流和相互提高。打开ansys/ls-dyna程序，在按钮下，分别定义表3-1中的参数。3.2 *vread命令读取数据从vb中得到的是一系列构成凸轮轮廓曲线的坐标值，在ansys中只有通过读取数据的方式获

19、得这些值，ansys中的参数模块对数组定义及数据的读取提供了非常方便快捷的方法，而且效率很高。此处用命令*dim定义数组，用*vread命令读取数据，大致方法是：根据数据源（此处vb计算得出的数据以文本文档的格式保存，如 shuju.txt）的数量及规格，通过*dim（或相应的gui方式，在按钮下进行）定义与之一一对应的数组，然后通过*vread命令读入，完整的apdl命令如下：*dim，xx,181,8 ！在ansys程序中，定义空白数组。*vread,xx(1,1),shuju,txt,ijk,8,181 ！从外部文件（vb得到的数据文件）中读取数据。(8f30.5)其中，(8f30.5)

20、必须得有，括弧中8代表数组列数，f代表浮点数，30代表数据源中每个数据的宽度，5代表小数位数，*vread不能直接在命令窗口中执行，只能通过创建宏或嵌套宏的方式执行。因为vb中的数据是以mm为单位，所以需要将其换成m,通过以下命令流实现：*dim,xx2,181,8*do,i,1,181*do,ii,1,8xx2(i,ii)=xx(i,ii)/1000*enddo*enddo 转换过程的结果示意图，图3-2，图3-3：图3-2 读取数据后的结果图3-3 数据换算后的结果至此，材料属性定义、参数定义、数据读取完成。第四章 建模及划分有限元网格平行分度凸轮机构由工作台、凸轮叶片、凸轮轴、转盘和滚子

21、组成，为了简化模型，降低计算时间，本文中不考虑工作台，把凸轮叶片和凸轮轴通过布尔运算成一个整体，把转盘和滚子通过布尔运算成一个整体。4.1 建立工作轮廓曲线读取完毕后，就可利用这组数据来构建凸轮的工作轮廓曲线，总的来说，ansys的建模功能是不强大的，但通过apdl的方式建模却可提高精度和建模效率。绘制凸轮的工作轮廓曲线主要用到的命令是*do（循环指令）和k（描关键点指令），bsplin（画样条曲线指令）三个指令的结合，具体如下：*do,j,0,hs-1,1 !此处表示循环181次k,1+j,xx2(1+j,1),-xx2(1+j,2) ! xx2为换算后的数组*enddoksel,1,hsb

22、splin,all ! 画样条曲线，构建与1号滚子接触的轮廓曲线，如图4-1， lplot 用同样的方法画与2号滚子接触的工作轮廓曲线图4-2，然后构建凸轮基圆，删去多余的线条即可得到凸轮工作轮廓曲线如图4-3： 图4-1 部分工作轮廓曲线 图4-2 部分工作轮廓曲线 图4-3 工作轮廓曲线画好轮廓曲线后，按钮，选中工作轮廓曲线，由线生成面，然后拉伸，结果分别如图4-4，图4-5。 图4-4由线生成面 图4-5 由面生成体4.2 凸轮安装及整个模型的建立凸轮的安装角为：由安装相位角cp=145.973，计算凸轮镜像角度，在第三象限中建立一个与y轴成该角度的工作平面，镜像凸轮即可得另一片凸轮，然

23、后绘制凸轮轴，转盘和滚子，最后通过布尔运算即得整个凸轮模型。由于刚建立好的平行分度凸轮机构不是在初始接触位置，所以还需要根据凸轮的参数，通过公式 d=180-cp-arcsin(rp·sin/)= 17.0213315 （2-13）计算凸轮需旋转的角度。（2-13）中，d为凸轮旋转的角度，cp为安装相位角，rp转盘节圆半径，为转盘的基准起始位置角，为凸轮的基准起始向径。算得凸轮旋转角度d后，将其旋转到初始接触位置，最终建立平行分度凸轮机构，如图4-6： 图4-6 ansys建立的平行分度凸轮机构4.3 定义材料属性在ansys/ls-dyna有限元分析中，选用solid164单元，该

24、单元是由8节点组成的单元，具有x、y、z三个方向的平移自由度，但solid164单元不具有转动自由度，为了模拟凸轮的旋转运动，在此引入具有6个自由度的shell163壳单元，对 shell163单元，单元为均匀厚度，定义节点 1处的壳厚为0.1，选择 s/rco-rotational hughes-liu面内多积分点改进型单元算法，以消除某种沙漏模态。在定义材料属性时，将材料3，材料4（即由shell163组成的单元）定义为刚性体，并约束刚性体x、y、z三个方向的平移自由度和rotx、roty绕两个轴的旋转自由度。这是因为被定义为刚性体的单元可以大大缩短有限元分析时间，在定义刚性体后，刚性体的

25、材料性质都集中在单元的质心上，在施加转动载荷时，直接施加在刚体的质心上，使刚体绕质心旋转，即将凸轮中心面和转盘中心面定义为刚体，并对其进行约束和施加载荷后，转盘和凸轮分别绕自己的中心轴线旋转。在按钮： 下，定义单元类型1，solid164；单元类型2，solid 164；单元类型3，shell163；单元类型4，shell 163。然后更改shell 163单元的关键字，选择 s/rco-rotational hughes-liu面内多积分点。在按钮下，定义实常数1、2，设置shell163单元节点1处的厚度为0.1。在按钮下，定义材料属性：材料号1为凸轮，材料号2为转盘，材料号3为凸轮内圈表

26、面，材料号4为转盘内圈表面，其中设置材料号3，4为刚体，并约束其x、y、z、rotx、roty方向的自由度。4.4 划分有限元网格在ansys/ls-dyna中，有限元模型的大小决定了计算时间，随着模型自由度的增加，计算时间加大，所以要尽量避免结构中的一个小单元而使得整体计算时间增长，为了满足足够的计算精度，同时又考虑到计算机时间，由于凸轮机构的不规则性导致了扫掠划分网格的困难，故采用部分自由划分网格的方式：接触区和整个凸轮采用扫掠划分网格，其他部分采用自由划分网格。划分网格具体步骤如下：在主菜单下，点击，在弹出的图4-7中，设置如图，点击ok后，选中凸轮和凸轮轴，点击，在弹出的图4-8中，选

27、择相应的材料号1，单元类型1，实常数1，单元所在坐标系，设置整体单元尺寸大小为0.006，扫略划分网格。图4-8 设置换分网格参数图4-7 选择凸轮实体 用同样的方法，选中转盘上在转动过程中可能与凸轮接触的四个滚子，选择相应的材料号2，单元类型2，实常数2，单元所在坐标系，扫略划分网格；设置整体单元尺寸大小为0.012，选中剩余的实体，选择相应的材料号2，单元类型2，实常数2，单元所在坐标系，划分网格；之后，划分凸轮刚体单元，选中凸轮轴内圈表面，选择相应的材料号3，单元类型3，实常数1，单元所在坐标系，划分网格；最后，划分转盘刚体单元，选中转盘内圈表面，选择相应的材料号4，单元类型4，实常数2

28、，单元所在坐标系，划分网格。最后获得完整的有限元模型如图4-9。图4-9 划分网格后的有限元模型第五章 施加载荷平行分度凸轮在实际工作中，凸轮作匀速转动，转盘作改进正弦加速度运动。考虑到凸轮对滚子的压力较大，冲击载荷也一直处在变化之中，所以，需要驱动电机的功率满足一定的要求，并避免电机的震动而引起的凸轮瞬时载荷过大，对凸轮造成损坏。5.1 创建part在ansys/ls-dyna中，不同的体是通过不同的part来界定的，即把具有相同的材料属性、单元类型和实常数的一组单元（或称为单元组）定义为具有唯一编号的part，根据凸轮转动的需要，共定义四个part,将前述的材料号1、2、3、4划分的四组单

29、元分别定义为part1、part2、part3、part4，其中，不同的part具有不同的属性。在按钮下，选择，创建四个part,如图5-1。 图5-1 创建四个part5.2 定义接触在ansys/ls-dyna中没有接触单元，而是通过在运动过程中可能接触的part，在他们之间建立接触，ansys/ls-dyna中提供了18种接触类型供用户选择，用户只需根据具体的工程需要来设置相关参数，满足实际要求。在定义接触时，选择接触面和目标面为凸轮part1和转盘part2，选用面-面（sts）接触，接触算法用罚函数法。并设置参数为：静摩擦系数为0.1，动摩擦系数为0.1，指数衰减系数为0，粘性系数为

30、0，临界粘性阻尼系数为0，接触时间在0-10000000s之间，接触刚度为0.1，其他均采用默认设置。在按钮、下，弹出图5-2的图，设置如图参数后，点击ok，在弹出的窗口中，选择part1,part2,即可定义接触。然后点击按钮，弹出如图5-3示图片，所有的接触控制参数均采用默认值。 图5-2 定义接触参数图5-3 设置接触控制参数5.3 施加载荷在实际工程中，凸轮在转动的过程中拨动转盘周期性的转动。凸轮以固定的角速度转动，受变化的驱动转矩，而转盘则以改进正弦加速运动，受变化的阻力转矩，在前面已经介绍如何约束和加载，并计算出了具体的数值。此处用上述定义的四个数组，分别对凸轮施加角速度和驱动转矩

31、，对转盘施加阻力转矩。施加载荷的方法如下：前面已经定义了四个数组，并读入了相应的数据值，所以此处直接调用数组即可，按钮，弹出如图5-4所示，在此框中首先对凸轮加载主转矩，具体设置如图，然后apply，用同样的方法加载凸轮驱动转矩和转盘阻力转矩。图5-4 给凸轮内圈表面（即part3,刚体）加载转速5.4 设置求解选项凸轮在0.6s内，完成一周的转动，在0.3s内完成一个分度，而凸轮与转盘的主要受力情况是在分度的0.3s内发生，理论上讲，计算时间应设置为略大于0.3s，因为，假如计算时间为0.3s,系统会自动计算到0.3s的后一两步，但为了减少整体计算时间，这里取0.3s，结果输出文件.rst的

32、输出步数为 100步，时间历程文件.his的输出步数为100步。在按钮，下，设置计算时间0.3s，时间步长0.9，在按钮下，选择输出格式为ansys和ls-dyna两种方式，设置输出节点应力，接触应力和节点能量，然后按钮写出k文件，为了获得节点的接触应力，把k文件里关键字*contact_node_to_surface中的spr、mpr控制参数均改为1，然后执行求解。第六章 结果分析ansys/ls-dyna中提供了多种结果分析方式，有通用后处理器/post1、时间历程处理器/post26、lsprepostd三种方式，本文采用前两种方式进行分析。为了获得较为直观的接触应力结果，首先观察动画，

33、然后找出最大应力发生的时间和区域，再通过时间历程处理器，进行精确的定点定位分析。6.1 通用后处理器/post1中滚子的受力分析为了提高凸轮分析效率，首先直接选择转盘为分析对象，显示其动画过程，按钮、在弹出的图6-1中，先如图设置，之后点击。图6-1 动画显示 设置框然后点击，此时会弹出如图6-2图，选中，单击ok。图6-2 设置显示结果动画单击ok后，会弹出如图6-3，按照该图进行设置，表明查看的是应力中的等效应力图，然后ok，等候片刻，即可进入动画显示。 图6-3 选择等效应力动画观察动画结果，在凸轮刚开始转动时，最大应力发生在凸轮和转盘中心曲面上，为14.4mpa，这是由于分别在凸轮和转

34、盘中心曲面上施加了转矩；在啮合的过程中，滚子上的应力变化规律为：转盘与凸轮在接触区的等效应力随着凸轮与滚子的啮入，应力会逐渐变化，凸轮先与第一个滚子啮合，再与第一、二个滚子同时接触，然后与第一、第二、第三个滚子同时接触，最后与第二、第三个滚子同时接触。第一个滚子上的最大应力发生在0.084s时，值为183 mpa，如图6-4；第二个滚子上的最大应力发生在0.141s时，值为205 mpa，如图6-5；第三个滚子上的最大应力发生在0.258s时刻，值为91.3mpa，如图6-6。图6-4 第一个滚子啮合时的最大应力云图 单位/pa图6-5 第二个滚子啮合时的应力最大云图 单位/pa图6-6 第三

35、个滚子啮合时应力最大时云图 单位/pa6.2 通用后处理器/post1中凸轮的受力分析用同样的方法，选中凸轮，单独对凸轮的应力变化进行动画观察。凸轮上的应力变化为：初始接触应力在凸轮的上叶片上，在0.06s时，为51.7 mpa，如图6-7，应力慢慢向凸轮叶片尖角处移动，在0.111s左右，凸轮与三个滚子同时啮合，如图6-8；在0.183s时，凸轮下叶片上出现最大应力，为417mpa，如图6-9，这是由于应力集中引起的瞬时峰值，对凸轮产生瞬时冲击载荷。如图6-7 0.06s初始啮合时的应力云图 单位/pa如图6-8 与三个滚子同时啮合时的应力云图 单位/pa如图6-9 凸轮凸轮下叶片最大应力时

36、的应力云图 单位/pa6.3 时间历程处理器/post26中节点的受力分析根据观察以上动画结果，查得，在第二个滚子上的最大应力节点为4098，其等效应力时间历程曲线如图6-10。图6-10 节点4098的等效应力曲线由图6-10知，4098节点在0.10s时，开始进入接触状态，在0.141s左右，被凸轮下叶片拨动，接触应力趋于最大值，为205mpa。凸轮下叶片上的最大应力节点为1723，其等效应力曲线如图6-11。图6-11 节点1723的等效应力曲线由图6-11可知，1723节点在0.183s左右，等效应力最大，值为417mpa，但时间很短，这说明在凸轮叶片上的1723节点，出现了应力集中，

37、但因为时间极为短暂，所以对凸轮影响不是很大，但是在实际工况中，应尽量避免这样的现象。综上所述，滚子上的最大应力发生在下面的一轮滚子上，在0.141s时刻，最大值为205mpa,而凸轮上的最大应力也是发生在下片叶片上，在0.183时刻，最大值为417mpa，因此，在凸轮设计时，应该刻意的提高下面滚子和叶片的轮廓曲线精度、安装精度，提高其质量性能。查材料参数表，合金结构钢的屈服强度为785mpa,远大于凸轮上的最大应力417 mpa，碳素工具钢的屈服强度为355mpa，大于滚子上的最大应力205 mpa，故此，该平行分度凸轮机构的设计满足实际要求。第七章 结论7.1 小结本文通过利用ansys/l

38、s-dyna有限元仿真分析软件，对凸轮的参数化建模和处理，分析了凸轮在实际工况中的受力情况，ansys/ls-dyna软件所获得的结果与理论计算比较相符，这充分说明了有限远分析软件的快速性和准确性。从ansys/ls-dyna分析的结果来看，滚子上的最大的接触应力发生在滚子的上下两脚处，而凸轮上的最大应力发生在凸轮叶片上，这为凸轮的强度设计和凸轮的改进提供了良好的依据。诚然，ansys/ls-dyna目前还处于不成熟阶段，有待各界人士的学习、利用、研究和改进。本文通过ansys/ls-dyna软件对平行分度凸轮的分析研究，得出四个创新点：（1）通过ansys/ls-dyna的*vread命令读

39、取由vb产生的平行分度凸轮轮廓曲线数据，实现了基于vb的平行分度凸轮轮廓的计算机辅助设计与基于ansys/ls-dyna的平行分度凸轮机构的应力分析的数据接口。（2）平行分度凸轮机构接触强度的薄弱环节在凸轮轮廓面，最大应力出现在凸轮轮廓面的棱边上。（3）用ansys/ls-dyna对平行分度凸轮机构进行应力分析，为平行分度凸轮机构的强度设计提供了理论参考。（4）采用apdl建模，实现了以更改凸轮基本参数而任意更改凸轮大小和以更换vb数据源而任意设计凸轮运动规律的模块化建模，大大提高了分析的效率。第八章 apdl方式建模及处理在ansys中，可以通过gui方式(即鼠标点击的方式)建模和apdl（

40、即命令流的方式）两种方式完成整个有限元的分析。本文在gui方式的基础上，对apdl也做了一定的研究，并给出了完整的命令流分析文件。注意：应将宏程序置于工作目录下。8.1宏程序shuju1.mac*dim,xx,181,8, !定义数组，181行，8列*vread,xx,shuju0,txt,jik,8,181 !读取数据，shuju0.txt，是由vb计算所得(8f30.5) !括弧及其内容必须得有，8数组列数，30数据宽度，5小数位数该程序主要用于定义空白数组，并读取vb中平行分度凸轮机构工作轮廓曲线的数据。8.2宏程序yundong.mac*dim,time,62,1 !定义数组，解释同上

41、*vread,time,time,txt,jik,1,62(1f8.7)*dim,zzhuansu,62,1*vread,zzhuansu,zzhuansu,txt,jik,1,62(1f12.9)*dim,zzhuanju,62,1*vread,zzhuanju,zzhuanju,txt,jik,1,62(1f12.4)*dim,czhuansu,62,1*vread,czhuansu,czhuansu,txt,jik,1,62(1f9.7)*dim,czhuanju,62,1*vread,czhuanju,czhuanju,txt,jik,1,62 (1f12.4)该程序主要用于定义空白数

42、组，并读取相应的数据。8.3命令流分析文件值得注意的是，由于在apdl分析过程中要读取数据源，为了分析的方便，本文将各数据源和各宏程序一律放在ansys的工作目录下。finish /clear/filname,tun0/title,haha/units,si!此下定义单元属性,考虑到凸轮适用于高速重载环境，所以凸轮选合金结构钢弹性模e取!200gpa，泊松比v取0.28,密度dens取7850kg/m3热!膨胀系数取12(10-6/摄氏度)；!从动件选碳素工具钢，弹性模量e取210gpa，泊松比v取0.25,密度dens取7850kg/m3/prep7et,1,164 !定义单元类型1，2，3

43、，4et,2,164 et,3,163 et,4,163 keyopt,3,1,7 !设置shell163单元，选择 s/rco-rotational hughes-liu面内多积分点 keyopt,4,1,7mp,ex,1,200e9 !定义1,2号材料mp,dens,1,7850mp,prxy,1,0.28mp,ex,2,210e9 !单元类型2 mp,dens,2,7850mp,prxy,2,0.25 mp,dens,3,7850 !定义3,4号材料mp,ex,3,200e9 mp,nuxy,3,0.28 edmp,rigi,3,7,4 !定义刚体，并约束x，y，z,rotx,roty方

44、向自由度mp,dens,4,7850mp,ex,4,210e9mp,nuxy,4,0.25 edmp,rigi,4,7,4 *set,_rc_set,1, r,1 rmodif,1,1, , ,0.1, , , , !定义实常数，节点 1处的壳厚为0.1rmodif,1,7, 0 *set,_rc_set,2, r,2 rmodif,2,1, , ,0.1, , , , rmodif,2,7, 0 !此下定义各参数,其中，从动件为改进正弦加速度规律,步长为1/prep7c=0.1 !中心距n=100 !凸轮转速 转/分 化为 n*/30弧度/秒i=4 !分度转位数h=2 !凸轮头数z=h*i

45、!滚子总数 r=0.010 !滚子半径 b=2*r !滚子厚度 a1=0.015 !转盘底厚a2=0.015 !转盘中厚a3=0.015 !转盘顶厚rp=0.046 !转盘节圆半径 rz1=rp+0.012 !转盘半径rz2=0.02 !转盘中心半径sd=0.002 !sd表示下片凸轮距底盘的距离 *afun,deg !改为度x=rp*cos(360/(2*i*h) !滚子参数y=rp*sin(360/(2*i*h) *afun,rad !改为弧度，系统默认ht=0.09 !凸轮总长hr1=0.035 !凸轮轴外半径hr=rz2 !凸轮轴内半径hl=0.015 ! hl为单片凸轮厚度hh=0.

46、035 !凸轮移动距离cf=180 !凸轮分度转角 amax=50.47863 ! 最大压力角cp=145.9573 !安装相位角 hs=181 !vb输出数据的行数 *afun,deg !此下，d为凸轮旋转的角度，用于旋转凸轮至初始接触位置 d=180-cp-asin(rp*sin(180/z)/(c*2+(rp)*2-2*c*rp*cos(180/z)*(1/2) *afun,rad /prep7shuju1 !读入数组数据 后附宏程序shuju1.mac内容*dim,xx2,181,8 !定义数组*do,i,1,181 !此循环命令用于将mm单位换算成m*do,ii,1,8xx2(i,i

47、i)=xx(i,ii)/1000*enddo*enddo k,10000,0,0,0 !描关键点，是辅助关键点k,10001,200/1000,0k,10002,-200/1000,0k,10003,0,200/1000k,10004,0,-200/1000k,10005,0,0,200/1000k,10006,0,0,-200/1000k,10009,-c,0,0l,10002,10001 !描线，是辅助线l,10004,10003l,10005,10006/pnum,line,1 !显示线*do,j,0,hs-1,1 !循环命令，描凸轮工作轮廓各坐标点k,1+j,xx2(1+j,1),-xx2(1+j,2)*enddoksel,1,hsbsplin