大多数情况下,ANSYS求解规模受限于机器配置(内存和外存)。即便选择方程求解器,如表1所示,求解规模总是有限度的自由度。为此,利用图1所示的空间悬臂结构(仅为试算而设),从直接划分网格到求解,可以大致了解你的机器配置所能求解的规模。
图1结构采用SOLID185单元,设置不同单元尺寸时,对应的自由度总数(方程数)不同,大致如下:
ESIZE=50,大约
ESIZE=40,大约140万自由度
ESIZE,30,大约320万自由度
ESIZE,20,大约1000万自由度
ESIZE,10,大约7600万自由度(网格划分和求解超花时间)
采用何种求解器,可以自己设置;若让程序选择,可采用“EQSLV,-1”设置。
1000万自由度规模时用PCG可轻松完成,结果如图2所示。
因此,当直接求解器不能求解时(缺省),可以试试其他求解器,或者让程序自动选择求解器,有可能会给你带来惊喜。
FINISH$/CLEAR$/PREP7
B=1000$H=2000$R=200$B1=600
H1=1200$L=10000$HD=1200
BLC4,,,B1/2,H1/2
CYL4,B1/2-R,H1/2-R,R,,,90
AOVLAP,ALL
ASEL,S,LOC,Y,H1/2-R,H1/2
ASEL,R,LOC,X,B1/2-R/2,B1/2
ADELE,ALL,,,1$ASEL,ALL
AADD,ALL$NUMMMP,ALL
BLC4,,,B/2,H/2$ASBA,2,1
WPOFF,B1/2-R$WPROTA,,,90
ASBW,ALL
WPOFF,,H1/2-R,R
WPROTA,,90$ASBW,ALL
WPOFF,,,-R-R/2$ASBW,ALL
KWPAVE,KP(B/2,H1/2 R/2,0)
WPROTA,,-45$ASBW,ALL
WPCSYS,-1
ARSYM,X,ALL$ARSYM,Y,ALL
NUMMRG,ALL$NUMCMP,ALL
AGEN,9,ALL,,,B
AGEN,1,ALL,,,B,,,,,1
ASEL,NONE
BLC4,,-H/2,B/2,H
BLC4,L-B/2,-H/2,B/2,H
WPROTA,,90$ASBW,ALL!
WPOFF,,,H1/2-R$ASBW,ALL
WPOFF,,,R R/2$ASBW,ALL
WPOFF,,,-(H1-R/2)$ASBW,ALL
WPOFF,,,-(R R/2)$ASBW,ALL
WPCSYS,-1$ASEL,ALL$NUMMRG,ALL
NUMCMP,ALL
K1=KP(0,0,0)$K,500,0,0,HD
L,K1,500$LINE1=_RETURN
VDRAG,ALL,,,,,,LINE1
LDELE,LINE1
ET,1,SOLID185
MP,EX,1,2.1E5$MP,PRXY,1,0.3
!修改尺寸设置改变自由度数
ESIZE,50!大约30万自由度
MSHKEY,1$VMESH,ALL
*get,nn,node,,count
nn=nn*3
ASEL,S,LOC,X,0$DA,ALL,ALL
ASEL,S,LOC,Y,H/2
SFA,ALL,1,PRES,1.0
ALLSEL,ALL
/SOLU
!选择求解器
!EQSLV,PCG!较SPARSE快
EQSLV,-1!自动选择求解器
SOLVE$FINISH
/POST1$PLNSOL,U,SUM
---------------------------------------------------------------------------------------------
版权声明:
原创文章,来源铁大王新敏,本文已经授权,欢迎分享,如需转载请联系作者。