matriceconductivte2d.asv

２Ｄ力学模型.有限单元法.巴黎六大实验教材.

%**********************************************************
% MatriceConductivite2D.m
%
%**********************************************************


function [K,Q]=MatriceConductivte2D(hx,hy,nx,ny,Nx,Ny)

%1) remplissage de la matrice avec des zeros.

K=sparse([]);
K(nx,nx)=0;

for i=1:nx & j=1:ny
    if j==1;
        K(i,j)=1
    elseif j==ny
        K(i,j)=1;
     if i==1 & j=2:(ny-1)
        K(i,j)=1;
        K(i+1,j)=-1;
     else if i==nx & j=2:(ny-1)
        K(i,j)=1;
        K(i-1,j)=-1;
      else K(i+1,j)=1;
           K(i,j)=-4;
           K(i-1,j)=1;
           K(i,j+1)=1;
           K(i,j-1)=1;
         end
 
         
         
%2) Creation du vecteur second membre
for i=1:nx & j=2:(ny-1)
Q(i,1)=T0;
Q(i,ny)=T1;
Q(1,j)=0;
Q(nx,j)=0;
end