关于ansys问题求助
大神们，刚学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用的不久，但以后将持续用下去。如果问题表述的不够清楚，还希望能回帖指正。谢谢！
我来回答
1你的密度的值不对，应该减小 3 个量级，如：1.5e-9、8.42e-9，等。
2你没有定义材料的比热，瞬态热分析的结果是不对的。
3solid87 是 10 节点四面体热单元，而 solid185 是 8 节点六面体结构单元，二者不能直接替换。
4从热分析转到结构分析，需要执行命令：ETCHG, TTS，相应菜单路径是：
Main Menu>Preprocessor>Element Type>Switch Elem Type
不知你是否执行过。如未执行，可以先看一下帮助文件中的有关内容。
因为对你的分网、加载、求解过程等不了解，只能先说这一些了。 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、你说的这一步已经做过了。
从你提供的一些材料数据来看，你使用的单位应该是： 长度 - mm， 时间 - 秒，弹性模量 (以及：压力、应力等) 的单位是 MPa。
在这样的情况下，质量的单位应该是 T (吨)，而不是 Kg。你给的密度值 1.5e-6， 其单位应该是 Kg/mm^3，需要转换为 T/mm^3，即相应的值应该是 1.5E-9。
这里涉及一个单位协调的问题，附件中有一份资料，你可以参考一下。
页:
[1]