MSC.Patran中体单元与壳单元的连接方法的探索.docVIP

MSC.Patran中体单元与壳单元的连接方法的探索* 陈继华 杜家政 隋允康 管昭 （北京工业大学数值模拟中心） 摘要：本文对用MSC.Patran建立体壳混合模型时，就怎样使体壳连接边的节点相互对应协调；以及怎样解决因为体壳单元自由度不同而使得两者之间用Equivalence（节点相等）不能固定的问题进行了些使用上的摸索；同时对该软件在计算应力集中问题的可靠性进行了检验，为用本软件建模提供了些可借鉴的方法。 Research of joint solid-element shell element in PATRAN Jihua Chen Jiazheng Du Yunkang Sui Zhao Guan (Numerical Simulation Center for Engineering Beijing Polytechnic University) Abstract：This paper is about how to make nodes of the mutual edges of the solids and the shells to assort with each other and how to solve the problem which is aroused by the difference of the free degrees of the solid and the shell so that “Equivalence” could not be used to solve it, when the compound model is founded by using MSC.Patran. Moreover, the paper makes some practices to test the software’s reliability of the calculational stress concentrate and gives some referenced methods for modeling with the software. 前言 在对一些工程实际问题建模时发现不能用单一的某种单元来处理，而是体、壳、梁或杆的混合模型。但由于这些单元的自由度上的差异造成点对点之间不能用Equivalence直接连接，比如体单元只有平动三个自由度，而壳单元有五个自由度，用Equivalence只能保证两者的连接没有相互移动，不能保证两者之间没有转动。本文就这个问题就怎样使用MPC （multitude point control）做了些实践性的摸索。 二、模型的建立和单元的划分 在建立单一的体、壳模型时，只要保证各个几何实体之间有公共的边、面或用Associate使之相互协调，以保证划分网格后用Equivalence命令来得到期望的有限元模型。但要建立体壳混合的有限元模型时，仅仅这样是不够的（对于B-rep体），因为这样并不能保证体和壳的公共边上的节点相互对应。通过反复的摸索发现可以在先保证体边和壳边协调的前提下，对体的表面及壳面先划分单元，然后再对体划分单元，这样体壳公共边上有位置相同的接点，如图1，然后删除体上的面单元或对这些单元不加属性。同理，当体的表面是比较复杂的曲面时，为了保证有限元模型与几何模型更接近也得用这种方法。 三、用MPC点连接体壳 建立如图2所视的体壳混合模型，体的长宽高分别为10、3、4，面的长为20，宽为3，两者的过度处用圆角连接，在壳的最右边施加与壳面垂直的均布力，其和为44，体的左端面固定。要实现体和壳的连接可以选用RSCCON Surf-Vol和Rigid之一。两者有所不同，前者dependent点是一项，independent点是两项，后者dependent点的项数不定，但independent点是一项；但两者每项的点数都必须相等，都必须在体上选至少两组点。前者必须使壳上的点作为dependent中的项，体上的两组点都是independent中的项；后者壳上的点可以在dependent中也可在independent中。 1. 交接面单元的不同划分的比较 这里首先对体壳实体模型的交接面（这个面一定是矩形）的网格划分进行试探，如下： 图3和图4的实体模型完全一样，但图3在交接面分了两行节点（没有第2行），而图4分了3行节点，如图1。然后都用RSCCON Surf-Vol方法建立MPC点，体上选用第1行和第3行的节点作为其independent项，第3行壳上的节点作为dependent项。计算结果如图3、4，其最大应力分别为1720和1660（不考虑单位），都在壳与体的连接处。但材料力学解为1716（图3是用科学记数法保留了3位有效数字