阅读说明：本技术 一种基于雪载荷作用下柔性体承载重量预测方法及系统 (method and system for predicting bearing weight of flexible body under snow load ) 是由 王单 王健 于 2019-10-28 设计创作，主要内容包括：本发明公开一种基于雪载荷作用下柔性体承载重量预测方法及系统,所述预测方法确定雪载荷在柔性体表面的分布模式；当所述分布模式为三角形堆垛,则根据第一载荷分布公式预测雪载荷作用下柔性体的承载重量；当所述分布模式为梯形堆垛,则根据第二载荷分布公式预测雪载荷作用下柔性体的承载重量,从而实现对雪灾的预防及检测；可见采用本发明中的方法进行预测,无需针对不同的柔性体进行不同实验,因此降低了实验方法的工作量。(The invention discloses flexible body bearing weight prediction methods and systems based on snow load, wherein the prediction method determines a distribution mode of snow load on the surface of a flexible body, when the distribution mode is triangular stacking, the bearing weight of the flexible body under the snow load is predicted according to a th load distribution formula, when the distribution mode is trapezoidal stacking, the bearing weight of the flexible body under the snow load is predicted according to a second load distribution formula, and therefore snow disaster prevention and detection are achieved.)

一种基于雪载荷作用下柔性体承载重量预测方法及系统

技术领域

本发明涉及柔性结构的承载能力测试技术领域，特别是涉及一种基于雪载荷作用下柔性体承载重量预测方法及系统。

背景技术

风中摇曳的柳树，水中摆动的水藻以及被雪压弯的树枝等都是大自然中常见的现象，这些柔性结构都蕴藏着一种智慧，即通过自身的结构变形降低了承受的外力载荷，从而提高自身的承载能力。研究柔性系统的阻力问题在生物、农业和森林业中有着重要的意义，如理解植物对于它们生存环境的适应性，预测和阻止农作物倒伏和强风把树连根拔起等现象。

柔性体的结构变形减小阻力对于工程中也十分有借鉴意义。工程中的结构大部分都是刚性的，刚性的结构在承受荷载时不发生显著的变形，但是相应承担的阻力也比较大。因此，在某些环境下，对变形要求不高时，可以考虑用柔性结构来降低承载阻力。目前，相关的工程例子有锥概念的风力发电机,微型飞行汽车的柔性机翼，拍动翼推进和烂泥污水处理过程中气体通过中空的柔性纤维转移等。

在雪灾严重的地区，树枝、房屋等往往被压得变形得很厉害，每年因为雪载荷造成的经济损失也十分显著。自然界的雪因为成分的复杂性，还会混有白霜、冰之类，直接研究有很大的局限性。目前，对于柔性体在风载与流体载荷作用下的承载力分析都已经有相关的理论工作，但是对于固体载荷作用下柔性体的承载作用，比如雪载荷作用下树枝的承载力分析，主要是基于实验的方法，研究雪的密度、凝聚力等参数对于承载力的影响，即针对不同的柔性体都需要进行不同实验，难以得到通用性的结论，因此实验方法工作量大。

发明内容

本发明的目的是提供一种基于雪载荷作用下柔性体承载重量预测方法及系统，预测在雪载荷作用下柔性体承受重量以及结构变形度，从而实现对雪灾的预防及检测。

为实现上述目的，本发明提供了一种基于雪载荷作用下柔性体承载重量预测方法，所述预测方法包括：

确定雪载荷在柔性体表面的分布模式；

当所述分布模式为三角形堆垛，则根据第一载荷分布公式预测雪载荷作用下柔性体的承载重量；

当所述分布模式为梯形堆垛，则根据第二载荷分布公式预测雪载荷作用下柔性体的承载重量。

可选的，所述预测方法还包括：

根据所述柔性体的承载重量确定柔性体的结构变形度。

可选的，所述根据所述柔性体的承载重量确定柔性体的结构变形度，具体公式为：

其中，W_f为柔性体的承载重量，ρ为雪载荷的密度，g为重力加速度，L为长度，α为雪载荷与柔性体表面作用的最大摩擦角。

可选的，所述当所述分布模式为三角形堆垛，则根据第一载荷分布公式预测雪载荷作用下柔性体的承载重量，具体包括：

根据第一载荷分布公式确定第一雪载荷分布高度函数；

根据所述第一雪载荷分布高度函数确定柔性体的承载重量。

