【求助】TO:wuhanwangtao: 怎样模拟正挤压的剪切？

ASKA

麻烦王老师给看一下，多谢了
  
finish
/cle
/NOPR   
/PMETH,OFF,0
KEYW,PR_SET,1  
KEYW,PR_STRUC,1  
KEYW,LSDYNA,1
  
!******************************
/UNITS,SI
!******************************
  
/plopts,info,0
!/view,,1,1,1
/VIEW, 1 ,,-1  
  
dw=0.001   !1 millimetre毫米 =0.03937 inch 英寸
  
C_PI=3.1415926  !PI的值
RHO=2800   !铝合金的密度
VTOT=0     !活塞的体积初始化
WT=0       !活塞的质量初始化
SURF=1E7   !压力
C_R=85*dw/2  !C_R为底面半径
C_F=SURF*C_PI*C_R*C_R       !C_F压力
           
ACAZ=0          !Z方向的加速度
  
/prep7  
ET,1,SHELL163   !铝合金和钢的二维有限元模型的et     
  
ET,2,solid164   !铝合金的solid164                    
ET,3,solid164   !钢的solid164    刚体   
  
MP,DENS,1,2800     !aluminum
MP,EX,1,72e9   
MP,NUXY,1,0.33
TB,BKIN,1,,,,   
TBDAT,1, 480e6                     
TBDAT,2, 27e9  
  
MP,DENS,2,7800
MP,EX,2,200e9
MP,NUXY,2,0.3
EDMP,RIGI,2,7,7  
  
wpro,,90.000000,               
  
!建模活塞
RECTNG,0,85/2000,20/1000+dw,31/1000+dw         
RECTNG,83/2000,85/2000,0/1000+dw,20/1000+dw
ADEL,ALL
LANG,1,5,-(90-5.5)
LFILLT,1,10,5/1000
LDELE,7,9,1
LDELE,2,6,4
L,3,6
  
LANG,3,10,90
LANG,2,11,90
LDIV,11
L,3,14
  
AL,4,1,7,6
AL,7,11,13,3
AL,12,9,2,13
AL,10,5,8,9
CM,HUOSAI,AREA
!分格活塞二维模型
  
TYPE, 1
mat,1
MSHAPE, 0, 2D  
MSHKEY, 1
ESIZE, ,5                       
AMESH, HUOSAI   
!建模分格活塞二维模型finish
   
K,1000, 0  
K,1001, 0, 0,5                     
TYPE, 2
mat,1
ESIZE, ,20                        
  
VROTAT, ALL, , , , , , 1000, 1001, 360,4  
CMSEL, S,HUOSAI
ACLEAR, ALL
ALLSEL
EPLOT  
  
!/eof
!1
  
!***********************
*DO,I,101,8100                    
  *GET,VCAI,ELEM,I,VOLU      
    VTOT=VTOT+VCAI      !总的体积
*enddo                       
   
WT=RHO*VTOT                 !质量WT
!**********************
ACAZ=C_F/WT          !Z方向的加速度
!/eof
!2
  
!建模制动面
wpro,,-90.000000,
k,1011,85.1*dw/2,0,0*dw
k,1012,108*dw/2,0,0*dw
  
k,1013,80*dw/2,0,-25*dw
k,1014,82*dw/2,0,-25*dw
k,1015,108*dw/2,0,-25*dw
  
k,1016,80*dw/2,0,-47*dw
k,1017,82*dw/2,0,-47*dw
k,1018,108*dw/2,0,-47*dw
  
/pnum,kp,1
/pnum,line,1
/pnum,area,1
gplot
  
l,1011,1012
l,1012,1015
l,1015,1014
l,1014,1011
  
l,1013,1014
l,1014,1017
l,1017,1016
l,1016,1013
  
l,1017,1018
l,1018,1015
  
!/eof
!3
  
al,82,83,84,85  !area65
al,84,91,90,87  !area66
al,86,87,88,89  !area67
  
asel,s,,,65   
asel,a,,,66
asel,a,,,67
  
CM,ZHIDONG,AREA
  
!/eof
!4
  
TYPE, 1
mat,2
ESYS,       0   
MSHAPE, 0, 2D  
MSHKEY, 1
ESIZE, ,10                       
AMESH, ZHIDONG
  
TYPE, 3
mat,2
ESIZE, ,20                        
  
