有限元法工具代码

利用有限元法实现刚度矩阵的集成与EIT正问题的实现。

%总刚度矩阵
function ASTIF=ASSEMBLE（NPOINNELEMYOUNGPOISSTHICKCOORDLNODS）
% 张成特定大小总刚矩阵，并赋初值 0
% ASTIF（1:2*NPOIN1:2*NPOIN）=0;
% % 组装总刚度矩阵
% for i=1:NELEM
%     % 计算单刚
%     [ESTIFSMATX]=FORMESTIF（iCOORDLNODSYOUNGPOISSTHICK）;
%     nt=LNODS（i:）; % 临时向量 nt用来记录当前单元的结点号
%     for j=1:3
%         for k=1:3
%             ASTIF（（nt（j）*2-1）:nt（j）*2（nt（k）*2-1）:nt（k）*2）=...
%             ASTIF（（nt（j）*2-1）:nt（j）*2（nt（k）*2-1）:nt（k）*2）+...
%             ESTIF（j*2-1:j*2k*2-1:k*2）;
%             % 利用结点位移编号对应关系将单元刚度分块叠加到总刚矩阵
%         end
%     end
% end
ASTIF（1:NPOIN1:NPOIN）=0;
% 组装总刚度矩阵
for i=1:NELEM
    % 计算单刚
    [ESTIFSMATX]=FORMESTIF（iCOORDLNODS）;
    nt=LNODS（i:）; % 临时向量 nt用来记录当前单元的结点号
    for j=1:3
        for k=1:3
            ASTIF（nt（j）nt（k））=...
            ASTIF（nt（j）nt（k））+...
            ESTIF（jk）;
            % 利用结点位移编号对应关系将单元刚度分块叠加到总刚矩阵
        end
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1051  2017-06-27 17:39  FEM\ASSEMBLE.m

     文件        492  2017-06-26 19:36  FEM\BOUNDARY.m

     文件        697  2017-06-27 17:51  FEM\FORMESTIF.m

     文件        264  2017-06-26 16:11  FEM\FORMLOAD.m

     文件        785  2017-06-27 10:08  FEM\gaus.m

     文件        690  2017-06-26 20:17  FEM\Ke.txt

     文件       1727  2017-07-16 16:10  FEM\main1.m

     文件        135  2017-06-27 17:36  FEM\tri_in.txt

     文件      21924  2017-06-29 19:16  FEM\tri_in1.txt

     文件        426  2017-06-26 16:13  FEM\WRITESTRESS.m

     目录          0  2017-07-16 15:57  FEM

----------- ---------  ---------- -----  ----

                28191                    11