可选的，所述当所述分布模式为梯形堆垛，则根据第二载荷分布公式预测雪载荷作用下柔性体的承载重量，具体包括：

根据第二载荷分布公式确定第二雪载荷分布高度函数；

根据所述第二雪载荷分布高度函数确定柔性体的承载重量。

可选的，所述第一载荷分布公式为：

其中，H(t)为第一雪载荷分布高度函数，α为雪载荷与柔性体表面作用的最大摩擦角，θ为变形后弧长s处梁截面法向与变形前截面法向的夹角，t＝s/L，s为弧长，L为长度，自重影响因子M＝m/ρL，ρ为雪载荷的密度，m为质量，柯西数CY＝ρgL⁴/B，g为重力加速度，B为柔度。

可选的，所述第二载荷分布公式为：

其中，H(u)为第二雪载荷分布高度函数，α为雪载荷与柔性体表面作用的最大摩擦角，θ为变形后弧长s处梁截面法向与变形前截面法向的夹角，t＝s/L，s为弧长，L为长度，自重影响因子M＝m/ρL，ρ为雪载荷的密度，m为质量，柯西数CY＝ρgL⁴/B，g为重力加速度，B为柔度，h_c为临界高度。

本发明还提供一种基于雪载荷作用下柔性体承载重量预测系统，所述预测系统包括：

分布模式确定模块，用于确定雪载荷在柔性体表面的分布模式；

第一承载重量确定模块，用于当所述分布模式为三角形堆垛，则根据第一载荷分布公式预测雪载荷作用下柔性体的承载重量；

第二承载重量确定模块，用于当所述分布模式为梯形堆垛，则根据第二载荷分布公式预测雪载荷作用下柔性体的承载重量。

可选的，所述第一承载重量确定模块，具体包括：

第一雪载荷分布高度函数确定单元，用于根据第一载荷分布公式确定第一雪载荷分布高度函数；

第一承载重量确定单元，用于根据所述第一雪载荷分布高度函数确定柔性体的承载重量。

可选的，所述第二承载重量确定模块，具体包括：

第二雪载荷分布高度函数确定单元，用于根据第二载荷分布公式确定第二雪载荷分布高度函数；

第二承载重量确定单元，用于根据所述第二雪载荷分布高度函数确定柔性体的承载重量。

根据本发明提供的具体实施例，本发明公开了以下技术效果：

本发明公开一种基于雪载荷作用下柔性体承载重量预测方法及系统，所述预测方法确定雪载荷在柔性体表面的分布模式；当所述分布模式为三角形堆垛，则根据第一载荷分布公式预测雪载荷作用下柔性体的承载重量；当所述分布模式为梯形堆垛，则根据第二载荷分布公式预测雪载荷作用下柔性体的承载重量，从而实现对雪灾的预防及检测；可见采用本发明中的方法进行预测，无需针对不同的柔性体进行不同实验，因此降低了实验方法的工作量。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例基于雪载荷作用下柔性体承载重量预测方法流程图；

图2为本发明实施例三角形堆垛与梯形堆垛示意图；

图3为本发明实施例最大摩擦角的确定示意图；

图4为本发明实施例柔性体受力示意图；

图5为本发明实施例基于雪载荷作用下柔性体承载重量预测系统结构图；

图6为本发明实施例三角形堆垛下梁的承受载荷与柯西数、自重因子之间的关系图；

图7为本发明实施例梯形堆垛下梁的承受载荷与柯西数、自重因子之间的关系图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明的目的是提供一种基于雪载荷作用下柔性体承载重量预测方法及系统，预测在雪载荷作用下柔性体承受重量以及结构变形度，从而实现对雪灾的预防及检测。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

图1为本发明实施例基于雪载荷作用下柔性体承载重量预测方法流程图，如图1所示，本发明提供一种基于雪载荷作用下柔性体承载重量预测方法，所述预测方法包括：

步骤S1：确定雪载荷在柔性体表面的分布模式；

步骤S2：当所述分布模式为三角形堆垛，则根据第一载荷分布公式预测雪载荷作用下柔性体的承载重量；

步骤S3：当所述分布模式为梯形堆垛，则根据第二载荷分布公式预测雪载荷作用下柔性体的承载重量。

本发明所述预测方法还包括：

步骤S4：根据所述柔性体的承载重量确定柔性体的结构变形度，具体公式为：

其中，W_f为柔性体的承载重量，ρ为雪载荷的密度，g为重力加速度，L为长度，α为雪载荷与柔性体表面作用的最大摩擦角。