VROTAT, ALL, , , , , , 1000, 1001, 360,4  
CMSEL, S,ZHIDONG
ACLEAR, ALL
ALLSEL
  
/type,1,8  
  
EPLOT
  
/eof
  
ESEL,S,MAT,,1   
NSLE,S   
CM,BOX,NODE  
  
ALLSEL,ALL   
ESEL,S,MAT,,2   
NSLE,S   
CM,TABLE,NODE   
  
ALLSEL,ALL
!EDPART,CREATE   
  
!生成接触ASSC         
EDCGEN,ASSC
  
!EDCONTACT,1.0,,,,,,2,,0
  
edclist,all  
  
finish
  
/solu
!修改了速度将公制改为英制  速度方向从y方向转为z方向   
EDVE,VELO,BOX,0,0,-200,0,0,0, , , , , ,  
  
*DIM,time,ARRAY,2,1,1, , ,   
*DIM,accc,ARRAY,2,1,1, , ,   
*SET,TIME(1,1,1) , 0
*SET,TIME(2,1,1) , 3.0e-4
*SET,ACCC(1,1,1) , ACAZ
*SET,ACCC(2,1,1) , ACAZ
  
   
EDLOAD,ADD,ACLZ,0,BOX,TIME,ACCC, 0, , , , ,  
  
TIME,3.0e-4,                !终止时间  
EDRST,50,         
EDHTIME,50,     
EDDUMP,1,   
EDENERGY,1,1,1,1
  
edhist,box
  
edout,all            
  
nsel,all
allsel,all
/replot
finish
  
!/eof
!         
/SOLU
/STAT,SOLU   
SOLVE   
FINISH   
  
/EXIT,ALL  

wuhanwangtao

我机子这周要算个项目
  
我周五帮你算算
看看问题在哪

ASKA

在此先谢谢了
  
敬候佳音！！

ASKA

我是算到第27步出错的
提示我CHECK MESSAGE
  
expanding   memory to     3835617
expanding   memory to     3835637
contracting memory to     3835637
contracting memory to     3835637
expanding   memory to     3835717
expanding   memory to     3946375
expanding   memory to     3946377
expanding   memory to     3950459
  
  storage needed for solution=   3950459
       0 t 0.0000E+00 dt 0.00E+00 flush i/o buffers
       0 t 0.0000E+00 dt 0.00E+00 flush i/o buffers
expanding   memory to     4280231
  
  The LS-DYNA time step size should not exceed     0.247E-07
  to avoid contact instabilities.  If the step size is
  bigger then scale the penalty of the offending surface.
  
       0 t 0.0000E+00 dt 0.00E+00 flush i/o buffers
       1 t 0.0000E+00 dt 3.14E-08 flush i/o buffers
       1 t 0.0000E+00 dt 3.14E-08 write ANSYS result file
     192 t 5.9890E-06 dt 3.14E-08 write ANSYS result file
     383 t 1.1982E-05 dt 3.14E-08 write ANSYS result file
     
是有一段错误提示，可是看不明白？
恳请王老师解释一下

ASKA

不对，拷贝错了
错误提示是下面一段：
  3994 t 1.2600E-04 dt 3.15E-08 write ANSYS result file
    4184 t 1.3198E-04 dt 3.13E-08 write ANSYS result file
    4376 t 1.3799E-04 dt 3.01E-08 write ANSYS result file
    4603 t 1.4398E-04 dt 2.45E-08 write ANSYS result file
    4919 t 1.4999E-04 dt 1.30E-08 write ANSYS result file
    5000 t 1.5098E-04 dt 1.06E-08 flush i/o buffers
  
  negative volume in brick element #    4100 cycle=    7743
    7743 t 1.5422E-04 dt 1.00E+06 write ANSYS result file
  
  E r r o r   t e r m i n a t i o n
  
  Storage required for explicit solution:      4280231
  
  T i m i n g   i n f o r m a t i o n
                        CPU(seconds)   %CPU  Clock(seconds) %Clock
  ----------------------------------------------------------------
  Initialization ....... 2.9000E+01    1.72     2.9000E+01    1.72
  Element processing ... 1.2140E+03   71.83     1.2071E+03   71.40
  Binary databases ..... 4.7000E+01    2.78     4.6125E+01    2.73
  ASCII database ....... 2.0000E+01    1.18     2.2250E+01    1.32
  Contact algorithm .... 1.9700E+02   11.66     2.0788E+02   12.30
  Contact entities ..... 0.0000E+00    0.00     0.0000E+00    0.00
  Rigid bodies ......... 1.8300E+02   10.83     1.7825E+02   10.54
  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            1.6900E+03  100.00     1.6906E+03  100.00
  
  Problem time       =    1.5422E-04

