压力载荷，结果是没有任何压力效果，为什么啊？

zlk2010

下面是命令流，用ansys/ls-dyna做的，直接复制命令流就行
初学者了呵呵，一个非常简单的过程，在板上某区域施加力，目的是产生塑性变形，但是结果是没有任何变形和压力梯度，就是算了和没算一样，搞不懂为什么，帮帮忙啊，应该不是很难，但是我搞了一整天了~~~

1. 
   
2. /clear
   
3. /filname,shuikongchongji
   
4. /title,chongji
   
5. /prep7
   
6. !国际单位制kg-m-s
   
7. /units,si
   
8. !单元定义
   
9. !*  
   
10. ET,1,SOLID164   
    
11. !*
    
12. !实体模型建立
    
13. cylind,0,0.0003,-0.0001,0,90
    
14. !划分网格
    
15. /PREP7  
    
16. !*  
    
17. LESIZE,ALL, , ,20, ,1, , ,1,
    
18. CM,_Y,VOLU  
    
19. VSEL, , , ,       1
    
20. CM,_Y1,VOLU
    
21. CHKMSH,'VOLU'   
    
22. CMSEL,S,_Y  
    
23. !*  
    
24. MSHAPE,0,3d
    
25. MSHKEY,1
    
26. VMESH,_Y1   
    
27. MSHKEY,0
    
28. !*  
    
29. CMDELE,_Y   
    
30. CMDELE,_Y1  
    
31. CMDELE,_Y2  
    
32. !*  
    
33. !材料力学参数
    
34. /PREP7
    
35. MPDE,DENS,1   
    
36. MP,DENS,1,7.8  
    
37. MPDE,EX,1
    
38. MP,EX,1,90000
    
39. MP,NUXY,1,0.36
    
40. TBDE,BISO,1   
    
41. TB,BISO,1,,,,
    
42. TBDAT,1,100   
    
43. TBDAT,2,80
    
44. 
    
45. !建立待加载的单元组件
    
46. allsel,all
    
47. csys,1
    
48. esel,s,live
    
49. nsel,r,loc,x,0,0.0001
    
50. nsel,r,loc,z,-0.000005,0
    
51. esln,r,1
    
52. CM,E,ELEM
    
53. !建立时间压力表
    
54. *DIM,time,ARRAY,5,1,1, , ,  
    
55. !*  
    
56. *SET,TIME(1,1,1) , 0.1  
    
57. *SET,TIME(2,1,1) , 0.2  
    
58. *SET,TIME(3,1,1) , 0.3  
    
59. *SET,TIME(4,1,1) , 0.4  
    
60. *SET,TIME(5,1,1) , 0.5  
    
61. *DIM,press,ARRAY,5,1,1, , ,
    
62. !*  
    
63. *SET,PRESS(1,1,1) , 100
    
64. *SET,PRESS(2,1,1) , 300
    
65. *SET,PRESS(3,1,1) , 500
    
66. *SET,PRESS(4,1,1) , 300
    
67. *SET,PRESS(5,1,1) , 100
    
68. !施加压力载荷
    
69. /solu
    
70. EDLOAD,ADD,PRES,6,E,PRESS,TIME, 0, , , , ,
    
71. !施加约束
    
72. allsel,all
    
73. /solu
    
74. FLST,2,1,5,ORDE,1   
    
75. FITEM,2,1   
    
76. !*  
    
77. /GO
    
78. DA,P51X,ALL,0
    
79. !求解时间
    
80. /solu
    
81. TIME,0.6,
    
82. allsel,all
    
83. solve
    
84. fini
    

复制代码

D_eagle

1# zlk2010
复制下来里面有很多乱码，建议直接上K文件。

xiaopiujy

1# zlk2010
复制下来里面有很多乱码，建议直接上K文件。
D_eagle 发表于 2011-7-6 08:06

d兄可以编辑楼主的帖子，把命令流复制下来！

wyc412721

我编辑了一下，可以复制代码

D_eagle

还是xiaopiujy兄想得周到~不明白lz的载荷曲线为什么是下图这样的，载荷曲线的横纵坐标搞反了吧。

zlk2010

哦哦  不好意思，的确是反了，但是改过来之后，好像还有问题

ywluck

建立好曲线数组之后，一般用画出来看一下：
edcurve,add,1,time,press
edcurve,plot,1

然后选择要施加压力的单元，举例如下：
nsel,s,loc,y,850
ESLN,S                  
cm,topelem,elem
/PSF,PRES,NORM,2,1,1      !以箭头方式显示面荷载
edload,add,pres,,topelem,,,0,1,100,,

zlk2010

为什么求解要十几个小时啊？？二维只需要两分钟！可是我是需要转隐式单元，必须用三维的，这个计算时间该怎么解决啊？一天不用干别的了

52orange

就这个问题干吗要用显示呢？隐式计算完全可以搞定啊。而且你在这个系统中建议加一定的阻尼，否则容易出现数值振荡而导致增加无谓的计算时间。

zlk2010

开始的时候是用隐式算，二维建模，发现网格畸变，严重失真，都没法看了。后来导入dyna 一下就成了，看起来还挺不错，然后变成三维建模，又有了更多问题，晕啊，受不了