swp2.m

有限元法计算平板受集中载荷的程序

function [point2,weight2]=swp2(nglx,ngly)

%-------------------------------------------------------------------
% 目的:
%     求二维高斯积分的积分点和权系数
%
% 变量：
%     nglx - x轴高斯积分点数
%     ngly - y轴高斯积分点数
%     point2 - 高斯积分点坐标
%     weight2 - 权系数
%-------------------------------------------------------------------

%  确定x,y轴最大的积分点数

   if nglx > ngly
      ngl=nglx;
   else
      ngl=ngly;
   end

%  初始化

   point2=zeros(ngl,2);
   weight2=zeros(ngl,2);

% 求出相应的积分点坐标和权系数

 [pointx,weightx]=swp1(nglx);    
 [pointy,weighty]=swp1(ngly);     

%  二维积分

 for intx=1:nglx                     
   point2(intx,1)=pointx(intx);
   weight2(intx,1)=weightx(intx);
 end

 for inty=1:ngly                    
   point2(inty,2)=pointy(inty);
   weight2(inty,2)=weighty(inty);
 end