part3浙江大学电气学院电磁场数值计算.pptxVIP

3.非线性有限元法(2D)1）牛顿-拉夫逊迭代法考察非线性方程式中系数k是u的函数，右端项 f 为常量。December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法若设则原问题求解就是求取由曲线g(u)与直线f的交点uDecember 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法为了求解上述方程，首先选取初始解u0，并将g(u)在u0点展开，有忽略高次项，有December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法这是一个线性代数方程组未知量 系数由此可得December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法是经过一次迭代得到的近似解。解得 后，可将其作为再次迭代的起点，继续用下式求取二次近似解December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法其过程可形象地表示为：December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法一般地，经过k次迭代都得到u(k)后求取u(k+1)的方程为经过适当次数迭代，上述方程的左端项接近于零，从而其解接近于收敛值。December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法一般地，经过k次迭代都得到u(k)后求取u(k+1)的方程为经过适当次数迭代，上述方程的左端项接近于零，从而其解接近于收敛值。December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法In summary构成了以函数的斜率为系数，以近似解的迭代差为变量，以余量为右端项的线性代数方程组，并由它求取每一次迭代的近似解。December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法In summary对于非线性代数方程组将上式左端项定义为向量{g}December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法Jacobion MatrixDecember 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法December 16, 20093.非线性有限元法(2D)1）牛顿-拉夫逊迭代法收敛判据December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式此时，偏微分方程为December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式对应的余量（弱式的）形式与线性问题一样，区别仅在于?与u的值有关December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式在一个单元内klh的计算同线性问题December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式在一个单元内December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式定义December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式A.线性媒质当媒质为线性时，此时[k]e与{ue}无关，因此有December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式A.线性媒质如果初始值 ，则与线性问题完全一致，一次即可解出December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质这时[k]e与{ue}有关。此时，对于有December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质在一个单元内，[k]e正比于{ue}。于是December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质令矢量磁位为磁通密度B，标量磁位为磁场强度HDecember 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质则December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质由磁化曲线得到December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质December 16, 20093.非线性有限元法(2D)2）在有限元法中的计算格式B. 非线性媒质December