ASME应力线性化分析和总结.docxVIP

实用文档 实用文档 文案大全 文案大全 第一部分 ASME应力线性化 路径的定义 MSC.Nastran for Windows/ASME/Stress Linearization的计算是基于美国ASME 锅炉和压力容器规范及中国JB4732-95《钢制压力容器-分析设计标准》中的规范研制开发的，适用于压力容器的分析设计 Stress Linearization模块借助于MSC.Nastran for Windows中建立直线的功能来定义应力的线性化路径， 定义方法较多, 详见 MSC.Nastran for Windows的使用手册 。图1所示的三种路径都是可取的 。 一般位置 (ii)部分在有限元网格外 iii) 通过节点或在单元的边上图1 线性化路径 N1 - N2 路径上的积分点数 n 可由用户自己定义（缺省值为50）， MSC.Nastran for Windows线性化模块只要求每条路径的积分点数不超过100个。 积分点位置用该点离路径起点的距离? 来表示，如图2所示。每个积分点的? 值由MSC.Na- stran for Windows线性化模块自动按公式 (1-1) 确定： (1-1) ? ? i ? 1 ? t i n ? 1 其中: ? i : 第 i 个积分点的值 n : 定义在路径上的积分点数 t : 路径长度 i : 积分点编号 i = 1,2,3,. n 图 2. 积分点位置 对轴对称问题,? MSC?.Nastran for Windows线性化模块还将进一步自动计算每个积分点 ?i 所对应的全局坐标 r z i, i , 如图 3 所示，以供线性化时用。 1. 2. 计算积分点应力 图 3. 一旦确定积分点位置, MSC.Nastran for Windows线性化模块便会自动按以下步骤获取每个积分点的应力张量。 下面以一个简单的模型为例, 详述积分点 应力张量的获取过程: 步骤一: 确定相关单元, 计算权重因子 如图 4 所示, 对路径上的每个积分点? i  , MSC.Nastran for Windows线性化模块将找到其 ? 所在的单元号; 根据该单元的映射关系, 积分点 i 被变换到相应的参数坐标中, 确定其参数 坐标；用该参数坐标及该单元的插值函数, 计算出该积分点对该单元各个节点 (例如Na, Nb, Nc, Nd) 的权重因子Wi (例如Wa, Wb, Wc, Wd)。 如果积分点恰好在二个单元的边或节点上, 则该积分点将属于最早搜索到的那个单元。 如果积分点不属于任何单元, 则MSC.Nastran for Windows线性化模块将显示提示信息, 并将积分点排除。 图 4 步骤二: 确定节点应力张量 若结果数据库中已直接含有节点应力张量 , 则直接读取即可 , 否则 MSC.Nastran for Windows线性化模块将按以下步骤自动外推出节点应力张量: 确定外围单元 如图 5 所示, 通过步骤一, 已找到路径所经过的单元 m1, m2 , m3 , m4 。 本步操作将找出有限元模型中与 m1? m4 相邻的所有单元, 即 m5? m12 。 外推节点应力 对每个找出的单元 (m1, m2, ....., m12), 根据各自的插值函数, 如图6所示, 用外推方法, 由单元积分点应力求得各节点上的应力, 并将该应力变换到统一的应力输出坐标中。 图5 图 6 均化节点应力 ? ? ik ?N K? ijk j ?N1 (2-1) k 其中: ? ik : 均化后第 k 号节点上的第 i 应力分量值 N k : 与节点 k 相连的单元数 ? ijk : 第 j 号单元在其 k 节点上的第 i 应力分量值 步骤三 内插求积分点应力 将均化后的节点应力与对应的权重因子相乘, 即可得到积分点上的各应力分量。 1. 3 应力线性化 定义线性化路径, 积分点数, 并求得所有积分点应力分量后, 如何进一步作线性化, 还根模型是否为轴对称 (axisymmetric)有关。 对于非轴对称模型, 用户可以选择一个直角坐标系来确定线性化应力分量 ,,? ? ? , , x y z , ? xy , ? yz , ?zx 具体方向。缺省值为全局坐标系 Coord 0。 对轴对称, 情况较非轴对称复杂： 随着半径的增大, 单位面积所对应的材料量增加, 致使中性面位置会沿半径向外移动一段距离Xf 。 必须将应力分量变换到截面坐标系 (Section coordinate) 中, 如图 7 所示, 使 x(或r)方向跟路径N1N2平行, y (或z)方向垂直于路径 N1N2。 图 7 中性面处几何的曲率半径对线性化有影响。 MSC.Nastran for Windows