连续介质动力学

从离散到连续

  考虑无限长的橡皮筋的纵向振动。首先离散处理,认为其由间隔为 aa、质量为 mm 的质量点以弹性系数为 kk 的轻弹簧连接而组成。设第 ii 个质量点相对平衡位置的偏离为 ηi\eta_i,则拉格朗日量写作:

L=i[12mη˙i212k(ηi+1ηi)2]=iaLiL = \sum_i \bigg[\frac{1}{2}m \dot{\eta}_i^2 - \frac{1}{2}k(\eta_{i+1} - \eta_i)^2\bigg] = \sum_i a L_i

其中

Li=12(ma)η˙i212ka(ηi+1ηia)2L_i = \frac{1}{2} \bigg(\frac{m}{a}\bigg) \dot{\eta}_i^2 - \frac{1}{2}ka\bigg(\frac{\eta_{i+1} - \eta_i}{a}\bigg)^2

之所以使用如上的分离,是因为可以认为 kaka 是杨氏模量 γ\gamma(根据胡克定理 F=γξF = \gamma \xiξ\xi 是延长量除以单位长度,这里取单位长度为 aa )。

然后使用平衡位置,而不是指标 ii 来重新描述离散系统:

ix,ηiη(x),ηi+1η(x+a)i \rightarrow x, \quad \eta_i \rightarrow \eta(x), \quad \eta_{i+1} \rightarrow \eta(x+a)

取连续极限:

ηi+1ηiaη(x+a)η(x)a=dηdx,aidx,maμ\frac{\eta_{i+1} - \eta_i}{a} \rightarrow \frac{\eta(x+a) - \eta(x)} {a} = \frac{\mathrm{d}\eta}{\mathrm{d}x}, \quad a\sum_i \rightarrow \int dx, \quad \frac{m}{a} \rightarrow \mu

则连续情形的拉格朗日量写作:

L=dx[12μη˙212γ(ηx)2]L = \int dx \bigg[\frac{1}{2}\mu \dot{\eta}^2 - \frac{1}{2}\gamma\bigg(\frac{\partial \eta}{\partial x}\bigg)^2\bigg]

如此,对于离散情形的欧拉-拉格朗日方程:

mη¨ik(ηi+1ηi)+k(ηiηi1)=0m\ddot{\eta}_i - k(\eta_{i+1} - \eta_i) + k(\eta_{i} - \eta_{i-1}) = 0

则有连续情形:

aμ2ηt2aγ2ηx2=0a\mu \frac{\partial^2 \eta}{\partial t^2} - a\gamma \frac{\partial^2 \eta}{\partial x^2} = 0

它是一个传播速度为 γ/μ\sqrt{\gamma/\mu} 的波动方程。下面我们来考虑直接从连续情形拉格朗日量得到它。

连续情形的最小作用量原理

我们需要为连续情形定义拉格朗日密度 hh。对于作用量 II

I[η]=h(η,η/x,η˙,x,t)dxdtI[\eta] = \int h(\eta, \partial\eta/\partial x, \dot{\eta}, x, t) \, \mathrm{d}x\mathrm{d}t

这表明 IIη\eta 的泛函。