板材冲压过程仿真

z_mj_2005

下面的APDL命令流，生产K文件，然后在LS-PREPOST打开，结果发现上模没动。APDL命令流如下：
fini
/clear
!--------------------------板材成形功能说明-------------------------!
!--1 模拟板材成形过程及应力、应变状态；
!--2
!-------------------------------------------------------------------!
!-------------------------- (1)设置工程属性-------------------------!
/FILNAME,MPF-6000-2500-16S-X
/UNITS,SI
/PREP7
/VIEW,,1,2,3
/PLOPTS,INFO,1
!------------------------- (2)建立实体模型--------------------------!
!* 建立钢板几何模型
K,1001,-1120,-0.5925E4,-500
K,1002,625,-0.5925E4,-500
K,1003,625,-0.5925E4,500
K,1004,-1120,-0.5925E4,500
L,1001,1002
L,1002,1003
L,1003,1004
L,1004,1001
AL,1,2,3,4
!* 导入模具几何模型
~SATIN,'MPF-6000-2500-16S','SAT','.\MPF\GModel\',SOLIDS,0    ! 以.sat格式导
                                                             ! 入几何模型
/NOPR   
/GO


!----------------- (3)定义单元类型、实常数与材料属性----------------!
! 设置单元类型
ET,1,SOLID164
ET,2,shell163
ET,3,SOLID164
KEYOPT,2,1,9          ! 9:全积分Belytschko-Tsay膜单元算法；10:B-W-C单元算法


! 定义材料模型
! 单元1
!R,1                  ! Solid164单元没有实常数
MP,DENS,1,2.28E-3
MP,EX,1,30E6
MP,NUXY,1,0.29
EDMP,rigid,1,7
! 单元2
R,2,5/6,3,0.016,0.016,0.016,0.016
!R,2,,,0.016          ! 定义钢板厚度为16mm  
MP,DENS,2,2.28E-3
MP,EX,2,30E6
MP,NUXY,2,0.29
TB,BISO,1,,,,
TBDATA,1,3.0e4
TBDATA,2,1.0e5
EDSHELL,20,-1,1,2,1,1
! 单元3
!R,1                  ! Solid164单元没有实常数
MP,DENS,3,2.28E-3
MP,EX,3,30E6
MP,NUXY,3,0.29
EDMP,rigid,3,6


! (4)划分网格，建立有限元模型
vsel,s,,,1,1,1        ! 选择体编号为1,此命令可以省去，因为TYPE已经选中单元1
!TYPE,1
MAT,1
ESIZE,125
MSHAPE,0,3D           ! 0：六面体；1:四面体
!VMESH,all     
VMESH,1
!vsel,s,,,1,3,1       ! 选择所用体


!ESIZE,125,0,
!MSHAPE,0,2D
!MSHKEY,1
!AMESH,6,6,1


! 划分板材网格   
Asel,all
Asel,s,,,1,1,1         ! 选择体号1,钢板
TYPE,2                 !
MAT,2                  !
REAL,2                 ! 激活常数类型2；否则显示"Nodal thikness is equal to zero"
ESIZE,25
MSHAPE,0,2D            ! 0：四边形；1:三角形
MSHKEY,1
AMESH,1


vsel,s,,,2,2,1     ! 选择上模
TYPE,3
MAT,3
ESIZE,125
MSHAPE,0,3D        ! 0：六面体；1:四面体
!SMRT,0.125
VMESH,all   
vsel,s,,,1,3,1     ! 选择所用体
allsel,all         ! 选择所对象：点、线、面、体、节点、单元


