第7讲——模拟技术.ppt

对称性 2. 轴对称结构 一维轴对称单元模拟柱壳 二维轴对称单元模拟三维实体 可用一维或二维单元来模拟三维实体。 对称性 3. 循环对称 uAn = uBn uAt = uBt 多点约束方程(MPC) 对称性 4. 重复对称 uAx = uBx 支撑的模拟 结构的支撑通常比较复杂，有限元法一般通过施加节点约束来模拟支撑。 梁的固定端支撑模拟 一端固支梁 在梁固支端中点施加完全约束 在梁固支端底部施加完全约束 完全固定支撑 支撑的模拟 固定支座支撑模拟 L L1 L1 L3 L2 d2 d1 kb kb 弹性支撑模拟 在工程实际问题中，通常会遇到弹性支撑，如图所示转子系统轴承支撑，并非刚性约束，可以简化为弹性支撑，当叠加单元刚度矩阵形成总体刚度矩阵后，在弹性支撑对应的节点位置加上相应的弹簧刚度Kb即可。 连接的模拟 不同单元类型之间存在连接问题，特别是节点自由度数不同的单元需要通过合适的方法或一组多点约束方程来模拟连接。 考虑汽轮机叶片和轮盘的连接，汽轮机叶片和轮盘通常完好的连接在一起，现建立有限模型来分析叶片—轮盘系统。 若轮盘和叶片均使用二维实体单元来模拟，如图中所示，执行共节点操作使得叶片和轮盘的交界面共用同一个节点就能保证完好的连接。 连接的模拟 若采用梁单元模拟叶片，采用二维实体模拟轮盘，这种模型能够有效减少节点数目，但需要采用适当的方法来模拟连接，其原因在于： 若采用由图中方法连接梁单元和二维实体单元，只能保证梁的节点和二维实体单元的节点平动位移的保持一致，而梁节点本身的转动自由度无法得到限制，会导致叶片相对于轮盘自由转动，显然这不符合实际工况，我们需要对梁的转动进行限制。 连接的模拟 解决二维实体网格上梁的转动的一种简单办法就是将梁单元延伸到实体里面至少两个节点，这就允许在梁和二维实体单元之间传递平动和转动变形。然而，这种方法人为地增加了连接区域的质量，影响动力分析的结果。 另一个方法就是建立多点约束方程： 应用位移元进行有限元分析时，未知的场函数是位移。而实际工程问题所需要的往往是应力的分布，特别是最大应力的位置和数值，为此，需要回到各个单元中，利用下列公式计算应力应变。 应变矩阵是有限元插值形函数对坐标进行求导后得到的矩阵，求导一次，插值多项式次数降低一次，所以位移元中应力和应变精度比位移精度低。应力应变解的近似性表现在： 单元内部一般不满足平衡方程 单元与单元交界面上应力应变一般不连续 在力的边界Sσ上一般不满足力的边界条件 应力应变近似解必然在精确解附近上下振荡 表明在单元内必存在一个最佳应力点，接近或等于精确解 应力计算 单元或节点应力平均 最简单的处理应力结果的方法是取相邻单元平均值或围绕节点各单元平均值/绕点平均法。 (1)取相邻单元平均值 该方法常用于3节点三角形单元中。由于应力近似解总是在精确解上下振荡，可以取相邻单元应力的平均值作为此2个单元合成的较大四边形单元形心处的应力。即： 取算术平均作为应力平均值 取面积加权平均作为应力平均值 平均应力 =(单元1应力＋单元2应力)/2 平均应力 =单元1应力×单元1 面积＋单元2应力×单元1 面积 单元1面积＋单元2面积 应力计算 (2)绕点平均 取围绕该节点周围的相关单元计算得到该节点应力的平均值。即： 应该注意这样得到的节点应力值是围绕该节点的有限区域内的应力平均值。这种处理方法对于等参元来说并不能从根本上改善节点应力精度差的问题。 m表示围绕节点i周围的全部单元 单元或节点应力平均 应力计算 总体应力磨平 由于位移元解得的应力场在全域是不连续的，我们可以用总体应力磨平的方法来改进计算结果，得到在全域连续的应力场。总体应力磨平的方法就是构造一个改进的应力解，使其在全域是连续的。具体做法：改进的应力解与有限元法求得得应力解之间应该满足加权最小二乘法原则，换言之，用有限元得到的不连续的有限个应力解值通过最小二乘拟合出一个连续的应力曲面。使得求解区域内总体应力光滑。其缺点是计算工作量十分庞大。 应力计算 单元应力磨平 为克服在整个求解域上进行应力总体磨平时存在的计算量大的困难，可以采用单元应力进行单元应力局部磨平。方法与应力总体磨平是一样的，只是此时拟合点为单元节点应力。因此，对于不同的单元，得到的拟合曲面不同，这时，可把相关单元求得的改进节点值再取平均作为最后的邻接节点上的应力值。 一般情况下，采用单元应力局部磨平处理可以得到较好的结果，而且计算量小，所以应用十分广泛。但是拟合误差明显比总体应力磨平要大。 应力计算 子域局部应力磨平及外推 鉴于单元局部应力磨平不能充分利用最佳应力点具有的高一阶精度，同时一般不能得到在单元交界面上的连续