下面对各个步骤进行详细论述：

步骤S1：确定雪载荷在柔性体表面的分布模式，如图2所示，本发明所述分布模式包括三角形堆垛和梯形堆垛，图2中的(a)为三角形堆垛，图2中的(b)为三角形堆垛中弧长s是经过截面质心的曲线，图2中的(c)为梯形堆垛，图2中的(d)为梯形堆垛中弧长s是经过截面质心的曲线。

步骤S2：当所述分布模式为三角形堆垛，则根据第一载荷分布公式预测雪载荷作用下柔性体的承载重量，具体包括：

步骤S21：根据第一载荷分布公式确定第一雪载荷分布高度函数，具体包括：

图4为本发明实施例柔性体受力示意图，如图4所示，弧长s是经过截面质心的曲线，θ为变形后弧长s处梁截面法向与变形前截面法向的夹角，h为变形后弧长s处梁截面的高度(厚度)，下面以三角形堆垛为例，推导梁的控制方程，具体推导过程如下：

建立大变形梁中的弯矩与转角之间的关系式：

其中，B为柔度，θ为变形后弧长s处梁截面法向与变形前截面法向的夹角，s为弧长；

对弯矩求微分有：

力的平衡关系：

X方向：

d(F_ssinθ)-d(F_Tcosθ)＝0 (3)

其中，F_T为截取微段上沿着弧长切向的力，F_s为截取微段上沿着弧长法线方向的力，d()为力在截取的微单元内的增量。

Y方向：

d(F_scosθ)+d(F_Tsinθ)＝ρghdx+mgdx (4)

其中，h为梁截面的高度，

几何关系为：

dx＝ds cosθ (5)

由(1)式和(2)式得：

F_s＝Bθ” (6)

从式(3)推出：

F_T＝F_stanθ+C＝Bθ”tanθ+C (7)

梁的自由端的边界条件为F_s＝F_T＝0，代入(7)式可以求出常数C＝0，则公式(7)整理得：

F_T＝Bθ”tanθ (8)

将(8)式代入(4)式，化简之后可以导出变形梁的平衡方程：

定义以下无量纲化坐标及无量纲量：

则将公式(10)代入公式(9)确定所述第一载荷分布公式为：

其中，H(t)为第一雪载荷分布高度函数，θ为变形后弧长s处梁截面法向与变形前截面法向的夹角，t＝s/L，s为弧长，L为长度，自重影响因子M＝m/ρL，ρ为雪载荷的密度，m为质量，柯西数CY＝ρgL⁴/B，g为重力加速度，B为柔度，α为雪载荷与柔性体表面作用的最大摩擦角，如图3所示，用小颗粒模拟雪载荷，由于内聚力作用，颗粒将在基底表面进行堆垛，能够堆积的最大角度记为最大摩擦角。

最后根据所述第一载荷分布公式确定第一雪载荷分布高度函数H(t)。

步骤S22：根据所述第一雪载荷分布高度函数确定柔性体的承载重量，具体公式为：

其中，t＝s/L，s为弧长，ρ为雪载荷的密度，g为重力加速度，L为长度，H(t)为第一雪载荷分布高度函数。

步骤S3：所述当所述分布模式为梯形堆垛，则根据第二载荷分布公式预测雪载荷作用下柔性体的承载重量，具体包括：

步骤S31：根据第二载荷分布公式确定第二雪载荷分布高度函数H(u)；具体确定过程与确定第二载荷分布公式相类似，具体不再一一赘述。所述第二载荷分布公式为：

其中，H(u)为第二雪载荷分布高度函数，α为雪载荷与柔性体表面作用的最大摩擦角，θ为变形后弧长s处梁截面法向与变形前截面法向的夹角，t＝s/L，s为弧长，L为长度，自重影响因子M＝m/ρL，ρ为雪载荷的密度，m为质量，柯西数CY＝ρgL⁴/B，g为重力加速度，B为柔度，h_c为临界高度，柔度B＝EI/W，临界高度h_c＝tanα·W/2L，E为弹性模量，I是梁截面的惯性矩，W为宽度。

弹性模量E＝mgx₀/(bh³/12)，其中，b为梁截面的宽度，h为梁截面的高度(厚度)，mgx₀为固定端的弯矩。

最后根据所述第二载荷分布公式确定第二雪载荷分布高度函数H(u)。

