基于LS—DYNA的斜齿轮动力学接触仿真分析.pdfVIP

第 3期 (总第 154期) 机 械 工 程 与 自动 化 No．3 2009年 6月 MECHANICAL ENGINEERING ＆ AUT0MATION Jun． 文章编号 ：l672—64l3(2009)O3—0009一O3 基于LS—DYNA 的斜齿轮动力学接触仿真分析 赵 鑫 ，梁义维 ，李珊珊 (1．太原理工大学 机械工程学院．山西 太原 O3OO24；2．兰州职业技术学院．甘肃 兰州 73OO7O) 摘要 ：某重型卡车变速箱在重载时出现输 出齿轮副折齿现象，而在进行静载试验时却有足够 的安全裕度，为此 进行计算机仿真以找 出导致该齿轮副破坏 的原因。首先利用 Pro／E分别建立在理想状态和变形状态下的实体 几何模型，然后将模型导入到ANsYS中．最后使用 ANSYS／Ls—DYNA求解器求解。结果表明。由变形造 成 的冲击应力显著。 关键词 ：斜齿轮副；ANSYS／Ls—DYNA；动态接触分析 ；仿真 中图分类号 ：THl32．41 文献标识码 ：A 1 概述 单元。SOIID164体单元，不需要设置实常数 ，采用缺 某重型卡车变速箱在重载时出现输出齿轮副折齿 省算法。对PIANE163面单元 ，设置实常数时需要定 现象，而在进行静载试验时却有足够的安全裕度，因此， 义 单 元 的均 匀厚 度 为 O．O1，并 选择 s／RC0一 需考虑动载冲击，而该冲击试验却难以实现，为此进行 ROTATI()NAIHUGHES—IIU 面内多积分点改进 计算机仿真以找出导致该齿轮副破坏的原因。事实上， 型单元算法，以消除某种沙漏模态。将斜齿轮副内圈定 在静载下由于啮合副的刚度会造成两轴不平行、中心距 义为刚性体．设置约束时只留两齿轮绕其心轴的转动。 变大等变形．最后变形导致啮合状态劣化。为此，本文 然后．设定网格尺寸并划分网格。为了提高计算的 基于 ANSYS／IS—DYNA 对该齿轮副在理想状态和 精度，又不至于计算时间太长。就需要控制好网格形状 变形状态下的冲击应力水平进行仿真研究。 与网格密度。这里采用六面体单元．对该斜齿轮副不同 2 建模 部位的网格密度进行有效控制．根据模型不同部位的重 齿轮副参数见表 1。齿轮的3D实体模型在Pro／E 要程度，选择合适的网格密度进行划分。渐开线圆柱斜 中建立，主要包括渐开线、螺旋线与过渡圆角的生成， 齿轮三维有限元网格生成的关键是斜齿轮端面网格的 再用混合扫描生成齿轮轮廓，生成齿轮 3D实体几何 划分，先划分斜齿轮端面的轮廓线，由于两接触面比较 模型；接着在 Pro／E中进行有变形和无变形的齿轮副 重要 ，因此两接触面的齿廓线应相对比较细，再采用 装配 ，其中变形为主动轮心轴 由一端面为支点向外扭 SwEEP来划分整个齿轮。并将斜齿轮的内圈面设置为 转 了0．3。。将 Pr0／E中装配好的模型通过 Pro／E与 PI，ANE163面单元，如图1所示 。 ANsYS之间的数据接 口导人ANSYS中。由于该齿轮 表 l 齿轮副参数 副只有两对齿进行啮合 ，为了减少在 ANSYS／IS—