st.m

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

       x=-0.3:0.018:0.3;y=-0.3:0.018:0.3;
       [X,Y]=meshgrid(x,y);
       Z=1./sqrt((X-0.2).^2.+Y.^2)-1./sqrt((X+0.2).^2+Y.^2);
       mesh(x,y,Z)
       hold on
       contour(x,y,Z,[-8,-5,-3,-2,1,-0.5,0.5,1,2,3,5,8],'k')
       hold off