如圖,我做的這個(gè),上下設(shè)置兩電極,中間是測試的材料,小的半圓是空氣,用的單元都是plane121,外圍的單元是110。
耦合電極上的節(jié)點(diǎn),設(shè)置上面電極加交流電壓,下面電極接地,求下面電極的電流。
我最后要的是材料的阻抗值,想的是用上面電極的電壓除以下面電極的電流得到。
程序如下,N1節(jié)點(diǎn)是下面電極上的一個(gè)節(jié)點(diǎn)。但是最后求得的電流實(shí)部虛部都是7.888609052E-31。
請問問題在哪??
解決之后更有重謝。!
FINISH
/CLEAR,START
/TITLE,SENSOR
/FILNAME,DETECTION,1
KEYW,MAGELC,1 !指定分析類型-電磁
/PREP7 !前處理器
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!變量定義
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
EW=0.05 !電極寬
EH=0.001 !電極厚
TW=0.100 !復(fù)合材料寬
TH=0.050 !復(fù)合材料厚
LEN=0.080 !空氣區(qū)域長度
V1=100
V0=0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!設(shè)定單元類型,建立幾何模型
/PNUM,AREA,1
ET,1,121
ET,2,121
ET,3,121
ET,4,110,1,1 !遠(yuǎn)場類型單元
EMUNIT,MKS
MP,RSVX,1,1E7 !空氣電阻率(歐.mm)
MP,RSVX,2,1.75E-5 !銅電極的電阻率(歐.mm)
MP,RSVX,3,1E5 !目標(biāo)板電阻率(歐.mm)
MP,PERX,1,1 !空氣的介電常數(shù)(F/mm)
MP,PERX,2,1E-5
MP,PERX,3,5 !目標(biāo)板介電常數(shù)(復(fù)合材料)
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 !空氣區(qū)域的建模
CYL4,0,0,1.2*LEN,0,,180 !遠(yuǎn)場區(qū)域的建模
AOVLAP,ALL
NUMCMP,AREA
NUMCMP,LINE
ASEL,S,AREA,,4 !復(fù)合材料檢測樣板
CM,YANGBAN,AREA
ASEL,S,AREA,,1 !電極1
CM,DIANJI1,AREA
ASEL,S,AREA,,3 !電極2
CM,DIANJI2,AREA
ASEL,S,AREA,,5 !空氣區(qū)域
CM,AIR,AREA
ASEL,S,AREA,,2 !遠(yuǎn)場
CM,YUANCHANG,AREA
ALLSEL,ALL
ALLSEL,ALL
DOF,CURR
ALLSEL,ALL
!!!!!!!!!!!!!!!!網(wǎng)格化
!!!!!!!!!復(fù)合材料試驗(yàn)樣板劃分
LESIZE,6,,,EW*500
LESIZE,18,,,(TW-EW)*250
LESIZE,19,,,(TW-EW)*250
LESIZE,4,,,EW*500
常數(shù)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
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!空氣區(qū)域劃分
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
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
!!!!!!!!!!!!!!遠(yuǎn)場區(qū)域劃分
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
!!!!!!創(chuàng)建遠(yuǎn)場標(biāo)志
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
![電場方面,求解電流時(shí)出錯(cuò)]()
小木蟲論壇 - 學(xué)術(shù)科研第一站.png |