电机学课程设计.docVIP

2010 级 《电机学》课程设计 组 员 张丰伟、彭永晶 彭鸿昌、邱天、李彦青 所在院系 电气与电子工程 班 级 电气1004班 日 期 2013年3月 作业评分 评阅人 设计题目 题目：铁心磁路计算 铁心磁路如图所示，磁路尺寸为：=3mm、=2mm、w=125mm、h=150mm、N1=2N2=100匝，铁芯宽度l=50mm、铁心厚度d=50mm。铁芯宽度和厚度均匀，忽略铁芯磁场边缘效应。 假设铁芯的磁导率为无穷大，（1）若I2=0，Φ1=6mWb，求I1（2）若I1=10A，I2=20A，求Φ1和Φ2。 若铁芯材料为DR510钢片（磁化曲线见教材），若I2=0，Φ1=6mWb，求I1和Φ1。 若I1=10A，I2=20A，编写程序，求Φ1和Φ2。 设计过程 由磁路基尔霍夫第一定律=0，可得，由于铁芯宽度和厚度均匀，截面积相等，则 由磁路基尔霍夫第二定律，得: 磁场关系: 则方程为： 假设铁芯的磁导率为无穷大，（1）若I2=0，Φ1=6mWb，求I1（2）若I1=10A，I2=20A，求Φ1和Φ。 铁芯的磁导率为无穷大，则铁芯磁路H=0，方程化简为 将I2=0，Φ1=6mWb带入上式方程，得I1=57.2958A 将I1=10A，I2=20A带入上式方程，得B1=0.418879T，B2=-0.628319T，则Φ1=B1ld=1.04720mWb，Φ2=B2ld=-1.57080mWb 若铁芯材料为DR510钢片（磁化曲线见教材），若I2=0，Φ1=6mWb，求I1和Φ1。 Φ1=6mWb,则B1=Φ1/(ld)=2.4T 采用编程法求解，由方程 可看作未知数为B2的一元非线性方程，通过二分法求解出B2，再计算出B3，然后由方程 计算出I1。 为提高计算精度，可使用拉格朗日插值法或曲线拟合法，对表中不能直接得到的数据进行近似计算。本题使用拉格朗日插值法。 拉格朗日插值法： 假设任意两个不同的xj都互不相同，那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为： 其中每个为拉格朗日基本多项式（或称插值基函数），其表达式为： 拉格朗日基本多项式的特点是在 上取值为1，在其它的点 上取值为。 计算结果：B1=2.4T ，B2=0.721749T ，B3=1.67825T Φ1=0.006Wb，Φ2=0.0018044Wb，Φ3=0b I1=57.2958A 改变程序，使收敛精度降低为10^-10，则运算精度降低，但循环次数减小，时间开销减下。 计算结果：Φ1=0.006Wb，Φ2=0.0018044Wb，Φ3=0b I1=57.2958A 若I1=10A，I2=20A，编写程序，求Φ1和Φ2。 编程用迭代法实现，计算时先用铁芯磁导率为无穷大情况时近似计算B1，B2，B3的值，将B3的值作为定值计算出H3的值带入方程（1）（2），用弦截法计算出B1和B2的值并计算出B3的值，与之前的值比较，若在允许的误差范围以内则计算结束输出计算值，否则将B3的值带入进行下一次循环，直到B3的精度达到要求。 计算步骤如下： 近似情况下求得 代入方程（1）（2），通过弦截法求得， 求得 与比较，若精度达到要求，则，，即为所求 若精度达不到要求，则重复2,3,4步骤，求得，直到精度足够即与之差足够小 通过得到的求 用近似计算结果充当第零次迭代值，因精确值与近似值偏差不是很大，可加快收敛，减小迭代次数，并可避免因步长太大或太小而不能收敛。 程序执行时可选择是否进行拉格朗日插值法或曲线拟合法，提高磁场强度的精确度。拉格朗日插值法计算量很大，时间消耗太大，故采用曲线拟合法。 因曲线变化范围较大，对整个曲线进行拟合比较困难，如使用6次多项式进行拟合，系数达到105数量级，R2为0.9995。因此，对曲线进行分段拟合： 图像如图所示，其中蓝点为磁化曲线上的点，不同颜色的曲线为拟合曲线。 磁化曲线（Ｂ－Ｈ） 磁化曲线（Ｈ－Ｂ） 流程图如下： 代码如下： #include iostream #include cstdlib #include cmath #include iomanip #include stdlib.h using namespace std; double delta1=0.003,delta2=0.002,w=0.125,h=0.15,n1=100,n2=50,l=0.05,d=0.05; double pi=3.1415