wuhanwangtao

你把 类存搞大点
  
用edstart  扩展内存试试
下面
Storage required for explicit solution: 4280231

ASKA

你的意思是加一块内存条？
还是加大C盘的虚拟内存的大小？
还是在initiative时加大default Memory Size 的大小？？
  
不好意思，我自己胡涂了

wuhanwangtao

我做了一下
  
活塞与挡面 模型图

wuhanwangtao

活塞单元

wuhanwangtao

最后一步，活塞的等效应力图

wuhanwangtao

最后一步，活塞的变形

wuhanwangtao

活塞受力动画

wuhanwangtao

阵个模型的动画

wuhanwangtao

finish  
/cle  
/NOPR  
/PMETH,OFF,0  
KEYW,PR_SET,1  
KEYW,PR_STRUC,1  
KEYW,LSDYNA,1  
  
!******************************  
/UNITS,SI  
!******************************  
  
/plopts,info,0  
!/view,,1,1,1  
/VIEW, 1 ,,-1  
  
dw=0.001 !1 millimetre毫米 =0.03937 inch 英寸  
  
C_PI=3.1415926 !PI的值  
RHO=2800 !铝合金的密度  
VTOT=0   !活塞的体积初始化  
WT=0 !活塞的质量初始化  
SURF=1E7 !压力  
C_R=85*dw/2 !C_R为底面半径  
C_F=SURF*C_PI*C_R*C_R !C_F压力  
       
ACAZ=0 !Z方向的加速度  
  
/prep7  
ET,1,SHELL163 !铝合金和钢的二维有限元模型的et  
  
ET,2,solid164 !铝合金的solid164  
ET,3,solid164 !钢的solid164 刚体  
  
MP,DENS,1,2800 !aluminum  
MP,EX,1,72e9  
MP,NUXY,1,0.33  
TB,BKIN,1,,,,  
TBDAT,1, 480e6  
TBDAT,2, 27e9  
  
MP,DENS,2,7800  
MP,EX,2,200e9  
MP,NUXY,2,0.3  
EDMP,RIGI,2,7,7  
  
wpro,,90.000000,  
  
!建模活塞  
RECTNG,0,85/2000,20/1000+dw,31/1000+dw  
RECTNG,83/2000,85/2000,0/1000+dw,20/1000+dw  
ADEL,ALL  
LANG,1,5,-(90-5.5)  
LFILLT,1,10,5/1000  
LDELE,7,9,1  
LDELE,2,6,4  
L,3,6  
  
LANG,3,10,90  
LANG,2,11,90  
LDIV,11  
L,3,14  
  
AL,4,1,7,6  
AL,7,11,13,3  
AL,12,9,2,13  
AL,10,5,8,9  
CM,HUOSAI,AREA  
!分格活塞二维模型  
  
TYPE, 1  
mat,1  
smrtsize,6
MSHKEY, 0
ESIZE, ,5  
AMESH, HUOSAI  
!建模分格活塞二维模型finish  
  
K,1000, 0  
K,1001, 0, 0,5  
TYPE, 2  
mat,1  
  
VROTAT, ALL, , , , , , 1000, 1001, 360,4  
CMSEL, S,HUOSAI  
ACLEAR, ALL  
ALLSEL  
EPLOT  
  
!/eof  
!1  
  
!***********************  
*DO,I,101,8100  
  *GET,VCAI,ELEM,I,VOLU  
    VTOT=VTOT+VCAI   !总的体积  
*enddo  
  
WT=RHO*VTOT !质量WT  
!**********************  
ACAZ=100!Z方向的加速度  
!/eof  
!2  
  
!建模制动面  
wpro,,-90.000000,  
k,1011,85.1*dw/2,0,0*dw  
k,1012,108*dw/2,0,0*dw  
  
k,1013,80*dw/2,0,-25*dw  
k,1014,82*dw/2,0,-25*dw  
k,1015,108*dw/2,0,-25*dw  
  
