非对称楔形重力坝受重力和液体压力的弹性力学解

袁学波, 袁江宏, 江晓禹, 李翔宇

袁学波, 袁江宏, 江晓禹等. 非对称楔形重力坝受重力和液体压力的弹性力学解. 力学与实践, xxxx, x(x): 1-8. DOI: 10.6052/1000-0879-24-437
引用本文: 袁学波, 袁江宏, 江晓禹等. 非对称楔形重力坝受重力和液体压力的弹性力学解. 力学与实践, xxxx, x(x): 1-8. DOI: 10.6052/1000-0879-24-437
Yuan Xuebo, Yuan Jianghong, Jiang Xiaoyu, et al. Elasticity solution for nonsymmetric wedge dams under gravity and hydrostatic pressure. Mechanics in Engineering, xxxx, x(x): 1-8. DOI: 10.6052/1000-0879-24-437
Citation: Yuan Xuebo, Yuan Jianghong, Jiang Xiaoyu, et al. Elasticity solution for nonsymmetric wedge dams under gravity and hydrostatic pressure. Mechanics in Engineering, xxxx, x(x): 1-8. DOI: 10.6052/1000-0879-24-437
袁学波, 袁江宏, 江晓禹等. 非对称楔形重力坝受重力和液体压力的弹性力学解. 力学与实践, xxxx, x(x): 1-8. CSTR: 32047.14.1000-0879-24-437
引用本文: 袁学波, 袁江宏, 江晓禹等. 非对称楔形重力坝受重力和液体压力的弹性力学解. 力学与实践, xxxx, x(x): 1-8. CSTR: 32047.14.1000-0879-24-437
Yuan Xuebo, Yuan Jianghong, Jiang Xiaoyu, et al. Elasticity solution for nonsymmetric wedge dams under gravity and hydrostatic pressure. Mechanics in Engineering, xxxx, x(x): 1-8. CSTR: 32047.14.1000-0879-24-437
Citation: Yuan Xuebo, Yuan Jianghong, Jiang Xiaoyu, et al. Elasticity solution for nonsymmetric wedge dams under gravity and hydrostatic pressure. Mechanics in Engineering, xxxx, x(x): 1-8. CSTR: 32047.14.1000-0879-24-437

非对称楔形重力坝受重力和液体压力的弹性力学解

基金项目: 西南交通大学2024年高水平育人课程教学改革项目(GK20240112),国家自然科学基金(12402114),四川省自然科学基金(2023NSFSC1301)和中央高校基本业务费项目(2682023CX053)资助。
详细信息
    作者简介:

    袁学波,讲师,从事柔性结构力学科学研究和弹性力学教学研究。E-mail:yuanxb@swjtu.edu.cn

    通讯作者:

    李翔宇,教授,研究方向为智能材料多场耦合力学。E-mail:zjuparis6@hotmail: .com

  • 中图分类号: O343

