mesh1.m

这是一个有关解二维弹性板的无网格EFG程序matlab源代码,是学习无网格法入门的必学程序!

function[x,numq]=mesh1(length,height,ndivl,ndivw)

numq=(ndivl+1)*(ndivw+1);
% SET UP NODAL COORDINATES
for i=1:(ndivl+1)
    for j=1:(ndivw+1)
        x(1,((ndivw+1)*(i-1)+j))=(length/ndivl)*(i-1);
        x(2,((ndivw+1)*(i-1)+j))=-(height(i)/ndivw)*(j-1)+height(i)/2;
    end
end