结构优化设计及matlab实现,基于ANSYS和MATLAB的结构优化设计

0引言结构静力分析是任何机械产品结构设计分析的重要方面。静力分析的目的是确定结构在静力条件下受指定载荷、温度和约束时的应力和位移分布,从而及时发现应力、应变危险点,为产品的结构优化设计进行指导。结构优化设计是将最优化理论应用于工程结构设计,在众多可行设计方案中寻找最适宜的设计方案。其中结构尺寸优化是在结构已经确定的情况下,通过调节参数(如板厚、梁截面尺寸等),达到优化结构性能的目的。ANSYS内部集成了结构优化设计的模块,可以直接利用有限元来辅助进行优化设计,但是需要学习APDL程序语言,对于一些简单的结构,其优化过程过于复杂。而MATLAB优化工具箱中包含了大量的优化求解函数,可以根据实际优化问题选择最适当的优化算法,简化了求解的过程;同时MATLAB语言简单易学,语法符合工程设计的习惯,具有很大的优越性。本文针对一般三维实体提出了一种分部优化的方法,其基本作法是:首先利用ANSYS对整体结构进行静力分析,然后根据已求得的各个零部件的受力状况建立结构优化设计的数学模型,最后利用MATLAB优化工具箱进行求解。这种优化方法不仅不需要编写繁琐的有限元程序,而且可以对优化结果进行验证以及采取进一步的优化措施,具有推广的价值。1三维实体优化设计的数学模型对于一般的三维实体,通常取各个零部件的板厚为设计变量,整体结构总质量最小为优化目标,考虑到强度、刚度等约束条件,其优化设计的数学模型通常为:设计变量:目标函数:1min()niiifXAx约束条件:h(X)0(k1,2,3,l)k==;g(X)0(j1,2,3,m)j=其中iA和ix分别为第i个零部件的截面积和厚度,为材料的密度,n为零部件的总个数。2计算实例分析结构如图1所示,上下为一个整体,材料为45号钢,其杨氏模量E200GPa,泊松比为0.3。其中面1受固定端约束,面2上施加了大小为40KN的均匀拉载荷,作用方向沿X正方向。现在要确定出该结构的变形和应力分布情况,同时求在满足强度条件下结构的质量最小设计。图1结构三维示意图【106】第35卷第10期2013-10(上)第35卷第10期2013-10(上)【107】2.1基于ANSYS的结构静力分析2.1.1定义单元类型和材料属性利用ANSYS8.0进行结构静力分析。选择高阶3维20节点固体结构单元(SOLID186),根据设计方案,材料属性为:杨氏弹性模量EPa11=210,泊松比0.3,密度33=7.8510kg/m,屈服强度MPas355。2.1.2建立模型采用自顶而下的建模方式。首先通过前处理中创建block模块分别绘制出两个长方体,再利用布尔加运算将两个长方体合并成一个整体,此时建模的工作完毕。2.1.3划分网格、加载并求解采用智能网格划分,在下立柱底面施加固定位移,上横梁右端面施加拉力载荷,大小p=40e3/0.02/0.03Pa,然后进行求解,求解后的模型如图2所示。图2求解后的模型2.1.4读取分析结果如图3所示,结构的最大位移为4.329mm,它出现在上横梁的最右端;同时从图中的虚线可以看出,结构的变形集中在上横梁的右侧,由于该处直接承受拉伸载荷,所以最容易产生变形。再如图4所示,结构的最大应力为1220MPa,它发生在下立柱的两侧面的底部,由于立柱与地面之间为固定端约束,此时立柱相当于一根悬臂梁,而该悬臂梁的截面面积较小,所以容易导致应力集中。图3结构的位移场分布等值线图4结构的等效应力场分布等值线图由以上的分析结果可知,在该工况下结构的最大应力和位移均超出了材料能够承受的范围。承受拉伸载荷的上横梁其每一部分的变形都在1mm以上