Learning Forum

[ty_jpn]

Dear everyone

 

When I run the following script, the analysis ends at the 8th loop．

I have the following outputs, what can I do about them?

*** WARNING *** CP = 2110.030 TIME= 05:22:15

Output entity is identical to the first input entity for VSBV command. No new entities created.

 *** ERROR *** CP = 2333.973 TIME= 05:25:55

There is at least 1 small equation solver pivot term (e.g., at the VOLT degree of freedom of node 15727). Please check for an insufficiently constrained model.

------------------------------------------------------------------

fini 
/cle 
 
C*******************************************************************************
C*** ループ外のパラメータ
C*******************************************************************************
height = 20*1e-6
width = 100*1e-6
separation = 10*1e-6

thickness_0 = (1*1e-6) ! 初期の電極の幅
d_thickness = 1e-6       ! 電極の幅の増分
n_thickness = 200          ! 調査する電極の幅のケース数

*dim,C_results,,n_thickness,4 ! キャパシタンス結果を保存する配列

eps0=8.854e-12       ! 真空の誘電率
d_dmn = 100*1e-6      ! 電極を超えるドメインの深さ
d_inf = 5*1e-6       ! 無限バウンダリ要素メッシュの深さ
esz = 10*1e-6         ! メッシュのサイズ
V = 3             ! 適用される電圧

C*******************************************************************************
C*** DOループ
C*******************************************************************************
*do,i,0,n_thickness-1    ! 電極の幅をループ

pars,all              ! file.parmにすべてのパラメータを保存

fini                  ! MAPDLデータベースをクリア（これによりパラメータも削除されます）
/cle

parr                  ! file.parmからMAPDLデータベースにパラメータを復元

thickness = thickness_0 + i*d_thickness ! i回目のループでの電極の幅

/title,THICKNESS = %thickness%

/vie,1,1,1,1
/sys,del file*.png

/PREP7

C*** 電気静止解析のための3D 8ノード要素を設定
ET,1,SOLID123

C*** 無限要素の定義
ET,2,111,2,1

C*** 材料特性（真空中の相対誘電率）
MP,PERX,1,1 ! 相対誘電率

C*** 最初のプレートを作成
BLOCK,0,width,0,height,0,-thickness
cm,electrode1_a,area

C*** 二つ目のプレートを作成
asel,none
BLOCK,0,width,0,height,separation,separation+thickness
cm,electrode2_a,area

C*** 周囲のドメイン
alls
cm,scrap1_v,volu
*get,xmin,kp,,mnloc,x
*get,xmax,kp,,mxloc,x
*get,ymin,kp,,mnloc,y
*get,ymax,kp,,mxloc,y
*get,zmin,kp,,mnloc,z
*get,zmax,kp,,mxloc,z
vsel,none
bloc,xmin-d_dmn,xmax+d_dmn,ymin-d_dmn,ymax+d_dmn,zmin-d_dmn,zmax+d_dmn
cm,scrap2_v,volu
alls
vsbv,scrap2_v,scrap1_v

C*** メッシュ
esiz,esz
vmes,all

C*** 電極の電圧
da,electrode1_a,volt
da,electrode2_a,volt,V

fini

C*** 解を求める
/solu
alls
save
solv
fini

C*** ポストプロセス
/post1
set ! MADLデータベースに結果を読み込む

etab,w,sene ! エネルギーに基づく
ssum
*get,Wnet,ssum,,item,w
C1 = 2*Wnet/(V**2)

nsel,s,d,volt ! 電荷と電圧に基づく
fsum
*get,q,fsum,,item,chrg
c2 = q/v

C_nominal=eps0*(height*width)/separation ! ノミナルフリンジフリーバリュー

/tla,-0.35,0.90,Capacitance (ANSYS, energy) = %C1% Fd ! 注釈
/tla,-0.35,0.85,Capacitance (ANSYS, Q/V) = %C2% Fd
/tla,-0.35,0.80,Capacitance (nominal fringe-free) = %C_nominal% Fd

alls ! キャップ付きセクションプロット
wpro,,,90
wpof,,,width/2
/type,1,7
/cpl,1

plns,volt
/sho,png $plns,volt $/sho,close $/wait,2
plve,ef,,,,vect,,on
/sho,png $plve,ef,,,,vect,,on $/sho,close $/wait,2

C_results(i+1,1)=thickness ! C_results配列に結果を保存
C_results(i+1,2)=C_nominal
C_results(i+1,3)=C1
C_results(i+1,4)=C2

*enddo  ! DOループの終了

C*******************************************************************************
C*** C_results配列をresults.txtテキストファイルに書き込む
C*******************************************************************************
*cfopen,results,txt
*vwrite
('           thickness             C_nominal       C_ANSYS (Energy)  C_ANSYS (Terminal Conditions)')
*vwrite,C_results(1,1),C_results(1,2),C_results(1,3),C_results(1,4)
(4X,E16.8,4X,E16.8,4X,E16.8,4X,E16.8)
*cfclose

pars,all
save

/eof

--------------------------------------------------

[ty_jpn]

Errors only occur in simulations when the thickness is 9 microns.