matriceconductivite2d.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,ny)=0;

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