缓冲器冲击试验模型出错

eastwood

我作的一缓冲器冲击试验,老有问题,请大家给我诊断一下:
下面是我的LOG文件:
/title, changchy new example
/filnam,testnew01
  
/prep7
/view,,1,2,3
/plopts,info,1
/pnum,mat,1
/num,1
ET,1,SOLID164   
!*   
KEYOPT,1,1,2
KEYOPT,1,5,0
!*   
!*   
  
*DIM,TABLE5,ARRAY,5,1,1, , ,
/INPUT,ansuitmp  
!*   
*CREATE,ansuitmp
*VREAD,TABLE5,'codata','dat',' ',1, , , , ,1,   
(5E10.3)
*END
/INPUT,ansuitmp  
!*
   
DENSCO=TABLE5(1,1,1)
GMAX=TABLE5(2,1,1)
G0=TABLE5(3,1,1)
K=TABLE5(4,1,1)
N=TABLE5(5,1,1)
  
DENSCO=TABLE5(1,1,1)
GMAX=TABLE5(2,1,1)
G0=TABLE5(3,1,1)
K=TABLE5(4,1,1)
N=TABLE5(5,1,1)
  
MP,DENS,1,DENSCO
TB,EVIS,1,,,,   
TBDAT,46,GMAX  
TBDAT,47,G0  
TBDAT,48,K
TBDAT,61,N
  
MP,DENS,2,7.85e-6
MP,EX,2,210  
MP,NUXY,2,0.3  
  
!EDMP,RIGI,3,4,0  
MP,DENS,3,7.85e-6
MP,EX,3,210  
MP,NUXY,3,0.3  
  
MP,DENS,4,7.85e-6
MP,EX,4,210  
MP,NUXY,4,0.3   
   
  
CYL4, , ,53, ,60, ,362  
CYL4, , ,52, , , ,362
BLC5, , ,20,20,362  
   
FLST,2,2,6,ORDE,2   
FITEM,2,2   
FITEM,2,-3   
VOVLAP,P51X
  
CM,_Y,VOLU   
VSEL, , , ,       1  
CM,_Y1,VOLU  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
VATT,       2, ,   1,       0   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
FLST,5,2,6,ORDE,2   
FITEM,5,3   
FITEM,5,-4   
CM,_Y,VOLU   
VSEL, , , ,P51X  
CM,_Y1,VOLU  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
VATT,       1, ,   1,       0   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
ESIZE,10,0,  
CM,_Y,VOLU   
VSEL, , , ,       3  
CM,_Y1,VOLU  
CHKMSH,'VOLU'   
CMSEL,S,_Y   
!*   
VSWEEP,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*  
ESIZE,20,0,   
CM,_Y,VOLU   
VSEL, , , ,       4  
CM,_Y1,VOLU  
CHKMSH,'VOLU'   
CMSEL,S,_Y   
!*   
VSWEEP,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*   
  
