无穷长棱柱的温度分布 | 分离变量法（六）| 偏微分方程（十八）

截面为矩形的无穷长棱柱，内部无热源。两相对侧面与外界绝热，另两相对侧面温度分别保持为$0^o$和与高度无关的稳恒分布，求此棱柱的温度分布。

解：可以看成二维问题。设矩形为$0\le x\le a,0\le y\le b$，温度分布为$u(x,y)$则有边值问题
$$\{\begin{array}{l}\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2}=0,0<x<a,0<y<b\\ u{|}_{x=0}=0,u{|}_{x=a}=f(y)\\ \frac{\partial u}{\partial y}{|}_{y=0}=0,\frac{\partial u}{\partial y}{|}_{y=b}=0\end{array}\text{\hspace{1em}(13)}$$
设$u(x,y)=X(x)Y(y)$，代入方程和齐次边界条件，分离得固有值问题
$$\{\begin{array}{l}{Y}^{\prime \prime }(y)+\lambda Y(y)=0\\ Y^{\prime }(0)=Y^{\prime }(b)=0\end{array}\text{\hspace{1em}(14)}$$
和关于$X(x)$的常微分方程
$$X^{\prime \prime }(x)-\lambda X(x)=0$$
对比S-L型方程（9），固有值问题（14）式中$q(y)\equiv 0$，且两端不出现第I、III类边界条件。由S-L定理，固有值$\lambda ={\omega }^2\ge 0$，具体可解得固有值
$${\lambda }_n=(\frac{n\pi }{b}{)}^2,n=0,1,2,\cdots ,$$
固有函数
$$Y_n(y)=cos\frac{n\pi }{b}y$$
和相应的
$$\begin{gathered} X_0(x)=C_0+D_0(x) \\ {X}_n(x)=C_{n}ch\frac{n\pi }{b}x+D_{n}sh\frac{n\pi }{b}x,n=1,2,\cdots \end{gathered}$$
将所有变量分离形状解叠加，令
$$u(t,x)=C_0+D_{0}x+\sum _{n=1}^{+\infty }(C_{n}ch\frac{n\pi }{b}x+D_{n}sh\frac{n\pi }{b}x)cos\frac{n\pi }{b}y$$
由（14）式中另一组非齐次边界条件
$$u{|}_{x=0}=C_0+\sum _{n=1}^{+\infty }{C}_{n}cos\frac{n\pi }{b}y=0$$
可推出$C_n=0,n=0,1,2,\cdots$;
$$u{|}_{x=a}=D_{0}a+\sum _{n=1}^{+\infty }{D}_{n}sh\frac{an\pi }{b}cos\frac{n\pi }{b}y=f(y)$$
这是$f(y)$在$[0,b]$上的Fourier余弦展开，由（11）式可得
$$\begin{gathered} D_0=\frac{1}{ab}{\int }_0^{b}f(y)dy \\ {D}_n=\frac{2}{bsh\frac{an\pi }{b}}{\int }_0^{b}f(y)cos\frac{n\pi }{b}ydy,n=1,2,\cdot \cdot \cdot \end{gathered}$$
最后得矩形域上边值问题（13）式的形式解
$$u(x,y)=\frac{x}{ab}{\int }_0^{b}f(\eta )d\eta +\sum _{n=1}^{+\infty }\frac{2}{b}(sh\frac{an\pi }{b}{)}^{-1}{\int }_0^{b}f(\eta )cos\frac{n\pi }{b}\eta d\eta sh\frac{n\pi }{b}xcos\frac{n\pi }{b}y$$