基于ansys软件的柴油机曲轴有限元分析.pdfVIP

基于 ANSYS 软件的柴油机曲轴有限元分析 1 引言 随着柴油机的不断强化，曲轴的工作条件愈加苛刻，保证曲轴的工作可靠性 至关重要，其设计是否可靠，对柴油机的使用寿命有很大影响，因此在研制过程 中需给予高度重视。由于曲轴的形状及其载荷比较复杂，对其采用经典力学的方 法进行结构分析往往有局限性。有限元法是根据变分原理求解数学物理问题的一 种数值计算方法，是分析各种结构问题的强有力的工具，使用有限元法可方便地 进行分析并为设计提供理论依据。 6105 高速柴油机曲轴为全支承式，总长 929mm，连杆颈直径为 63mm，主 轴颈直径为 70mm，分别在第 1、第 3、第4、第 6 连杆轴颈上设有非轴线对称 平衡块。本文采用 ANSYS 有限元软件，在静应力计算部分采用整体曲轴模型进 行有限元分析和模态分析，并就单元选择及网格划分对应力的影响做了分析比较。 2 整体曲轴有限元模型的建立 传统的曲轴分析，国内外多采用单拐或 1 /2 单拐模型。但这种方法还不能反 映整体曲轴内部应力场的分布状态，有些学者也采用连续梁理论对曲轴进行分析 计算，但把象曲轴这样复杂的结构简化成连续梁，计算结果显然是不准确的。因 此，为了较为准确地计算曲轴强度和了解曲轴内部的应力状况，本文采用曲轴三 维整体模型，对曲轴进行静强度和刚度的有限元分析。 2.1 有限元网格的划分 由于曲轴结构复杂，利用有限元软件进行建模时很难保证与图纸上的曲轴结 构完全一致，因此建模时必须简化。为了减少应力集中，曲轴上不同截面的结合 处都有半径不一的倒角，如果在建模时考虑这些倒角和油孔，则会使有限元的网 格非常密集，这就大大地增加了模型的单元数量，花费大量的求解时间，而且生 成的网格形状也不理想，降低了求解精度，因此在整体曲轴建模时仅考虑主轴颈、 曲轴轴颈与曲拐连接处的过渡圆角。 采用 ANSYS 有限元软件，根据曲轴的结构特点，结合有限元分析软件中所 提供的单元类型，选择 10 节点的四面体单元 Solsd92。在建模时，首先采用较 疏的网格实体模型进行网格划分，然后将连杆轴颈和曲轴主轴颈的圆角处进行网 格细划。剖分后形成的曲轴有限元网格见图 1。图中的曲轴模型共有 114622 个 节点，71982 个单元。 2.2 载荷状况的确定 曲轴在工作时承受缸内的气体压力、往复和旋转质量惯性力的作用，根据已 给定的发动机参数，通过发动机动力学计算，求得此发动机连杆轴颈载荷的最大 值，及随后曲轴再转过 120 度, 240 度, 360 度,480 度,600 度时连杆轴颈载荷 的数值，见表 1。 由于曲轴主要是因弯曲而破坏的，因此对曲轴受到飞轮处的扭转力可暂不考 虑，为简便起见，可假设对发火的气缸，当活塞处于上止点位置时连杆轴颈载荷 达到最大值 Pa，对六缸发动机只需考虑各个气缸分别处于压缩终了活塞在上止 点位置时的受力状况即可，发动机各缸的发火顺序为 1—5—3—6—2—4，可将 曲轴的受力状况用图形表示出来，其中三缸、四缸发火时曲轴受力状况见图 2。 2.3 曲轴的边界条件 2.3.1 力边界条件 根据传统的方法及有限宽度轴颈油膜压力应力分布规律，并忽略油孔处压力 峰值突变的影响，假定力边界条件为:载荷沿曲柄销轴向均匀分布，沿曲柄销径 向 120 度角范围内按余弦规律分布，如图 3 所示。 对于余弦规律的分布曲线方程可通过面积积分方式求得，这里求得余弦曲线 方程为 式中:F—对应于不同的 l 作用于曲柄销上的力; F 均布—假设力为均布时，平均分布于每个节点上的力; l—力作用点对应曲轴中心对称面的弧长。 根据以上所得力的分布公式，可求得各个离散单元的受力，然后将其分配到各个 节点上，从而得到等效节点边界力。考虑曲轴惯性力的影响，输入以额定转速 2500n/min 时的角速度 261. 7rad/s，有限元程序会自动将惯性力加在每一个节 点上。 2.3.2 支撑边界条件 将主轴承对曲轴的支撑视为弹性支座，设弹簧刚度为 K，认为 K 值在曲轴纵 向对称面内沿主轴颈均布，对于图 1 的曲轴有限元模型，可视 K 均分在曲轴纵向 对称面内主轴颈中截面左右的两个对称点上，本文中我们取 K 值为 3X1.0E+11N/m，这一刚度接近于主轴承的实际刚度。 在进行有限元分析时，为模拟曲轴的全支撑情况，约束弹簧对主轴颈的径向 位移;为模拟止推轴承的作用，可将曲轴左端面靠近轴心的对称四个节点的 Z 向 位移为 0，以防止曲轴沿轴向产生刚体位