[求助] 关于扭矩的问题

dqpi-wht

我前一阵做过一个橡胶圆盘扭转的问题（半径254mm，高度25.4mm,一端固定，另一个端面想到固定端面旋转10度，且让圆盘的高度不变，求轴向反力及端面扭矩，最后再把轴向反力与扭矩和理论值进行比较，选用单元186，材料用neo－hookean，c1＝0.690mp，即剪切模量为1.379mp，认为是不可压缩材料），做这道题的目的只是为了和理论结果对照一下，因为涉及到加扭矩，所以问过很多人如何加扭矩，我得到的答复有很多种：
1 转化成一对力偶，由于橡胶是软材料，我想如果转化成一对力偶可能不太合适，所以没这么做
2 把扭矩转化成面上剪切力，我自己不知道如何转化，结果也是没做成
3 在国外有限元网站上也问过，得到的回答是用mechanical toolbar加，因为在mechanical toolbar上可以直接加扭矩，我对ansys也不是学的很好，去mechanical toolbar上试过一次，好像只能算线性的，但是橡胶是非线性的，不知道能不能加了扭矩在回到ansys里继续算，我自己后来也没做了
4 有个老外也建议用表面效应单元，加剪切力，和上面转化成面上的剪切力一样吧
5 给ansys公司的技术人员发邮件请教，他说用刚性接触面加，但是结果扭矩和理论值接近，但是轴向力却很小，
6 建立刚性面直接加扭矩或转角，把一个端面建成刚性面，把载荷加在一个质点单元上，这个我做过了，扭矩算的很理论值很接近，但是得到轴向力比理论值几乎大了一倍
7 直接把端面节点的坐标系改为柱坐标，把转角加到节点上，算出的扭矩和轴力都要小于理论值，小很多，不明白究竟为什么，也不知道这个算法和方法6有什么区别，为什么结果却不一样
上面这些算法都要让轴向高度不变
哪位大侠有空的时候可以帮俺算算，上面这些方法都有那些毛病，哪种方法是最适合这个问题的，请指点
这里还有另一个问题就是在为什么我加转角位移时用degree和rad，同样的角度，为什么得到的结果不一样，ansys技术人员说，要用rad不能用degree，为什么啊

还要补充一下：
通过上面的方法，至少有扭矩和理论值接近了，为什么轴向反力总是不对，不知道是不是我自己查看的方法不对，请大家指点。
求解完成后，我选中施加扭矩的那个端面上的所用节点，用菜单
General Postproc >nodal calcs > Total Force Sum
可以看到一些节点力和扭矩的和，不知道这样做对不对？

freefree

一般都是用MASS偶合的方式做，不过你说的轴向力问题确实没有关注过，请高手门给讲讲吧