步骤S32：根据所述第二雪载荷分布高度函数确定柔性体的承载重量，具体公式为：

其中，u＝1-t，t＝s/L，s为弧长，ρ为雪载荷的密度，g为重力加速度，L为长度，H(u)为第二雪载荷分布高度函数。

图5为本发明实施例基于雪载荷作用下柔性体承载重量预测系统结构图，如图5所示，本发明还提供一种基于雪载荷作用下柔性体承载重量预测系统，所述预测系统包括：

分布模式确定模块1，用于确定雪载荷在柔性体表面的分布模式；

第一承载重量确定模块2，用于当所述分布模式为三角形堆垛，则根据第一载荷分布公式预测雪载荷作用下柔性体的承载重量；

第二承载重量确定模块3，用于当所述分布模式为梯形堆垛，则根据第二载荷分布公式预测雪载荷作用下柔性体的承载重量。

所述预测系统还包括：

结构变形度确定模块，用于根据所述柔性体的承载重量确定柔性体的结构变形度，具体公式为：

其中，W_f为柔性体的承载重量，ρ为雪载荷的密度，g为重力加速度，L为长度，α为雪载荷与柔性体表面作用的最大摩擦角。

所述第一承载重量确定模块2，具体包括：

第一雪载荷分布高度函数确定单元，用于根据第一载荷分布公式确定第一雪载荷分布高度函数；所述第一载荷分布公式为：

其中，H(t)为第一雪载荷分布高度函数，α为雪载荷与柔性体表面作用的最大摩擦角，θ为变形后弧长s处梁截面法向与变形前截面法向的夹角，t＝s/L，s为弧长，L为长度，自重影响因子M＝m/ρL，ρ为雪载荷的密度，m为质量，柯西数CY＝ρgL⁴/B，g为重力加速度，B为柔度。

第一承载重量确定单元，用于根据所述第一雪载荷分布高度函数确定柔性体的承载重量，具体公式为：

其中，t＝s/L，s为弧长，ρ为雪载荷的密度，g为重力加速度，L为长度，H(t)为第一雪载荷分布高度函数。

所述第二承载重量确定模块3，具体包括：

第二雪载荷分布高度函数确定单元，用于根据第二载荷分布公式确定第二雪载荷分布高度函数，所述第二载荷分布公式为：

其中，H(u)为第二雪载荷分布高度函数，α为雪载荷与柔性体表面作用的最大摩擦角，θ为变形后弧长s处梁截面法向与变形前截面法向的夹角，t＝s/L，s为弧长，L为长度，自重影响因子M＝m/ρL，ρ为雪载荷的密度，m为质量，柯西数CY＝ρgL⁴/B，g为重力加速度，柔度B＝EI/W，临界高度h_c＝tanα·W/2L，E为弹性模量，I是梁截面的惯性矩，W为宽度。

第二承载重量确定单元，用于根据所述第二雪载荷分布高度函数确定柔性体的承载重量，具体公式为：

其中，u＝1-t，t＝s/L，s为弧长，ρ为雪载荷的密度，g为重力加速度，L为长度，H(u)为第二雪载荷分布高度函数。

图6为本发明实施例三角形堆垛下梁的承受载荷与柯西数、自重因子之间的关系图；如图6所示，图6中(a)表示当α＝30°和M＝0.01时，结构变形度与柯西数之间的关系；图6中(b)表示当α＝30°和M＝0.01、CY＝0.01梁的变形，此时梁基本上是刚性的；图6中(c)表示当α＝30°和M＝0.01、CY＝100梁的变形，此时梁有一定程度的变形；图6中(d)表示当α分别等于30°，45°和60°并且M＝0.01时

随CY的变化曲线；图6中(e)表示当α等于30°，M分别等于0.1,0.01和0.001时

随CY变化的曲线。

图7为本发明实施例梯形堆垛下梁的承受载荷与柯西数、自重因子之间的关系图，如图7所示，图7中(a)表示当α＝30°和M＝0.01时，结构变形度与柯西数之间的关系；图7中(b)表示当α＝30°和M＝0.01、CY＝0.01梁的变形，此时梁基本上是刚性的；图7中(c)表示当α＝30°和M＝0.01、CY＝100梁的变形，此时梁有一定程度的变形；图7中(d)表示当α分别等于30°，45°和60°并且M＝0.01时随CY的变化曲线；图7中(e)表示当α等于30°，M分别等于0.01,0.05和0.001时

随CY变化的曲线。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。