September 15, 2023 at 5:10 amYuuhi TakabayashiSubscriber
Maxwell cannot be used due to licensing issues.
I created the program below to simulate the effect of fringe electric fields on parallel plate capacitors in MAPDL, but the simulation does not work well./CLEAR/PREP7! プレートの寸法（μm単位をmm単位に変換）width = 0.1height = 0.02thickness = 0.003separation = 0.01! 3D 8ノード要素の静電解析用に設定ET,1,SOLID122! 無限要素を定義ET,2,111,2,1! 材料特性（真空の相対誘電率）MP,PERX,1,8.854e-12! 最初のプレートを作成BLOCK,0,width,0,height,0,thickness! 二番目のプレートを作成BLOCK,0,width,0,height,separation,separation+thickness! メッシュ設定ESIZE,0.005 ! 要素サイズを設定! すべてのボリュームにメッシュを適用ALLSEL,ALLVMESH,ALL/SOLU! 静電解析を指定ANTYPE,STAT! 境界条件ALLSELCSYS,0! 最初のプレートの上面に0Vを適用ASEL,S,LOC,Z,thicknessNSLA,S,1D,ALL,VOLT,0! 二番目のプレートの下面に1Vを適用ALLSELASEL,S,LOC,Z,separationNSLA,S,1D,ALL,VOLT,1! 解を求めるSOLVEFINISH/POST1PLNSOL,VOLT ! ノード解（電圧）をプロット! 電場エネルギーを抽出して容量を計算/OUT, CAPACITANCE_CALCULATION, , TXT*GET, W, ENER, , STAT, ENEV = 1C = 2*W/(V**2)*STATUS, C/OUTPlease tell me what needs to be improved.thank you
September 15, 2023 at 11:14 amAshish KumarForum Moderator
Dear everyone, Maxwell cannot be used due to licensing issues. I created the program below to simulate the effect of fringe electric fields on parallel plate capacitors in MAPDL, but the simulation does not work well. /CLEAR /PREP7 ! Plate dimensions (convert from μm to mm) width = 0.1 height = 0.02 thickness = 0.003 separation = 0.01 ! Configured for electrostatic analysis of 3D 8-node elements ET,1,SOLID122 ! Define infinite element ET,2,111,2,1 ! Material properties (relative dielectric constant in vacuum) MP,PERX,1,8.854e-12 ! Create your first plate BLOCK,0,width,0,height,0,thickness !Create second plate BLOCK,0,width,0,height,separation,separation+thickness ! Mesh settings ESIZE,0.005 ! Set element size !Apply mesh to all volumes ALLSEL,ALL VMESH,ALL /SOLU ! Specify electrostatic analysis ANTYPE,STAT ! boundary condition ALLSEL CSYS,0 !Apply 0V to the top of the first plate ASEL,S,LOC,Z,thickness NSLA,S,1 D,ALL,VOLT,0 !Apply 1V to the bottom of the second plate ALLSEL ASEL,S,LOC,Z,separation NSLA,S,1 D,ALL,VOLT,1 ! Find the solution SOLVE FINISH /POST1 PLNSOL,VOLT ! Plot node solution (voltage) ! Calculate capacity by extracting electric field energy /OUT, CAPACITANCE_CALCULATION, , TXT *GET, W, ENER, , STAT, ENE V = 1 C = 2*W/(V**2) *STATUS, C /OUT Please tell me what needs to be improved. thank you
September 15, 2023 at 5:45 pmBill BulatAnsys Employee
Please try the listing below. In electrostatic analyses, one models/meshes the domain enclosing the electrodes (not the electrodes themselves). The electrodes are usually represented by voids in a mesh of the domain that encloses them. I didn't bother to add the INFIN111 elements (it's a little tedious - pardon me for getting cheap on you). You can increase the value of parameter d_dmn in successive runs until the value of the calculated capacitance stops changing significantly.
If memory serves, INFIN111 will not account for additional electrostatic energy of the infinite domain that extends beyond the boundaries of the SOLID122/123 mesh. So if you do use them, my second approach for assessing the capacitance would be, I believe, the more accurate way (divide the charge reaction on one of the electrodes by the applied voltage difference).
Billfini/cle/vie,1,1,1,1/sys,del file*.png/PREP7C*** Plate dimensions (we'll assume this is in meters)eps0=8.854e-12 ! FREE SPACE PERMITTIVITYwidth = 0.1height = 0.02thickness = 0.003separation = 0.002 !0.010d_dmn = 0.010 ! DEPTH SURROUNDING DOMAIN EXTENDS BEYOND ELECTRODESd_inf = 0.005 ! DEPTH OF INFINITE BOUNDARY ELEMENT MESHesz = 0.003 ! MESH SIZEV = 1 ! APPLIED VOLTAGEC*** Configured for electrostatic analysis of 3D 8-node elementsET,1,SOLID123C*** Define infinite elementET,2,111,2,1C*** Material properties (relative dielectric constant in vacuum)!MP,PERX,1,8.854e-12 ! THIS IS ABSOLUTE PERMITTIVITY IN MKS, WHICH WORKS OK, BUT...MP,PERX,1,1 ! WE GENERALLY USE RELATIVE PERMITTIVITY INSTEADC*** Create your first plateBLOCK,0,width,0,height,0,-thicknesscm,electrode1_a,areaC*** Create second plateasel,noneBLOCK,0,width,0,height,separation,separation+thicknesscm,electrode2_a,areaC*** SURROUNDING DOMAINallscm,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,zvsel,nonebloc,xmin-d_dmn,xmax+d_dmn,ymin-d_dmn,ymax+d_dmn,zmin-d_dmn,zmax+d_dmncm,scrap2_v,voluallsvsbv,scrap2_v,scrap1_vC*** MESHesiz,eszvmes,allC*** ELECTRODE VOLTAGEda,electrode1_a,voltda,electrode2_a,volt,VfiniC*** SOLVE/soluallssavesolvfiniC*** POST PROCESS/post1set ! READ RESULTS INTO MADL DATABASEetab,w,sene ! BASED ON ENERGYssum*get,Wnet,ssum,,item,wC1 = 2*Wnet/(V**2)nsel,s,d,volt ! BASED ON CHARGE & VOLTAGEfsum*get,q,fsum,,item,chrgc2 = q/vC_nominal=eps0*(height*width)/separation ! NOMINAL FRINGE FREE VALUE/tla,-0.35,0.90,Capacitance (ANSYS, energy) = %C1% Fd ! ANNOTATION/tla,-0.35,0.85,Capacitance (ANSYS, Q/V) = %C2% Fd/tla,-0.35,0.80,Capacitance (nominal fringe-free) = %C_nominal% Fdalls ! CAPPED SECTION PLOTSwpro,,,90wpof,,,width/2/type,1,7/cpl,1plns,volt/sho,png $plns,volt $/sho,close $/wait,2plve,ef,,,,vect,,on/sho,png $plve,ef,,,,vect,,on $/sho,close $/wait,2/eof
September 17, 2023 at 5:19 pmYuuhi TakabayashiSubscriber
Thank you for your answer.
I tried it right away and was able to successfully analyze it. It never occurred to me to mesh the area instead of the electrode.
I have an additional question: What kind of command operation should I use to extract the electric field leaking from the electrode?
Where is the capacitance worth when there is a fringe?
September 17, 2023 at 3:15 pmYuuhi TakabayashiSubscriber
September 17, 2023 at 4:30 pmYuuhi TakabayashiSubscriber
Where is the capacitance worth when there is a fringe?
- You must be logged in to reply to this topic.
Boost Ansys Fluent Simulations with AWS
Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
- Solver Pivot Warning in Beam Element Model
- Saving & sharing of Working project files in .wbpz format
- Understanding Force Convergence Solution Output
- User manual
- An Unknown error occurred during solution. Check the Solver Output…..
- What is the difference between bonded contact region and fixed joint
- The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
- whether have the difference between using contact and target bodies
- Defining rigid body and contact
- Colors and Mesh Display
© 2023 Copyright ANSYS, Inc. All rights reserved.