【求助】厚壁封头的冲压问题！

bizarrEaster

各位大大好，小弟刚学ANSYS一个礼拜，毕业设计中要求做厚壁封头的冲压，模型如下图：
模型中上、下模具采用mesh200单元分网，然后扫掠成solid164单元（根据书的冲压例子做的），而板料直接用solid164单元分网。

没有人教，有很多问题搞不太懂，请教各位大大：
1. 由于模型是轴对称的，我只建立了如上图所示的1/4模型是否可行？
2. 不知道我的分网方式对吗？听说板料要用自适应网格，可是好像只适合壳单元...
3. 上、下模具采用刚体，在材料属性已经约束自由度了，那还需不需要在实体模型上约束？
4. 我只在板料上的两个截面设置了轴对称约束，是否就不用加别的约束了？
5. 我把位移-时间数组加到上模具上没有显示箭头，奇怪？
6. 进行求解时马上就弹出了solution is done，可是没有结果...


本来我的毕设题目是做一个大尺寸厚壁封头数据库系统的，老师说太简单了，要我在软件里添加有限元分析的功能，而资料只有上届学长写的论文（他用abaqus分析），所以很多材料的信息也不懂，所以当中的材料属性都是参考书籍然后暂时定下来的。请各位大大帮帮忙，谢谢。

顺便也把我的数据库文件打包（在附件里）及把命令流文件贴上，请大大看看是不是我的模型建立有问题，还是单元的选择不对...

刚开始学，问了很幼稚的题目非常抱歉...
fini
/cle
/filname,thesis

/prep7
*set,h,130            !毛胚厚度
*set,mr,1550        !毛胚半径
*set,lr,1153        !拉环半径
*set,cir,850        !冲头内径(半径)
*set,cor,1010        !冲头外径(半径)
*set,cl,2000        !冲头身长

!属性定义
et,1,164            
et,2,200,6

mp,ex,1,200e9
mp,nuxy,1,.27
mp,dens,1,7800e-9
tb,plaw,,,,1
tbdata,1,310e6
tbdata,2,763e6
tbdata,4,40.0
tbdata,5,5.0
tbdata,6,.75

mp,ex,2,200e9
mp,nuxy,2,.27
mp,dens,2,7800e-9
edmp,rigid,2,6,7

mp,ex,3,200e9
mp,nuxy,3,.27
mp,dens,3,7800e-9
edmp,rigid,3,7,7

!几何实体建模
!毛胚
k,1,0,0,0
k,2,mr,0,0
k,3,mr,h,0
k,4,0,h,0
a,1,2,3,4

vrotat,all,,,,,,1,4,90
vsel,all
mat,1
real,1
type,1
vmesh,all

!冲头
k,,0,h+1,0
k,,1010,h+cor+1,0
k,,1010,h+cor+1+cl,0
k,,0,h+1+160,0
k,,850,1141,0
k,,850,1141+cl,0
k,,0,1140,0
larc,7,8,13,1010
larc,10,11,13,850
lstr,8,9
lstr,11,12
lstr,7,10
lstr,9,12
al,10,11,12,13,14,15

type,2
mat,2
asel,s,AREA,,6,,,0
amesh,all
extopt,aclear,1
extopt,esize,20
type,1
vrotat,all,,,,,,7,10,90

!拉环
k,,lr,-101,0
k,,lr,-371,0
k,,lr+200,-371,0
k,,lr+200,-1,0
k,,lr+100,-1,0
k,,lr+100,-101,0
k,,0,-5,0
larc,18,22,23,100
lstr,18,19
lstr,19,20
lstr,20,21
lstr,21,22
al,25,26,27,28,29

type,2
mat,3
asel,s,AREA,,13,,,0
amesh,all
extopt,aclear,1
extopt,esize,20
type,1
vrotat,all,,,,,,1,24,-90

!轴对称约束
esel,s,mat,1

!创建PART号
allsel
edpart,create

!建立载荷数组
*dim,time,array,15,1
time(1,1)=0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9
time(11,1)=1,1.1,1.2,1.3,1.4

*dim,punch,array,15,1
punch(1,1)=0,-200,-300,-400,-500,-600,-700,-800,-900,-1000
punch(11,1)=-1100,-1200,-1300,-1400,-1500

sunnygu

我试着回答一下。局限于我的水平

1. 由于模型是轴对称的，我只建立了如上图所示的1/4模型是否可行？
可以。
2. 不知道我的分网方式对吗？听说板料要用自适应网格，可是好像只适合壳单元...
因为冲压属于大变形，会引起网络畸形，因此采用自适应网络，自适应网络也适应于实体单元。

3. 上、下模具采用刚体，在材料属性已经约束自由度了，那还需不需要在实体模型上约束？

不明白材料属性设置约束。实体模型还是需要加边界条件吧。MSC 软件就不需要定义。

4. 我只在板料上的两个截面设置了轴对称约束，是否就不用加别的约束了？

你的模型不需要其他约束了吧。

5. 我把位移-时间数组加到上模具上没有显示箭头，奇怪？

没检查你的网络，不能回答。

6. 进行求解时马上就弹出了solution is done，可是没有结果...

