有限元分析薄板挠度(附C程序).docx

1问题描述 某周边简支非均匀的矩形（或圆形）板在均布载荷作用下挠度过大。结合实际，提出集中改进设计方案，并进行对比分析。 2.问题分析 不均匀板有两种主要的情况，结构不均匀和材料不均匀，结构不均匀是指板的厚度不是常量，材料不均匀体现在板的弹性模量和泊松比是变化的。另外，有的板可以是以上两种情况的混合情形。 不均匀板与均匀板的有限元问题有哪些差别呢？下面从均匀板问题推导出非均匀板有限元问题的解决方法。 2.1应力应变 先以结构不均匀板为例来讨论。假设一矩形板长为2，宽为2，厚度沿，不均匀，由一函数描述，但仍然符合薄板假设。对于均匀板，显然是一个常数。设挠度为，则板内应变向量可以表示为 应力应变关系为 弯矩扭矩矩阵 这里就体现出不均匀板和均匀板的区别了。积分完毕后，可以得到 其中薄板的弯曲系数矩阵 是关于薄板总体坐标的函数，所以对各个分单元都是不同的。 各单元的弯曲系数矩阵可以采用单元中心处的代替。那么就可以得出一系列的弯曲系数矩阵。如果单元划分得足够细，是可以代替真实解的。 2.2单元分析 可以将板分为边长为0.25的矩形小单元，每一个单元都是一样的。对于任何一个单元的节点，都有3项独立的位移 位移模式 形状函数矩阵是一个的行向量 其中 单元刚度矩阵 很明显，积分式中包含了弹性系数矩阵，而不同单元的弹性系数矩阵是不同的，所以，即便单元划分相同，得到的单元刚度矩阵也不同。对于均匀板，相同形式的单元，刚度矩阵是相同的。均匀与非均匀的差别，完全体现在弹性系数矩阵上。但是非均匀板的一些结果可以间接地运用。 矩形单元四节点单元刚度矩阵是一个规律性很强的对称矩阵。矩阵中待求的独立元素只有21个。 其中  如果以单元中心点的参数代替单元的参数，那么非均匀的薄板单元刚度矩阵与均匀板差别不大，即各个单元的刚度矩阵都需要利用单元中心的参数来代替。这个结论对材料不均匀板也是实用的。下面运用这种近似来求解非均匀薄板问题。 3问题求解 求解的模型仍然是边长为2的正方形薄板，材料的弹性模量，泊松比。 图1 薄板单元 3.1 均匀板 将板分为均匀的正方形单元，等，可以得到精度不等的数值解。 厚度为0.035的薄板，在载荷作用下可以计算出板的挠度。可以看出挠度是完全对称的，并且最大挠度出现在板的几何中心上。 图2 均匀板均布载荷下变形（整体变形和中线变形） 3.2 结构非均匀板 通过一个均匀板的计算，可以看出对称均匀板受到均布载荷后，挠度也是完全对称的。但是对于非均匀板，受到均布载荷挠度如何分布，最大挠度在什么地方，是最关键的问题。只有找到最大挠度的位置，以及最大挠度是多少才能采取相应的解决办法。 下面解决一个结构非均匀薄板。设板在x方向上厚度从0.02线性变化到0.05（平均厚度为0.035，与上一均匀板厚度相同），厚度关于中面对称，其余参数和前面的均匀板相同。结构如图。 图3 不均匀板模型 同样，将板分为（64个）正方形单元，在均布载荷作用下，板的形变已经完全不一样了。最大挠度不再出现在几何中心处。很显然，挠度最大处向板较薄的一侧推移了。值得注意的是，这种不均匀薄板在薄弱环节挠度往往增大得很快，平均厚度相同的薄板受到相同载荷，不均匀板的挠度往往大得多。 图4 不均匀板在均布载荷下变形 但是新的问题出现了。实际情况确实是挠度最大处向板较薄的方向推移了，可是目前的网格还是比较粗，所有的位移量都是针对节点计算的，并不能精确反映出现的位置。所以，要想较精确地得到具体位置，网格需要加密。 但针对现有的精度，我们可以提出几种解决大挠度的方案。 1）在最大挠度处添加支撑（固支或简支） 如果在最大挠度处添加支撑，那么该处的挠度就已固定。简支时，转角自由；固支时，转角为0或者给定。假如限制最大挠度处挠度为0.1。下面分别就固支和简支情况给出数值解。 图5 薄弱环节采用简支 图6 薄弱环节采用固支 2）在最大挠度出添加固支肋条 在最大挠度出使挠度为0，可以得到下面的情形。可以通过不断改变添加的位置找到最合适的方法。显然，图中不是最好的方法。 图7 薄弱处添加肋条 3）施加外力 例如给薄板最大挠度出施加一集中载荷500000N ，方向与均布载荷相反。可以看到最大挠度（0.12）也有明显减小。 图8 薄弱处施加外力 4）改变材质（材质分布） 不均匀薄板本来就可以是材质的不均匀。当改变材料性质可行时，那么这也是一种好的方法。下面通过改变材料的性质改变薄板的形变。 首先应该清楚，在薄板较薄的位置，结构的刚度较小。如果结构不能改变，那么需要适当地调整材料的性质。采用不同的材料，或者对材料进行不同的处理。那么，下面我们让薄板较薄的半边弹性模量为275Gpa，看看会有什么发生。 图9 改变材质（分布） 通过对比发现，薄板的最大挠度随着较薄区域的刚度的增加，最大挠度明显减小；如果全域