hy1966 发表于 2012-7-24 18:15
1、D,ALL,ALL !施加固定约束
这也包括电压约束为0,显然不正确
2、MP,RSVX,2,0 ...
谢谢您的指点,我又把所有的节点又都耦合了一遍,但是计算时显示“There are elements with negative and positive electric charge reaction solutions for the volt degree of freedom in the model. The analysis is not valid.”,请您帮我看一下是哪里出了问题呢
以下是我的命令流
/prep7
*SET,tpzt,0.005
*SET,tcem,0.02
*SET,r,0.05
ET,1,PLANE223,1001
MP,PERX,1,635
MP,PERY,1,730
MP,PERZ,1,730
TB,PIEZ,1,,,1
TBDATA,18,5.25E-10
TBDATA,11,5.25E-10
TBDATA,7,-1.065e-10
TBDATA,4,-1.065e-10
TBDATA,1,3E-10
TB,ANEL,1,1,,1
TBTEMP,0
TBDATA,1,1.61E-11,-5.52E-12,-5.52E-12
TBDATA,7,1.24E-11,-3.98E-12
TBDATA,12,1.24E-11
TBDATA,16,3.91E-11
TBDATA,19,3.27E-11
TBDATA,21,3.91E-11
MP,PERX,2,1e10
TB,ANEL,2,1,21,1
TBTEMP,0
TBDATA,,7.28E-012,-2.185E-012,-2.185E-012,0,0,0
TBDATA,,1.6E-012,0,0,0,0,2E-012
TBDATA,,0,0,0,1.894E-011,0,0
TBDATA,,0,0,0,,,
CSYS,1
CYL4,0,0,r+tpzt,0,r,30,l
CYL 4,0,0,r+tpzt+tcem,0, r+tpzt,30
CYL 4,0,0, r+2*tpzt+tcem,0, r+tpzt+tcem,30
CYL 4,0,0, r+2*tpzt+2*tcem,0, r+2*tpzt+tcem,30
CYL 4,0,0, r+3*tpzt+2*tcem,0, r+2*tpzt+2*tcem,30
CYL 4,0,0, r+3*tpzt+3*tcem,0, r+3*tpzt+2*tcem,30
CYL 4,0,0, r+4*tpzt+3*tcem,0, r+3*tpzt+3*tcem,30
AGLUE,1,2,3,4,5,6,7
ASEL,S,AREA,,1
AATT,1,,1
ASEL,S,AREA,,9
AATT,1,,1
ASEL,S,AREA,,11
AATT,1,,1
ASEL,S,AREA,,13
AATT,1,,1
ASEL,S,AREA,,8
AATT,2,,1
ASEL,S,AREA,,10
AATT,2,,1
ASEL,S,AREA,,12
AATT,2,,1
ASEL,ALL
LESIZE,3,,3 !以线3为标准线建立角度间隔为3度的线尺寸
MSHKEY,1 !设置映射网格划分
MSHAPE,0,2D
SHPP,OFF !关闭单元形状检查以允许生成边缘实体网格
AMESH,ALL
NUMMRG,KP !关键点编号压缩合并
CSYS,1
NSEL,S,LOC, Y,0
D,all,UX,0
D,all,UY,0 !施加固定约束
LSEL,S,LINE,,3
NSLL,S,1
CM,BOT,NODE
D,BOT,VOLT,0
LSEL,S,LINE,,25
NSLL,S,1
CM,TOP,NODE
CP,2,VOLT,ALL
*GET,N2,NODE,0,NUM,MIN
D,N2,VOLT,100
LSEL,S, , , 1
NSLL,S,1
CM,pzt0,NODE
CP,3,VOLT,ALL
LSEL,S, , , 5
NSLL,S,1
CM,pzt1,NODE
CP,4,VOLT,ALL
LSEL,S, , , 9
NSLL,S,1
CM,pzt2,NODE
CP,5,VOLT,ALL
LSEL,S, , , 13
NSLL,S,1
CM,pzt3,NODE
CP,6,VOLT,ALL
LSEL,S, , , 17
NSLL,S,1
CM,pzt4,NODE
CP,7,VOLT,ALL
LSEL,S, , , 21
NSLL,S,1
CM,pzt5,NODE
CP,8,VOLT,ALL
NSEL,S,LOC,X,0.0525
CP,9,VOLT,ALL
NSEL,S,LOC,X,0.065
CP,10,VOLT,ALL
NSEL,S,LOC,X,0.0775
CP,11,VOLT,ALL
NSEL,S,LOC,X,0.090
CP,12,VOLT,ALL
NSEL,S,LOC,X,0.1025
CP,13,VOLT,ALL
NSEL,S,LOC,X,0.115
CP,14,VOLT,ALL
NSEL,S,LOC,X,0.1275
CP,15,VOLT,ALL
NSEL,S,LOC,X,0.05125
CP,16,VOLT,ALL
NSEL,S,LOC,X,0.05375
CP,17,VOLT,ALL
NSEL,S,LOC,X,0.06
CP,18,VOLT,ALL
NSEL,S,LOC,X,0.070
CP,19,VOLT,ALL
NSEL,S,LOC,X,0.07625
CP,20,VOLT,ALL
NSEL,S,LOC,X,0.07875
CP,21,VOLT,ALL
NSEL,S,LOC,X,0.085
CP,22,VOLT,ALL
NSEL,S,LOC,X,0.095
CP,23,VOLT,ALL
NSEL,S,LOC,X,0.10125
CP,24,VOLT,ALL
NSEL,S,LOC,X,0.10375
CP,25,VOLT,ALL
NSEL,S,LOC,X,0.110
CP,26,VOLT,ALL
NSEL,S,LOC,X,0.120
CP,27,VOLT,ALL
NSEL,S,LOC,X,0.12625
CP,28,VOLT,ALL
NSEL,S,LOC,X,0.12875
CP,29,VOLT,ALL
/SOLU
ANTYPE,STATIC
NSEL,ALL
SOLVE !求解
|