盒体板材厚度优化.docx

盒体板材厚度优化问题概述 优化的模型为一个盒体。分别有上下左右四个表面组成，每个表面的厚度的初始值为2mm。盒子的外尺寸为100*100*1000。 盒子的材料为16Mn，杨氏弹性模量为2.06E5MPa。密度为7.85g/cm^2。材料的强度为345MPa。 盒体的受载状况：在盒体长度方向的中央，受到一个200kgf的载荷；而对盒子体的两端进行支撑。模型建立基础模型建立：（图1） 第一次建模对于载荷加载不熟悉，建立的模型比较粗糙。首先，进行几何建模。在空间中分别设置（0,0,0），（1000,0,0），（1000,0,100），（0,0,100），（0,100,0），（1000,100,0），（1000,100,100），（0,100,100），（500,100,0），（500,100,100）。然后，生成边界线。之后，再生成上下左右前后表面。重点是要将中间的（500,100,0），（500,100,100）连线指定为硬线。这样几何模型就建立完毕了。之后，进行材料和单元属性的设置。材料的设置如图2，单元属性的设置如图3。图2图3然后，进行网格划分。先设定网格的密度，在模型的下侧两条边界上设定mesh seed，每边上的单元个数设定为20个。先对模型的上表面进行网格划分，需要注意的是由于有硬线的存在，Mesher应选为Paver。同样，之后再对模型的前后左右以及下表面进行网格划分之后，实现加载。在硬线上平均分布有四个节点，每个加载500N的力。两端中一端六个自由度全部约束，另外一端的x方向的自由度放松。有关建模的附加方法： 在两边的约束点，使用弹簧的话能过有效的降低应力集中。 在中间的加载力，使用一个面的分布力加载也能够有效的简单载荷力处的应力集中。 重点是使用MPC建立自由度之间的约束，消除结构的奇异点。只留下y方向的自由度，其余方向的节点自由度固结掉。在设置是自由度是可以多选的。有限元分析模型分析过程常见错误分析： 由于在软件操作过程中，软件自身在内存中会存储大量的数据。如果在短时间给操作了好几步，就会造成数据bdf生成出错。比如说我遇到的错误就是在dbf文件中，载荷力标签没有加上。所以在生成Nastran提交文件后，可以查看一下bdf文件，看一下所有标签是否都已经出现。查看Nastran解算结果中的错误：对应的模型： 查看.06文件，这个文件中的第一个错误点是，由于Spring单元和三角形单元是不相容的。在设定Spring单元时，Spring单元的节点自由度只有一个，而三角形单元的节点自由度有六个。那么就出现了GRID POINT SINGULARITY TABLE中所示的奇异点。造成矩阵无法计算。优化分析结果： 由于只是练习，模型使用的是不加Spring并且使用集中载荷的有限元模型。优化结果如下： 重量的单位为T（吨）；板厚的单位为mm（毫米） 优化结果为：重量1.722*10^-3（吨），四个板厚分别为0.1,0.1,0.1,1.5mm。