qtd.m

四面体有限元程序,单元刚度的计算及总刚度矩征集成

h=0.041;x=-0.3:h:0.3;y=-0.3:h:0.3;
[X,Y]=meshgrid(x,y);
Z=1./sqrt((X-0.2).^2.+Y.^2)-1./sqrt((X+0.2).^2+Y.^2);
[PX,PY]=gradient(-Z,h);
figure(1)
contour(x,y,Z,[-12,-8,-5,-3,-2,-1,-0.5,-0.1,0.1,0.5,1,2,3,5,8,12],'b')
hold on
quiver(X,Y,PX,PY,'k')