24小时热门版块排行榜

返回列表

lemony1314

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 980.4
散金: 120
帖子: 73
在线: 76.5小时
虫号: 1790071
注册: 2012-05-02
性别: GG
专业: 光学和光电子学

[求助] 电场方面，求解电流时出错已有1人参与

如图，我做的这个，上下设置两电极，中间是测试的材料，小的半圆是空气，用的单元都是plane121，外围的单元是110。
耦合电极上的节点，设置上面电极加交流电压，下面电极接地，求下面电极的电流。
我最后要的是材料的阻抗值，想的是用上面电极的电压除以下面电极的电流得到。
程序如下，N1节点是下面电极上的一个节点。但是最后求得的电流实部虚部都是7.888609052E-31。
请问问题在哪？？
解决之后更有重谢！！！
FINISH
/CLEAR,START
/TITLE,SENSOR
/FILNAME,DETECTION,1
KEYW,MAGELC,1 !指定分析类型-电磁

/PREP7       !前处理器
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!变量定义
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

EW=0.05 !电极宽
EH=0.001 !电极厚

TW=0.100 !复合材料宽
TH=0.050 !复合材料厚

LEN=0.080 !空气区域长度

V1=100
V0=0

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!设定单元类型，建立几何模型
/PNUM,AREA,1
ET,1,121
ET,2,121
ET,3,121
ET,4,110,1,1 !远场类型单元

EMUNIT,MKS
MP,RSVX,1,1E7 !空气电阻率(欧.mm)
MP,RSVX,2,1.75E-5 !铜电极的电阻率(欧.mm)
MP,RSVX,3,1E5 !目标板电阻率(欧.mm)

MP,PERX,1,1 !空气的介电常数(F/mm)
MP,PERX,2,1E-5
MP,PERX,3,5 !目标板介电常数（复合材料）

RECTING,-TW/2,TW/2,EH,TH+EH
RECTING,-EW/2,EW/2,0,EH
RECTING,-EW/2,EW/2,EH+TH,TH+2*EH
CYL4,0,0,LEN,0,,180 !空气区域的建模
CYL4,0,0,1.2*LEN,0,,180 !远场区域的建模

AOVLAP,ALL
NUMCMP,AREA
NUMCMP,LINE

ASEL,S,AREA,,4             !复合材料检测样板
CM,YANGBAN,AREA

ASEL,S,AREA,,1             !电极1
CM,DIANJI1,AREA

ASEL,S,AREA,,3             !电极2
CM,DIANJI2,AREA

ASEL,S,AREA,,5             !空气区域
CM,AIR,AREA

ASEL,S,AREA,,2             !远场
CM,YUANCHANG,AREA

ALLSEL,ALL
ALLSEL,ALL

DOF,CURR

ALLSEL,ALL
!!!!!!!!!!!!!!!!网格化
!!!!!!!!!复合材料试验样板划分

LESIZE,6,,,EW*500
LESIZE,18,,,(TW-EW)*250
LESIZE,19,,,(TW-EW)*250
LESIZE,4,,,EW*500
常数LESIZE,16,,,(TW-EW)*250
LESIZE,17,,,(TW-EW)*250
LESIZE,2,,,TH*500
LESIZE,1,,,TH*500

LSEL,S,LINE,,6
LSEL,A,LINE,,18
LSEL,A,LINE,,19

LCCAT,ALL
ALLSEL,ALL

CMSEL,S,YANGBAN
MSHAPE,0
MSHKEY,0
AATT,3
TYPE,3
AMESH,ALL
ALLSEL,ALL

!!!!!!!!!!!!!电极划分
CMSEL,S,DIANJI1
MSHAPE,0
MSHKEY,1
AATT,2
TYPE,1
AMESH,ALL
ALLSEL,ALL

LSEL,S,LINE,,14
LSEL,A,LINE,,15
LCCAT,ALL
ALLSEL,ALL
CMSEL,S,DIANJI2
MSHAPE,0
MSHKEY,1
AATT,2
TYPE,2
AMESH,ALL
ALLSEL,ALL

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!空气区域划分
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
CSYS,1
LSEL,S,LOC,X,LEN
LESIZE,ALL,,,40
CMSEL,S,AIR
!AESIZE,AIR,30
TYPE,3
AATT,1
MSHAPE,1
MSHKEY,0
AMESH,ALL
ALLSEL,ALL

!!!!!!!!!!!!!!远场区域划分
CSYS,0
LSEL,S,LOC,X,1.1*LEN
LSEL,A,LOC,Y,1.1*LEN
LESIZE,ALL,,,1
CMSEL,S,YUANCHANG
TYPE,4
MSHAPE,0
MSHKEY,1
AMESH,ALL
ALLSEL,ALL

!!!!!!创建远场标志
ALLSEL,ALL
CSYS,1
NSEL,S,LOC,X,1.2*LEN
SF,ALL,INF
D,ALL,VOLT,V0

CSYS,0
ALLSEL,ALL

!!!!耦合电极自由度
ESEL,S,TYPE,,2
NSLE,S
CP,1,CURR,ALL
ALLSEL,ALL

N1=NODE(EW/2,EH,0)
ESEL,S,TYPE,,1
NSLE,S
CP,2,VOLT,ALL
ALLSEL,ALL

CMSEL,S,DIANJI2
DA,ALL,VOLT,V0
ALLSEL,ALL
FINISH

/SOLU
ANTYPE,HARM
CMSEL,S,DIANJI1
DA,ALL,VOLT,V1
ALLSEL,ALL
HARFRQ,10000
!NSUBST,10
KBC,1
OUTPR,BASIC,ALL
OUTRES,ALL,ALL,
SOLVE
SAVE,RESULTT.db
SAVE
FINISH

/POST1

SET,1
*GET,CURREAL,NODE,N1,CURR
SET,1,1,,1
*GET,CURIMAG,NODE,N1,CURR
*STATUS
FINISH

小木虫论坛 - 学术科研第一站.png

回复此楼