微带中电压和电场分布的有限元仿真（Matlab实现）

荧光Matlab

已于 2024-10-20 20:18:11 修改

阅读量446

点赞数 5

文章标签： matlab 开发语言

于 2024-10-20 19:59:32 首次发布

本文链接：https://blog.csdn.net/Matlab88888/article/details/143097632

版权

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

💥1 概述

该程序的目标是求解微带线内部区域的稳态电压分布。然后，利用该图获得该区域中的电场分布，并使用矢量场图（Quiver 图）。微带线在坐标系中由 x 表示的 0 到 7 个单位以及 y 表示的 0 到 2 个单位组成。该区域的所有边上都有四个电容器极板。该区域由在 y 处从 0 到 1 个单位的相对介电常数为 1 的材料和在 y 处从 1 到 2 个单位的相对介电常数为 9.6 的材料组成。微带线沿 z 方向无限延伸。区域的边被赋予狄利克雷边界条件。水平微带线放置在区域的中心。微带线用 1 伏特的电压激励，周围所有的极板都保持在 0 伏特。

📚2 运行结果

主函数部分代码：

clear all;
clc;

%Dimensions of the simulation grid in x (xdim) and y (ydim) directions
xdim=7;
ydim=2;

%Area of square elements in the grid
area_sq_elem=0.0625^2;

%Step size in the grid in cartesian directions viz. side of the square
%element
step_size=area_sq_elem^0.5;

%No of square elements
sq_elems=(xdim*ydim)/area_sq_elem;

%Dimension in x direction in steps
dim1=xdim/step_size;

%Dimension in y direction in steps
dim2=ydim/step_size;

%No of triangular elements
no_of_elems=sq_elems*2;

%Total no of points for voltage distribution
total_no_pts=(dim1+1)*(dim2+1);

%grid_pts matrix consisting of numbering of all the points in the domain in
%a page read fashion from left to right and then top to bottom along with
%the x and y coordinates of the points in the adjacent columns of the
%matrix
grid_pts=zeros(total_no_pts,3);

%look_up_table matrix consisting of numbering of all the triangular
%elements in a page read fashion just like before along with the cardinal
%numbers of the three points making the triangular elements (obtained from
%the grid_pts matrix) in the adjacent columns of the matrix
look_up_table=zeros(no_of_elems,4);

%Populating the first columns of the above two matrices with serial numbers
i1=1:1:no_of_elems;
look_up_table(:,1)=i1';
i1=1:1:total_no_pts;
grid_pts(:,1)=i1';



% Polpulating the x coordinate and y coordinate columns of the grid_pts
% matrix with the appropriate values from the knowledge of the grid
% geometry
for i1=1:1:total_no_pts
    if rem(i1,dim1+1)==0
        grid_pts(i1,3)=((i1/(dim1+1))-1)*step_size;
    else
        grid_pts(i1,3)=(floor(i1/(dim1+1)))*step_size;
    end
    if rem(i1,dim1+1)==0
        grid_pts(i1,2)=dim1*step_size;
    else
        grid_pts(i1,2)=(i1-(dim1+1)*floor(i1/(dim1+1))-1)*step_size;
    end
end