外文翻译模板

本科毕业设计（论文）外文翻译题目插座面板注塑模具设计姓名陈小勇专业机械设计制造及其自动化学号 201233467指导教师李玉龙郑州科技学院机械工程学院二○一六年三月参数化建模滚珠丝杠主轴达塔拉•米莫塔•格罗•光费尔摘 要产品开发过程的数值优化可以成功地应用于产品设计的早期阶段。

在滚珠丝杠驱动器很常见的情况下，动态现象大多数根据滚珠丝杠本身的几何形状而定。

轴向和扭转刚度相同的丝杠，最大速度和加速度不仅取决于伺服电机，也取决于丝杆直径，凹槽斜率和球半径。

此外联轴器的设计参数影响使优化变得更加困难。

为了捕捉这些影响，有效的数据（通常是有限元或MBS ）模型是必要的。

在这项工作中，一个新的更准确和有效的计算滚珠丝杠主轴轴向和扭转刚度被提出。

我们分析得到描绘的丝杠几何参数对大多数刚度的依赖关系的参数方程。

此外，我们增加一个确定函数的分析模型，从而提高了准确性。

在许多例子帮助下，所提出的分析模型针对有限元模型和目录数据进行了验证。

1 绪论滚珠丝杠主轴的轴向和扭转刚度中对滚珠丝杠驱动器动态特性起着重要作用，因为它基本上决定了滚珠丝杠驱动器的第一个和第二个特征值。

当用有限元建模时，滚珠丝杠驱动器的螺纹通常被忽略并且一些平均直径被用来建立简化的滚珠丝杆模型。

因此，关键是得到最接近的平均直径。

在大多数关于前人建模与仿真下，滚珠丝杆传动建模集中在滚珠丝刚螺母和滚珠丝杠主轴部件。

Jarosch 比较了不同类型的滚珠丝杠，但考虑到主轴简化为圆柱体，直径等于主轴外径，从而忽视了削减主轴螺纹。

随着了解的实际轴向uz k 和单位长度的螺杆扭转刚度zk ϕ，平均直径可以被计算为E k d uzuz m π4,=（1）4,32G k d zz m πϕϕ= （2）杨氏模量和剪切模量分别为E 和G 。

平均直径总是比主轴外径小。

对于每个刚度我们得到两个不同的平均直径。

这取决于每个应用的平均直径的最好选择。

这两个直径也可以做到线性组合。

一般滚珠丝杠制造商提供轴向刚度数据，但没有扭转刚度。

基于这个原因我们使用有限元法（FEM ）来计算两者滚珠丝杠主轴轴向和扭转刚度。

使用完全参数化的有限元计算模型，我们也可以不用滚珠丝杠主轴刚度。

此外，由于参数范围是不离散的，我们可以结合使用滚珠丝杠驱动器参数优化模型。

这种方法的困难是如何有效地计算轴向和扭转刚度。

有些作品提供了计算扭曲梁的性能计算方法，但只适用于抗弯刚度或弯曲特征频率，它的作用对滚珠丝杠驱动器是不太重要的。

2 详细的参数化有限元模型图1是我们描述的三维滚珠丝杠生成方法。

在有限元软件ANSYS 宏的帮助下，这个过程是全自动的。

我们的滚珠丝杠几何模型是参数化的，因此可以生成任意几何形状。

这几何模型描述是以下六个参数：主轴直径d1，主轴直径d2，滚珠槽半径Rs ，主轴间距Ph ，主轴长度Ls 和螺纹nT 的数量。

由于生产商不提供滚珠槽半径的数据，但提供滚珠直径Dw ，我们使用了振荡确定滚珠槽半径的关系：54.0=wsD r .这样的计算模型刚度可以非常准确的，但也很昂贵。

为了减少自由度数提高精度，我们划分滚珠丝杠在核心筒（0.9d2）和带螺纹圆柱。

建成线性，弹性和杨氏的弹性模量2m N 210.109E =与泊松比ν=0.3的模型。

