电磁场生医-实验3要点.doc

《电磁场数值》实验报告一、实验目的 掌握二维电流场有限元分析的方法，编制相应程序，包括有限元系数矩阵的生成，边界条件的处理，方程组的求解等，培养解决实际电磁场问题的能力。 二、实验条件 硬件：微型计算机。 软件：MATLAB 6.0软件。 三、实验内容 1、针对二维电流场，根据变分原理或加权余量法推导其有限元方程形式。 2、编制其有限元分析程序，并进行求解。 、实验步骤 1、预习：推导生物医学电磁场中二维稳态电流场的有限元离散方程。 2、根据场域剖分的结果，利用三角形单元的形状函数求解系数矩阵。 3、处理第二类和第一类边界条件。 4、求解有限元线性方程组。 实验结果 nsx=3;nsy=1; %材料种类 nsx1=[4,4,4];nsy1=8; %剖分为几行几列 x1=[0,3,8,12];y1=[0,8]; %交界面坐标 num_nodex=sum(nsx1)+1; %行节点数 num_nodey=sum(nsy1)+1; total_node=num_nodex*num_nodey; %总节点数 total_element=(num_nodex-1)*(num_nodey-1)*2; %总单元数 rou1=[2,3,1]