CM,_Y,VOLU   
VSEL, , , ,       1  
CM,_Y1,VOLU  
CHKMSH,'VOLU'   
CMSEL,S,_Y   
!*   
VSWEEP,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*   
!********************************************************************************
!*   
CYL4, , ,52, , , ,100   
FLST,3,1,6,ORDE,1   
FITEM,3,2   
VGEN, ,P51X, , , , ,363, , ,1   
BLC5, , ,20,20,100   
FLST,3,1,6,ORDE,1   
FITEM,3,5   
VGEN, ,P51X, , , , ,363, , ,1   
/REPLOT  
FLST,2,2,6,ORDE,2   
FITEM,2,2   
FITEM,2,5   
VOVLAP,P51X  
/DIST, 1 ,1.371742,1
/REP,FAST   
vlist, all   
FLST,5,2,6,ORDE,2   
FITEM,5,5   
FITEM,5,-6   
CM,_Y,VOLU   
VSEL, , , ,P51X  
CM,_Y1,VOLU  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
VATT,       3, ,   1,       0   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
ESIZE,10,0,  
CM,_Y,VOLU   
VSEL, , , ,       5  
CM,_Y1,VOLU  
CHKMSH,'VOLU'   
CMSEL,S,_Y   
!*   
VSWEEP,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*   
ESIZE,20,0,  
CM,_Y,VOLU   
VSEL, , , ,       6  
CM,_Y1,VOLU  
CHKMSH,'VOLU'   
CMSEL,S,_Y   
!*   
VSWEEP,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*   
!********************************************************************
BLC5, , ,200,200,100
FLST,3,1,6,ORDE,1   
FITEM,3,2   
VGEN, ,P51X, , , , ,470, , ,1   
/VIEW, 1 ,,1
/ANG, 1  
/REP,FAST   
/AUTO, 1
/REP
/VIEW, 1 ,,-1   
/ANG, 1  
/REP,FAST   
CM,_Y,VOLU   
VSEL, , , ,       2  
CM,_Y1,VOLU  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
VATT,       4, ,   1,       0   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
ESIZE,50,0,  
CM,_Y,VOLU   
VSEL, , , ,       2  
CM,_Y1,VOLU  
CHKMSH,'VOLU'   
CMSEL,S,_Y   
!*   
VSWEEP,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*   
!*******************************************************************************
ESEL,S,MAT,,4   
/REPLOT  
NSLE,S   
/REPLOT  
NPLOT   
CM,TOP,NODE  
ALLSEL,ALL   
EPLOT   
EDPART,CREATE   
!*   
!EDCGEN,ASSC, , ,0,0,0,0,0,0, , , ,0,10000000,0,0
!edcgen,assc  
!EDCONTACT,0.1,0,2,0,1,1,1,4,0   
!EDSP,OFF,' ',' ',1,  
!*   
ASEL,S,LOC,Z,0   
/REPLOT  
APLOT   
APLOT   
/USER,  1   
/VIEW,  1, -0.103713418686E-01, -0.979168504745    ,  0.202784305564
/ANG,   1,   2.94834066170   
/REPLO   
FLST,2,3,5,ORDE,3   
FITEM,2,1   
FITEM,2,11   
FITEM,2,17   
!*   
/GO  
DA,P51X,UZ,  
ALLSEL,ALL   
EPLOT   
/AUTO, 1
/REP
!***********************************************************************
*DIM,TIME1,ARRAY,3,1,1, , ,  
!*   
*SET,TIME1(1,1,1) , 0.0  
*SET,TIME1(2,1,1) , 1  
*SET,TIME1(3,1,1) , 2
*DIM,AY1,ARRAY,3,1,1, , ,   
!*   
*SET,AY1(1,1,1) , -9800
*SET,AY1(2,1,1) , -9800
*SET,AY1(3,1,1) , -9800
  
FINISH   
/SOL
!*   
EDVE,VELO,TOP,0,0,-1000,0,0,0, , , , , ,
EDLOAD,ADD,AZ,0,TOP,TIME1,AY1, 0, , , , ,   
EDENERGY,1,1,1,1
EDBVIS,1.5,0.06,
EDHGLS,0.1,  
TIME,0.03,   
EDRST,100,   
EDHTIME,1000,   
EDDUMP,1,   
!*   
!*   
!*   
EDOUT,SPCFORC   
!*********************************************************************************
FINISH   
/PREP7   
!*   
ASEL,S,LOC,Z,463
APLOT   
NSLA,S,1
CM,BOTTON,NODE   
ALLSEL,ALL   
ASEL,S,LOC,Z,470
/REPLOT  
NSLA,S,1
/REPLOT  
CM,TOP1,NODE
ALLSEL,ALL   
EPLOT   
!*   
EDCGEN,ASTS,TOP1,BOTTON,0,0,0,0,0, , , , ,0,10000000,0,0
FINISH   
/SOL
/STATUS,SOLU
!SOLVE   
/PREP7   
ESEL,S,MAT,,1,3,1   
VPLOT   
EPLOT   
NSLE,S   
CM,BODY,NODE
!*   
EDCGEN,ASSC,BODY, ,0.1,0,0,0,0,0.1, , , ,0,10000000,0,0
ALLSEL,ALL   
EPLOT   
EDMP,HGLS,1,5,0.1,1.5,0.06, , ,  
EDMP,HGLS,2,5,0.1,1.5,0.06, , ,  
EDMP,HGLS,3,5,0.1,1.5,0.06, , ,  
EDMP,HGLS,4,5,0.1,1.5,0.06, , ,  
  
