FWD模拟

kuresu

我想模拟FWD的落锤撞击铺面板，但是contact的部份一直弄不好，msg档都显示不收敛或是overclosure，想请问有人可以弄个input给我参考参考吗??拜托拜托…

zhangmeng

overclosure是很正常的，是在计算过程中还是初始增量就有问题
是显式还是隐式

kuresu

这是我的input..一直不知道那里有错..

                               
登录/注册后可看大图

+Input file
**
**To compare with the static central load
**use *Dload to compare with the exact solution                               
**
*********************************
***** prepare for the model *****
*********************************
**
*NODE
100000,0,0,0
100082,410,0,0
182000,0,410,0
182082,410,410,0
*NGEN,NSET=LEFT
100000,182000,1000
*NGEN,NSET=RIGHT
100082,182082,1000
*NFILL,NSET=DECK_1
LEFT,RIGHT,82,1
*NCOPY,CHANGE NUMBER=100000,OLD SET=DECK_1,SHIFT,NEW SET=DECK_2
0,0,-5
0., 0., 0., 0., 0., 1., 0.
*NCOPY,CHANGE NUMBER=200000,OLD SET=DECK_1,SHIFT,NEW SET=DECK_3
0,0,-10
0., 0., 0., 0., 0., 1., 0.
*NODE,NSET=B1
5001,200,200,5
5002,205,200,5
5003,210,200,5
5004,210,205,5
5005,210,210,5
5006,205,210,5
5007,200,210,5
5008,200,205,5
5009,205,205,5
6001,200,200,10
6002,205,200,10
6003,210,200,10
6004,210,205,10
6005,210,210,10
6006,205,210,10
6007,200,210,10
6008,200,205,10
6009,205,205,10
7001,200,200,15
7002,205,200,15
7003,210,200,15
7004,210,205,15
7005,210,210,15
7006,205,210,15
7007,200,210,15
7008,200,205,15
7009,205,205,15
*NODE,NSET=A1
500001,200,200,0
500002,210,200,0
500003,210,210,0
500004,200,210,0
**
**************************************
***** this is a pavement element *****
**************************************
**
*ELEMENT,TYPE=C3D20R
1,300000,300002,302002,302000,100000,100002,102002,102000,300001,301002,
302001,301000,100001,101002,102001,101000,200000,200002,202002,202000
*ELGEN,ELSET=PAVEMENT
1,41,2,1,41,2000,41,1,0,0
*ELEMENT,TYPE=C3D8R,ELSET=BLOCK
6001,5001,5002,5009,5008,6001,6002,6009,6008
6002,5002,5003,5004,5009,6002,6003,6004,6009
6003,5009,5004,5005,5006,6009,6004,6005,6006
6004,5008,5009,5006,5007,6008,6009,6006,6007
6005,6001,6002,6009,6008,7001,7002,7009,7008
6006,6002,6003,6004,6009,7002,7003,7004,7009
6007,6009,6004,6005,6006,7009,7004,7005,7006
6008,6008,6009,6006,6007,7008,7009,7006,7007
*ELEMENT,TYPE=SPRING2,ELSET=RUBER
3001,5001,500001
3002,5003,500002
3003,5005,500003
3004,5007,500004
*SPRING,ELSET=RUBER
3,3
788000000
*ELEMENT,TYPE=DASHPOT2,ELSET=DAMPING
4001,5001,500001
4002,5003,500002
4003,5005,500003
4004,5007,500004
*DASHPOT,ELSET=DAMPING
3,3
500
*ELEMENT,TYPE=S4R,ELSET=CONTACTFACE
60002,500001,500002,500003,500004
**
*******************************
***** Material definition *****
*******************************
**
*SOLID SECTION,ELSET=PAVEMENT,MATERIAL=PAVEMENT
*FOUNDATION
PAVEMENT,F2,100
*MATERIAL,NAME=PAVEMENT
*ELASTIC
4E+06,0.15
*DENSITY
2410.2
*SOLID SECTION,ELSET=BLOCK,MATERIAL=BLOCK
*MATERIAL,NAME=BLOCK
*ELASTIC
2E+11,0.27
*DENSITY
7830
*SHELL SECTION,ELSET=CONTACTFACE,MATERIAL=CONTACTFACE
0.001
*MATERIAL,NAME=CONTACTFACE
*ELASTIC
4E+06,0.33
*DENSITY
8450
**
*******************************
***** Boundary conditions *****
*******************************
**
*NSET,NSET=X1,GENERATE
100000,100082,1
200000,200082,1
300000,300082,1
*NSET,NSET=X2,GENERATE
182000,182082,1
282000,282082,1
382000,382082,1
*NSET,NSET=Y1,GENERATE
100000,182000,1000
200000,282000,1000
300000,382000,1000
*NSET,NSET=Y2,GENERATE
100082,182082,1000
200082,282082,1000
300082,382082,1000
*BOUNDARY
500001,1,2,0
500002,1,2,0
500003,1,2,0
500004,1,2,0

*****************************
***** Surface & Contact *****
*****************************
*SURFACE,NAME=CONTACTFACE,TYPE=NODE
A1,0.001
*SURFACE,NAME=PAVEMENTFACE
PAVEMENT,S2
*CONTACT PAIR,INTERACTION=IMPACT
CONTACTFACE,PAVEMENTFACE
*SURFACE INTERACTION,NAME=IMPACT
*********************
***** LOAD CASE *****
*********************
*INITIAL CONDITIONS,TYPE=VELOCITY
B1,3,-2
*STEP,INC=5000
*DYNAMIC
0.1,10
*DLOAD,OP=NEW
BLOCK,GRAV,10,0,0,-1
*********************************************************
*RESTART,WRITE
*EL PRINT,POSITION=AVERAGED AT NODES
*NODE PRINT
U
*END STEP

[ 本帖最后由 kuresu 于 2007-3-27 09:21 编辑 ]

zhangmeng

你可以改用显式动力学看一下运动过程

kuresu

想请问一下有办法看到落锤在整个运动周期中产生的力吗??我只会看位移…想看落锤的力，不是应力喔，是落锤落下打到地面整个时间产生的作用力历时…

zhangmeng

用刚体，输出参考点的力

kuresu

我是用standard去做的，我只会standard…standard不是不能设刚体吗? 而且我就是不知道怎么让一个node或element输出力   Orz

zhangmeng

谁说的！刚体可以用在显式和隐式中
只要消除刚体位移就可以了

kuresu

喔…抱歉，那请问要怎么输出参考点的力呢?指令是什么? 谢谢指教

zhangmeng

设定好了刚体，直接就可以从参考点上得到了

kuresu

我的问题是…我不知道input里的指令要用什么才能输出作用力。ABAQUS Viewer里的选项我也只找的到位移和应力等，就是找不到可以看作用力历时的选项

zhangmeng

在history output里设定RF不就行了吗？

kuresu

RF不是反力吗?

zhangmeng

作用力和反作用力是大小相等方向相反的
哥哥，你也太狠了！
非要刨根问底