[讨论]APDL命令流挑战赛一

[windows251]

编写生成如下图形的apdl命令流，参数自设，只要形近就行。不同的人因为对apdl命令的掌握程度和理解不同会写出不一样的代码。大家都出来晒晒！！
(下面给出了一些命令流，以供参考)

[windows251]

fini
/clea
/filn,wheel
/title,wheel-csys
R=40
r1=20
r=8
nn=5
/prep7
/view,,1,1,1
csys,0
cylind,0,40,0,20,,
cm,cyc,volu
wpro,-90,,
wpoff,r1*sin(180/nn),r1*cos(180/nn),
cswpla,11,1,,
cylind,0,8,0,20,,
csys,1
vgen,nn,2,,,,360/nn,,,0
cm,hole,volu
csel,all
vsbv,cyc,hole
csys,0
save,all
finish

[windows251]

/PREP7
CYL4,0,0,5, , , ,3
CYL4,3,0,1, , , ,3
CSYS,1
VGEN,8,2, , , ,45, , ,1
ALLSEL,ALL
VSEL,U, , , 1
cm,delete,volu
ALLSEL,ALL
VSBV,1,delete
CSYS,0
save,all
finish

[windows251]

“自底向上”建模，点-—线——面——体
fini
/clear
/filn,wheel
/title,wheel
/prep7
k,1,0,0,0
k,2,0,0,12
k,3,1,0,0
circle,1,20,2,3,,,
al,all
*get,numkp,kp,,num,max
wpoff,0,10,
k,numkp+1,0,10,0
k,numkp+2,0,10,2
k,numkp+3,1,10,0
circle,numkp+1,4,numkp+2,numkp+3,,,
al,5,6,7,8
CSYS,1
aGEN,6,2, , , ,60, , ,1
cm,hole,area
asba,1,hole,,delete,delete
csys,0
lstr,1,2
*get,numline,line,,num,max
VDRAG,8,,,,,,numline,,,,,
save,all
fini

[dnping]

虽然不懂，但是还是支持版主的工作……

[windows251]

finish
/CLEAR,NOSTART
/prep7
CYL4,0,0,0.5,36,0.2,0
CYL4,0.3,0.1,0.07,360,0,0
ASBA,1,2
VOFFST,3,0.1,,
CSYS,1
VGEN,10,1,,,,36,,,1
CYL4,0,0,0.2,,,,0.1
VGLUE,all
save,all
finish

[windows251]

fini
/cle
/title,test
/prep7
n=5 !!!孔的个数
cyl4,,,10,,,,4
*do,i,1,n
wpoffs,6,,
CYLIND,2, ,0,4,0,360
*get,vnum,volu,,num,max
vsel,u,,,vnum
cm,v_all,volu
vsel,all
vsbv,v_all,vnum
wpoffs,-6,,
wprota,360/n,,
numcmp,volu
*enddo
fini

[windows251]

/prep7
R1=5
R2=0.8
cyl4,,,R1,,,360,1
*afun,deg
cyl4,2.5,0,R2,,,360,2
cyl4,-2.5,0,R2,,,360,2
csys,4
wprota,22.5
vsymm,y,2,3
wprota,22.5
vsymm,y,2,3
wprota,22.5
vsymm,y,2,3
csys,0
vsel,s,,,2,9,,1
cm,del2,volu
allsel,all
vsbv,1,del2
cmdele,del2
finish

[windows251]

fini
/clea
/prep7
pi=acos(-1)
R=10
r1=6
r2=2
cylind,R,,0,3
*do,i,0,4
cyl4,r1*cos(i*2*pi/5),r1*sin(i*2*pi/5),r2,,,,3
vsbv,1,2,1,delete,delete
numcmp,volu
*enddo

[windows251]

相对来说，第五个和第七个的模型与要求比较符合。
第五个命令流的通用性比较强，可以通过修改n值来建立不同的模型。

[wenlihui20062189]

太喜欢楼主了！我又学到好多东西！以后一定常来逛逛！

第七个我改了一下不知对不对，望各位大侠指教：
fini
/clea
/prep7
pi=acos(-1)
R=10
r1=6
r2=2
n=5 !孔个数
m=3 !板厚度
cylind,R,,0,m
*do,i,0,n-1
cyl4,r1*cos(i*2*pi/n),r1*sin(i*2*pi/n),r2,,,,3
vsbv,1,2,1,delete,delete !这里没看懂
numcmp,volu
*enddo