!****************************************************************************
FINISH   
/SOL
/STATUS,SOLU
SOLVE   
wpro,,90.000000,
/GRAPHICS,POWER  
!*   
/TYPE,1,5   
/CPLANE,1   
/SHADE,1,1   
/REPLOT  
!*   
!*   
/DSCALE,1,1.0   
/REPLOT  
!*   
FINISH   
/POST1   
SET,LAST
!*   
PLNS,S,EQV   
ANDATA,0.5, ,0,0,0,1,0,1
!*

mehere

很乱，你应该删除多余的命令，将剩余的参数化，要不别人没法看。

htbbzzg

!*  由于缺少  codata.dat  数据文件    (5E10.3)，无法试算，只能对 log 文件作一些简单分析，感觉有问题之处写在下面
!*  1  后面定义了材料 2-4，从密度数据来看，应该是使用了 kg-mm-s 单位制，材料为钢材，但其中弹性模量似乎不对，在 kg-mm-s 单位制中弹性模量的单位
!*     应该是 kPa，因此钢的弹性模量应该是 210 (GPa) -> 210000 MPa -〉 210E6 kPa， 你给为 210，其单位应该是 GPa，数值差了1.0e6。
!*     由于弹性模量数值给得太小了，计算时会造成变形过大，使结构扭曲，出现负体积错误等。
!*  2  定义接触的语句
!*     EDCGEN,ASSC,BODY, ,0.1,0,0,0,0,0.1, , , ,0,10000000,0,0  
!*    似乎不对，使用 ASSC 方式时，不要指定任何接触实体，默认为全部实体，故该语句应改为：
!*    EDCGEN,ASSC, , ,0.1,0,0,0,0,0.1, , , ,0,10000000,0,0  
!*  关键问题应该是第一条。
  
!  缓冲器冲击试验         
!  color=black]  
!  我作的一缓冲器冲击试验,老有问题,请大家给我诊断一下:  
!  下面是我的LOG文件:  
/title, changchy new example  
/filnam,testnew01  
  
/prep7  
/view,,1,2,3  
/plopts,info,1  
/pnum,mat,1  
/num,1  
ET,1,SOLID164  
!*  
KEYOPT,1,1,2  
KEYOPT,1,5,0  
!*  
!*  
  
*DIM,TABLE5,ARRAY,5,1,1, , ,  
/INPUT,ansuitmp  
!*  
*CREATE,ansuitmp  
*VREAD,TABLE5,'codata','dat',' ',1, , , , ,1,  
(5E10.3)  
*END  
/INPUT,ansuitmp  
!*  
   
DENSCO=TABLE5(1,1,1)  
GMAX=TABLE5(2,1,1)  
G0=TABLE5(3,1,1)  
K=TABLE5(4,1,1)  
N=TABLE5(5,1,1)  
  
DENSCO=TABLE5(1,1,1)  
GMAX=TABLE5(2,1,1)  
G0=TABLE5(3,1,1)  
K=TABLE5(4,1,1)  
N=TABLE5(5,1,1)  
  
MP,DENS,1,DENSCO  
TB,EVIS,1,,,,  
TBDAT,46,GMAX  
TBDAT,47,G0  
TBDAT,48,K  
TBDAT,61,N  
  
!*  以下定义了材料 2-4，从密度数据来看，应该是使用了 kg-mm-s 单位制，材料为钢材，但其中弹性模量似乎不对，在 kg-mm-s 单位制中钢的弹性模量一般为 210000 MPa
MP,DENS,2,7.85e-6  
MP,EX,2,210                !!   应改为   210e6,  下同
MP,NUXY,2,0.3  
  