ELASTICITY SOLUTION FOR NONSYMMETRIC WEDGE DAMS UNDER GRAVITY AND HYDROSTATIC PRESSURE

  • 摘要:

    弹性力学教材给出了楔形重力坝的应力场理论解,但假定坝体受液压边界与重力方向平行且下端无限延伸,这与实际工程重力坝不符。本文针对更加贴近实际工程的非对称楔形重力坝,采用叠加原理和半逆解法推导了应力场的理论解,并与有限元分析结果进行了对比。理论解与有限元数值结果在远离坝底区域吻合,说明了理论解的准确性。本文研究结果不仅对重力坝的工程设计具有一定的指导意义,还为弹性力学提供了具体的教学实践。

    Abstract:

    The textbook of Elasticity presents a theoretical solution to the stress filed of a wedge dam. This model assumes that the boundary subjected to hydrostatic pressure is parallel to the gravity and that the lower end extends infinitely, which does not align with engineering design. In the current work, we establish a mechanics model for a nonsymmetric wedge dam that is more practically relevant, and derive the analytical solution to the stress field by virtue of the principle of superposition and the semi-inverse method. Finite element analysis is performed to obtain numerical results for the stress field, demonstrating the accuracy of the theoretical solution in regions far from the dam base. This work provides valuable insights for the engineering design of gravity dams and contributes a concrete example to the theory of elasticity.

  • 在弹性力学的经典教材[1,2]中,平面问题的直角坐标解答章节介绍了楔形重力坝受重力和液体压力作用的弹性力学问题,通过纯三次多项式应力函数推导了应力场的理论解,即Levy解。此模型对重力坝做了3个简化假定:(1)坝体承受液体压力的边界与重力方向平行;(2)坝体下端无限延伸,可自由变形;(3)坝体上端为尖顶且自由。需要指出的是,这些假设在实际工程中往往并不完全成立。在教学实践中,引导学生意识到简化假设的局限性,并鼓励学生思考更加贴近工程设计的力学模型,对于提升教学质量以及培养学生的理论建模能力、批判思维能力和理论与实践相结合的思维能力具有重要的促进作用。

    刘章军等[3]和周浩等[4]摒弃假定(3),按照教材中楔形体受重力和液体压力以及楔形体在楔顶受力的基本解答,采用叠加法对受压面铅垂的矩形变截面重力坝求解了应力场。类似地,肖先志等[5]摒弃假定(3),针对直角梯形截面重力坝,结合楔形体基本解答与叠加法构造了直角坐标应力函数,推导了应力场的解析表达。在实际工程设计中,大部分重力坝承受液体压力的边界是倾斜的,而非平行于重力方向;此外,坝体是有限高的,坝底受到地基约束,而非自由延伸。然而,笔者在调研中尚未发现理论文献摒弃假定(1),也未见文献定量分析弹性力学解对有限高重力坝的适用性。

    本文不再遵循假定(1),考虑一个非对称楔形体受重力和液体压力的弹性力学问题,采用半逆解法与叠加原理推导了平面应力场的理论解。采用ABAQUS通用软件建立有限高重力坝的平面应变有限元模型,有限元分析结果与理论解在远离坝底区域吻合。本文工作不仅能为重力坝的工程设计提供一定指导,还能深化弹性力学教学内容,兼具工程应用价值与教学实践意义。

    图1为非对称楔形重力坝受重力和液体压力的力学模型。重力坝视为很长的等截面柱体,其变形视为平面应变问题。直角坐标系(O-xy)的原点取为坝顶Ox轴正方向为水平向右,y轴正方向为垂直向下。楔面OAOBy轴正方向的夹角分别为αβ。左楔面OB承受液体压力,液面与坝顶平齐。重力平行于y轴正方向,重力加速度为g。坝体和液体的密度分别为ργ。坝体下端视为无限长,可自由变形。

    图  1  非对称楔形体受重力和液体压力的弹性力学模型与子问题分解

    由叠加原理,楔形重力坝受重力和液体压力的变形问题可分解为两个子问题。下面通过半逆解法推导两个子问题的应力场理论解。

    在极坐标系(O-rθ)中,用量纲分析方法寻求该问题的应力函数ϕ1。描述问题的载荷与几何参数包括:γgαβr和θ,其中γg表征了作用于楔面OB液体压力的大小,量纲为[力]/[长度]-3αβ和θ为无量纲量;r的量纲为[长度]。显然,应力分量是γgαβr和θ的函数表达,而应力分量的量纲为[力]/[长度]-2。从量纲上看,应力分量的形式应为γgrF(α, β, θ),其中F为无量纲函数。应力函数在r的幂次上比应力分量高两次,因此应力函数可假设为

    $$ {\phi _1} = {r^3}f(\theta ) $$ (1)

    式中,f(θ)是包含γgα和β的有量纲函数。应力函数必须满足相容方程$ \Big( \dfrac{{{\partial ^2}}}{{\partial {r^2}}} + \dfrac{1}{r}\dfrac{\partial }{{\partial r}} + \dfrac{1}{{{r^2}}} \cdot \dfrac{{{\partial ^2}}}{{\partial {\theta ^2}}} \Big)^2{\phi _1} = 0 $,将式(1)代入其中,可求得

    $$ f(\theta ) = A\cos \theta + B\sin \theta + C\cos 3\theta + D\sin 3\theta $$ (2)

    因此,应力函数为

    $$ {\phi _1} = {r^3}\left( {A\cos \theta + B\sin \theta + C\cos 3\theta + D\sin 3\theta } \right) $$ (3)

    进而,各应力分量的表达式求解为

    $$ \left. \begin{split} &\sigma _r^{(1)} = - 2r\big( - A\cos \theta - B\sin \theta + 3C\cos 3\theta + \\ &\qquad 3D\sin 3\theta \big) \\ & \sigma _\theta ^{(1)} = 6r\big( A\cos \theta + B\sin \theta + C\cos 3\theta +\\ &\qquad D\sin 3\theta \big) \\ &\tau _{r\theta }^{(1)} = - 2r\big( - A\sin \theta + B\cos \theta - 3C\sin 3\theta + \\ &\qquad 3D\cos 3\theta \big) \end{split} \right\} $$ (4)

    右楔面OAθ = π/2−α)为自由边界,左楔面OBθ = π/2)承受液体压力,故有应力边界条件

    $$ \left. {\begin{array}{*{20}{l}} {{{\left( {\sigma _\theta ^{(1)}} \right)}_{\theta = \pi /2 - \alpha }} = 0,}&{{{\left( {\tau _{r\theta }^{(1)}} \right)}_{\theta = \pi /2 - \alpha }} = 0} \\ {{{\left( {\sigma _\theta ^{(1)}} \right)}_{\theta = \pi /2 + \beta }} = - \gamma gr\cos \beta ,}&{{{\left( {\tau _{r\theta }^{(1)}} \right)}_{\theta = \pi /2 + \beta }} = 0} \end{array}} \right\} $$ (5)

    将应力分量表达式代入边界条件,4个待定参数求解为

    $$ \left. {\begin{array}{*{20}{l}} {A = - \dfrac{1}{8}\gamma g\dfrac{{\sin (4\alpha + 4\beta ) + \sin (4\alpha + 2\beta ) + 2\sin (2\alpha + 4\beta ) + 2\sin (2\alpha + 2\beta ) - 2\sin (2\beta )}}{{\cos (4\alpha + 4\beta ) - 4\cos (2\alpha + 2\beta ) + 3}}} \\ {B = - \dfrac{1}{8}\gamma g\dfrac{{\cos (4\alpha + 4\beta ) + \cos (4\alpha + 2\beta ) - 2\cos (2\alpha + 4\beta ) - 2\cos (2\alpha + 2\beta ) - \cos (2\beta ) + 1}}{{\cos (4\alpha + 4\beta ) - 4\cos (2\alpha + 2\beta ) + 3}}} \\ {C = - \dfrac{1}{{24}}\gamma g\dfrac{{\sin (4\alpha + 2\beta ) + 2\sin (2\alpha - 2\beta ) + \sin (4\alpha ) + 2\sin (2\alpha ) + 3\sin (4\beta ) + 3\sin (2\beta )}}{{\cos (4\alpha + 4\beta ) - 4\cos (2\alpha + 2\beta ) + 3}}} \\ {D = - \dfrac{1}{{24}}\gamma g\dfrac{{\cos (4\alpha + 2\beta ) + 2\cos (2\alpha - 2\beta ) + \cos (4\alpha ) + 2\cos (2\alpha ) - 3\cos (4\beta ) - 3\cos (2\beta )}}{{\cos (4\alpha + 4\beta ) - 4\cos (2\alpha + 2\beta ) + 3}}} \end{array}} \right\} $$ (6)

    将上述4个参数代回式(4),即可得子问题1在极坐标下的应力场理论解。由转换关系$ x = r\cos \theta $和$ y = r\sin \theta $,直角坐标下应力分量可表达为

    $$ \left. {\begin{array}{*{20}{l}} {\sigma _x^{(1)} = \dfrac{{\sigma _r^{(1)} + \sigma _\theta ^{(1)}}}{2} + \dfrac{{\sigma _r^{(1)} - \sigma _\theta ^{(1)}}}{2}\cos 2\theta - \tau _{r\theta }^{(1)}\sin 2\theta } \\ {\sigma _y^{(1)} = \dfrac{{\sigma _r^{(1)} + \sigma _\theta ^{(1)}}}{2} - \dfrac{{\sigma _r^{(1)} - \sigma _\theta ^{(1)}}}{2}\cos 2\theta + \tau _{r\theta }^{(1)}\sin 2\theta } \\ {\tau _{xy}^{(1)} = \dfrac{{\sigma _r^{(1)} - \sigma _\theta ^{(1)}}}{2}\sin 2\theta + \tau _{r\theta }^{(1)}\cos 2\theta } \end{array}} \right\} $$ (7)

    参考弹性力学教材中楔形体问题的直角坐标解答,由量纲分析将应力函数取为纯三次多项式,即

    $$ {\phi _2} = a{x^3} + b{x^2}y + cx{y^2} + d{y^3} $$ (8)

    坝体受x方向和y方向的体力分量分别为$ 0 $和$ \rho g $(常体力),因此应力分量表达为

    $$ \left. {\begin{array}{*{20}{l}} {\sigma _x^{(2)} = 2cx + 6dy} \\ {\sigma _y^{(2)} = 6ax + 2by - \rho gy} \\ {\tau _{xy}^{(2)} = - 2bx - 2cy} \end{array}} \right\} $$ (9)

    右楔面OA($ y = x\cot \alpha $)和左楔面OB($ y = - x\cot \beta $)均为自由边界,故应力边界条件为

    $$ \left. {\begin{array}{*{20}{l}} {\cos \alpha {{\left( {\sigma _x^{(2)}} \right)}_{y = x\cot \alpha }} - \sin \alpha {{\left( {\tau _{xy}^{(2)}} \right)}_{y = x\cot \alpha }} = 0} \\ {\cos \alpha {{\left( {\tau _{xy}^{(2)}} \right)}_{y = x\cot \alpha }} - \sin \alpha {{\left( {\sigma _y^{(2)}} \right)}_{y = x\cot \alpha }} = 0} \end{array}} \right\} $$ (10)
    $$ \left. {\begin{array}{*{20}{l}} { - \cos \beta {{\left( {\sigma _x^{(2)}} \right)}_{y = - x\cot \beta }} - \sin \beta {{\left( {\tau _{xy}^{(2)}} \right)}_{y = - x\cot \beta }} = 0} \\ { - \cos \beta {{\left( {\tau _{xy}^{(2)}} \right)}_{y = - x\cot \beta }} - \sin \beta {{\left( {\sigma _y^{(2)}} \right)}_{y = - x\cot \beta }} = 0} \end{array}} \right\} $$ (11)

    将应力分量表达式代入边界条件,4个待定参数求解为

    $$ \left. {\begin{array}{*{20}{l}} {a = \dfrac{1}{{12}}\rho g\dfrac{{2\sin (2\alpha + 4\beta ) - \sin (4\alpha + 2\beta ) + 2\sin (2\alpha - 2\beta ) - \sin 4\alpha + \sin 2\alpha + \sin 4\beta - \sin 2\beta }}{{\cos (4\alpha + 4\beta ) - 4\cos (2\alpha + 2\beta ) + 3}}} \\ {b = \dfrac{1}{4}\rho g\dfrac{{\cos (4\alpha + 2\beta ) + 2\cos (2\alpha + 2\beta ) - 2\cos (2\alpha + 2\beta ) - \cos 4\alpha - \cos 4\beta + 1}}{{\cos (4\alpha + 4\beta ) - 4\cos (2\alpha + 2\beta ) + 3}}} \\ {c = \dfrac{1}{4}\rho g\dfrac{{\sin (2\alpha + 4\beta ) - \sin (4\alpha + 2\beta ) - 2\sin (2\alpha - 2\beta ) + \sin 4\alpha + \sin 2\alpha - \sin 4\beta - \sin 2\beta }}{{\cos (4\alpha + 4\beta ) - 4\cos (2\alpha + 2\beta ) + 3}}} \\ \begin{split} d =& \dfrac{1}{{12}}\rho g\big[ \cos (4\alpha + 4\beta ) - 2\cos (4\alpha + 2\beta ) - 2\cos (2\alpha + 4\beta ) - 4\sin 2\alpha \sin 2\beta + \\ &\cos 4\alpha + 2\cos 2\alpha + \cos 4\beta + 2\cos 2\beta - 3 \big]/\left[{{\cos (4\alpha + 4\beta ) - 4\cos (2\alpha + 2\beta ) + 3}}\right] \end{split} \end{array}} \right\} $$ (12)

    将上述4个参数代回式(9),即得到子问题2在直角坐标下的应力场理论解。

    基于子问题1与子问题2的应力场理论解,楔形重力坝在重力与液体压力作用下的应力分量最终可通过叠加原理表示为

    $$ {\sigma _x} = \sigma _x^{(1)} + \sigma _x^{(2)},\;\;{\sigma _y} = \sigma _y^{(1)} + \sigma _y^{(2)},\;\;{\tau _{xy}} = \tau _{xy}^{(1)} + \tau _{xy}^{(2)} $$ (13)

    若左楔面OB为铅垂面(即β = 0),式(6)中各参数简化为

    $$\left. \begin{split} A =& - \dfrac{1}{4}\gamma g{\cot ^3}\alpha ,\;B = - \dfrac{1}{8}\gamma g + \dfrac{1}{8}\gamma g{\cot ^2}\alpha \\ C = &- \dfrac{1}{{12}}\gamma g{\cot ^3}\alpha ,\;D = \dfrac{1}{{24}}\gamma g + \dfrac{1}{8}\gamma g{\cot ^2}\alpha \end{split}\right\}$$ (14)

    由极坐标应力解(式(4))与坐标变换公式(式(7)),子问题1在直角坐标系下的应力分量表达为

    $$ \left. {\begin{array}{*{20}{l}} {\sigma _x^{(1)} = - \gamma gy} \\ {\sigma _y^{(1)} = - 2\gamma gx{{\cot }^3}\alpha + \gamma gy{{\cot }^2}\alpha } \\ {\tau _{xy}^{(1)} = - \gamma gx{{\cot }^2}\alpha } \end{array}} \right\} $$ (15)

    式(12)中各参数简化为

    $$ a = \dfrac{1}{6}\rho g\cot \alpha ,\;b = c = d = 0 $$ (16)

    子问题2在直角坐标系下的应力分量(式(9))表达为

    $$ \sigma _x^{(2)} = 0,\;\sigma _y^{(2)} = \rho gx\cot \alpha - \rho gy,\;\tau _{xy}^{(2)} = 0 $$ (17)

    叠加式(15)与式(17),应力场理论解在β = 0情况下退化为

    $$ \left. \begin{split} &{{\sigma _x} = \sigma _x^{(1)} + \sigma _x^{(2)} = - \gamma gy} \\ & {\sigma _y} = \sigma _y^{(1)} + \sigma _y^{(2)} = (\rho g\cot \alpha - 2\gamma g{{\cot }^3}\alpha )x +\\ &\qquad (\gamma g{{\cot }^2}\alpha - \rho g)y \\ & {{\tau _{xy}} = \tau _{xy}^{(1)} + \tau _{xy}^{(2)} = - \gamma gx{{\cot }^2}\alpha } \end{split} \right\} $$ (18)

    该式与弹性力学教材中的经典表达式(即Levy解)一致。

    通过有限元分析对重力坝进行数值模拟。坝体材料的杨氏模量和泊松比分别设置为E = 30.0 GPa和μ = 0.2,坝体材料的密度为ρ = 2500 kg/m3,液体密度为γ = 1000 kg/m3,重力加速度为g = 9.8 m/s2。坝体的铅垂高度H取为200 m。采用ABAQUS软件建立平面有限元模型(图2(a)),坝体被离散成精细的四节点平面应变单元(CPE4R)。离散单元的最小尺寸取为0.1 m,足够保证有限元分析的网格收敛性。坝体底部的所有自由度均被约束。采用DLOAD子程序对楔面OB施加线性分布的液体压力。设定5条水平路径(y = 20,50,100,150和180 m),提取各应力分量的有限元分析值。

    图  2  含铅垂楔面的重力坝示意图以及各应力分量沿5条水平路径的理论与有限元结果对比

    对于α = 15o的楔形重力坝,图2(b)和图2(c)展示了3个应力分量(σxσyτxy)沿5条水平路径的分布情况,离散点为有限元分析结果,实线为式(18)的理论计算结果。正应力σx沿x方向没有变化且为负值,正应力σy沿x方向单调减少,剪应力τxy沿x方向从零逐渐线性减小。可见,当水平路径离坝体底部相对较远时(如y ≤ 150 m),应力分量σxσyτxy的理论解与有限元结果非常吻合。然而,在接近坝体底部区域(如y = 180 m),理论解与有限元结果的偏差非常明显。

    本文理论模型假定坝体无限高且底部可自由变形,而有限元分析模型是有限高的且底部完全约束。理论解与有限元分析结果在靠近坝体底部区域的偏差来源于底部边界条件的差异。对于α = 15o的楔形重力坝,图3展示了左楔面不同位置处(y = 50,100和180 m)的正应力(σy)OB随重力坝高度H的变化关系。对于3个不同的左楔面位置,当重力坝高度H较小时,(σy)OB随高度H减小,这是因为底部边界效应减弱。当重力坝高度H足够大时,(σy)OB与高度H几乎无关,且与理论解吻合。

    图  3  含铅垂楔面的重力坝不同位置处左楔面应力(σy)OB随重力坝高度H的变化关系

    工程重力坝的坝体通常由脆性的混凝土材料制成,其抗拉能力远小于抗压能力。在重力坝的结构设计中,应使坝体处于压应力状态。根据Levy解,右楔面OA正应力(σy)OA= −γgycot2α为负,左楔面OB正应力(σy)OB= (γgcot2αρg)y的正负依赖于γgcot2αρg的相对大小。如图4(a)所示,给定坝体材料与液体的密度(ρ = 2500 kg/m3γ = 1000 kg/m3),增加角度α可降低正应力σy的水平。图4(b)展示了(σy)OB随角度α的变化。理论上,若角度α大于安息角arctan(γ/ρ)1/2 ≈ 32o,(σy)OB为负值,σy沿x方向始终小于零。有限元分析预测的安息角与理论计算结果接近。由此可见,基于坝体无限长假定的Levy解对于重力坝设计具有一定的指导意义。但是,坝底附近正应力σy的理论解低于有限元结果(图3),因此重力坝的整体安全设计仍需借助有限元分析等数值模拟手段。

    图  4  不同角度α下应力σy沿水平路径y = 50 m的分布及左楔面应力(σy)OB随角度α的变化关系

    在重力坝的实际工程设计中,承受液体压力的左楔面OB往往并非为铅垂面(β > 0),如图5(a)所示。此处,针对α = 15oβ = 5o对应的非对称楔形重力坝,图5(b)和图5(c)展示了3个应力分量沿着5条水平路径的分布情况。该楔形重力坝的应力分布情况与Levy解存在较大差别。正应力σx沿x方向近似线性增大,且始终保持为负值。剪应力τxy在左楔面上不再为零,沿x方向近似线性减小,且也始终保持为负值。相比于α = 15oβ = 0的情况,该重力坝的正应力σy水平较低且分布更加集中。与前面结果类似,理论解与有限元结果在远离坝底区域非常吻合,但在接近坝底区域具有明显偏差。

    图  5  含非铅垂楔面的重力坝示意图以及各应力分量沿5条水平路径的理论与有限元结果对比

    图6(a)展示了不同角度β下正应力σy沿水平路径y = 50 m的分布情况,此处固定α = 15o。可见,增加左楔面的倾斜角度β可降低正应力σy的水平,从而提升重力坝的安全性。按照理论计算,当β大于临界值25o时,左楔面正应力(σy)OB为负值,如图6(b)所示。该理论临界值与有限元结果非常吻合。

    图  6  不同角度β下应力σy沿水平路径y = 50 m的分布及左楔面应力(σy)OB随角度β的变化关系

    由上述结果可知,本文针对非对称楔形重力坝建立的弹性力学理论解具有较高的准确性。在工程重力坝的优化设计实践中,本文理论解有望用于提供一组良好的几何设计初值,从而有效地减少优化迭代次数,提升设计效率。

    本文的理论模型和应力场求解综合了叠加原理、基于量纲分析的半逆解法、应力坐标变换公式、直角坐标解答以及极坐标解答等重要的弹性力学教学内容,推导难度比教材中楔形体基本问题稍大,适合用于思维能力扩展。在教学中,笔者将本文模型应用于本科生课程“弹性力学”,具体设计如下。

    (1)“平面问题的直角坐标解答”的“楔形体受重力和液体压力”章节

    首先,笔者讲解Levy解后,向学生展示我国三峡大坝、礼亨水库等典型重力坝的实景图和截面图,启发引导学生思考工程重力坝与教材模型的差异,引出Levy解的3个简化假定。

    然后,展示有限高重力坝的有限元分析模型和仿真结果,通过对比分析Levy解的正确性。

    最后,引导学生建立摒弃简化假定的弹性力学模型。摒弃假定(1)的力学模型即为本文模型。通过叠加原理,将问题分解为两个子问题。针对子问题1,说明采用直角坐标解答的推导难度,为后续的极坐标解答留下伏笔。针对子问题2,采用教材中的纯三次多项式推导应力场的理论解。

    (2)“平面问题的极坐标解答”的“楔形体在楔顶或楔面受力”章节

    首先,笔者在完成“楔面受均布载荷作用的问题”的讲解后,引入本文的子问题2(即楔面受非均匀、线性分布载荷),采用量纲分析和半逆解法推导应力场的理论解。

    然后,叠加子问题1和子问题2,给出非对称楔形体受重力与液体压力作用下应力场的理论解,展示有限元分析与理论解的对比。简要分析本文模型对工程重力坝设计的指导意义。

    最后,引导学生同时摒弃3个简化假定,建立具有非直角梯形重力坝的力学模型,作为课后探索。

    笔者采用启发式教学法和案例教学法,通过上述流程将本文内容融入到弹性力学的教学设计中。通过课堂讲解与案例分析,学生直接接触到实际工程中的弹性力学问题,从而将课堂理论与工程实践相结合,激发了学生的学习热情,并获得了积极的课堂反馈。读者可以进一步采用分组讨论、翻转课堂等形式实施本文内容的教学,鼓励学生自主进行理论建模、求解与分析,提升其批判性思维和实践能力。

    针对非对称楔形重力坝受重力和液体压力的弹性力学问题,本文采用半逆解法和叠加原理推导了应力场的理论解,并与有限元分析进行了对比。本文所得理论解与有限元分析结果在远离坝底的区域非常吻合,显示出较高的准确性。该理论解对重力坝的工程设计具有一定的理论指导意义和实际应用价值。本文的理论求解过程融合了叠加原理、半逆解法、直角坐标解答与极坐标解答等诸多弹性力学内容,不仅深化了对教材中经典问题的探讨,还能启发学生思考经典模型的局限性,从而有助于提升教学质量。同时摒弃3个简化假定的非直角梯形重力坝更加符合工程设计,其在重力、液体压力以及坝顶载荷作用下的弹性力学解值得进一步探索。

  • 图  1   非对称楔形体受重力和液体压力的弹性力学模型与子问题分解

    图  2   含铅垂楔面的重力坝示意图以及各应力分量沿5条水平路径的理论与有限元结果对比

    图  3   含铅垂楔面的重力坝不同位置处左楔面应力(σy)OB随重力坝高度H的变化关系

    图  4   不同角度α下应力σy沿水平路径y = 50 m的分布及左楔面应力(σy)OB随角度α的变化关系

    图  5   含非铅垂楔面的重力坝示意图以及各应力分量沿5条水平路径的理论与有限元结果对比

    图  6   不同角度β下应力σy沿水平路径y = 50 m的分布及左楔面应力(σy)OB随角度β的变化关系

  • [1] 徐芝纶. 弹性力学, 第5版. 北京: 高等教育出版社, 2016

    Xu Zhilun. Elasticity, 5th edn. Beijing: Higher Education Press, 2016 (in Chinese)

    [2] 王光钦, 丁桂宝, 杨杰. 弹性力学, 第3版. 北京: 清华大学出版社, 2015

    Wang Guangqin, Ding Guibao, Yang Jie. Elasticity, 3rd edn. Beijing: Tsinghua University Press, 2015 (in Chinese)

    [3] 刘章军, 叶永, 周宜红等. 用楔形体解答求解矩形变截面梁及其适用范围. 力学与实践, 2012, 34(2): 71-74
    [4] 周浩, 王博, 徐建国. 利用弹性理论求混凝土重力坝应力分布. 人民黄河, 2012, 34(6): 120-124

    Zhou Hao, Wang Bo, Xu Jianguo. Solving the stress distribution of concrete gravity dam on the basis of elastic theory. Yellow River, 2012, 34(6): 120-124 (in Chinese)

    [5] 肖先志, 肖奉英, 黄模佳. 梯形重力坝的应力函数解及其有限元验证. 力学与实践, 2023, 45(1): 113-118

    Xiao Xianzhi, Xiao Fengying, Huang Mojia. Stress function solution for trapezoidal gravity dam and its finite element verification. Mechanics in Engineering, 2023, 45(1): 113-118 (in Chinese)

图(6)
计量
  • 文章访问数:  46
  • HTML全文浏览量:  4
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-10-29
  • 录用日期:  2025-01-12
  • 网络出版日期:  2025-01-16

目录

/

返回文章
返回