【分享】用LS-DYNA的切断解析

ansyscai

                /prep7
  et,1,162
  keyopt,1,3,2
  mp,dens,1,7.865e-6
  mp,ex,1,200e6
  mp,nuxy,1,0.27
  tb,plaw,1,,,1
  tbdata,1,310e3
  tbdata,2,768e3
  tbdata,3,0
  tbdata,4,40.0
  tbdata,5,5.0
  tbdata,6,0.75
  edmp,rigid,2,6,7
  mp,dens,2,7.865e-6
  mp,ex,2,200e3
  mp,nuxy,2,0.27
  edmp,rigid,3,6,7
  mp,dens,3,7.865e-6
  mp,ex,3,200e3
  mp,nuxy,3,0.27
  edmp,rigid,4,7,7
  mp,dens,4,7.865e-6
  mp,ex,4,200e3
  mp,nuxy,4,0.27
  edmp,hgls,1,5
  rectng,5,9.5,0.01,2.01,  
  rectng,9.5,10.5,0.01,2.01,   
  rectng,10.5,15,0.01,2.01,   
  rectng,8,9.99,3,5,   
  rectng,10,12,2.1,4.1,   
  rectng,10,15,0,-2,   
  /pnum,area,1
  aplot
  aglue,1,2,3
  asel,s,,,1
  asel,a,,,7
  asel,a,,,8
  aatt,1
  asel,s,,,4
  aatt,2
  asel,s,,,5
  aatt,3
  asel,s,,,6
  aatt,4
  asel,all
  /pnum,line,1
  /pnum,area,0
  lplot
  lsel,s,,,27
  lsel,a,,,3
  lsel,a,,,28
  lesize,all,,,20,5
  lsel,s,,,1
  lesize,all,,,20,0.2
  lsel,s,,,25
  lsel,a,,,26
  lesize,all,,,11
  lsel,s,,,2
  lsel,a,,,4
  lsel,a,,,6
  lsel,a,,,10
  lesize,all,,,20
  lsel,all
  esize,,1
  mshape,0
  mshkey,1
  amesh,all
  /pnum,mat,1
  /num,1
  edcgen,ass2d,,,0.6
  edpart,create
  partsel,'plot'
  fini
  /solu
  *dim,mtime,,2
  *dim,move,,2
  mtime(1)=0,0.01
  move(1)=0,-3
  *dim,ftime,,3
  *dim,force,,3
  ftime(1)=0,0.001,1
  force(1)=0,-500e3,-500e3
  edload,add,rbuy,,2,mtime,move
  edload,add,rbfy,,3,ftime,force
  time,0.01
  edrst,50
  edhtime,100
  edout,glstat
  save
  solv
  fini
  /post1
  eplot
  /user,1
  /dscale,,1
  set,last
  pldisp
  avprin,,0.5
  ples,eppl,eqv
  andyna,0.1,,,,,1

foransys

希望ansyscai对程序做一个简要的说明。针对程序中的关键技术做个简评。
以便给推荐到电子期刊发表。

ansyscai

（用LS-DYNA的切断解析―補充説明）目的：但求共同進歩
在ANSYS/LS-DYNA、作為材料特性有能設定破損応変値的材料模型。通過設定此破損応変値、当相当塑性応変Von Mises到達破損応変値時其要素将不被計算。根拠此原理可進行破壊模型的仿真。本例是根拠此机能進行切断解析的。解析模型以及条件等如以下。
本解析把此作為２次元平面応変問題処理。被加工材料為鋼相当的２直線近似移動硬化則模型、沖刀、固定台、工件的圧件為鋼相当的剛体。
解析模型的基本単位為、mm、kg、sec。（関連単位為、mN、kPa。）
要素形式： &nbspLANE162Keyopt(3)=2 平面応変
材料1：被切削工件  ２直線近似移動硬化則模型
　材料常数  楊氏模量  200×106 kPa
  箔松比  0.27
  質量密度  7.865×10-6 kg/mm3
  降伏点  200×106 kPa
  接線係数  200×106 kPa
  応変速度（c）  40.0 1/sec
  応変速度（p）  5
  破損応変値  0.75
材料2：沖刀  剛体
　材料常数  楊氏模量  200×106 kPa
  箔松比  0.27
  質量密度  7.865×10-6 kg/mm3
  約束  唯UY自由
