maple计算地基均布载荷下岩土内应力分布（正确）

问题介绍：http://blog.sina.com.cn/s/blog_8dd135170102vzfo.html

-----------------------------------------------------------------------------

 

restart;
a := -10;
b := 10;
q := 1000000;
`σx` := proc (x, y) options operator, arrow; -q*(arctan((y+b)/x)-arctan((y-b)/x)+x*(y+b)/(x^2+(y+b)^2)-x*(y-b)/(x^2+(y-b)^2))/Pi end proc;

with(plots), with(plottools);
plot3d(`σx`(x, y), x = 0 .. 120, y = -40 .. 40);

 

-----------------------------------------------------------------------------

 restart

a := -10:

 b := 10:

q := 1000000:

`σx` := proc (x, y) options operator, arrow; -q*(arctan((y+b)/x)-arctan((y-b)/x)+x*(y+b)/(x^2+(y+b)^2)-x*(y-b)/(x^2+(y-b)^2))/Pi end proc

with(plots), with(plottools)

plot3d(`σx`(x, y), x = 0 .. 120, y = -40 .. 40)