球磨机回转体的强度有限元分析.docxVIP

球 磨 机 回 转 体 的 强 度 有 限 元 分 析

单元荷载列阵和整个结构的节点荷载。由于作用在单元上的外载荷也都按静力等效的原则移置到相应节点上，所以每个单元只受到节点力的作用，单元的节点力列阵表示为{F}勹

(2) 单元分析J46)[47)[48J

有限元法的单元分析就是以单元的节点位移列阵{li丫为基础，结合具体的单元类型分别确定出该单元的位移分虽列阵{f}、应变分摄列阵｛c}、应力分量列阵｛叶和节点力列阵｛F丫的过程。对本文选取的三维空间单元solid45和solid92,单元位移分瘟列阵、应变分揽列阵和应力分贵列阵分别表示为：

{u,v,wrT

,{{c心 心 约在心},＇｛，气，6），心 房 g，心｝

根据选定的位移模式，可以先推导出用节点位移表示单元内任意一点位移的关系式，其矩阵形式为：

{J}=[N]{li} 0 -l)

其中[N]为形函数矩阵，它的元素是位宜坐标的函数。

对千线性弹性空间力学问题，由儿何方程表示应变与位移的关系为：

｛叶＝［a]{J}=[a][N｛]

砰＝［B］伈?｝

(3-2)

[B]=[o)[N] 0-3)

式中：［B］－－－－－变－转应换矩阵；

同 －－－－微－分－算子矩阵：在空间间题一中[a)为

。6 0 0

。

勿改 a 。

勿

。。a

。[a]＝|8 6

。

。匈灰

。

oz|

(3-4)

6 6

oz 句

8 。a

az 改

将应变矩阵 (3-2)代入物理方程{a}=[D｛]

叶中，得：

回 ＝［D)[B｛)砰＝［s］｛5}＇ （3-5)

其中[S]=[D](B]为应变转换矩阵，［D]为弹性矩阵，在空间问题中：

应的等效节点荷载，即F=F户E +F60 +

分别为：

，F，，。这些等效 节点载荷的具体表达式

Ff＝俨fdV

Fr=J炉TdS

尺＝j印。dV

凡＝j炉DG。dV

(3-10)

n:为了分析整体结构的应力、

n:

平衡状态下，单元的总势能 等于其形变势能U与外力势能尸之和。设整体结构被离散化为ne个单元，则离散后结构的总势能nP 为：

rrp=frr＇广钞立V 0- 11)

?I e=I ??I

将单元分析中求得的各个单元节点位移列阵{6}＇代入上式，可得：

＝芷－固[K｛]外－卧·｝r

＝芷－固[K｛]

外－卧

·｝r {R}＇ （3-12)

由最小势能原理可知：

e=12 ．?I

叮p

6{6} =0

0 -13)

将式 (3-II）和式 (3-12)代入式(3-13）中，可得：

n nc

I[K]飞｝＝区[Pr

=I ??I

0 -14)

式(3-l4)充分说明了结构的整体刚度矩阵与整体节点载荷列阵可以分别由组成结构的所有单元刚度矩阵和单元节点载荷列阵叠加得到。通过求解方程式(3-14)可以得到节点位移，然后由单元的应变转换矩阵、应力转化矩阵等即可求出单元的应力和应变。

回转体有限元模型前处理

前处理是指构建分析所用的有限元模型。在ANSYS中，有限元模型的建立又可分为直接法和间接法。直接法是根据结构的几何外形直接创建节点和单元得到有限元模型。而对千节点及单元数目较多的具有复杂几何外形的结构系统，则采用先创建几何模型，再对该模型进行实体网格划分，以建立有限元模型的间接

法。本文所要分析的球磨机回转体虽然几何结构规则，但所需节点单元数虽较多，所以在建模过程中使用间接法，其前处理过程主要包括创建实体模型、定义单元属性、划分有限元网格和修正模型等几项内容。

由千分析需要，所建立的球磨机回转体实体模型包括进、出料端端盖和装有大齿圆的简体部。为了合理利用计算时间，达到最佳计算效率，在保证分析结果准确性的前提下，对球磨机回转体模型进行了如下简化：

(I)本章的主要分析目的为确定回转体的最大｝莎力值和危险点位置，因此在建模过程中去掉了螺纹孔、很小的倒角等不影响此结果的要素；

(2)由千衬板和中空轴内部的衬套等属于磨损件，无法确定实际尺寸，所以在建模过程中并未创建其实体模型，只蒂在施加载荷时将其重舟按等效密度加到相应位置上，实现重力补偿；

(3)端盖与中空轴为焊接连接，在建模过程中将其视为连续的整体，不考

虑焊缝对分析结果的影响；

(4) 将大齿圆简化为分度圆直径相同的圆环实体结构，定义其材料特性时使用等效密度；

(5)模型通过两端的中空轴与轴承部相配合，由千轴承部 具