function?BF=BFI（XPYPXVYVIWJWKWBXBYXNUNTGES）
%?XP?YP?節(jié)點(diǎn)坐標(biāo)
%?XV?YV?構(gòu)成三角形小區(qū)域的節(jié)點(diǎn)坐標(biāo)
%?IW?JW?KW?頂點(diǎn)編號
%?BX?BY?體力分量給定值
%?XNU?是泊松比
%?NT?三角形小區(qū)域數(shù)
%?GE?是剪切模量
%?S?三角形小區(qū)域面積
%?用7點(diǎn)高斯積分公式

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

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件??????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