!EDMP,RIGI,3,4,0  
MP,DENS,3,7.85e-6  
MP,EX,3,210  
MP,NUXY,3,0.3  
  
MP,DENS,4,7.85e-6  
MP,EX,4,210  
MP,NUXY,4,0.3  
   
  
CYL4, , ,53, ,60, ,362  
CYL4, , ,52, , , ,362  
BLC5, , ,20,20,362  
  
FLST,2,2,6,ORDE,2  
FITEM,2,2  
FITEM,2,-3  
VOVLAP,P51X  
  
CM,_Y,VOLU  
VSEL, , , , 1  
CM,_Y1,VOLU  
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1  
VATT, 2, , 1, 0  
CMSEL,S,_Y  
CMDELE,_Y  
CMDELE,_Y1  
!*  
FLST,5,2,6,ORDE,2  
FITEM,5,3  
FITEM,5,-4  
CM,_Y,VOLU  
VSEL, , , ,P51X  
CM,_Y1,VOLU  
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1  
VATT, 1, , 1, 0  
CMSEL,S,_Y  
CMDELE,_Y  
CMDELE,_Y1  
!*  
ESIZE,10,0,  
CM,_Y,VOLU  
VSEL, , , , 3  
CM,_Y1,VOLU  
CHKMSH,'VOLU'  
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y  
CMDELE,_Y1  
CMDELE,_Y2  
!*  
ESIZE,20,0,  
CM,_Y,VOLU  
VSEL, , , , 4  
CM,_Y1,VOLU  
CHKMSH,'VOLU'  
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y  
CMDELE,_Y1  
CMDELE,_Y2  
!*  
  
CM,_Y,VOLU  
VSEL, , , , 1  
CM,_Y1,VOLU  
CHKMSH,'VOLU'  
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y  
CMDELE,_Y1  
CMDELE,_Y2  
!*  
!********************************************************************************  
!*  
CYL4, , ,52, , , ,100  
FLST,3,1,6,ORDE,1  
FITEM,3,2  
VGEN, ,P51X, , , , ,363, , ,1  
BLC5, , ,20,20,100  
FLST,3,1,6,ORDE,1  
FITEM,3,5  
VGEN, ,P51X, , , , ,363, , ,1  
/REPLOT  
FLST,2,2,6,ORDE,2  
FITEM,2,2  
FITEM,2,5  
VOVLAP,P51X  
/DIST, 1 ,1.371742,1  
/REP,FAST  
vlist, all  
FLST,5,2,6,ORDE,2  
FITEM,5,5  
FITEM,5,-6  
CM,_Y,VOLU  
VSEL, , , ,P51X  
CM,_Y1,VOLU  
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1  
VATT, 3, , 1, 0  
CMSEL,S,_Y  
CMDELE,_Y  
CMDELE,_Y1  
!*  
ESIZE,10,0,  
CM,_Y,VOLU  
VSEL, , , , 5  
CM,_Y1,VOLU  
CHKMSH,'VOLU'  
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y  
CMDELE,_Y1  
CMDELE,_Y2  
!*  
ESIZE,20,0,  
CM,_Y,VOLU  
VSEL, , , , 6  
CM,_Y1,VOLU  
CHKMSH,'VOLU'  
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y  
CMDELE,_Y1  
CMDELE,_Y2  
!*  
!********************************************************************  
BLC5, , ,200,200,100  
FLST,3,1,6,ORDE,1  
FITEM,3,2  
VGEN, ,P51X, , , , ,470, , ,1  
/VIEW, 1 ,,1  
/ANG, 1  
/REP,FAST  
/AUTO, 1  
/REP  
/VIEW, 1 ,,-1  
/ANG, 1  
/REP,FAST  
CM,_Y,VOLU  
VSEL, , , , 2  
CM,_Y1,VOLU  
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1  
VATT, 4, , 1, 0  
CMSEL,S,_Y  
CMDELE,_Y  
CMDELE,_Y1  
!*  
ESIZE,50,0,  
CM,_Y,VOLU  
VSEL, , , , 2  
CM,_Y1,VOLU  
CHKMSH,'VOLU'  
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y  
CMDELE,_Y1  
CMDELE,_Y2  
!*  
!*******************************************************************************  
ESEL,S,MAT,,4  
/REPLOT  
NSLE,S  
/REPLOT  
NPLOT  
CM,TOP,NODE  
ALLSEL,ALL  
EPLOT  
EDPART,CREATE  
!*  
!EDCGEN,ASSC, , ,0,0,0,0,0,0, , , ,0,10000000,0,0  
!edcgen,assc  
!EDCONTACT,0.1,0,2,0,1,1,1,4,0  
!EDSP,OFF,' ',' ',1,  
!*  
ASEL,S,LOC,Z,0  
/REPLOT  
APLOT  
APLOT  
/USER, 1  
/VIEW, 1, -0.103713418686E-01, -0.979168504745 , 0.202784305564  
/ANG, 1, 2.94834066170  
/REPLO  
FLST,2,3,5,ORDE,3  
FITEM,2,1  
FITEM,2,11  
FITEM,2,17  
!*  
/GO  
DA,P51X,UZ,  
ALLSEL,ALL  
EPLOT  
/AUTO, 1  
/REP  
!***********************************************************************  
*DIM,TIME1,ARRAY,3,1,1, , ,  
!*  
*SET,TIME1(1,1,1) , 0.0  
*SET,TIME1(2,1,1) , 1  
*SET,TIME1(3,1,1) , 2  
*DIM,AY1,ARRAY,3,1,1, , ,  
!*  
*SET,AY1(1,1,1) , -9800  
*SET,AY1(2,1,1) , -9800  
*SET,AY1(3,1,1) , -9800  
  