材料3：工件圧件  剛体
　材料常数  楊氏模量  200×106 kPa
  箔松比  0.27
  質量密度  7.865×10-6 kg/mm3
  約束  唯UY自由
材料4：固定台  剛体
　材料常数  楊氏模量  200×106 kPa
  箔松比  0.27
  質量密度  7.865×10-6 kg/mm3
  約束  全約束
/prep7
! LS-DYNA用２維要素「PLANE162」定義
et,1,162
keyopt,1,3,2
!材料模型的定義
!材料１：具有破損応変値的２直線近似移動硬化則塑性模型、
!材料２～４：定義為剛体。（定義材料２後、複製到３，４
!有関剛体的約束条件、在定義材料模型時進行
mp,dens,1,7.865e-6
mp,ex,1,200e6
mp,nuxy,1,0.27
tb,plaw,1,,,1
tbdata,1,310e3
tbdata,2,768e3
tbdata,3,0
tbdata,4,40.0
tbdata,5,5.0
tbdata,6,0.75
  
edmp,rigid,2,6,7
mp,dens,2,7.865e-6
mp,ex,2,200e3
mp,nuxy,2,0.27
  
edmp,rigid,3,6,7
mp,dens,3,7.865e-6
mp,ex,3,200e3
mp,nuxy,3,0.27
  
edmp,rigid,4,7,7
mp,dens,4,7.865e-6
mp,ex,4,200e3
mp,nuxy,4,0.27
  
edmp,hgls,1,5 !hourglass control定義
!建模
!細分切断部分的mesh
  rectng,5,9.5,0.01,2.01,  
  rectng,9.5,10.5,0.01,2.01,   
  rectng,10.5,15,0.01,2.01,   
  
  rectng,8,9.99,3,5,   
  rectng,10,12,2.1,4.1,   
  rectng,10,15,0,-2,   
  
  /pnum,area,1
  aplot
  
  aglue,1,2,3
asel,s,,,1
asel,a,,,7
asel,a,,,8
aatt,1
asel,s,,,4
aatt,2
asel,s,,,5
aatt,3
asel,s,,,6
aatt,4
asel,all
/pnum,line,1
/pnum,area,0
  lplot
lsel,s,,,27
lsel,a,,,3
lsel,a,,,28
lesize,all,,,20,5
lsel,s,,,1
lesize,all,,,20,0.2
lsel,s,,,25
lsel,a,,,26
lesize,all,,,11
  
lsel,s,,,2
lsel,a,,,4
lsel,a,,,6
lsel,a,,,10
lesize,all,,,20
lsel,all
  
esize,,1
mshape,0
mshkey,1
  
amesh,all
/pnum,mat,1
/num,1
edcgen,ass2d,,,0.6　!接触的定義、摩擦係数0.6
edpart,create！為定義荷重、建part
partsel,'plot'
fini
/solu
  *dim,mtime,,2　!為定義荷重　定義行列
  *dim,move,,2　!為定義荷重　定義行列
mtime(1)=0,0.01
move(1)=0,-3
  
*dim,ftime,,3!為定義荷重　定義行列
*dim,force,,3!為定義荷重　定義行列
  ftime(1)=0,0.001,1
  force(1)=0,-500e3,-500e3
  
  edload,add,rbuy,,2,mtime,move！荷重定義
  edload,add,rbfy,,3,ftime,force　　!荷重定義
！解析終了時間、輸出的設定
  time,0.01！終了時間
!定義文件輸出周期
  edrst,50      ! 輸出step数
  edhtime,100　 !輸出step数
  edout,glstat
  
save
solv
fini
!結果表示
/post1
eplot
/user,1
/dscale,,1
set,last
pldisp
avprin,,0.5
ples,eppl,eqv
  
andyna,0.1,,,,,1! animation

lgh790723

不错。十分感谢！！！

zlj

鼓励原创! 顶一下

liuzongfahjx

楼主好人呀，向你致敬

taojunhust

呵呵，好强啊，我在学习几年看能否达到这个水平。呵呵。

asdzhy2005

呵呵，学习学习

lyju

强啊学习学习

zax77

很不错.请教高手,这个模型可以用来分析如剪刀切断钢筋的过程吗?尤其是对于剪刀的边角等问题?