- 积分
- 3
- 注册时间
- 2011-9-30
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2012-6-29 21:01:46
|
显示全部楼层
来自 陕西西安
zwqbill 发表于 2012-6-29 21:00
有子结构,帮助文档e卷的8.23例子。。就是没有软件设置过程。。
可以参考这个帖子,http://forums.simwe. ...
不让发连接。。
论坛现在传不上图片,过几天再试.抱歉.不过应该也可以看懂的. 关于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% 以上的计算时间,模型越大,效果越名显,可以看以一下结果对比图.【转载CAENET】 |
|