关于ansys问题求助

world20110

大神们，刚学ansys没多久，需要用ansys Mechanical研究anand本构方程,因为只需要试一下ansys使用anand本构进行热-结构耦合分析的整个流程，验证其可行性即可，因此建立建立一个简单模型，对于具体数据的准确性也没有深究，是从一些文献里面找到的。使用的ansys软件为12.0版本。其命令流如下所示： /prep7 ET,1,SOLID87     !单元类型 MP,DENS,1,1.5e-6    ! pcb  密度 MP,EX,1,11e6        !弹性模量 MP,KXX,1,15000      !导热系数 mp,nuxy,1,0.28      !泊松比 mp,alpx,1,1.5e-5    !热膨胀系数 MP,DENS,2,8.42e-6        ! 焊料 MP,EX,2,31.5e6 MP,KXX,2,57000 mp,nuxy,2,0.4 mp,alpx,2,2.4e-5 TB,ANAND,2                  ! VISCOPLASTIC MODEL BY ANAND TBDATA,1,12.41E3 TBDATA,2,9400 TBDATA,3,4.1E6 TBDATA,4,1.5 TBDATA,5,0.303 TBDATA,6,1379E6 TBDATA,7,13.79E6 TBDATA,8,0.07 TBDATA,9,1.3 block,0,10,0,20,0,5 block,2,8,2,18,5,10 模型如下： 网格划分采用自由网格划分，设置为6和8. 问题：1、在温度分布计算阶段，能够的到完成温度分布的计算，得到温度分布图。然后进行后续的热循环应变分析的时候，却在current   LS的时候自动停止，查看错误文件，显示错误为： *** ERROR ***     SUPPRESSED MESSAGE    CP =      13.812   TIME= 19:24:09 Rate-dependent plasticity: flow function for Anand not defined for       element 60 and material 2.         2、在查阅一些文献资料后得知anand对单元类型有一定的要求。查看自带帮助文档将其中一个单元类型改为solid185. 但在温度分布计算时，选择分析类型的时候，会和前面使用solid87不同，但是也能够得到温度分布图。   后面进行热循环应变分析的时候，却发现无法找到rth文件。            希望各位前辈不吝赐教，本人ansys用的不久，但以后将持续用下去。如果问题表述的不够清楚，还希望能回帖指正。谢谢！

我来回答

htbbzzg

1  你的密度的值不对，应该减小 3 个量级，如：1.5e-9、8.42e-9，等。
2  你没有定义材料的比热，瞬态热分析的结果是不对的。
3  solid87 是 10 节点四面体热单元，而 solid185 是 8 节点六面体结构单元，二者不能直接替换。
4  从热分析转到结构分析，需要执行命令：ETCHG, TTS，相应菜单路径是：
    Main Menu>Preprocessor>Element Type>Switch Elem Type
    不知你是否执行过。如未执行，可以先看一下帮助文件中的有关内容。
因为对你的分网、加载、求解过程等不了解，只能先说这一些了。

world20110

htbbzzg 发表于 2014-4-4 09:09
1  你的密度的值不对，应该减小 3 个量级，如：1.5e-9、8.42e-9，等。
2  你没有定义材料的比热，瞬态热分 ...

您好！谢谢您的回复。
1、我的密度值一直都是用的这个啊，以前没有出现问题，只是在现在将本构方程换为anand之后才出现现在的问题。另外为什么会与密度值有关，请教一下？呵呵···
2、我使用的是热稳态分析，再进行与结构分析的耦合，因此也就没有使用到比热容。
3、可能我没有表述清楚，我没有将solid87直接替换为solid185，因为在ansys帮助文档里看到anand本构对单元类型有要求，所以将要使用anand本构方程的体的单元类型改为了solid185，其它体的单元类型没变。主要是改了之后，进行稳态热分析，能够得到温度分布，却不能得到rth文件，很纳闷？
4、你说的这一步已经做过了。

htbbzzg

    从你提供的一些材料数据来看，你使用的单位应该是： 长度 - mm， 时间 - 秒，弹性模量 (以及：压力、应力等) 的单位是 MPa。
    在这样的情况下，质量的单位应该是 T (吨)，而不是 Kg。你给的密度值 1.5e-6， 其单位应该是 Kg/mm^3，需要转换为 T/mm^3，即相应的值应该是 1.5E-9。
    这里涉及一个单位协调的问题，附件中有一份资料，你可以参考一下。    有限元中的材料性能单位.doc (149.5 KB, 下载次数: 15)

2014-4-4 15:48 上传

点击文件名下载附件