matlab使用有限元方法求解偏微分方程

大量的有限元法求解偏微分方程的实例程序，程序里有详细的解释语句

%----------------------------------------------------------------------------
% Example
%   A barrel vault has radius r=25 ft subtended angle of 80 degrees
%   length 50 f and thickness 3in. 
%   Elastic modulus is 3 msi Poisson‘s ratio is 0 and
%   weight of the shell is 90 lb per square ft. Two curved edges are
%   supported by rigid diaphragm and the other two edges are free. 
%   Solve using 4 by 4 elements of a quarter of the shell.         
%
% Variable descriptions                                                      
%   k = element matrix in the local axes                                            
%   ke = element matrix in the global axes
%   kb = element matrix for bending stiffness
%   ks = element matrix for shear stiffness
%   km = element matrix for membrane stiffness
%   f = element vector
%   kk = system matrix                                             
%   ff = system vector                                                 
%   disp = system nodal displacement vector
%   gcoord = coordinate values of each node
%   nodes = nodal connectivity of each element
%   index = a vector containing system dofs associated with each element     
%   pointb = matrix containing sampling points for bending term
%   weightb = matrix containing weighting coefficients for bending term
%   points = matrix containing sampling points for shear term
%   weights = matrix containing weighting coefficients for shear term
%   bcdof = a vector containing dofs associated with boundary conditions     
%   bcval = a vector containing boundary condition values associated with    
%           the dofs in ‘bcdof‘                                              
%   kinmtsb = matrix for kinematic equation for bending
%   matmtsb = matrix for material property for bending
%   kinmtss = matrix for kinematic equation for shear
%   matmtss = matrix for material property for shear
%   kinmtsm = matrix for kinematic equation for membrane
%   matmtsm = matrix for material property for membrane
%   tr3d = transformation matrix from local to global axes
%
%----------------------------------------------------------------------------            

%------------------------------------
%  input data for control parameters
%------------------------------------

clear
nel=16;                   % number of elements
nnel=4;                  % number of nodes per element
ndof=6;                  % number of dofs per node
nnode=25;                 % total number of nodes in system
sdof=nnode*ndof;         % total system dofs  
edof=nnel*ndof;          % degrees of freedom per element
emodule=3e6;            % elastic modulus
poisson=0.0;             % Poisson‘s ratio
t=3;                   % plate thickness
nglxb=2; nglyb=2;        % 2x2 Gauss-Legendre quadrature for bending 
nglb=nglxb*nglyb;        % number of sampling points per element for bending
nglxs=1; nglys=1;        % 1x1 Gauss-Legendre quadrature for shear 
ngl

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

     文件      12089  2000-05-12 06:54  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\Ex10101.m

     文件      15466  2000-05-12 09:18  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\Ex10102.m

     文件       8737  2000-05-29 08:21  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\Ex1091.m

     文件       8753  2000-05-29 08:22  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\Ex1092.m

     文件       1971  2000-05-29 08:22  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX11102.M

     文件       2089  2000-05-29 08:22  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX11111.M

     文件       2549  1999-09-13 11:48  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\Ex1281.m

     文件       4956  1999-09-13 11:29  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\Ex1291.m

     文件       4980  1999-09-13 11:29  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\Ex1292.m

     文件       4854  1999-09-13 11:29  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\Ex1293.m

     文件       4849  1999-09-13 11:29  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\Ex1294.m

     文件       4493  2007-05-29 14:31  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX351.asv

     文件       4442  2000-05-29 08:23  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX351.M

     文件       4643  2000-05-29 08:23  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX352.M

     文件       4362  2000-05-29 08:23  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX353.M

     文件       5359  2000-05-29 08:23  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX5101.M

     文件       7262  1995-03-05 09:19  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX5111.M

     文件       7624  1995-03-05 13:13  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX5112.M

     文件       7204  1995-03-05 10:17  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX5113.M

     文件       9032  1999-03-09 04:42  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX5114.M

     文件       8671  1995-03-05 22:07  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX5115.M

     文件       5214  2000-05-29 08:24  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX5121.M

     文件       8119  2000-05-29 08:24  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX591.M

     文件       7516  2000-05-29 08:24  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX592.M

     文件       1501  2000-05-29 08:24  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX651.M

     文件       1843  2000-05-29 08:25  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX652.M

     文件       2059  2000-05-29 08:25  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX653.M

     文件       3364  2000-05-29 08:25  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX661.M

     文件       9687  1996-04-01 04:57  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX662.M

     文件       6798  2000-05-29 08:25  fem_matlab大量的有限元法求解偏微分方程的实例程序（程序有解释）\fem_matlab\fem_matlab\EX741.M

............此处省略93个文件信息