Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.
Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.
Varying capacitance due to thickness reduction
Posted 19 apr 2011, 07:34 GMT-4 Low-Frequency Electromagnetics, Results & Visualization, Structural Mechanics Version 4.0a 12 Replies
Please login with a confirmed email address before reporting spam
I am simulating a capacitor with a compliant dielectric and compliant electrodes in COMSOL 3.5a, using the Structural mechanics and electrostatics modes. I am using 2D axial symmetry
I wanted to know whether it is possible to plot the change in capacitance which occurs due to structural deformation of the capacitor, e.g. increase in surface area and consequent reduction in thickness.
Any suggestions would be helpful.
Thank you,
Aseem
Please login with a confirmed email address before reporting spam
check the lumped port/terminal entries, and check the doc, there is an example I beleive in the Model library
--
Good luck
Ivar
Please login with a confirmed email address before reporting spam
C11_emes variable plots the capacitance only when the voltage is being applied.
I was wondering if I could also plot capacitance which occurs because of structural deformation as well.
Because, capacitance can also be calculated as C = (epsilon)*(Area_dielectric / Thickness_dielectric).
Thanks
-Aseem
Please login with a confirmed email address before reporting spam
try it out you have the integration operator there, you can also calculate via the energy
--
Good luck
Ivar
Please login with a confirmed email address before reporting spam
No. The energy method does not work either.
Is it possible to track the change in distance between two points or nodes in COMSOL as the geometry gets deformed? Then I might be able to use the above mentioned formula and calculate change in capacitance as a function of the thickness.
Any ideas?
Thanks
Aseem
Please login with a confirmed email address before reporting spam
i think this can be done with parametric sweep in 4.0 and 4.0a but i dont know weather the provision is available in 3.5a. you can declare any quantity as a parameter and under study you can mention the range of values this parameter can take. this variation can be discrete or continuous. if it is discrete u should have those values in a text file and you should load them. to view results you should define boundary probe or domain probe depending on weather you want results on a boundary or in a domain.
regards
deepak
Please login with a confirmed email address before reporting spam
I too have the same issue of deriving the varying capacitance in capacitive pressure sensor. Thing is i am not sure as to how to enter the formulae. Where can i get the syntax of such expressions.please help.if i get that i can work out something.
Thanks a lot
Please login with a confirmed email address before reporting spam
mpheval(model,’w’) %for gradient of deformation in zZ direction for example
But with the capacitance, let say es.C11 the mpheval(model,’es.C11’) doesn’t work.
Is there any other way to obtain this value?
Any hint or discuss will be very useful for me. Thank you so much
LOAN
Please login with a confirmed email address before reporting spam
Try by applying a small potential difference across the electrodes and then plotting the C11_emes component.
Worked for me.
I hope this helps.
Please login with a confirmed email address before reporting spam
My problem is similar to RV Coll's. However, what i'm looking for is deriving capacitance by the variation in the gap of the dielectric b/w the plates.
I'm applying pressure on the top plate>plate deforms by displacing into the gap>gap spacing reduces>capacitance changes because of the formula on paper
C=epsilon*A/g
However, to integrate capacitance with the mechanics part i'm facing difficulties.
I'm attaching my model. In the 'variables' field under 'Definitions' is the formula
epsilon0/(2[um]+w)
Now, 2[um] is the gap between the plates and epsilon0 is defined in 'parameters'. What i'm not sure is the 'w'. w must be the deflection of the top plate. I do not know if it is being derived from the results 'solid.disp'. Nor am i getting proper results when i change the gap (like 10[um] instead of 2[um])
If some one could help me figure that out....i shall,forever, be grateful.
Cheers!
Abhi
Attachments:
Please login with a confirmed email address before reporting spam
a few comments:
In V4 there are many physical constants aready defined, epsilon0 among others, but you should not forget to add the suffix "_const" to the name, this tag tells COMSOL to use its internal values, hence use "epsilon0_const"
u,v,w are the 3D dependent variables solved for in "structural", so indeed w is the displacement along "z".
solid.disp = sqrt(u^2+v^2+w^2) the total displacement vector norm (check it in the equations view, turn on via preferences)
So for me the formula is correct, and too get to C = epsilonr*epsilon0*Area/Distance = 1.854E-15 F from my estimation for a 1.5[um] separation epsilonr=1and 10[um] radius disc capacitor with "0" deformation or pressure applied.
So what is so wrong ?
Do not forget that the displacements/deformations are scaled to be "visible" as the membrane is only deflecting at most 42 nm for 5 bars pressure difference
--
Good luck
Ivar
Please login with a confirmed email address before reporting spam
Please login with a confirmed email address before reporting spam
V and p are they dependent variables or one defined with a parameter ?
many cases are possible, if you define V as a parameter (or a parametric or "continuation" solver sequence) and caclulate the resulting p than you can plot p versus v in the Derived values section which creates a table and plot out the table in a line graph
--
Good luck
Ivar
Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.
Suggested Content
- KNOWLEDGE BASE Solving electromagnetic heating models with high-frequency loads
- FORUM Varying thickness in 2D cantilever model
- BLOG Changing the Dimensions of a Model Using Shape Optimization
- KNOWLEDGE BASE Solving Wave-Type Problems with Step Changes in the Loads
- BLOG How to Model Electrochemical Resistance and Capacitance