均匀穿孔铝板在均布荷载作用下的简单计算探究.docVIP

均匀穿孔铝板在均布荷载作用下的简单计算探究 　　摘 要：通过分析找到均匀开孔和未开孔铝板之间对应的关系，从而可以根据未开孔铝板的结果推测出均匀开孔铝板的结果。实现方法：利用有限元软件建模分析一块未开孔的铝板，得到其最大的变形及应力。然后建立相同规格具有不同孔隙率的一系列均匀穿孔铝板，在相同条件下分别对其进行分析得到各自对应最大的变形及应力。在整个对比分析过程中，仅有孔隙率发生改变，因而可以归纳出孔隙率与铝板最大的变形及应力之间的关系。 　　关键词：穿孔铝板；孔隙率；均匀载荷；计算 　　引言 　　穿孔铝板在实际工程中运用越来越广泛，在建模分析时，需要按照实际穿孔情况利用有限元软件建模分析，由于开孔较多，需要花费较多时间进行建模分析，并且当板块尺寸过大时，会产生过多的单元数量，出现软件内存不足而无法计算的情况。为了提高计算效率，简化穿孔铝板的计算，文章进行试验分析得出孔隙率与铝板最大的变形及应力之间的关系。希望对相关工程有参考价值。 　　1 建模分析得到初始数据 　　1.1 单元选择 　　因为最后得到的关系是由孔隙率（开孔的总面积/铝板的总面积）确定的，所以在建模分析过程中，除了孔隙率一个变量外，其他条件的都必须保持一致。实际工程中开孔都比较小，为了保证孔边的精度，每个圆弧设置16个节点，足以保证圆孔的精度。同时为了保证开孔板和未开孔板边界约束条件完全一致，就要求孔边单元外沿节点间距和数目与未开孔单元一致，因此选用strand7软件中如图1所示单元组对比分析。 　　1.2 模型设定 　　为了便于统计单元总数，设定板块的长宽尺寸为1000mm×1000mm，板厚3mm。将板块划分为100个单元组。没开孔的板块中每个单元组包含4×4个4节点单元，开孔的板块中每个单元组包含8个8节点单元。这样能够保证开孔和没开孔的铝板边界上的节点数目及位置保持一致。采用软件默认铝板材料的属性即可满足要求。假定均布荷载为1kPa，并且沿板四周约束每个节点X，Y，Z方向的平动自由度。 　　1.3 软件结果统计 　　通过调整倒角率来改变孔径直径，进而达到改变孔隙率的要求。经过strand7模拟实验得到不同孔隙率开孔铝板最大的变形与对应未开孔铝板变形之比（dr），及不同孔隙率开孔铝板最大的应力与对应未开孔铝板应力之比（str）。 　　为了便于观察随着孔隙率变化，铝板最大变形及最大应力的变化趋势，建立对应的分布图。 　　根据图2中分布情况，分别采用线性和多项式方式，利用mathcad对其进行拟合。拟合结果如图3。 　　由拟合结果可知，线性拟合最大的残差为0.019。拟合度为0.996，非常接近1，表明对其采用该线性拟合得到的线性函数非常接近试验数据。该线性函数为Line（x）=0.517x+1.014。多项式拟合最大的残差为0.134，拟合度为0.97，比较接近1，表明该拟合的函数比较接近试验数据。该拟合函数为： 　　2 结果分析、推广与验证 　　（1）因为约束条件，荷载形式不一致，铝板的变形及应力的计算式就不一致，所以，该结果仅适合于四周约束X，Y，Z三个方向平动自由度且均布荷载的情况。 　　（2）因为在试验过程中，仅仅变化了孔隙率，其他的变量没有改变，同时最后的结果是基于未开孔时的非线性结果乘以一个相应的系数而得到的，所以该结果适用于不同的厚度，风压，长宽尺寸。 　　下面演示其计算方式。首先，计算出相同规格未开孔铝板的最大变形d0和最大应力？滓0，比如试验中的4.681mm和10.3013MPa。接着计算开孔铝板与未开孔铝板变形的比值Line（x），以孔隙率为0.5为例，Line（0.5）=1.2725。那么孔隙率为0.5的铝板变形为d0×Line（0.5），即1.2725×4.681mm=5.957mm。试验中的结果是5.9083mm，误差为0.0487mm。该穿孔铝板的最大应力为？滓0×dxs（0.5）=21.85MPa。试验中的结果为22.447 　　8MPa，误差为0.5978MPa。结果满足工程运用要求。 　　以尺寸为1400mm×1600mm，厚度为2.5mm铝板在1.5kPa均布风压作用下结果来验证分析。经过建模计算分析，该铝板最大变形为9.9mm，最大应力为14.5MPa。当相同规格铝板孔隙率为0.2时，在相同边界条件及荷载作用下，最大变形为11.5439mm，应力为35.6MPa。利用推导公式计算最大变形为9.9×Line（0.2）=11.062mm，比较接近建模计算结果11.54。满足工程要求。最大应力为14.5×dxs（0.2）=29.5MPa。相对建模结果35.6MPa，存在偏差。偏差率为17%。对于该计算结果，假设穿孔铝板材质为3003-H24，根据《铝合金结构设计规范》，该铝板强度设