个人收集的一个例程,不过由于我没做过这个方面的,所以没用过,据说里边有一定的错误,你可以参考一下基本思路吧
ANSYS做的轮胎模型
- fini
- /cle
- /prep7
- /title,wheel-3
- !_____________________________________定义单元类型__________________________________________________
- et,1,45 !用来给土层,轮辋以及基层来分网格
- et,2,181
- !____________________________________________________________________________________________________
- !____________________________________定义轮胎材料参数_________________________________________________
- !___________________________定义超弹材料,用9常数模型(单位为MPa)___________________
- mp,ex,1,8000e6
- mp,prxy,1,0.4
- !____________________________定义正交弹性材料(单位为MPa)___________________________
- mp,ex,2,1.03e10
- mp,prxy,2,0.33
- mp,ex,3,2.1e11
- mp,prxy,3,0.3
- save,material4,db
- /prep7
- !_______________________________________定义道路几何参数___________________________________________
- areaload=0.557800e6 !单位mpa
- arearadiu=0.1066 !单位m
- top=0.5493
- thick1=0.04+top !单位m
- thick2=0.06
- thick3=0.07
- thick4=0.25
- thick5=0.27
- thick6=3.5
- distanceofaxie=1.5*arearadiu !两轮中心距
- rightb=1.50
- leftb=-1.50
- frontb=1.50
- backb=-1.50
- depthb=-(thick1+thick2+thick3+thick4+thick5+thick6)
- !___________________________________________________________________________________________________
- !______________________________________定义道路材料参数_____________________________________________
- ex1=1200e6 !单位MPa
- prxy1=0.35
- dens1=2.509e-03 !单位kg/m3
- ex2=1000e6
- prxy2=0.35
- dens2=2.447e-03
- ex3=1000e6
- prxy3=0.35
- dens3=2.432e-03
- ex4=1400e6
- prxy4=0.2
- dens4=2e-03
- ex5=800e6
- prxy5=0.2
- dens5=1.800e-03
- ex6=35e6
- prxy6=0.4
- dens6=1.900e-03
- !___________________________________________________________________________________________________
- !________________________________定义道路材料参数__________________________________________________
- mp,ex,4, ex1 !定义道路第一层材料参数
- mp,PRXY,4,prxy1
- mp,dens,4,dens1
- mp,ex,5, ex2
- mp,PRXY,5,prxy2
- mp,dens,5,dens2
- mp,ex,6, ex3
- mp,PRXY,6,prxy3
- mp,dens,6,dens3
- mp,ex,7, ex4
- mp,PRXY,7,prxy4
- mp,dens,7,dens4
- mp,ex,8, ex5
- mp,PRXY,8,prxy5
- mp,dens,8,dens5
- mp,ex,9, ex6
- mp,PRXY,9,prxy6
- mp,dens,9,dens6
- !__________________________________________________________________________________________________
- /prep7
- !______________________________建立轮胎模型_________________________________________________________
- k,1,0,0,0
- k,2,0,0,0.03
- k,3,0,0.01,0
- k,4,0.297,0,0
- k,5,0.297,0,0.2
- k,6,0.297,0,-0.2
- k,7,0.297,0,0.01
- k,8,0.297,0,-0.01
- k,9,2,0,0.01
- k,10,0,0,-0.005
- k,11,0,0,0.005
- k,12,0.4028,0,0
- k,13,0.5425,0,0
- k,14,0.4028,0.01,0
- circle,1,0.297,2,4,360,1
- l,5,6
- LGEN,2,4,,1,0.01
- l,6,19
- l,5,18
- lsel,s,,,4,7,1
- al,all
- lsel,all
- VDRAG,1,,,,,,1,2,3
- vsel,s,,,1,3,1
- vadd,all !4
- torus,0.4028,0.1362,0.1422 !1
- torus,0.4028,0.1422,0.1465 !2
- vsbv,2,4,sepo,dele,keep !3,5
- vdele,3
- vsbv,1,4,sepo,dele,keep !2,3
- vdele,2
- wpoffs,,, 0.110836
- WPSTYLE,,,,,,,,1
- wpstyle,0.05,0.1,-1,1,0.003,0,0,,5
- vsbw,4,speo,dele
- vdele,2
- wpoffs,,, -0.110836*2
- vsbw,1,speo,dele !轮辋宽0.11085
- vdele,4
- l,11,10 !40
- circle,11,0.297,2,7,360,1
- lsel,s,,,41,43,2
- lsel,a,,,56
- al,all !11
- lsel,all
- vdrag,11,,,,,,40 !1
- vsel,s,,,1,2
- vadd,all !4
- vsel,all
- vglue,3,4,5 !1内胎,2外胎,6轮辋
- !分网格
- HPTCREATE,area,11,, coord,0,0,0.005
- HPTCREATE,area,73,, coord,0,0,-0.005
- vsel,s,,,6
- type,1
- mat,3
- mshape,1,3D
- mshkey,0
- esize,0.04
- !要粗化
- vmesh,all
- vsel,all !1336,1178
- csys,0
- !给外轮胎分网格
- !HPTCREATE,area,21,, coord,0,0.5493,0
- vsel,s,,,2
- type,2
- mat,1
- mshape,1,3D !定义单元形状
- mshkey,0
- esize,0.02 !定义为自由网格
- !可能要局部粗化
- vmesh,all
- vsel,all
- !给轮胎束带层分网格
- vsel,s,,,1
- type,2
- mat,2
- mshape,1,3D
- mshkey,0
- esize,0.02
- !可能要局部粗化
- vmesh,all
- !在内胎的内面上施加面荷载
- !13,14,67-70
- asel,s,,,13,14
- asel,a,,,67,70,1
- nsla,s
- SF,all,pres,0.63e6 !压力单位为Pa
- csys,0
- save,finimodel,db
- /prep7
- !____________________________________________________定义路面_______________________________________________
- block,rightb,leftb,-top,-thick1,frontb,backb !3
- block,rightb,leftb,-thick1,-(thick1+thick2),frontb,backb !4
- block,rightb,leftb,-(thick1+thick2),-(thick1+thick2+thick3),frontb,backb !5
- block,rightb,leftb,-(thick1+thick2+thick3),-(thick1+thick2+thick3+thick4),frontb,backb!7
- block,rightb,leftb,-(thick1+thick2+thick3+thick4),-(thick1+thick2+thick3+thick4+thick5),frontb,backb!8
- block,rightb,leftb,-(thick1+thick2+thick3+thick4+thick5),depthb,frontb,backb!9
- vglue,3,4,5,7,8,9 !3,10,11,12,13,14
- vsel,s,,,3
- type,1
- mat,4
- mshape,1,3D !定义单元形状
- mshkey,0
- esize,0.0.5 !定义为自由网格
- vmesh,all
- vsel,s,,,10
- type,1
- mat,5
- mshape,1,3D !定义单元形状
- mshkey,0
- esize,0.15 !定义为自由网格
- vmesh,all
- vsel,s,,,11
- type,1
- mat,6
- mshape,1,3D !定义单元形状
- mshkey,0
- esize,0.18 !定义为自由网
- vmesh,all
- vsel,s,,,12
- type,1
- mat,7
- mshape,1,3D !定义单元形状
- mshkey,0
- esize,0.2 !定义为自由网格
- vmesh,all
- vsel,s,,,13
- type,1
- mat,8
- mshape,1,3D !定义单元形状
- mshkey,0
- esize,0.3 !定义为自由网格
- vmesh,all
- vsel,s,,,14
- type,1
- mat,9
- mshape,1,3D !定义单元形状
- mshkey,0
- esize,0.5 !定义为自由网格
- vmesh,all
- nsel,all
- vsel,all
- asel,all
- fini
- /solu
- !__________________________________给道路加约束_______________________________________________
- csys,0
- nsel,all
- nsel,r,loc,y,depthb
- d,all,all,0
- nsel,all
- nsel,r,loc,x,rightb
- nsel,a,loc,x,leftb
- d,all,ux,0
- nsel,all
- nsel,r,loc,z,frontb
- nsel,a,loc,z,backb
- d,all,uz,0
- save,pavemodel,db
- fini
- /prep7
- !______________________________加接触对_________________________________________
- r,1001
- et,11,targe170
- et,12,conta174
- KEYOPT, 12,2,1 !定义使用罚函数
- KEYOPT, 12,12,1 !定义为无限大,忽略滑动
- KEYOPT, 12,7,3 !对下一个子步预测一个最小的时间增量
- asel,s,,,22
- asel,a,,,100,101
- nsla,s
- !nsel,r,loc,y,-0.4,-0.5493,0.0001
- cm,_contact12,node
- type,12
- real,1001
- esurf,,buttom
- asel,s,,,9 !目标面,选择较为粗的网格
- nsla,s
- !nsel,r,loc,x,-0.8,0.8,0.0001
- !nsel,r,loc,z,-0.7,0.7,0.0001
- cm,_target12,node
- type,11
- real,1001
- esurf,,top
- save,contact,db
- fini
- fini
- /solu
- nsel,s,,,206
- nsel,a,,,412
- d,all,all,0
- !d,all,uy,0
- vsel,s,,,1,2
- nsel,r,loc,x,0
- d,all,ux,0
- vsel,s,,,1,2
- vsel,a,,,6
- nsel,r,loc,z,0
- d,all,uz,0
- antype,static
- asel,all
- nsel,all
- time,1
- solve
- 但是其接触部分的UY,UX方向位移明显错误(问题1)
- 当我把材料一改为超弹材料
- mp,ex,1,9000
- mp,prxy,1,0.49
- MPTEMP,,,,,,,,
- MPTEMP,1,0
- MPDATA,DENS,1,,1150
- TB,mooney,1,1,9,MOON
- TBTEMP,0
- TBDATA,,0.770,,0.0671,,,0.0955
- TBDATA,,,,,,,
- MPTEMP,,,,,,,,
- MPTEMP,1,0
复制代码
|