请问有搞有限元 计算的么

Quark的虎牙

用ABAQUS搞有限元计算，可是由于是新手，想找人讨论一些问题。
$送花$ $送花$

klaus1

me, whats up.
wangjiasheng@hotmail.com
my msn

joytuhh

我也用到一点，不过很少，可以一起交流一下。
joytuhh@hotmail.com

eisenstange

可以考虑建学习小组了

康熙大帝

你想问啥??

Quark的虎牙

原帖由 康熙大帝 于 2007-6-11 15:45 发表
你想问啥??

$m9$ 这个问题还真不好回答，就是计算中会碰到些问题，希望有经验的能指点一下 $m12$

康熙大帝

不明白你到底要知道啥

Quark的虎牙

原帖由 康熙大帝 于 2007-6-11 22:04 发表
不明白你到底要知道啥

接触上的一些问题。看来您是高手，指点一下把。
我做了刚体和变形体的接触问题模型，给刚体定义了下降位移，可是等计算完了，并没有位移出现，也没有应力产生，请问是什么原因呢

康熙大帝

那你刚体和变形体的月约束咋做的？

Quark的虎牙

原帖由 康熙大帝 于 2007-6-12 12:04 发表
那你刚体和变形体的月约束咋做的？

变形体底面4点，ENCASTRE
刚体参考点 Z 方向给位移

康熙大帝

End ABAQUS/Standard Analysis
ABAQUS JOB ******** COMPLETED
你看看。log，有这些文字否。
如果没有，看。dat里有没有错误提示。
如果没有，你按plot deformed shape试试

[ 本帖最后由 康熙大帝 于 2007-6-12 14:18 编辑 ]

Quark的虎牙

原帖由 康熙大帝 于 2007-6-12 15:10 发表
End ABAQUS/Standard Analysis
ABAQUS JOB ******** COMPLETED
你看看。log，有这些文字否。
如果没有，看。dat里有没有错误提示。
如果没有，你按plot deformed shape试试

我发现跟刚体的约束有关，如果我这样定义刚体

U1=U2=UR1=UR2=UR3=0, U3=-2E-5
There is zero MOMENT everywhere in the model based on the default criterion. please check the value of the average MOMENT during the current iteration to verify that the MOMENT is small enough to be treated as zero. if not, please use the solution controls to reset the criterion for zero MOMENT.
应力分布跟期望相似

但是如果我这样定义刚体

U1,U2,UR1,UR2,UR3不定义，U3=-2E-5

The system matrix has 1 negative eigenvalues.

Solver problem. Zero pivot when processing D.O.F. 4 of 1 nodes. The nodes have been identified in node set WarnNodeSolvProbZeroPiv_4_2_2_1_1.

Solver problem. Zero pivot when processing D.O.F. 5 of 1 nodes. The nodes have been identified in node set WarnNodeSolvProbZeroPiv_5_2_2_1_1.

Solver problem. Zero pivot when processing D.O.F. 6 of 1 nodes. The nodes have been identified in node set WarnNodeSolvProbZeroPiv_6_2_2_1_1.

没有应力没有变形，计算后什么都没发生

这是为什么呢。
而且如何解决情况一下的warning呢

$送花$

康熙大帝

你的长度单位是mm？
然后位移-2E-5？？:o
你的刚体和变形体是接触的？

Quark的虎牙

原帖由 康熙大帝 于 2007-6-12 17:14 发表
你的长度单位是mm？
然后位移-2E-5？？:o
你的刚体和变形体是接触的？

长度单位是米
刚体和变形体一开始是分开的，然后接触，然后压痕
$郁闷$ 有什么不对么
ABAQUS默认的单位不是米么，我看说明是默认国际单位制阿

康熙大帝

abaqus没有单位，单位完全按照你的input file给的数据。
你的两个体之间初始的距离是多少？ 方向对吗？
你不妨用mm当单位试试,你给的值也许小于缺省的TOLERANCE

Quark的虎牙

我是用的CAE，不会写INP文件，怕写不好。
看了下生成的INP文件，上面没有单位，如果是单独写INP的话，是可以定义单位的。但是我的问题是我要跟实验结果想比较，所以是需要单位的，如果在CAE中无法实现单位一说的话，就比较麻烦。我是不是可以在CAE生成的INP文件中加入对单位的定义，然后打开，运行。没这样做过，不知道行不行$郁闷$

老虾

ABAQUS对单位的要求就是整体一致.比如使用米,牛,那么应力就是帕,如果使用毫米,牛,应力就是兆帕.工程上常用的是mm,MPa单位,在涉及质量和能量时,单位换算比较麻烦.
rigid body的位移是加在RP上,RP有6个DOF,这6个DOF必须全部加以约束,否则就有rigid motion,所以第二种BC定义肯定错的,zero pivot就是欠越束导致的rigid motion.
第一种BC定义,应该是由于过小的位移.2e-5,也就是20微米.大多数情况下,远小于单元大小,不足以发生接触.
我不知道你用的是什么单元,俺用continuum element的时候,两个part在接触之前也有zero force的warning,你现在是zero moment,看起来比较类似.

至于单元,如果换成mm制,整个模型需要重新修改,看你这个模型是否需要继续用下去,还有就是涉及密度,热等方面的时候,材料属性切记切记切记要小心.