--- title: "从边界应力，推到有限厚度冰壳的Love数" description: "求解同密度液体上的不可压缩壳，加入自引力，再用同一位移场推导偏心率潮汐的总功率。" canonical_url: "https://sandboxphysics.com/zh/stories/ocean-worlds/derive-a-finite-thickness-love-response/" language: "zh-CN" page_type: "story" alternate_languages: en: "https://sandboxphysics.com/stories/ocean-worlds/derive-a-finite-thickness-love-response/" last_modified: "2026-10-28" --- [从材料记忆走到潮汐加热与海洋世界](https://sandboxphysics.com/zh/stories/ocean-worlds/) AC12 · L2 # 从边界应力，推到有限厚度冰壳的Love数 Sandbox Physics · 2026-10-28 壳厚不仅改变整体刚度，还会改变冰在内部怎样运动。我们可以直接求解有限厚度边界问题，不必把厚壳当成薄膜。四个弹性边界条件给出几何刚度，再借助黏弹性对应原理得到复Love响应，以及符号能够独立检验的热功率。 ## 开始前需要什么 先读Maxwell章节。本章使用线性弹性、二阶球谐、复振幅和基础面积积分，并逐项规定归一约定。 ## 先选择一个能完整算清的内部 天体半径取五百千米，恒定密度为每立方米一千千克。固态壳与无黏液体内部具有相同密度，均不可压缩。壳层的Maxwell性质均匀，液体不承受剪应力；潮汐变化足够慢，内部力平衡忽略惯性。 同密度在这里是一条实质性假设。内部边界移动时，没有密度反差，也就不额外产生界面引力质量。诱导重力势全部来自外表面的移动，二阶潮汐的重力Love数恰好是径向Love数的五分之三。密度分层会打破这个固定比例。 液体半径与总半径之比，等于一减去壳厚与总半径之比。液体半径归零时，天体成为全固态；模型真正包含这一极限，没有在中心留下假想液体点。这里没有岩核、可压缩性、海洋共振、旋转驱动的内部动力学，也不让冰温度随深度改变。 $$ g=\frac{4\pi G\rho R}{3},\quad x=1-\frac dR,\quad k_2=\frac35h_2 $$ 保存例子的表面重力为每平方秒零点一三九七八六米。 对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 ## 先解形状，再求整体振幅 先将二阶角向图样与径向变化分开。把表面径向函数归一为一，内部径向函数写成四个幂次的线性组合。不可压缩条件确定切向函数，弹性平衡方程给出压力项；剩下四个系数，正好由边界条件确定。下式径向坐标已经除以总半径。 外表面的切向牵引为零。液体界面同样不承受切向牵引；扣掉共同的液体静水载荷后，弹性法向牵引扰动也为零。加上表面归一条件，共有四条线性方程，实验直接解这个系统。全固态极限中，中心正则性排除负幂项，再由两个外边界条件得到独立的正则解。 外表面无量纲法向牵引定义几何刚度。薄壳时，它近似为二十四倍壳厚除以十一倍半径；实心球时为十九除以五。默认五十千米冰壳得到零点二三二六九一。表面的切向位移也来自同一组系数，因此可见壳边界和能量积分使用的是同一份力学解。 $$ U(s)=a s+b s^3+c s^{-2}+d_0 s^{-4},\quad V(s)=\frac{sU^\prime+2U}{6},\quad p(s)=7b s^2+2c s^{-3} $$ 最后一个系数与真实壳厚用不同符号；压力与牵引采用省略了位移、模量公共因子的归一方式。 $$ T_t=V^\prime+\frac{U-V}{s},\quad T_r=-p+2U^\prime,\quad U(1)=1,\quad T_t(1)=T_t(x)=T_r(x)=0 $$ 四个条件适用于液体半径非零时；实心球的正则解另为径向坐标一次项八倍减三次项三倍，再除以五。 $$ z_h=\frac{24}{5}\frac{19-75x^3+112x^5-75x^7+19x^{10}}{24+40x^3-45x^7-19x^{10}} $$ 先从边界条件数值求出位移，再用论文闭式表达式独立检查几何刚度。 对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 ## 让自引力与弹性一起抵抗载荷 径向位移会改变重力势和静水压力。表面法向平衡中，向内的静水恢复作用与施加潮汐势、诱导势以及额外弹性牵引相平衡。把整个平衡式除以密度与外加势的乘积，就得到径向Love数、负一、负重力Love数和几何弹性阻力之和。 代入同密度Love数关系后求解，流体自引力恢复系数成为五分之二。再把弹性剪切模量替换为复Maxwell模量，这就是线性对应步骤。实恢复系数保持为正，被动材料的正虚部刚度则产生负虚部位移响应。 零频率的Maxwell天体有充分时间松弛剪应力，径向Love数趋于二点五，重力Love数趋于一点五。理想弹性天体则保留有限刚度，只是响应没有虚部。这是两种不同的取极限方式，流体极限与弹性零耗散极限需要分别检验。 $$ h_2-1-k_2+z_h\frac{\mu^*}{\rho gR}h_2=0\quad\Longrightarrow\quad h_2=\left[\frac25+z_h\frac{\mu^*}{\rho gR}\right]^{-1},\qquad k_2=\frac35h_2 $$ 密度、重力与半径的乘积具有压力单位，整个分母无量纲。 对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 ## 把默认冰壳从头算到尾 取五十千米壳厚、零点零零五偏心率、三点五吉帕弹性模量和十的十四点五次方帕秒黏度。液体半径占比为零点九，四个牵引方程给出几何刚度零点二三二六九一；天体的密度、重力和半径相乘，得到约六千九百八十九点三一万帕的压力尺度。 公转频率下，材料章节得到的复模量约为三点一八四三九加一点零零二五一虚数吉帕。先除以压力尺度，再乘几何刚度，最后加流体恢复系数五分之二，复分母约为十一点零零一六加三点三三七六零虚数。给分子分母乘复共轭求逆，再乘五分之三，就得到重力响应。 表面势振幅随方向变化，在球面上将其模长取最大，得到下式。乘径向Love数模长并除以重力，便得到三点零八三一九米；相位滞后为十六点八七六四度。最后，完整功公式给出一百点零零九吉瓦。独立体积耗散检验使用内部应变场，而不是把这个最终标量功率重复计算一遍。 $$ h_2=\frac1{11.0016+3.33760i}\simeq0.0832352-0.0252513i,\qquad k_2\simeq0.0499411-0.0151508i $$ 压力尺度、几何刚度与模量，都进入这份明确的数值代入。 $$ |\widehat W|_{\max}=\frac32\sqrt{\frac{32}{7}}\frac{GM_pR^2}{a^3}e,\qquad |\widehat u_r|_{\max}=\frac{|h_2|}{g}|\widehat W|_{\max} $$ 这里在所有表面方向中选取最大的谐波径向振幅；场景则在指定相位显示真实位移。 对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 ## 距离变化和方向摆动，都要保留 取天体固定第一轴朝向行星的平均方向，第三轴沿自转轴。平均同步旋转扣掉平均公转角。偏心率一阶中，距离变化改变潮汐强度，光学天平动改变它的方向；下面的复对称无迹张量同时包含两种作用，稳态潮汐已扣除。 用单位径向向量在张量两侧收缩，就得到表面势的复振幅。径向响应等于Love数乘这个势，再除以重力。壳内部的径向和切向函数来自牵引边界解，共同组成位移向量。液体内部采用正则、不可压缩势流延拓，径向运动连续，允许切向滑移。 轨道画面使用开普勒椭圆，施加潮汐保留偏心率首阶；两者的差别从偏心率二阶开始。允许的偏心率很小，不能因为画了一条更扁的椭圆，就认为内部模型已经包含非线性轨道谐波。扩展到大偏心率，需要同时升级载荷和内部响应。 $$ \widehat W(\boldsymbol n)=\boldsymbol n^{\mathsf T}\widehat{\boldsymbol S}\boldsymbol n,\qquad \widehat{\boldsymbol S}=\frac{GM_pR^2}{2a^3}e\begin{pmatrix}6&-6i&0\\-6i&-3&0\\0&0&-3\end{pmatrix} $$ 虚数非对角项表示正时间虚指数约定下的光学天平动。 $$ \widehat{\boldsymbol u}(r,\boldsymbol n)=\frac{h_2}{g}\left[U(r/R)\widehat W\boldsymbol n+V(r/R)\boldsymbol\nabla_\Omega\widehat W\right] $$ 外表面径向函数为一；液体延拓与半径成正比，在中心保持正则。 对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 ## 把功积起来，而不是涂一张热图 不可压缩条件把体积中的潮汐功转成表面势与径向速度的积分。两个真实谐波相乘后取平均，等于一个复振幅与另一个共轭振幅之积实部的一半。代入位移响应，便留下Love数的负虚部；被动滞后因此得到正功率。 角向积分必须使用完整张量。对称无迹张量对应的势振幅模平方积分，等于八倍圆周率除以十五，再乘张量范数平方。对角项和光学天平动项合起来，才给出最终功率中的二十一除以二。若漏掉方向摆动，这个系数就会改变。 保存的默认例子中，外部做功得到约一百点零零九吉瓦。另一路体积计算对位移做笛卡尔数值微分，组成对称应变张量，再积分冰壳中的黏性剪切损耗，结果与外部功一致。两种计算都没有把画面的蓝白材料颜色解释成温度。 $$ \mathcal P=\rho R^2\int\left\langle W\dot u_r\right\rangle d\Omega=-\frac{\rho R^2\omega}{2g}\operatorname{Im}h_2\int|\widehat W|^2d\Omega $$ 符号直接来自已经规定的谐波约定和径向位移响应。 $$ \int|\widehat W|^2d\Omega=\frac{84\pi}{5}\left(\frac{GM_pR^2}{a^3}\right)^2e^2,\quad \mathcal P=-\frac{21}{2}\frac{GM_p^2R^5}{a^6}\omega e^2\operatorname{Im}k_2 $$ 真实的同步偏心率潮汐频率为轨道平均运动；诊断用谐波探测在改变频率时保持空间载荷振幅固定。 对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 · Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025) 保存的数值例子、单位与方法 合成不可压缩、同密度天体。施潮天体质量为五乘十的二十六次方千克，轨道半长轴三十万千米。四个弹性牵引条件确定壳位移，复Maxwell对应确定滞后与发热；这些数值是预测，只有明确标为拟合的行使用观测式合成噪声。 | 半径、壳厚与偏心率 | $R=500\,\mathrm{km},\quad d=50\,\mathrm{km},\quad e=0.00500$ | | --- | --- | | 两层密度与表面重力 | $\rho=1000\,\mathrm{kg\,m^{-3}},\quad g=0.139786\,\mathrm{m\,s^{-2}}$ | | 有限壳几何刚度 | $z_h=0.232691$ | | 径向与重力响应 | $h_2=0.0832352-0.0252513i,\quad k_2=0.0499411-0.0151508i$ | | 最大径向振幅与总功率 | $u_{r,\max}=3.08319\,\mathrm{m},\quad\mathcal P=100.009\,\mathrm{GW}$ | 场景扣掉稳态潮汐并放大周期径向形状。每种材料与壳厚分别在参考偏心率下归一显示增益，比较真实振幅请使用米制读数。外部做功和体积耗散的独立积分核对总功率。 ## [移动边界，检验边界条件](https://sandboxphysics.com/zh/labs/tidal-locking/) 1. 从薄壳拖到厚壳，解释液体半径怎样改变，并把真实米制位移与显示增益一起比较。 2. 先预测偏心率加倍后功率变为四倍，再核对读数，最后切到圆轨道预设。 ## 练习与解释性答案 1. 全固态极限为什么需要单独的正则解，不能把负径向幂直接代到中心？ 展开解释 负幂在零半径发散，必须先令对应系数归零。剩下的正则解给出十九除以五的几何刚度，并恢复均匀不可压缩实心球的Love响应。 2. 保留轨道尺度与材料，把偏心率加倍，位移和功率分别怎样变？ 展开解释 首阶潮汐势和位移都加倍，功率包含载荷与速度之积，因此变为四倍。同一种材料的显示增益以固定参考偏心率设定，所以画面周期形变也会加倍。 3. 真实卫星有更密的岩核，还能强制使用重力与径向Love数的固定五分之三比例吗？ 展开解释 不能。密度反差会增加内部移动质量，改变引力边界问题；必须先解分层结构，再比较两种观测。固定比例属于本同密度例子的假设。 [上一篇：材料的记忆，怎样把潮汐循环变成热](https://sandboxphysics.com/zh/stories/ocean-worlds/from-a-maxwell-material-to-tidal-lag/) [下一篇：一个复Love数，究竟能识别哪些内部信息](https://sandboxphysics.com/zh/stories/ocean-worlds/what-a-complex-love-number-can-identify/) ## 来源与计算范围 1. [Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1](https://arxiv.org/html/1410.4735v1) 固定为2014年10月17日预印本，Icarus正式发表DOI10.1016/j.icarus.2014.10.027。附录A式105–116给出有限厚度均匀不可压缩壳；附录C讨论Maxwell对应原理，附录E推导功与功率，附录G规定偏心率和光学天平动潮汐。图8–9、12–13展示近似与热模型依赖。核读相关正文和附录；本站的同密度Maxwell例子比真实卫星简单得多。 2. [Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)](https://www.nature.com/articles/s41586-025-09818-x) 正式版2025年12月17日，DOI10.1038/s41586-025-09818-x，CC BY 4.0。核读出版方存档全文PMC12711566及Methods。图1给出重力响应估计，图2比较有条件的内部后验，图3是推断内部的示意图。图1误差棒为三倍标准差，图2椭圆为一倍标准差。原文正虚部表示其耗散约定；本站用正时间虚指数，被动Love响应的虚部为负。