球磨机实际结构中盖板、筋板、支板是焊接在轴上的，有限元分析中分别采用了壳单元与体单元来模拟，为保证交界面上的位移协调，顺利传递力及力矩等，就必须使单元之间联结。
    对于同种单元的不同类型有很多种处理方法，如进行简单的形函数插值等。但对于不同类型的单元，由于单元节点型式与配置的不一致，简单的插值则不能胜任。如壳单元与体单元，壳单元节点具有六个自由度，体单元节点只有三个平动自由度。进行简单的节点自由度耦合仅仅能耦合相同的三个平动自由度，那么，壳单元的转动自由度则会无约束，从而导致求解错误，或者根本求解不出来。
    解决这一问题通常采用多点约束方程法，即构建连接区域节点的约束方程。这一方法的基本思想是保证交界面上位移的协调性，恰当地提出不同类型单元在交界面上的协调条件。Ansys软件中将这种方法与先进接触计算功能的线性接触算法结合，由程序内部定义多点约束方程(MPC)，如图3-6，能极方便地完成实体一壳一梁模型的混合“装配”。图中的影响距离表示实体单元内的节点，凡与壳单元任意边界的距离在影响距离之内的都用来构建约束方程。
    采用这种方法，在球磨机的壳与体连接区域，壳上部分定义为接触单元175，轴上部分定义为目标单元170，实现了壳单元与体单元顺利连接。
    较后，得到了陶瓷球磨机的离散化有限元模型如图3—7所示。其中壳单元有11008个，实体单元有7822个，接触单元有1080个，目标单元有1608个，总共单元数21518个，节点数为46861个。