轮齿结构优化设计
/batch
*set,r1,725
*set,r2,795
*set,b,35
*set,r3,r2+b
*set,ang,360/28
*set,ang1,ang/2
*set,h,175
*set,h1,75
*set,du,0.02*h
*set,sigma,3*324
/prep7
et,1,solid92
mptemp,,,,,,,,
mptemp,1,0
mpdata,ex,1,,2.06e5
mpdata,prxy,1,,0.28
cyl4,,,r1,0,r2,ang,h
cyl4,,,0,0,r1,ang,h
cyl4,,,r2,ang1,ang,h
cyl4,,,r2,ang1,r3,ang,h1
n,1,0,0,0
n,2,0,0,h
csys,1
n,3,r2,ang1,h1
n,4,r3,ang1,h1
csys,0
vadd,1,2,3
lesize,22,,,2,,,,,1
lesize,27,,,2,,,,,1
lesize,29,,,2,,,,,1
lesize,39,,,2,,,,,1
esize,25,0,
mshkey,0
mshkey,0
mshape,1,3d
vmesh,4
nummrg,all,,,,low
numcmp,all
/view,1,1,1,1
/triad,ltop
/replot
fini
/solu
csys,1
asel,s,loc,y,0
asel,a,loc,y,ang
da,all,symm
allsel,all
asel,s,loc,x,r2,r3
asel,r,loc,z,h1
da,all,uz,
allsel,all
asel,s,loc,z,0
asel,r,loc,x,0,r1
sfa,all,1,pres,7.7
allsel,all
csys,0
solve
fini
/post1
etable,evol,volu,
ssum
*get,vtot,ssum,,item,evol
*set,wt,28*7.85*vtot*1.0e-6
path,aa,2,30,20,
ppath,1,1
ppath,2,2
pdef,,s,eqv,avg
pdef,,u,z,avg
*get,s1,path,,max,seqv
*get,uz,path,,max,uz
path,bb,2,30,20,
ppath,1,3
ppath,2,4
pdef,,s,eqv,avg
*get,s2,path,,max,seqv
eplot
save
lgwrite,gate,lgw,,comment
fini
/opt
opanl,'gate','lgw'
opvar,h,dv,140,175,,
opvar,h1,dv,55,90,,
opvar,b,dv,30,50,,
opvar,s1,sv,,204,,
opvar,s2,sv,,sigma,,
opvar,uz,sv,,du,,
opsave,'gate','opt',''
opvar,wt,obj,,,50,
optype,firs
opfrst,12,,,
save
opexe
oplist,all,,0
/gropt,axnsc,1.5,
plvaropt,h
/image,save,h,bmp
plvaropt,wt
/image,save,wt,bmp
/image,save,wt,bmp
plvaropt,h1,b
/image,save,hb,bmp
plvaropt,s1,s2
/image,save,s12,bmp
finish

希望大家多发例子 共同学习