图.1滚珠丝杠主轴建模与ANSYS为了计算滚珠丝杠轴向和扭转刚度，我们需要滚珠丝杠的轴向力和扭矩一个滚珠丝杠应用在两个不同的静态载荷步。

滚珠丝杠的另一端必须在同一方向的限制，以防止刚体运动。

同时滚珠丝杠两端应该能够在径向方向自由地转动。

表面上制约两个单一节点（TARGE170元素），我们运用这些制约因素和力。

制约和节点力分布在滚珠丝杠的一端，只有轴向和切向方向（CONTA 174）通过接触节点，见图.2。

3 目录数据比较为了验证我们的模型，我们产生的40种不同的Bosch-Rexroth 滚珠丝杠主轴的模型。

轴向刚度仿真结果与由Bosch-Rexroth 提供的目录资料相比，提供了数据。

作为对比的基准，我们使用的轴向刚度的分析方程，DINISO3408-4符合，主轴长度为：()()4cos 20απw uz D d E k -=（3）()()4cos 20απw uz D d E k -=公称直径d0、滚珠和滚珠凹槽之间的接触角α。

分析数值和产品目录显示，没有依赖主轴螺距。

但是，主轴螺距确定了有限元模型。

此外，如果我们区分产品目录和分析数值之间有限元和分析数值之间的偏差，我们可以看到产品目录的数据显示的最小的偏差。

这可以解释产品目录值仅仅是圆形的分析值。

图.2应用约束和滚珠丝杠轴轴向力4 滚珠丝杠分析模型我们计划推出一滚珠丝杠传动的分析模型，捕捉螺纹的作用。

拥有一个精确的分析方法来计算滚珠丝杠刚度的分析方程，并且优先考虑有效率的参数模型，而不是有限元模型，因为那需要大量的迭代。

4.1 主要思想对于一个无扭曲统一长度的等截面A 和惯性力矩I 的主轴，轴向和扭转刚度的分析方程为：EA k uz = （4）GIk z =ϕ （5）横截面面积为A ，极惯性矩为I 。

滚珠丝杠有固定截面，这截面沿着滚珠丝杠螺纹角为ϕ，根据纵坐标Z 和主轴间距Ph0确定：02h P zπϕ=. （6）螺纹数量的影响丝杠的刚度。

在我们的分析方法中，我们假设丝杠轴径A 和扭转刚度I 正比和，并且和Ph0相乘：()0h uz P EAf k = （7）()0h z P G I f k =ϕ （8）从公式7和8我们可以看出，我们分析推导滚珠丝杠主轴的刚度的两个步骤： - 分析推导A 和I - 数值计算确定f （Ph0）4.2 分析推导截面特性滚珠丝杠螺纹横截面面积的范围为：()πϕϕ20,21≤≤=d r c （9）另一条曲线（钟形曲线图.3）的形状取决于参数 d2，Rs 和Ph0。

在我们的推导中，我们认为在大多数情况下，螺纹的轮廓用哥特式形状。

这种简化只导致滚珠丝杠刚度的小偏差。

为了得到一个钟形解析函数，在图.3中我们削减扭曲的圆柱体。

图.3确定滚珠丝杠横截面积曲线平面垂直于丝杠的纵向轴。

忽略了扭曲的圆柱体曲率我们就可以简化，并且生成的段可近似于椭圆形：()1c o s 222220222=⎪⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛-'+⎪⎭⎫⎝⎛'p s h s s r P y r r d x αϕπ （10）螺纹斜率Ap ：()s h p r d P 220+=πα （11）切割曲线的结果由切割直线Y0=0切割椭圆得到：()()22022c o s 12ϕπαϕ⎪⎪⎭⎫⎝⎛--⎪⎭⎫⎝⎛+=s p h s s k r P r r d r （12）表示为：s r d k +=221s r k =2()sP h r P k πα2c o s 03=（13）我们简化公式12：()223211ϕϕk k k r k --= （14）为了计算丝杠横截面面积，我们需要圆（9）和钟形曲线（14）的交点（r0| U0）。

