大学生方程式三球销半轴ANSYS静力学强度仿真校核

一、半轴极限工况计算

FSC赛车半轴所承受的最大扭矩可能出现在以下两种工况：

(1)工况一：在水平路面上以最大驱动力矩起步。

通过计算的得出赛车起步的地面提供的最大驱动力矩Mo和一档起步发动机提供的最大力矩M1：

M₀ = Fx × r = 1191.52 N·m

M₁ = Tmax × i₀ × i₁ × ig × η = 1145.03 N·m

由计算结果可知 M₁ < M₀，即在起步瞬间，发动机提供的扭矩小于地面所能提供的最大驱动力矩不会产生滑转。因此，取最大驱动力矩 M₀ = 1145.03 N·m 作为设计依据。对于德雷克斯勒差速器,当左右半轴输出转速相同时，其扭矩均匀分配到左右半轴。由于水平赛道起步时，左右半轴可视为转速相等，所以半轴所承受的扭矩 M= 1/ 2 M0，即 M为 572.515N.m。

(2)工况二：在过弯时，变速器处于二档，发动机以最大扭矩输出，差速器锁死，此时单侧半轴将承受较大扭矩。通过比较两种工况下半轴所需承受的最大扭矩，取较大值作为半轴设计的理论条件。

通过计算的得出赛车起步的地面提供的最大驱动力矩Mo和一档起步发动机提供的最大力矩M1：

M₀ = Fx × r = 1191.52 N·m

M₁ = Tmax × i₀ × i₁ × ig × η = 801.52N·m

由上可得 M2<M0'，即发动机提供的最大力矩全部用于驱动车轮。 对于德雷克斯勒差速器，赛车前行时选用的Ramp angles 为 40°，其对应的锁紧扭矩百分比ca=60%。可得单边半轴所承受的最大扭矩为 M'=0.6M0，即 M'=480.912N.m。

综上，以 M 的值 572.515N.m 作为半轴的设计理论条件。

二.对半轴的理论直径进行计算

(1)对于韧性材料 d 切应力计算公式为：

τ=0.577[δs]

[δs]=δs/n

其中，材料选用 Tc4，其屈服强度δs 取 1050MPa，安全系数 n 取为2，可得 Tc4 的许用切应力为[τ]≥292。425MPa.

(2)半轴的直径：

d≥

可得 d≥21.52mm

(3)因为小号三球销的小径为 22.6mm 大径为24.6mm，因此我们选用直径为 22mm的轴。

三.对半轴进行静态力学仿真

(1)选择仿真类型

选择静态结构(将其拖拽到项目原理图处)

(2)导入材料

1.打开工程数据添加材料

2.选择钛合金(TC4)并进行添加

(3)导入catia三维建模图形

1.导入三维图

2.进入spaceclaim几何结构

将半轴上的三球销简化成如图所示:

(4)打开模型

(5)赋予材料给半轴和三球销

(6)对不同物件接触定义参数

1.对接触选择种类

2.定义接触参数

工程常用摩擦系数区间 0.1~0.2，取 0.15 贴合真实装配状态。

高斯点检测是 ANSYS 接触求解精度更高的算法，适配三球销曲面接触

法向刚度控制接触面之间 “穿透量” 的软硬程度，因数越小，接触越柔和：不会大幅影响应力结果：仅弱化数值震荡，不改变稳态接触应力分布,是多曲面摩擦接触的标准调参手段。

(7)进行网格划分

1.选择尺寸调整网格划分

2.定义网格单位

选择全部几何体,单位尺寸改成3mm(可以对接触面高应力处进行网格细分,同时也可以进行网络无关性验证)

(8)对物体进行约束

1.对分析设置进行更改

开启弱弹簧的目的:三球销半轴如果是靠接触传递载荷（比如万向节球销和球笼的接触、花键配合接触），在求解初始步，接触还没完全建立贴合时，模型会存在轻微的刚体漂移自由度，很容易报「矩阵奇异、存在刚体位移」的错误，求解直接终止,把模型的微小刚体位移锁住,开弱弹簧可以大幅降低收敛报错概率.

开启大挠曲的目的:开启大挠曲后，求解器会每一步迭代都更新节点坐标和刚度矩阵，考虑结构变形后几何形状变化带来的刚度变化，结果更贴近真实工况,开启大变形才能准确模拟接触位置的应力分布，不会低估应力集中,减小偏差。

2.一端添加圆柱面支撑,并将轴向改成自由

3.对另外一端添加力矩

(9)仿真求解

1.添加仿真结果

2.检查仿真结果

总变形

等效应力

安全系数