MARC子结构应用实例(原创,)

CAE联盟( http://forums.caenet.cn,) NaCl原创,转载请注明出处.
论坛现在传不上图片,过几天再试.抱歉.不过应该也可以看懂的.
关于MARC中子结构的应用,帮助文档里有三个例子,比较典型的是齿轮分析和两个圆柱体接触分析的
例子,齿轮的例子中指出,完整模型求解为1002S,利用子结构方法求解,总时间为264S,效率得高近
3/4,可见效果较好,本人最近想利用子结构解决问题,苦于没有参考书籍及文献,文献中多以ANSYS的
子结构进行分析,而MARC的帮助文档中没有一个完整的操作设置过程,在网上搜到的也都是求助,着
实为难了一番,研究数日,小有收获,与大家分享.费话不多说,上例子。
参照两个圆柱接触的例子，做了两个长方体接触的模型，见图1，下端固定,上端X向固定,Y向施加向下位移,
其它完整模型分析不详述，求解共计32S。
子结构求解步骤：
1，将每个变形体分为两部分，分别为子结构部分和超单元部分，如图2。
2，求解超单元部分的刚度矩阵。将完整模型文件另存一个,如sub1,只保留上部变形体的超单元部
分，其余单元删除，在超单元上边界和下边界（与子结构相连的节点）添加边界条件DOF-SET
NODES，固定其X，Y两个方向的位移。命令：boundary conditions-〉structural-〉more-〉dof-
set nodes-〉displacement x/y。
3，超单元设置。在loadcase中，将类型设置为superelement。命令：loadcase-〉type-〉
superelement/
在proerties中，选中dof-set nodes，并在matrix name中设置一个矩阵名称（很重要，后面要用
到）。
4，设置DMIG输出，在JOBS中的DMIG OUT菜单中，选中stiffness matrix，在弹出菜单中，依次选
中element stiffness/global stiffness 及frequency中的once，coord system中的global，
elements中的all elements。
5，保存模型进行求解，计算时间约5S，正常退出号为3022。
6,相同的方法,把完整模型另存sub2,将下部变形体的超单元部分刚度矩阵求出,注意matrix name中
设置一个矩阵名称要换成不同于上部的,在文件夹中会分别生成两个刚度矩阵的存储文件
XXX_XXX_dmigst_0000
7,两个超单元的刚度矩阵求出后,就开始计算关键部位了.再把完整模型另存为sub,将定义为超单元
的两个部分删掉,注意将带有边界条件的节点保留.
8.在JOBS菜单中,进入include file菜单,在下面的model section栏里选add file,分别将上面生成
的两个刚度矩阵存储文件导入.OK确认.
9.进入input file text菜单，依次输入如下4行命令：
k2gg，XXX（第一个刚度矩阵名称），1，
include，XXX_XXX_dmigst_0000（第一个刚度矩阵存储文件名称）
k2gg，XXX（第二个刚度矩阵名称），1，
include，XXX_XXX_dmigst_0000（第二个刚度矩阵存储文件名称）
10.接下来就可以求解了,时间约为5S,这样求解刚度矩阵各5S,主程序计算5S,共计约15S,节省了50%
以上的计算时间,模型越大,效果越名显,可以看以一下结果对比图.