由于交叉点是在圆内，很明显，210d r =（15）角c ϕ等于（9）和（14）：221130211⎪⎪⎪⎪⎭⎫⎝⎛--=k d k k ϕ （16）丝杠截面面积可分为Ac 和Ak 区域，如图.4。

Ac 包括一个涵盖角2(0ϕπ-)的圆形机构，这样我们得到：()02144ϕ-=d A c （17）由曲线（14）包涵在0和U0之间（由于相对于X 轴对称）可以计算AK 区域：()⎰⎰-==-022232121212ϕϕϕϕϕd k k kd r A K K （18）经过一番计算，我们得到：()⎰⎰⎰---+=000022321232220222112ϕϕϕϕϕϕϕϕd k k k d k k d k k A K （19）第一和第二积分可以很容易地使用标准的积分公式计算。

第三个积分较为复杂，所以我们用数学符号求解：()tc o n s a ax a ax x ax tan sin 121122+⎪⎪⎭⎫⎝⎛+-=-⎰（20）图.4包围的丝杠截面面积使用公式19和20，我们得到A K ：()()⎪⎪⎭⎫ ⎝⎛+---+=303202302130232202221s i n 13k k a k k k kk k k A k ϕϕϕϕϕ （21） 类似于横截面积，我们划分极惯性矩在两附加的部分：Ic 和Ik 区域。

第一部分是圆的转动惯量：()04132ϕπ-=d I c （22） 对其余区域Ak 转动惯量，可以计算为：()⎰⎰--==0004223214121ϕϕϕϕϕd k k kd r I k （23）公式24类似18它可以得出同一组的其他公式19和20。

经过一番计算，我们得到的极惯性矩：()()()()0324103sin 32323323122120231200132232312230232221242233220222164241ϕπϕϕϕϕϕϕ-+⎪⎪⎪⎭⎫⎝⎛--+-+-+⎪⎭⎫⎝⎛--+++=d k a k k k k k k k k k k k k k k k k k k k k I （24）我们比较公式，得到22和25横截面的所有的参数设置数值计算值。

提出了错误的百分比，如图. 5。

对于最大的错误在于以下0.16％，这是微不足道的，而我的最大错误是较高的，但仍低于0.52％。

4.3 主轴间距因素对于一个给定d1，d2，rs 和NT 的滚珠丝杆，螺距影响的刚度可以分为两类：截面的影响和扭转量的影响。

截面的解析表示为公式22和25。

扭曲影响的正式表达Ph0通过公式7和8得到。

函数f （Ph0）是未知的但我们期望的特定值f （0）=1和f （∞）= 1转换为无螺距的影响，而对一些螺距的中间点的影响将达到最大。

Fisher 分布符合这些条件的最小的参数。

鉴于这些考虑，我们提出以下螺距函数：()110100+⎪⎪⎭⎫ ⎝⎛+⎪⎪⎭⎫⎝⎛=nh mh h b d P d P a P f （25）图.5 解析和数值计算滚珠丝杠第一部分计算面积A 和极惯性I 的误差百分比这类似于Fisher 的分布。

这个函数的主要缺点是系数的M ，N ，A 和B 是在转弯的最高点P h0max |fmax 。

我们给予不同的10％和300％d1的螺距。

在参数化有限元模型中，我们可以计算出f （Ph0）的不同值，如图. 6。

有限元计算结果得出了最相近的，反过来又得到M ，N 的系数a 和b 的值：()1m a x 10-⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫⎝⎛⎪⎭⎫ ⎝⎛-=f m n bn m n bn a d P m m n b m nh （26）图6 间距影响函数和其数值拟和我们使用曲线拟合来确定m，n。

得到的最佳拟合为m = 297.89, n= 299.32, a = 0.10 and b =-0.09，见图.6。