基于薄壳有限元法的带旋转自由度膜元研究-材料加工工程专业论文.docxVIP

摘要橡塑模具制品一般具有厚度方向尺寸较面内尺寸很小且形状复杂等特点，所 摘要 橡塑模具制品一般具有厚度方向尺寸较面内尺寸很小且形状复杂等特点，所 以在对制品进行翘曲变形分析时，就可采用薄壳有限元法。目前用于薄壳翘曲变 形分析的单元主要有平板壳元、退化壳元和多变量壳元，其中平板壳元应用最为 广泛，也最为成熟。所谓的平板壳元就是用许多小平板元构成的多面体去逼近真 实的曲面，许多著名软件如ABAQUS、ADINA、ANSYS、MARC、NASTRAN、 SAP的壳体分析都主要应用平板壳单元。但是平板壳元也有自身的问题：膜刚度 近似差，膜板结合为壳单元和一维单元与二维单元连接时困难。 解决上述问题的最好办法无疑是在膜元中引入垂直于该膜单元平面绕单元 平面法线旋转的旋转自由度。关于此类单元的构造，Allman采用的插值方法得 到广泛的采用，具有里程碑的意义。目前，引入旋转自由度后的膜元都不同程度 的提高了膜单元的计算精度，但是在一些带旋转自由度的膜元中存在一多余的零 能模式，几乎所有此类膜元所引入的旋转自由度都不具有连续介质力学中关于旋 转自由度的定义，同时鲜有对旋转自由度自身的正确性进行考核。因此在解决一 维与二维连接问题时结果很不理想、甚至是没法解决，并且有些单元推导过程过 于繁琐、不实用。 为了给注塑模CAE中的制件翘曲分析提供一个可靠的单元，基于以上带旋 转自由度膜元存在的问题，本文在生成平板壳元时的工作主要集中在对带旋转自 由度膜元的研究上，总结如下： 1．在6节点三角形单元的基础上，通过对切向位移一次插值和法向位移二 次插值消除了边中间节点，成功地引入了旋转自由度。同时引入了一个 待定系数a，并由连续介质力学中关于转动分量的定义确定其取值，同 时给出了一个附加方程，进而得到一个简洁有效的带有旋转自由度的平 面三角形单元，记为OAT3。单元推导简单，列式简洁。具体算例表明， 此单元无多余的零能模式，仍保持了原有Allman单元的精度，引入的 旋转自由度具有连续介质力学关于转动分量的定义，且首次对旋转自由 度自身的正确性进行了校核。 2．在10节点三角形单元的基础上，通过对切向位移一次插值和法向位移 三次插值，成功地消除了边中间节点和单元中心节点，推导得到另一个 带旋转自由度的3节点三角形单元，记为AT4。算例表明，此单元无多 余的零能模式，精度较无旋转自由度的膜元有所提高，不足之处是单元 较刚。 3． 生成了基于离散Kirchhoff理论的板弯曲单元和空间一维杆单元。具体 算例证明了以单元OAT3为膜元的平板壳元具有良好的计算精度，很好 的解决了膜板结合为壳单元和一维单元与二维单元连接时的问题。 关键词：旋转自由度； 零能模式； 膜元； 三角形单元； 有限元； 平板壳元； 薄壳； III AbstractMost Abstract Most injection—molded parts are complicated thin—shell like structure，SO the thin shell finite element method is usually used to analyze the warpage and deformation of the parts．Generally，we have several types of thin shell elements to be used：flat-shell element，degenerated shell element and multi—variable shell element，of which the flat-shell element is adopted most widely due to its simplicity and is also selected by most famous finite element softwares，such as ABAQUS，ADINA，ANSYS，MARC， NASTRAN and SAP．However，the fiat-shell has its own disadvantages：The precision of membrane element is not perfect and it is difficult to connect one—dimensional element and two．dimensional element． To resolve the problems，some scholars thought that the drilling degree of freedom should be added to membran