几个电磁场问题的计算机模拟.ppt

电力线之美 P 王泉 问题的提出 电场线是美的 问题的提出 电场线是美的 电场线是怎么绘制出来的？ 电场线方程（1） 电场线方程（1） 电场线方程（1） 电场线方程（1） 电场线方程（1） 电场线方程（2） 电场线方程（2） 电场线方程（2） 电场线方程（2） 程序流程图(Start) 程序流程图(End) 一些电场线（线性四极子） 一些电场线 一些电场线 一些电场线 一些电场线 一些电场线（平面电三极子） 一些电场线 一些电场线（电偶极子） 一些电场线（平面电四极子） 一些电场线 一些电场线 一些电场线（任意情形） 等势曲线 一些等势曲线 一些等势曲线 一些等势曲线 一些等势曲线 总结 数值解法的强大威力(“笨” 但是 有效) 谢谢! Generated By Mathematica 4.1 令： PositiveCharge = {2 + 1i, 2}; NegativeCharge = {-2 + 0i, -1}; Generated By Mathematica 4.1 令： PositiveCharge = {2 + 0i, 1}; NegativeCharge = {-2 + 0i, 3}; * * 电场线的定义：电场所在空间中的一系列想象的曲线，曲线上每一点的切线方向都与该点的场强方向一致。 定义 数学化 抽象 电荷 (Q1,Q2…Qn) 库仑定律 叠加原理 空间各点的场强（大小，方向） 孤立点电荷的电场线方程 ： Solve it! 复杂一些的情况 ： Solve it ? ？？？ 数值方法（欧拉折线法） 基本思想： 寻求精确解的折线逼近 流程： 点A (x,y) 以A为起始点，A的切线方向为方向，以s为步长作线段 A的切线方向 例： 线性电四极子 （位于 N(-L,0), O(0,0), P(L,0)，带电量分别为 + Q,-2Q, + Q的电荷系统 ） 思想：欧拉折线法 初始条件：在距离电荷很近的地方取点 电荷初始化 转动角θ= 0.01 开始 测试点A=电荷P; φ = θ;点集Points[]={P的坐标};测试点的初位移l=0.01 上一测试点PreviousPoint=测试点A 测试点A=(PreviousPointX+Step*Cosφ ,PreviousPointY+Step*Sinφ ) 利用库仑定律计算测试点A的受力F(用复数表示) φ =F的幅角; Points[]+={A的坐标} 侧试点A是否满足条件？ YES NO NO YES θ 2π ? Points[]中的点以线段连接成为场线; Graphic[]+=场线 注： 1.测试点A带一个单位正电荷 2.φ 代表测试点A要走的方向 3.Points[]为存储A所经过的点的数组 4.Graphic[]为存储场线的数组 5.测试点A满足的条件:不能离电荷太近 不能离远点太远 θ += π/12 结束 打印出Graphic[]中的所有场线 令： PositiveCharge = {2 + 0i, 1}; NegativeCharge = {-2 + 0i, 1}; MiddleCharge = {0 + 0i, -2}; {坐标，带电量} Generated By Mathematica 4.1 令： PositiveCharge = {2 + 0i, 1}; NegativeCharge = {-2 + 0i, 1}; MiddleCharge = {0 + 0i, 0}; Generated By Mathematica 4.1 令： PositiveCharge = {2 + 0i, 1}; NegativeCharge = {-2 + 0i, -1}; MiddleCharge = {0 + 0i, 0}; Generated By Mathematica 4.1 令： PositiveCharge = {2 + 0i, 3}; NegativeCharge = {-2 + 0i, -1}; MiddleCharge = {0 + 0i, 0}; Generated By Mathematica 4.1 令： PositiveCharge = {2 + 3i, 5}; NegativeCharge = {-2 + 0i, -1}; MiddleCharge = {0 + 0i, 0}; Generated By Mathematica 4.1 令： PositiveCharge = {2 + 0i, 1}; NegativeCharge = {-2 + 0i, 1}; MiddleCharge = {0 + 3.4641i, -2 };