600g-t土工离心机吊篮力学分析与强度设计.docxVIP

? ? 600g-t土工离心机吊篮力学分析与强度设计 ? ? 杨玉明，刘远东，王沪毅 （中国工程物理研究院 总体工程研究所，四川 绵阳 621900） 土工离心机是用来模拟重力场，进行土工建筑物缩模试验研究的重要设备。土工离心模型试验是将1/n的模型置于离心机中，在n倍g离心加速度的空间进行试验。由于惯性力与重力绝对等效，且高加速度不会改变工程材料的性质，从而使模型和原型的应力、应变相等，变形相似，塑性区发展和破坏过程相似，以获取原型变形破坏机理和破坏过程的物理模拟试验技术[1—2]。土工离心机由传动系统、转臂系统、吊篮几大部件构成。吊篮是试验件的承载结构，试验件放置于模型箱中，模型箱安装在吊篮底部平台上，吊篮上端通过销轴连接于转臂的末端。在工作状态下，吊篮自身所受的重力产生离心力，离心机转速越高，离心力越大。模型箱及试验件也将产生很大的离心力，该离心力作用于吊篮平台上，因此，吊篮是关键的承力系统。这要求吊篮是高比强度、高比刚度的结构体系，在其强度设计和力学分析中，通常需对吊篮进行优化设计，文献[3—6]针对多个型号的离心机吊篮结构整体开展了优化设计，得到了满意的结果。 文中的土工离心机的容量为600g-t，有效旋转半径为5 m，有效旋转半径处（模型试验中心）的最大加速度达到200g，最大转速达到19.8 r/s，是目前国内较大的土工离心机。在对离心机吊篮初始设计方案进行力学分析的基础上，对吊篮开展了优化设计。因吊篮由多个部件构成，对吊篮作整体优化分析约束条件太多，不易收敛到满意的结果。不同于文献[3—6]，文中对结构的各个部件开展优化分析，在对吊篮结构作整体计算。通过优化设计，得到了满足强度要求的吊篮结构参数。 1 计算方法 1.1 接触算法 吊篮部件之间的连接为接触配合连接，吊耳与上端的向心关节轴承和下端的销轴均为接触配合，底部平台与销轴之间也为接触配合，在结构的力学响应分析中需要采用接触算法。 接触问题属于高度复杂的非线性问题，表现在接触的区域随时间不断变化，约束条件非常复杂，接触界面的摩擦条件也十分复杂。分析复杂的接触非线性力学问题，通常采用Lagrange乘子法和罚方法[7—11]。 Lagrange乘子法和罚方法是把约束条件乘以Lagrange乘子λ或罚参数α，并引入离散结构系统的变分公式，从而处理复杂的约束条件。Lagrange乘子作为独立的变量，即在平衡方程中增加了场变量，从而导致计算量增加。罚函数法是引入罚参数以消去平衡方程中的某一场变量，其困难在于罚参数的选择，会出现收敛性问题。 在接触分析中，通过位移约束来阻止网格的交叠，而且还要控制接触面网格的相互滑动。接触模型把可能发生相互接触的体（面或线）分为接触体和目标体，在增量方法的某步迭代中，当接触体上的节点k进入目标体的目标单元内时，若发生交叠，则通过返回节点k到目标单元的边界上。当接触体上的节点k是在滑动接触中，则位移约束应仅施加法向约束。约束方程为： 对于Lagrange乘子法，增量平衡方程为： 将式（1）的约束方程引入罚方法的平衡方程，其增量迭代方程为： ANSYS软件具有较强的接触分析功能，支持点-点接触、点-线接触、点-面接触、点-体接触、面-面接触、面-体接触，通过调节法向接触刚度、初始接触关闭条件、切向罚刚度等实常数，可以提高收敛速度以及计算精度。文中采用Ansys软件增进的拉格朗日方法——Augmented Lagrangian Method来分析结构的接触非线性。该方法可寻找到精确的拉格朗日乘子，对罚函数进行修正迭代[12]。 1.2 优化分析方法 优化设计包括设计变量、目标函数和约束条件三个要素。目标函数表述为： 式中：x1，x2，…，xn为设计变量。 约束条件为： 或： G（x）是不等式约束，H（x）是等式约束。目标函数是优化所要求最优设计性能，一般以结构的某要素作为目标函数，如质量最小或体积最小等。文中以Ansys软件来作结构的优化分析。 2 吊篮有限元分析及优化设计 2.1 吊篮结构 离心机吊篮由2个吊耳、底部平台以及销轴构成，如图1所示。其中，底部平台是由4根梁、加强筋板、底板和顶板焊接在一起的整体，底部平台的顶板上放置试验用模型箱。销轴用于连接底部平台与吊耳，吊篮的上端通过向心关节轴承与转臂系统的拉力带叉耳相连接。吊耳与销轴采用高强度钢34CrNi3Mo，底部平台采用Q690。离心机运转时，离心机主轴带动转臂旋转，转臂带动吊篮旋转，吊篮及模型箱受到离心力的作用。 图1 吊篮结构Fig.1 Structure of nacelle 2.2 有限元分析结果 吊耳上端孔与向心关节轴承为接触配合，连接吊耳与底部平台的销轴与吊耳下端孔和梁的孔也为接触配合，在这些接触配合部位建立面-面接触单元。由于结构和载荷的对称性，取1/4吊篮来