mesh1.m

来自「这是一个有关解二维弹性板的无网格EFG程序matlab源代码,是学习无网格法入门」· M 代码 · 共 10 行

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