我不知道是不是你没加载结果。

在workbench里面我给出一个例子，比你的要复杂的多，直接是模板。建议你参考。如果是为了论文，不想深究的话，我的例子比你的要复杂的多。

你的这个过程几乎所有软件都有帮助里介绍，很简单。

你用的ls-dyna

bizarrEaster

非常感谢楼上的回答，我知道我做的其实只是皮毛，
我这里的图书馆很少ANSYS的书，
我参照的都大都是网上的例子，而冲压例子大都是薄壁，用的是shell163单元
而我要做的是厚壁，要观察板料的厚度变化的，我选用solid164单元是否有错？虽然我用ls-dyna已经省略了很多单元了，可是我还是不太懂单元该怎么选取...

为上下模具设置刚体材料属性时，我已限制了上模具的除y方向以外的自由度，而下模具限制了所有自由度。

至于加载我确实是加载了，我列表显示也可以看到，可是solve后很快就solution is done了，我进行求解时出现两个警告对话框，一个是a check of your model data produced 1 warning，另一个是a check of your load data produced 1 warning.

再请教各位大大谁有厚壁冲压的例子吗？我找到的都是薄壁的，还是其实没有区别？请原谅我的无知，谢谢...

以下是我的messag文件内容：

contracting memory to           1 beginning of keyword reader
expanding   memory to       21002  
expanding   memory to       46356  
expanding   memory to      127374  
expanding   memory to      157374  
expanding   memory to      217374  
expanding   memory to      238578  
expanding   memory to      259089  
expanding   memory to      303572  
expanding   memory to      342727  

Memory required to process keyword     :       342727
contracting memory to           1 end of keyword reader

LS-DYNA will perform a structural only analysis

LS-DYNA user input                                                      
                         ls-dyna ls971s.50638 R4.2    date 02/04/2009
expanding   memory to        2625 material, node and element data
expanding   memory to        2705 material properties
expanding   memory to       27865 element data
expanding   memory to       34155 gradient of damage data
expanding   memory to       64270 user control for contact interfaces
expanding   memory to       64274 for boundary cond. vectorization
expanding   memory to       66251 arbitrary element numbering input
expanding   memory to      122005 process element material data
expanding   memory to      124011 load curve data
contracting memory to      122723 regularize time history data
contracting memory to      122723 rigidwall definitions
expanding   memory to      125898 velocity b.c
expanding   memory to      129168 nodal constraint sets
expanding   memory to      152356 contact interfaces
expanding   memory to      159630 implicit contact 1
contracting memory to      153356 contact segment data 1
expanding   memory to      178775 contact segment data 3
expanding   memory to      240216 contact definition
expanding   memory to      249960 rigid body
expanding   memory to      266293 contact friction
expanding   memory to      266338 joints
contracting memory to      266338 shell solid interfaces
contracting memory to      266338 accelerometers
expanding   memory to      266438 rigid bodies for disp. termination
expanding   memory to      266558 cyclic symmetry
contracting memory to      266558 transformation to move nodes
expanding   memory to      266561 thermal expansion 0
expanding   memory to      266564 flexible body 2
contracting memory to      266564 implicit joints

E r r o r   t e r m i n a t i o n


s t o r a g e   a l l o c a t i o n   

Memory required to complete solution   :       342727
Additional dynamically allocated memory:        88821
                                   Total:       431548

T i m i n g   i n f o r m a t i o n
                        CPU(seconds)   %CPU  Clock(seconds) %Clock
  ----------------------------------------------------------------
  Initialization ....... 0.0000E+00    0.00     0.0000E+00    0.00
  Element processing ... 0.0000E+00    0.00     0.0000E+00    0.00
  Binary databases ..... 0.0000E+00    0.00     0.0000E+00    0.00
  ASCII database ....... 0.0000E+00    0.00     0.0000E+00    0.00
  Contact algorithm .... 0.0000E+00    0.00     0.0000E+00    0.00
  Contact entities ..... 0.0000E+00    0.00     0.0000E+00    0.00
  Rigid bodies ......... 0.0000E+00    0.00     0.0000E+00    0.00
  Implicit Nonlinear ... 0.0000E+00    0.00     0.0000E+00    0.00
  Implicit Lin. Alg. ... 0.0000E+00    0.00     0.0000E+00    0.00
  ----------------------------------------------------------------
  T o t a l s            0.0000E+00    0.00     0.0000E+00    0.00

Problem time       =    0.0000E+00
Problem cycle      =         0
Total CPU time     =         0 seconds (   0 hours  0 minutes  0 seconds)
CPU time per zone cycle  =         0 nanoseconds
Clock time per zone cycle=         0 nanoseconds

Number of CPU's    1
NLQ used/max      96/   96
Start time   04/15/2010 07:20:30  
End time     04/15/2010 07:20:31  
Elapsed time       1 seconds(  0 hours  0 min.  1 sec.) for       0 cycles

E r r o r   t e r m i n a t i o n

一笑泯恩仇

楼主你刚学一个星期，就写出这么长的程序，而且做出相对比较复杂的问题，膜拜。想当初，我花了两个月的时间才做出一个焊接的模拟过程，哈哈