3-D Poisson's Equation using PDE Toolbox on Solidworks part

2 views (last 30 days)
I am attempting to solve a simple paralle plate capacitor problem in order to test my understanding of using the PDE Toolbox in conjunction with Solidworks files. What I am attemtping to set up is a simple Parallel Plate made from a Solidworks part that is a simple cylinder with a top and bottom plate. The top plate has a voltage of 1,000 while the rest of the faces have a voltage of 0. If my interpertation of the PDE Toolbox documentation is correct, this is reflected in the stated boundary conditions and the source-term f.
The issue arises when I attempt to exam the results.nodalSolution to find the electric field and I am still getting only values of 1,000 at the top plate; nothing at the bottom. I believe this is also indicated in my attempt to show the 3-D colormap; its all blue. Attached below is my code and the associated model. Thank you for the help and and if any further clarification is needed just ask.
model = createpde(3)
importGeometry(model,'Parallel Plate Capacitor.STL')
figure(1)
pdegplot(model,'FaceLabels','on','FaceAlpha',0.2) % creates simple 3D model
axis equal
applyBoundaryCondition(model,'dirichlet','Face',1:12,'u', [0;0;0]) % try with no BC
applyBoundaryCondition(model,'dirichlet','Face',14,'u',[1000;1000;1000]) % BC for top plate
applyBoundaryCondition(model,'dirichlet','Face',13,'u',[0;0;0]) % BC for bottom plate
generateMesh(model)
pdeplot3D(model)
f = [1000;1000;1000] % source charge (constant)
CA = specifyCoefficients(model,'m',0,...
'd',0,...
'c',1,...
'a',0,...
'f', f)
results = solvepde(model);
v = results.NodalSolution
pdeplot3D(model,'ColorMapData',v(:,3), 'FaceAlpha', 0.1)

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!