FINISH  
/SOL  
!*  
EDVE,VELO,TOP,0,0,-1000,0,0,0, , , , , ,  
EDLOAD,ADD,AZ,0,TOP,TIME1,AY1, 0, , , , ,  
EDENERGY,1,1,1,1  
EDBVIS,1.5,0.06,  
EDHGLS,0.1,  
TIME,0.03,  
EDRST,100,  
EDHTIME,1000,  
EDDUMP,1,  
!*  
!*  
!*  
EDOUT,SPCFORC  
!*********************************************************************************  
FINISH  
/PREP7  
!*  
ASEL,S,LOC,Z,463  
APLOT  
NSLA,S,1  
CM,BOTTON,NODE  
ALLSEL,ALL  
ASEL,S,LOC,Z,470  
/REPLOT  
NSLA,S,1  
/REPLOT  
CM,TOP1,NODE  
ALLSEL,ALL  
EPLOT  
!*  
EDCGEN,ASTS,TOP1,BOTTON,0,0,0,0,0, , , , ,0,10000000,0,0  
FINISH  
/SOL  
/STATUS,SOLU  
!SOLVE  
/PREP7  
ESEL,S,MAT,,1,3,1  
VPLOT  
EPLOT  
NSLE,S  
CM,BODY,NODE  
!*  
EDCGEN,ASSC,BODY, ,0.1,0,0,0,0,0.1, , , ,0,10000000,0,0  
ALLSEL,ALL  
EPLOT  
EDMP,HGLS,1,5,0.1,1.5,0.06, , ,  
EDMP,HGLS,2,5,0.1,1.5,0.06, , ,  
EDMP,HGLS,3,5,0.1,1.5,0.06, , ,  
EDMP,HGLS,4,5,0.1,1.5,0.06, , ,  
  
!****************************************************************************  
FINISH  
/SOL  
/STATUS,SOLU  
!  SOLVE  
wpro,,90.000000,  
/GRAPHICS,POWER  
!*  
/TYPE,1,5  
/CPLANE,1  
/SHADE,1,1  
/REPLOT  
!*  
!*  
/DSCALE,1,1.0  
/REPLOT  
!*  
FINISH  
/POST1  
SET,LAST  
!*  
PLNS,S,EQV  
ANDATA,0.5, ,0,0,0,1,0,1  
!*