二维场的数值计算实验报告

上传人:ba****u6 文档编号:168579622 上传时间:2022-11-11 格式:DOCX 页数:5 大小:168.52KB
收藏 版权申诉 举报 下载
二维场的数值计算实验报告_第1页
第1页 / 共5页
二维场的数值计算实验报告_第2页
第2页 / 共5页
二维场的数值计算实验报告_第3页
第3页 / 共5页
资源描述:

《二维场的数值计算实验报告》由会员分享,可在线阅读,更多相关《二维场的数值计算实验报告(5页珍藏版)》请在装配图网上搜索。

1、二维场的数值计算实验报告一实验目的通过采用两种不同的计算方法,上机计算二维场的电位分布,了解数值法与 解析法的优缺点和区别,掌握超松弛迭代方法。二实验内容1根据分离变量法求得的二维场的用级数表示的解,编写程序,上机计算二 维场中各点的电位;2根据有限差分法原理,编写超松弛迭代法程序,上机计算二维场中各点的 电位。三实验题目如图所示横截面为矩形的无限长直导体槽,内填空气。已I IIT IT(0, y ) = V (a, y ) = V(x,0 ) = V (x, b ) = V34o1. 根据分离变量法求得的二维场的用级数表示的解,编写程序, 上机计算二维场中各点的电位;2. 根据有限差分法原理

2、,编写超松弛迭代法程序,上机计算二维场中各点的电 位。四.实验设计过程1 分离变量法:在利用解析法求解边界问题时,一般都是利用麦克斯韦方程,在已知边界 值的情况下,结合条件球场电位函数。结果是一个复杂的级数函数表达式。虽然 不能具体体现场的分布,但可以清楚地算出场中的每一个数值,非常的精确。根 据上述题目,电位分布与z轴无关,并且槽内点电荷密度P二0,由此可以知道 其满足的拉普拉斯方程及边界条件如下:d 2 d 2+ dx 2dy 2(x, y)二 0 ,y=0(x, y )= Vy=b3通过对第一个方程进行分离变量法得到方程的通解为:山,y )=(C sin k x + C cos k x)

3、D sinh k y + D cosh k y)1x2x1y2y再将后四个方程边界条件代入求得电位函数:n为奇数1n兀.n兀smhy sm -n sinh 巴 b a aa2.有限差分法在题目中,有限差分法与分离变量法的最主要区别是对于(1)中的拉普拉斯方程进行 近似计算,将场域划分为网格节点,将网格节点电位作为未知数的差分近似替代 该点的偏导数,从而把待求得偏微分方程的边值问题转化为一组相应的差分方程 问题。最后利用超松弛迭代法计算出数值。由于点电荷密度为0则拉普拉斯方程的有限差分方程形式为具体步骤如下:1. 网格设置2. 区域设置3. 设置边界条件4. 应用模式设置选择(静电学)5方程参数

4、设定五实验程序先假设V、V、V、V的数值分别为0, 0, 100V, 0, a和b的数值为10和5cm。12341分离变量法: cleara=0;X,Y=meshgrid(0:0.05:2, 2:0.05:4);for n=1b=4*100/pi; a=b*sinh(n*pi*Y/2).*sin(n*pi*X/2)./n*sinh(n*pi*4/2); a=a+a;enda=a;figurecontour(a,20);2有限差分法:clearhx=40;hy=40;vl=ones(hy, hx);vl(hy,:)=ones(l, hx)*100;vl(l,:)=zeros(l, hx);for

5、 i=l:hyvl(i, 1)=0;v1(i, hx)=0;endv2=v1; maxt=1; t=0;k=0while(maxt1e-9)k=k+1;maxt=0;for i=2: hy-1for j=2:hx-1v2(i, j)=(v1(i, j+1)+v1(i+1, j)+v2(i-1, j)+v2(i, j-1)/4; t=abs(v2(i, j)-v1(i, j);if(tmaxt)maxt=t;endendendv1=v2;endcontour(v2, 20)hold onx=1:1:hx;y=1:1:hy;Gx, Gy=gradient(- v2, 10, 10); quiver

6、(xx, yy, Gx, Gy, 8, r) hold off六实验结果1分离变量法实验图形七实验结果分析比较和总结对比上述两种方法求解思路以及所得到的实验图形可知,虽然它们的解题思路不同,但是在通过程序进行画图后,他们的电位分布图是一样的,都无一例外 地真实准确地展示了静态场的电位分布。但是,他们也有不同的地方,分离变量 法求得的是函数,然后通过离散来展示场的分布;而有限差分法得到的解本身是 离散的,通过编程不断迭代来得到数值以接近分离变量法得到的函数。求解静态场分布的方法是多样的,但最终借助于MATLAB这个程序使我们 的结果趋于一致,可见我们求解问题并不一定只局限于一个方法,多样的方法最 终都会走向一个结果,我们需要综合这些方法,在其中摸索发现,以期有更好的 、p r r.前进。

展开阅读全文
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

copyright@ 2023-2025  zhuangpeitu.com 装配图网版权所有   联系电话:18123376007

备案号:ICP2024067431-1 川公网安备51140202000466号


本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!