边界元法的MATLAB程序

利用MATLAB编写的边界元法程序。

function BF=BFI（XPYPXVYVIWJWKWBXBYXNUNTGES）
% XP YP 节点坐标
% XV YV 构成三角形小区域的节点坐标
% IW JW KW 顶点编号
% BX BY 体力分量给定值
% XNU 是泊松比
% NT 三角形小区域数
% GE 是剪切模量
% S 三角形小区域面积
% 用7点高斯积分公式

GT（1:）=[0.333333330.059715870.470142060.470142060.797426990.101286510.10128651];
GT（2:）=[0.333333330.470142060.059715870.470142060.101286510.797426990.10128651];
GT（3:）=[0.333333330.470142060.470142060.059715870.101256510.101286510.79742699];
WT=[0.2250.132394160.132394160.132394160.125929180.125929180.12592918];
BF=zeros（12）;
DE=4*pi*（1-XNU）;
for l=1:NT
    B=zeros（22）;
    for k=1:7
        XK=XV（IW（l））*GT（1k）+XV（JW（l））*GT（2k）+XV（KW（l））*GT（3k）;
        YK=YV（IW（l））*GT（1k）+YV（JW（l））*GT（2k）+YV（KW（l））*GT（3k）;
        ra=sqrt（（XP-XK）^2+（YP-YK）^2）;
        r1=（XK-XP）/ra;
        r2=（YK-YP）/ra;
        B（11）=B（11）+S（l）*（（3-4*XNU）*log（1/ra）+r1^2）*WT（k）/（2*DE*GE）;
        B（12）=B（12）+S（l）*r1*r2*WT（k）/（2*DE*GE）;
        B（22）=B（22）+S（l）*（（3-4*XNU）*log（1/ra）+r2^2）*WT（k）/（2*DE*GE）;
    end
    BF（1）=BF（1）+B（11）*BX（l）+B（12）*BY（l）;
    BF（2）=BF（2）+B（12）*BX（l）+B（22）*BY（l）;% 矩阵相乘
end
end

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

     文件      15343  2007-12-26 12:48  2dBEM\mymainfunction.m

     文件       1212  2007-12-20 14:00  2dBEM\BFI.m

     文件       2139  2007-12-26 09:46  2dBEM\DSKIJ.m

     文件       1438  2007-12-21 10:47  2dBEM\DTKIJ.m

     文件        374  2007-12-23 15:05  2dBEM\HGII.m

     文件       1415  2007-12-25 19:07  2dBEM\HGIJ.m

     文件         79  2007-12-11 15:17  2dBEM\BC_codeX.txt

     文件         79  2007-12-11 15:40  2dBEM\BC_codeY.txt

     文件        408  2007-12-22 12:36  2dBEM\BC_X.txt

     文件        400  2007-12-11 15:42  2dBEM\BC_Y.txt

     文件       1705  2007-12-26 12:48  2dBEM\elementcoordinate.txt

     文件        403  2007-12-26 12:48  2dBEM\internalcoordinate.txt

     文件       3308  2007-12-26 12:48  2dBEM\BoundaryConditions.txt

     文件       1186  2007-12-26 12:48  2dBEM\interCoorUP.txt

     文件       1747  2007-12-26 12:48  2dBEM\pointcoordinate.txt

     文件       1747  2007-12-26 12:48  2dBEM\pointp.txt

     文件       1747  2007-12-26 12:48  2dBEM\pointu.txt

     文件        463  2007-12-11 14:40  2dBEM\inputdata_x.txt

     文件        460  2007-12-11 14:50  2dBEM\inputdata_y.txt

     目录          0  2007-12-26 12:49  2dBEM

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

                35653                    20