k,1016,80*dw/2,0,-47*dw  
k,1017,82*dw/2,0,-47*dw  
k,1018,108*dw/2,0,-47*dw  
  
/pnum,kp,1  
/pnum,line,1  
/pnum,area,1  
gplot  
  
l,1011,1012  
l,1012,1015  
l,1015,1014  
l,1014,1011  
  
l,1013,1014  
l,1014,1017  
l,1017,1016  
l,1016,1013  
  
l,1017,1018  
l,1018,1015  
  
!/eof  
!3  
  
al,82,83,84,85 !area65  
al,84,91,90,87 !area66  
al,86,87,88,89 !area67  
  
asel,s,,,65  
asel,a,,,66  
asel,a,,,67  
  
CM,ZHIDONG,AREA  
  
!/eof  
!4  
  
TYPE, 1  
mat,2  
smrtsize,6
  
MSHKEY, 0
AMESH, ZHIDONG  
  
TYPE, 3  
mat,2  
  
VROTAT, ALL, , , , , , 1000, 1001, 360,4  
CMSEL, S,ZHIDONG  
ACLEAR, ALL  
ALLSEL  
  
/type,1,8  
  
EPLOT  
  
/eof  
  
ESEL,S,MAT,,1  
NSLE,S  
CM,BOX,NODE  
  
ALLSEL,ALL  
ESEL,S,MAT,,2  
NSLE,S  
CM,TABLE,NODE  
  
ALLSEL,ALL  
!EDPART,CREATE  
  
!生成接触ASSC  
EDCGEN,ASSC  
  
edclist,all  
  
finish  
  
/solu  
!修改了速度将公制改为英制 速度方向从y方向转为z方向  
EDVE,VELO,BOX,0,0,-200,0,0,0, , , , , ,  
  
*DIM,time,ARRAY,2,1,1, , ,  
*DIM,accc,ARRAY,2,1,1, , ,  
*SET,TIME(1,1,1) , 0  
*SET,TIME(2,1,1) , 3.0e-4  
*SET,ACCC(1,1,1) , ACAZ  
*SET,ACCC(2,1,1) , ACAZ  
  
    
EDLOAD,ADD,ACLZ,0,BOX,TIME,ACCC, 0, , , , ,  
  
TIME,3.0e-4, !终止时间  
EDRST,50,         
EDHTIME,50,      
EDDUMP,1,  
EDENERGY,1,1,1,1  
  
edhist,box  
  
edout,all  
  
nsel,all  
allsel,all  
/replot  
finish  
  
!/eof  
!  
/SOLU  
/STAT,SOLU  
SOLVE  
FINISH  
  
/EXIT,ALL

wuhanwangtao

应该是你网格化分有问题：）））））））））））））））））））

ASKA

谢谢
我真幸运，遇到一个热心人，感激之情溢于言表。
  
拜托你用bottom显示
然后/type,all,8  Q-zbuff 显示
再在post1里  animate一下von mise
看看有剪切么？
我看了你做的整个的动画，看上去好像活塞还没有接触到制动台面就停止了，
而我们的实际工作时是有一点剪切的
原因在哪里呢？
  
三番五次打搅你，实在是不好意思，不过大家都是做这个东西，互相帮助，共同进步嘛（我脸皮太厚，明明是自己打扰人，还在找接口）
  
再次表示我衷心的谢意！！
  
把剪切进行到底！！！
  
真像一部电影的名字了《没完没了》！呵呵

ASKA

wuhanwangtao wrote:
修改接触的高级设置
  
EDCONTACT，中的PENCHK
  
你把命令流给我
让我看看
  
我想应该是参数有问题

  
你说的EDCONTACT，中的PENCHK
是怎么设定的？
  
我觉着EDCGEN,ASSC中没有体现出摩擦，是不是该加上？
  
你有没有遇到像我那样的错误（但是仍然可以看到一部分结果），？
  
奇怪了，以前做分析，如果在分析时出现error，肯定歇菜了，一步结果也不会
看到，这个分析虽然出现error，但是还可以读出一部分substep的结果，正是搞不懂，你知道这是怎么回事么？

wuhanwangtao

最后一步的sxy

wuhanwangtao

最后一步的sxz

wuhanwangtao

最后一步的syz