工程电磁场教案国家精品课华北电力学院崔翔第2章第二部分.docx
- 文档编号:10097199
- 上传时间:2023-02-08
- 格式:DOCX
- 页数:34
- 大小:262.85KB
工程电磁场教案国家精品课华北电力学院崔翔第2章第二部分.docx
《工程电磁场教案国家精品课华北电力学院崔翔第2章第二部分.docx》由会员分享,可在线阅读,更多相关《工程电磁场教案国家精品课华北电力学院崔翔第2章第二部分.docx(34页珍藏版)》请在冰豆网上搜索。
工程电磁场教案国家精品课华北电力学院崔翔第2章第二部分
2.4电介质中的电场
1.电位移矢量
由高斯定理,得
整理得▽∙(ε0E+P)=ρ
定义电位移矢量:
D=ε0E+P=ε0(1+χe)E=εE
其中,ε=ε0(1+χe)=εrε0,εr=ε/ε0=(1+χe)
2.介电常数
上式分别给出了介质的介电常数和相对介电常数。
从而电介质中电场问题可简洁地归结为场量D、E或位函数ϕ的定解问题。
例1:
同轴电缆其长度L远大于截面半径,已知内、外导体半径分别为a和b。
其间充满介电常数为ε的介质,将该电缆的内外导体与直流电压源U0相联接。
试求:
(1)介质中的电场强度E;
(2)介质中Emax位于哪里?
其值多大?
[解]:
(1)设内、外导体沿轴线方向线电荷密度分别为+τ和-τ。
由应用高斯定理,得
即
所以
(a<ρ
由因为
则
得
(a<ρ
(2)最大场强位于内导体表面(ρ=a),其值为
3.边界条件
介质分界面上的边界条件:
跨越分界面的一狭小的矩形回路l如图所示,且令∆l2→0而∆l1足够地短。
求电场强度在l上的环量,有
即E1t=E2t或en⨯(E2-E1)=0
上式表明,在介质分界面上电场强度的切向分量是连续的。
跨越分界面的一个扁平圆柱体S如图所示,令两个底面∆S足够小且平行于分界面,圆柱面高度∆l→0。
求电位移矢量在圆柱面的通量,有
式中分界面上法线方向单位矢量en规定为由介质1指向介质2,σ是分界面上可能存在的自由电荷面密度。
从而得
D2n-D1n=σ或en⋅(D2-D1)=σ
一般两种介质分界面上不存在自由电荷(σ=0),此时有
D1n=D2n或en⋅(D2-D1)=0
上式表明,在介质分界面上电位移矢量的法向分量是连续的。
对于两种线性且各向同性介质,应用上述边界条件,得
E1sinα1=E2sinα2,ε1E1cosα1=ε2E2cosα2
两式相除,得
上式综合表述了场量在介质分界面上遵循的物理规律,称为静电场的折射定律。
导体表面上的边界条件:
设导体为媒质1、导体外介质为媒质2,并考虑到导体内部电场强度和电位移矢量均为零且其电荷只能分布在导体表面,得
E1t=E2t=0,D2n-D1n=D2n=σ
式中,σ是导体表面的电荷面密度。
上式说明在导体表面相邻处的电场强度E和电位移D都垂直于导体表面,且电位移的量值等于该点的电荷面密度(需注意en是导体表面的外法线单位矢量)。
一般写为
Et=0或en⨯E=0;Dn=σ或en⋅D=σ
4.边界条件的电位表达
介质分界面:
由于介质分界面上E1t=E2t,显然可以得出
ϕ1=ϕ2
即电位连续在介质分界面上是连续的。
又由于D2n-D1n=σ和
最后可以得出,边界条件的电位表示为
ϕ1=ϕ2,
导体表面上的边界条件:
ϕ=C,
式中,C是由所论静电场导体系统决定的常数。
例2:
图示平行板电容器,其极板间介质由两种绝缘材料组成,介质的分界面与极板平行。
设电容器外施电压为U0,试求:
(1)两绝缘材料中的电场强度;
(2)极板上的电荷面密度。
[解]:
(1)在电压U0下,并应用分界面的边界条件,得
,
(2)极板A上的电荷面密度为
极板B上的电荷面密度为σ'=-D2n=-ε2E2=-σ
讨论:
本例中,设εr2>εr1,则E1>E2。
在实际中,如果因制造工艺上的不完善性,使极板与绝缘材料间留有一空气层,设绝缘材料的相对介电常数为εr2,则空气层中电场强度E1将为绝缘材料中电场强度E2的εr2倍,这很容易由于空气层被击穿而导致电容器的损坏。
2.5边值问题
1.泛定方程
由∇∙D=ρ、D=εE和E=-∇ϕ,得
∇∙D=∇∙εE=-∇∙ε∇ϕ=ρ=〉∇∙ε∇ϕ=-ρ
对于均匀介质ε为常数,得
∇2ϕ=-ρ/ε
上式称为电位ϕ的泊松方程,式中
,称为拉普拉斯算子,在直角坐标系中∇2ϕ
对于场中无自由电荷分布(ρ=0)的区域,泊松方程退化为拉普拉斯方程,即
∇2ϕ=0
2.边界条件
第一类边界条件(狄利赫莱条件):
场域边界S上的电位分布已知,即
式中rb为相应边界点的位置矢量。
它与泛定方程构成第一类边值问题。
第二类边界条件(纽曼条件):
场域边界S上电位的法向导数分布已知,即
当f2(rb)取零时,称为第二类齐次边界条件。
它与泛定方程构成第二类边值问题。
第三类边界条件(混合条件):
场域边界S上电位及其法向导数的线性组合已知,即
它与泛定方程构成第三类边值问题。
无限远边界条件:
对于电荷分布在有限域的无边界电场问题,在无限远处有
即电位ϕ在无限远处趋于零,ϕ(r)|r→∞=0
介质分界面条件:
当场域中存在多种媒质时,还必须引入不同介质分界面上的边界条件,常称为辅助的边界条件。
静电场边值问题:
就是在给定的边界条件下,求解满足泊松方程或拉普拉斯方程的电位函数。
3.直接积分法
对于一些具有对称结构的静电场问题,电位函数仅是一个坐标变量的函数。
静电场边值问题可归结为常微分方程的定解问题。
这时可以直接积分求解电位函数。
例1:
图示二块半无限大导电平板构成夹角为α的电极系统。
设板间电压为U0,试求导电平板间电场。
[解]:
本例为平行平面场问题,选极坐标系进行分析。
显然电位仅是变量φ的函数,可以写出如下的第一类边值问题:
将泛定方程直接积分二次,得通解为
ϕ=C1φ+C2
由给定的两个边界条件,得
,C2=0
所以
例2:
求真空中球状分布电荷所产生的空间电场强度和电位分布,设电荷体密度为
[解]:
设球状电荷分布内、外的电位分别为ϕ1和ϕ2,显然,ϕ1满足泊松方程,ϕ2满足拉普拉斯方程。
由于电荷分布的球对称性,选球坐标系,有
(0 (r≥a) 边界条件为 ; ; 可解得ϕ1和ϕ2的通解为 代入边界条件,得 C1=0,C4=0,C2= ,C3= 最终得电位函数的解为 (r≤a) (r≥a) 利用球坐标系中的梯度表达式,求得 (r≤a) (r≥a) 可见,以上结果与应用高斯定理求得的结果完全一致。 4.分离变量法 基本思路: 当待求电位函数是二个或三个坐标变量的函数时,分离变量法是直接求解偏微分方程定解问题的一种经典方法。 对于拉普拉斯方程对应的边值问题,其步骤是: 首先,结合场域边界形状,选用适当的坐标系;其次,设待求电位函数由两个或三个各自仅含一个坐标变量的函数乘积组成,并代入拉普拉斯方程,借助于“分离”常数,将拉普拉斯方程转换为两个或三个常微分方程;第三,解这些常微分方程并以给定的定解条件决定其中的待定常数和函数后,即可解得待求的电位函数。 一般而言,当场域边界和某一正交曲线坐标系的坐标面相吻合时,分离变量法往往是一种简便而有效的方法。 直角坐标系中的平行平面场问题: 设电位函数为ϕ(x,y),满足拉普拉斯方程: 设电位函数有分离变量形式,即 ϕ(x,y)=X(x)Y(y) 代入拉普拉斯方程,整理得 显然,上式两边在x和y取任意值时恒成立,即等式两边应该恒为同一常数。 记该常数(常称为分离常数)为λ,这样,上式即转化为两个常微分方程 , 式中,分离常数λ可取0、mn2>0和-mn2<0,可分别得出如下三种形式的解,即 当λ=0时X(x)=A10+A20x;Y(y)=B10+B20y 当λ=mn2>0时X(x)=A1nchmnx+A2nshmnx;Y(y)=B1ncosmny+B2nsinmny 当λ=-mn2<0时X(x)=A1n'cosmnx+A2n'sinmnx;Y(y)=B1n'chmny+B2n'shmny 当mn取不同值时,上述解的线性组合便构成了拉普拉斯方程的通解,即 最后,可根据给定的定解条件,通过傅里叶级数展开方法,确定各个待定常数。 例3: 长直接地金属槽的横截面如图所示,其侧壁与底面电位均为零,顶盖电位为ϕ0。 求槽内电位分布。 [解]: 依题意,本问题为第一类边值问题,即 由于电位函数在x方向具有周期性、在y方向具有单调性,得A1n=0和A2n=0。 通解为 由边界条件,当x=0和y=0时,ϕ=0,得 A10=0,A1n'=0,B10=0,B1n'=0 即 又因为当x=a时,ϕ=0,得 C0=0, (n=1,2,3,⋯⋯) 故得 最后,当y=b时,ϕ=ϕ0,代入上式,有 作傅里叶正弦级数展开, 积分,得 又上式,得 本问题的电位函数解答为 本问题的等位线的分布如图虚线所示。 圆柱坐标系中的平行平面场问题: 设电位函数为ϕ(ρ,φ),满足拉普拉斯方程: 令电位函数为ϕ(ρ,φ)=R(ρ)Q(φ),代入上式,并理得 式中n2为分离常数,上式转化为下列两个常微分方程: 当n=0时R(ρ)=A10+A20lnρ;Q(φ)=B10+B20φ 当n≠0时R(ρ)=A1nρn+A2nρ-n;Q(φ)=B1ncosnφ+B2nsinnφ 得电位函数的通解为 由给定的边界条件,即可确定上式中的各个待定常数,最终得到待求的电位函数。 例4: 一个横截面半径为a,介电常数为ε1的长直介质圆柱体放置在均匀的外电中(场强为E0,方向与介质圆柱的轴线相垂直),介质外的介电常数为ε2,如图所示。 求圆柱体放入后,场域中的电位和电场强度。 [解]: 采用圆柱坐标系,且令z轴与圆柱轴重合,外电场方向与x轴同向,如图所示。 分别以ϕ1和ϕ2表示圆柱内外的电位函数。 首先,确定定解条件。 选坐标原点为电位参考点,即ϕ1=0,ρ=0 因而,均匀外电场E0=E0i对应的电位函数为 ϕ0=-xE0=-E0ρcosφ 显然,当ρ→∞时介质圆柱体产生的极化电场应当消失,在ρ→∞处的电位应与均匀外电场对应的电位ϕ0相一致,即 ϕ2=ϕ0=-E0ρcosφ,ρ→∞ 在圆柱表面ρ=a处,介质分界面的边界条件为 ϕ1=ϕ2 由本例图示可以看出,电场分布关于x轴对称,即ϕ1,2(ρ,φ)=ϕ1,2(ρ,-φ),这意味着特解Q(φ)是偶函数,所以,B20=B2n=0。 另外,根据场的对称性可以推知,y轴是电位等于零的等位线,即ϕ(ρ,±π/2)=0,也就是 A10=A20=0,B1ncos(±nπ/2)=0 可知n取奇整数。 然而,对ρ→∞的电位分析可知,B1n=0(n=3,5,…)。 综上讨论,得 ϕ1=(C1ρ+D1ρ-1)cosφ,ρ≤a ϕ2=(C2ρ+D2ρ-1)cosφ,ρ≥a 由ρ=0,ϕ1=0可得D1=0。 又由ρ→∞,ϕ2=-E0ρcosφ得C2=-E0。 最后利用介质分界面的边界条件,有 C1acosφ=(-E0a+D2/a)cosφ ε1C1cosφ=ε2(-E0-D2/a2)cosφ 解得 , 电位函数的解答为 ,ρ≤a ,ρ≥a 圆柱介质内外的电场强度为 可见,圆柱体内为均匀电场,且与外加均匀电场方向一致,如图画出了当ε2>ε1时的电场分布图。 值得注意的是,此时E1>E0,这表明若电介质内部有细长的空气泡是,则气泡内的电场强度增强,可能导致击穿绝缘损坏。 5.唯一性定理 本段将证明满足给定边界条件的泊松方程或拉普拉斯方程的解是唯一的,这也称作静电场的唯一性定理。 如图所示为充满均匀介质和置有n个导体的场域。 场域空间V的边界为S1,S2,…,Sn及外边界面S0。 设V中存在两个电位函数ϕ1和ϕ2,对于给定第一类或第二类边界条件,均满足泊松方程,即 , 令ϕd=ϕ1-ϕ2,因此 ∇2ϕd=0 利用格林公式 令ϕ=ψ=ϕd,代入上式得 如图所示,上式场域V的边界面S=S0+S1+S2+…+Sn。 如果所设的这两个不同的电位函数的解答ϕ1和ϕ2,在全部边界面上都应有相同的第一类边界条件或第二类边界条件,则它们在相应边界面Si上的差值 或 。 代入上式,有 这说明,场域V内ϕd的梯度处处为零,即V内所有场点上的ϕd值与其在各导体表面S1、S2、…、Sn上的值是相同的。 对于第一类边值问题,由于在导体表面上已知ϕd=0,所以整个场域内必有ϕd=0,由此得证ϕ1=ϕ2,即解唯一。 对于第二类边值问题而言,即已知各导体表面上的面电荷分布,此时ϕd=C,即电位ϕ1和ϕ2之间可能相差一个常数,但采用相同的电位参考点将导致C=0,所以解仍是唯一的。 静电场唯一性定理的重要意义在于,求解静电场问题时,不论采用哪一种解法,只要在场域内满足相同的偏微分方程、在边界上满足相同的给定边界条件,就可确信其解答是正确的。 2.6镜象法 镜象法的实质是以一个或多个位于场域边界外虚设的镜象(等效)电荷替代实际边界上未知的较为复杂的电荷分布,将原来具有边界的非均匀空间变换成无限大均匀媒质的空间,从而使计算过程得以简化。 根据唯一性定理可知,这些等效电荷的引入必须维持原问题边界条件不变,以保证原场域中的静电场分布不变。 通常这些等效电荷位于镜象位置,故称镜象电荷,由此构成的分析方法即称为镜象法。 1.对无限大接地导电平面的镜像 点电荷情况: 设有一点电荷q位于距无限大接地导电平面上方h处,其周围介质的介电常数为ε,如图所示。 显然,电位函数在场域内满足如下边值问题 ∇2ϕ=0(除去点电荷所在点) 边界条件为ϕ|y=0=0 可以设想,在场域边界外引入一个与点电荷q呈镜象对称的点电荷q'=-q,并将原来的导体场域由介电常数为ε的介质所替换。 这样,原场域边界面(z=0)上的边界条件ϕ=0保持不变,而对应的边值问题被简化为同一均匀介质ε空间内两个点电荷的电场计算问题。 根据唯一性定理可知,其解答有效区域仅限于图示上半部分介质场域。 应用镜像法,得待求电位为 无限大接地导电平面上的感应电荷的面密度分布为 式中负号表示感应电荷与点电荷q的极性相反。 对感应电荷作面积分,得 上式表明镜象电荷q'确实等效了无限大接地导电平面上的全部感应电荷。 此外,上述方法很容易推广到图示的由半无限大导电平面形成的劈形边界且其夹角为π的整数分之一的情况。 如图所示夹角为π/3的导电劈可以引入5个镜象电荷,以保证劈形边界电位为零的边界条件。 线电荷情况: 图示线电荷τ及其镜像电荷如图示。 由高斯定理得P点的电场强度为 图线电荷的镜象 (a)线电荷对无限大接地平面 (b)线电荷的镜象 -τ ∙ +τ ∙ ∙ o y P(x,y) 介质ε0 导体 x b ρ1 D D (-b,0) y P(x,y) ε0 x +τ ρ1 ε0 (b,0) ρ2 eρ2 eρ1 o 现任取Q点为电位参考点,则P点电位为 设在无限大接地导电平面上,即ρ1=ρ2时,ϕ=0,即电位参考点Q应选在接地导电平面上,所以C=0。 由上式,场中任意点电位为 由上式,可以进一步获得其等位线分布。 按等位线定义有ρ2/ρ1=K,平方得 整理,得 显然,上式为直角坐标系中圆的方程。 所以在xoy平面上,等位线分布是如图虚线所示的一簇圆。 对应于某一给定的K值,圆心坐标是 ,圆半径是 。 对于每个等位圆轨迹而言,圆半径a、圆心到原点的距离h和线电荷至原点的距离b三者间关系为 h2=a2+b2 亦即,a2=h2-b2=(h+b)(h-b) 这表明,两线电荷(±τ)位置对每个等位圆的圆心来说,满足圆的几何上反演的关系。 此外,当P点位于y轴右侧时,因ρ2/ρ1=K>1,ϕP皆为正值;当P点位于y轴左侧时,则ϕP皆为负值。 2.对无限大介质平面的镜像 (a)无限大介质平面上的点电荷 q ε2 en ε1 D1 h D2 q' (b)上半空间电场计算的镜象 θ ︒ E1" E1t P r θ E1' D1 h (c)下半空间电场计算的镜象 E2 ε2 ε2 θ E2t P r q" E2n ︒ D2 h 图无限大介质平面镜像 点电荷情况: 对于图示无限大介质平面上的点电荷边值问题也可采用镜象法。 上下半无限空间中的电场是由点电荷q及其分界面上的束缚电荷共同产生的。 对于介质为ε1的上半空间的电场计算,其分界面上的束缚电荷可归结为在均匀介质ε1中镜像点电荷q';对于介质为ε2的下半空间的电场计算,其分界面上的束缚电荷可归结为在均匀介质ε2中点电荷q"-q。 镜象电荷q'和q"的量值,可以通过分界面上的边界条件确定如下。 对于分界面上任意点P,由其上的边界条件E1t=E2t和D1n=D2n,得 解得 对于线电荷τ与无限大介质平面系统的电场,可类比推得。 3.电轴法 两半径相同的圆柱导体电场: 基于线电荷对无限大接地导电平面的镜像分析,我们可以进而讨论两同半径、带有等量异号电荷的平行长直圆柱导体间的电场问题。 此时,尽管圆柱导体表面电荷面密度不是常量,但沿轴向单位长电荷分布(线密度τ)是相同的,圆柱导体表面为等位面。 若设想圆柱导体表面与线电荷对应的等位面重合,即可以用等效线电荷计算圆柱导体外的电场分布,该线电荷就是圆柱导体表面电荷的等效电荷,如 (a)同半径两线输电线系统 (b)电轴法图示 y x a a h o -τ +τ D h b b 图电轴法 x a a h h o y -τ +τ D d 图所示。 为表述方便,成这个线电荷为圆柱导体的等效电轴,这种方法称为电轴法。 设圆柱导体半径为a,间距为2h,电轴间距为2b。 三者之间的关系为 例1: 半径为a的传输线平行于地面,传输线轴心对地高度为h,对地电位为U0,如图所示。 试求: (1)大地上方传输线的电场; (2)场域最大电场场强的位置及其数值。 [解]: (1)首先,由电轴法确定电轴的位置,得 大地上方任意场点P处的电位为 由传输线表面点A的电位U0,得 => 大地上方任意场点P处的电位为 (2)显然,最大场强将出现在导线相距地面最近处,即点A处,有 两半径不同的圆柱导体电场: 如图所示,设两平行长直圆柱导体半径分别为a1和a2,对于图(a)其轴心距d=h1+h2,对于图(b)其轴心距为d=h2-h1(设a2>a1)。 可以应用电轴法计算这两种情况的电场问题。 其关键问题仍然是确定等效电轴的位置。 显然 h12=b2+a12,h22=b2+a22,d=h2±h1 已知a1、a2和d,联立求解上式三个方程得 x d y a2 a1 h1 o -τ (a)平行传输线的电轴法图示 h2 b b +τ x y a2 a1 h1 o -τ (b)偏心同轴电缆的电轴法图示 h2 b b +τ 图半径不同圆柱导体的电轴法 , , 4.对导体球的镜像 导体球接地情况: 如图所示,设导体球半径为a,点电荷q至球心距为d。 设等效导体球表面感应电荷的镜象电荷为-q'且位于球内的球心与点电荷的连线上,其到球心的距离为b。 在导体球表面上任取一点P,得 => 整理,得[q2(a2+b2)-q'2(a2+d2)]+2a(q'2d-q2b)cosθ=0 上式对于任意的θ值恒成立,故有 解得 可以看出,点电荷q和其镜象电荷-q'的位置,满足球反演的几何关系。 根据q及-q'即可方便地计算点电荷在接地导体球外的电场分布。 可以证明,接地导体球面上感应电荷的总量等于-q'。 导体球不接地情况: 此时如导体球原不带净电荷,即呈中性,为使导体球表面上等电位,除引入镜象电荷-q'外,还应在原导体球的球心处再引入一个镜象电荷q"=q'。 同理,对呈电性的不接地导体球和位于导体球腔内的点电荷的电场计算问题,也可以应用镜像法进行计算。 例2: 图示为半径为a的接地导体球壳外置有一沿直径方向的线段电荷,线段的一端距球心为d。 求导体球壳上总的感应电荷。 [解]: 应用点电荷对接地导体球的镜像,有 元电荷为τdt,元电荷的位置为d+t;镜像元电荷为τ’dx=-aτdt/(d+t),镜像元电荷的位置为x+a2/d=a2/(d+t)。 所以,导体球壳上总的感应电荷为 图线段电荷的镜像 2.7电容与部分电容 电容或部分电容是导体系统的重要的集总电气参数,在是电网络中电容元件的重要参数,也是导体系统静电场的集总体现。 一般而言,需要借助于电场分析来计算。 1.两导体的电容 一般两导体电容的计算过程为: 给定两导体携带的电荷±q计算其电场分布和其间电位差U或给定两导体间电位差U,通过计算其电场分布和其携带的电荷±q,最后按定义计算电容C=q/U。 例1: 两半径为a、轴心距为d的平行长直圆柱导体构成一对均匀传输线,试求其单位长电容。 [解]: 应用电轴法,令h=d/2。 首先确定电轴位置 ,基于电轴法的分析结果,两导体表面最近距离对应的点A1(h-a,0)和点A2(-h+a,0)的电位差为 从而,均匀传输线的单位长度电容 通常有h>>a,此时b≈h,故 此外,对于h>>a的情况,也可以采用高斯定理计算。 设均匀传输线单位长线电荷密度为τ,则两导体轴心连线上距带正电荷导体x处的电场强度为 两导体间的电位差为 显然,有上式计算的电容与电轴法获得的结果相同。 从本例电容表达式可以看出,电容与导体之间施加的电压或携带的电荷量无关,只与导体的形状、相互位置和电介质有关,是导体系统自身固有电气参数。 2.部分电容 对于多导体需要引入部分电容概念。 静电独立系统: 系统的电场分布只与系统内各带电导体的形状、相互位置和电介质的分布有关,而与系统外的带电导体无关,并且所有电位移通量全部从系统内的带电导体发出又全部终止于系统内的带电导体。 现考察由(n+1)个导体组成的静电独立系统。 令各导体按0-n顺序编号,其相应的带电量分别为q0,q1,…,qk,…,qn。 由定义,知 q0+q1+…+qk+…+qn=0 选0号导体为电位参考点,即ϕ0=0,应用叠加原理,可得各个导体电位与各个导体上电荷的关系为 写成矩阵形式,为 {ϕ}=[α]{q} 式中,系数αij称为电位系数,其涵义不难从以下定义式得到理解 式中,αii称为自有电位系数,αij(i≠j)称为互有电位函数。 显然,电位系数只与导体的形状、相互位置以及电介质的介电常数有关。 当给出各个导体的电位时,有前式,得 {q}=[α]-1{ϕ}=[β]{ϕ} 或 式中,系数βij称为感应系数,与电位系数之间的关系为 式中,∆是[α]行列式,Aji是相应的代数余子式。 βii称为自有感应系数,βijj(i≠j)称为互有感应系数,即 显然,感应系数也只
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 工程 电磁场 教案 国家 精品 华北电力 学院 崔翔第 第二 部分