求解时程序跳出?????

浏览:168448
代码如下:
------------------------------------------------------
!/BACH,LIST
/TITLE,CAPACITANCE OF TWO LONG RECTNG ABOVE A GROUND PLANE
/FILNAME,CAPACITANCE_SENSOR_2,1
KEYW,MAGELC,1
/PREP7 !前处理器
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! 变量定义,类型及单元材料定义
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
PW=50 !电极宽(mm)
G=10 !间隙宽
PH=4 !电极厚度
R0=500 !空气区域半径
BW=500 !基板参数(宽,高)
BH=10
TW=300 !目标板参数(宽,高,距离)
TH=10
TD=10
V1=10
V0=0
*DIM,CON,CHAR,12
CON(1)='CON1','CON2','CON3','CON4','CON5','CON6','CON7','CON8','MUBIAO','BIANJIE','JIBAN','AIR' !对模型各部分面进行元件定义
*DIM,CON_TEMP,CHAR,8
CON_TEMP(1)='COND1','COND2','COND3','COND4','COND5','COND6','COND7','COND8' !对各极板节点进行定义
/PNUM,AREA,1 !以面积方式按序号排列显示
ET,1,121 !单元类型定义
ET,2,110,1,1 !远场类型单元
ET,3,99 !复合材料类型单元
EMUNIT,EPZRO,8.854E-6 !介电常数uMKSV级
MP,PERX,1,1 !空气的介电常数
MP,PERX,2,3.45 !基板的介电常数
MP,PERX,3,3 !目标板参数(复合材料)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! 模型的建立
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
RECTNG,0,PW,BH,PH+BH !极板建模(宽PW,高PH)
AGEN,8,ALL,,1,G+PW,,,,,0 !复制8个同样的极板
AGEN,1,ALL,,,-3.5*G-4*PW,,,,,1 !向左移动极板,使中点位于原点处
RECTNG,-BW/2,BW/2,0,BH !基板建模(宽BW,高BH)
RECTNG,-TW/2,TW/2,TD+BH,TH+BH+TD ! 复合材料板建模(宽TW,高TH,距基板TD)
CYL4,0,0,R0,0,,180 !空气区域建模(半径为R的半圆)
CYL4,0,0,2*R0,0,,180 !远场区域的建模(半径为2R的半圆)
ALLSEL,ALL
AOVLAP,ALL !各面积重叠
NUMCMP,AREA !面积号压缩
*DO,I,1,12 !模型各部分面积元件的定义,以便后续操作
ASEL,S,AREA,,I
CM,CON(I),AREA
*ENDDO
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! 下面为网格划分
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
CMSEL,S,JIBAN
AESIZE,JIBAN,5 !对所选择的面进行单元尺寸的设置,5为单元尺寸值
MSHAPE,1 !0时为四边形剖分,1时为三角行剖分
MSHKEY,0 !0时为自由网格剖分,1时为映射式剖分
AATT,2 !指定所选面的材料属性
TYPE,1 !指定所选面的单元属性
AMESH,ALL !网格划分
MSHAPE,0
MSHKEY,1
*DO,I,1,8 !对极板进行网格划分,使用DO,ENDDO 命令四边形映射剖分
CMSEL,S,CON(I)
AESIZE,CON(I),5
AMESH,ALL
NSLA,S,0 !选择面上所依附的节点,0时仅在所选择面内的节点,1时为面上的所有节点
CM,CON_TEMP(I),NODE !将生成的节点定义为元件,以便后续计算
*ENDDO
CMSEL,S,MUBIAO !复合材料板网格划分
SMRTSIZ,7 !采用智能网格剖分,剖分细化程度为7(1--7)
MSHAPE,1 !材料为3,三角形自由网格剖分
MSHKEY,0
AATT,3
!TYPE,3
AMESH,ALL
CSYS,1 !坐标转换为极坐标系,0为笛卡尔坐标系,1为极坐标,2为球坐标,3为Y轴球坐标
LSEL,S,LOC,X,R0 !选取半径为R0的圆弧
LESIZE,ALL,,,15 !将所选圆弧划分为15等份
!SMRTSIZ,OFF
CMSEL,S,AIR !对空气区域划分网格
AESIZE,AIR,30 !对所选择的面进行单元尺寸的设置,30为单元尺寸值
TYPE,1 !单元类型为1,三角形自由网格剖分
MSHAPE,1
MSHKEY,0
AMESH,ALL
CSYS,0 !转换为笛卡尔坐标系
LSEL,S,LOC,X,1.5*R0 !选择远场区域的直线
LSEL,A,LOC,Y,1.5*R0
LESIZE,ALL,,,1
CMSEL,S,BIANJIE
TYPE,2
SMRTSIZ,OFF !单元类型为2,四边形映射网格剖分
MSHAPE,0
MSHKEY,1
AMESH,ALL
NSEL,S,LOC,Y,0
CM,COND9,NODE !创建地边界元件
ALLSEL,ALL
CSYS,1
NSEL,S,LOC,X,2*R0 !设置远场标志
SF,ALL,INF
CSYS,0
ALLSEL,ALL
FINISH
/SOLU
CMATRIX,1,'COND',9,0
FINISH
/PREP7
CMSEL,S,COND1
D,ALL,VOLT,V1
CMSEL,S,COND2
CMSEL,A,COND3
CMSEL,A,COND4
CMSEL,A,COND5
CMSEL,A,COND6
CMSEL,A,COND7
CMSEL,A,COND8
CMSEL,A,COND9
D,ALL,VOLT,V0
FINISH
/SOLU
ANTYPE,static,new
EQSLV,JCG
FINISH
邀请回答 我来回答

当前暂无回答

回答可获赠 200金币

没解决?试试专家一对一服务

换一批
    App下载
    技术邻APP
    工程师必备
    • 项目客服
    • 培训客服
    • 平台客服

    TOP