!--------------------- (6)定义接触-------------------！
ESEL,S,TYPE,,1              ! 激活单元类型为1的单元
CM,eDie,ELEM                ! 将已经激活的单元编成组，单元组名称为eDie
NSLE,S                      ! 激活单元组eDie中的所有节点
CM,nDie,NODE                ! 将已经激活的节点编成组，节点组名称为nDie
ALLSEL,ALL
EDPART,CREATE,1,eDie        ! 创建组元，编号为1
ESEL,S,TYPE,,2              ! 激活单元类型为2的单元
CM,eGangban,ELEM            ! 将已经激活的单元编成组件，单元组名称为eGangban
NSLE,S                      ! 激活单元组件eGangban中的所有节点
CM,nGangban,NODE            ! 将已经激活的节点编成组件，节点组名称为nGangban
ALLSEL,ALL
!EDPART,CREATE,2,nGangban    ! 创建组元，编号为2
EDPART,CREATE
!EDCGEN,NTS,2,1,0.15,0.125,0,0,0,,,,,0,10000000
EDCGEN,ASSC, , ,0,0,0,0,0, , , , ,0,10000000,0,0 ! Specifies contact
                                                 ! parameters for an explicit     
                             ! dynamics analysis.
!EDCONTACT,0.1,0,2,0,1,1,1,4,0
!EDCGEN,NTS,2,3,0.15,0.125,0,0,0,,,,,0,10000000
!EDCONTACT,0.1,0,2,0,1,1,1,4,0
!--------------- (5)定义初始条件与边界条件------------------！
!ESEL,S,TYPE,,3    ! 激活单元类型为3的单元
ESEL,S,MAT,,3
CM,ePunch,ELEM     ! 将已经激活的单元编成组件，单元组件名称为ePunch
NSLE,S             ! 激活单元组ePunch中的所有节点
CM,nPunch,NODE     ! 将已经激活的节点编成组件，节点组件名称为nPunch
! 定义位移载荷
*DIM,CTIME,ARRAY,10,1,1,,,
*DIM,DISP,ARRAY,10,1,1,,,
*SET,CTIME(2,1,1),0.001
*SET,CTIME(3,1,1),0.002
*SET,CTIME(4,1,1),0.003
*SET,CTIME(5,1,1),0.004
*SET,CTIME(6,1,1),0.005
*SET,CTIME(7,1,1),0.006
*SET,CTIME(8,1,1),0.007
*SET,CTIME(9,1,1),0.008
*SET,CTIME(10,1,1),0.010     ! 时间如果超过程序终止运行时间，会没反应。
*SET,DISP(2,1,1),-0.010
*SET,DISP(3,1,1),-0.020
*SET,DISP(4,1,1),-0.030
*SET,DISP(5,1,1),-0.040
*SET,DISP(6,1,1),-0.050
*SET,DISP(7,1,1),-0.060
*SET,DISP(8,1,1),-0.070
*SET,CTIME(9,1,1),-0.080  
*SET,CTIME(10,1,1),-0.08341
ALLSEL,ALL                   ! 将全部的节点、单元都激活
EDPART,CREATE,3,ePunch       ! 创建组元，编号为3
EDLOAD,ADD,RBUY,0,3,CTIME,DISP,0,,,,,
!EDVEL,VGEN,nPunch,0,100,0,0,0,0  ! define initial vecolicity


! 约束板材一端面
DL,4,,ALL


! 约束凹模底座
ALLSEL,ALL          ! 将全部的节点、单元都激活
FLST,2,1,5,ORDE,1   
FITEM,2,6   
/GO
DA,P51X,ALL,0  
DTRAN               ! 将实体约束转换为节点约束
ALLSEL,ALL          ! 将全部的节点、单元都激活


!-------------------- (7)定义求解参数-------------------！
TIME,0.015
EDRST,100
EDHTIME,100
EDOUT,NCFORC
EDOUT,RCFORC
EDOUT,MATSUM
ALLSEL,ALL
SAVE
FINISH
!--------------- (8)进入求解其，输出K文件---------------！
/SOLU
EDWRITE,LSDYNA,'MPF-6000-2500-16S-X','K',''  ! 输出K文件

z_mj_2005

麻烦那位高手指点一下！