[发明专利]一种基于直接微分法的输电导线响应敏感性分析方法

权利要求书

1.一种基于直接微分法的输电导线响应敏感性分析方法，其特征是包括以下步骤：

步骤S1，建立基于DDM直接微分法的输电导线宏观结构响应敏感性计算方法

直接微分法计算过程是在每一时步计算收敛后，对离散后的结构动力平衡方程直接求导，得到公式：

M(θ)u··(t,θ)+C(θ)u·(t,θ)+R(u(t,θ),θ)=F(t,θ)---(1);]]>

式中t为时间，θ为敏感性参数，u(t,θ)为有限元宏观响应，即节点位移，M表示质量，C表示阻尼矩阵，R(u(t,θ),θ)为与应力历史有关的结构反力，F(t,θ)为外力；公式中上标点表示对时间求一阶导，两点表示二阶；

对公式(1)使用Newmark-β法对时间离散得离散后的动力平衡方程：

Ψ(un+1)=F~n+1-[1β(Δt)2Mun+1+αβ(Δt)Cun+1+R(un+1)]=0---(2);]]>

其中公式(2)中Δt表示时间增量，表示n+1时步时的外力,α和β是NewMark积分常数；

取α＝1/2，β＝1/4,即为无条件稳定时，Ψ(u_n+1)为残值；

F~n+1=Fn+1+M[1β(Δt)2un+1β(Δt)u·n-(1-12β)u··n]+C[αβ(Δt)un-(1-αβ)u·n-(Δt)(1-α2β)u··n]---(3);]]>

公式(3)对θ求导，得敏感性分析基本方程：

[1β(Δt)2M+αβ(Δt)C+(KTstat)n+1]dun+1dθ=-(1β(Δt)2dMdθ+αβ(Δt)dCdθ)un+1-∂R(un+1(θ),θ)∂θ|un+1+dF~n+1dθ---(4);]]>

公式(4)中代表第n+1时步的刚度矩阵，为宏观位移响应对于θ的敏感性或者梯度，为需要求解的量；

为结构内力敏感性，表示节点位移固定；

结构内力敏感性计算方法如下：

∂(R(un+1(θ)),θ)∂θ|un+1=Ae=1Nel{∫Ω(e)BT∂(σ(ϵn+1(θ)),θ)∂θ|un+1dΩ(e)}---(5);]]>

其中表示对从e＝1到e＝Nel时的各项集成，B表示应变和应力的关系，B^T表示B的转置；

敏感性计算框架中有两个步奏：第一步是求固定位移的内力梯度第二步为求位移不固定条件下导数此项为下一时步计算的基础；

求解公式(4)即可得到宏观位移响应的敏感性；

步骤S2，建立基于DDM直接微分法的输电导线细观结构响应敏感性计算方法

公式(5)中即为输电导线细观分层应力的敏感性；由于输电导线分层应力是从每根输电导线丝集成得到的；输电导线从内到外的第i层每一根输电导线丝的内力为：

p_i＝σ_i(ε_i-ε_iT)A_i(6)；

式中，σ_i(ε_i-ε_iT)指第i层输电导线丝的应力σ_i是第i层输电导线应变ε_i和温度应变ε_iT的函数，A_i表示第i层导线面积；

将这些力在绞线轴向投影叠加后得到输电导线截面的内力表达式：

P(ϵ0)=Σi=1npinicosβi---(7);]]>

式中，P表示宏观下绞线截面上受力，n_i表示第i层输电导线丝的数量,β_i表示第i层绞线的旋转角；

因此输电导线的等效截面应力为：

σ=P(ϵ0)A---(8);]]>

将上述公式(6)-(8)对θ求导：

∂pi∂θ=∂σi∂θ(ϵi-ϵiT)Ai+∂(ϵi-ϵiT)∂θσiAi+∂Ai∂θσi(ϵi-ϵiT)---(9);]]>

∂P(ϵ0)∂θ=Σi=1n(∂pi∂θnicosβi+∂cosβi∂θpini)---(10);]]>

∂σ∂θ=1(A)2(∂P(ϵ0)∂θA-∂A∂θP(ϵ0))---(11);]]>

即得到细观输电导线响应的计算公式；

步骤S3，考虑塑性和屈服行为的输电导线非线性响应敏感性计算方法

当输电导线丝受拉达到屈服极限后，输电导线丝进入非线性阶段，输电导线丝非线性力学行为用一维von Mises J₂塑性材料模型模拟：

弹性预测

Δλ^Trial＝0 (12)；

(ϵn+1p)Trial=ϵnp---(13);]]>

αn+1Trial=αn---(14);]]>

(ϵ‾n+1p)Trial=ϵ‾np---(15);]]>

σn+1Trial=E(ϵn+1-ϵnp)---(16);]]>

σy,n+1Trial=σy,n---(17);]]>

DDM方法分析弹性截断敏感性

∂∂θ(Δλ)=0---(18);]]>

∂ϵn+1p∂θ=∂ϵnp∂θ---(19);]]>

∂αn+1Trial∂θ=∂αn∂θ---(20);]]>

∂ϵ‾n+1p∂θ=∂ϵ‾np∂θ---(21);]]>

∂σn+1∂θ=E(∂ϵn+1∂θ-∂ϵnp∂θ)+∂E∂θ(ϵn+1-ϵnp)---(22);]]>

∂σy,n+1∂θ=∂σy,n∂θ---(23);]]>

塑性纠正阶段

Δλ=|σn+1Trial-αn+1Trial|-σy,nE+Hiso+Hkin---(24);]]>

σn+1=σn+1Trial-E×Δλ×sgn(σn+1-αn+1)---(25);]]>

αn+1=αn+1Trial+Hkin×Δλ×sgn(σn+1-αn+1)---(26);]]>

σ_y,n+1＝σ_y,n+H_iso×Δλ (27)；

∂(Δλ)∂θ=(E+Hiso+Hkin)×((∂σn+1Trial∂θ-∂αn+1Trial∂θ)×nn+1-∂σy,n∂θ)(E+Hiso+Hkin)2-(∂E∂θ+∂Hiso∂θ+∂Hkin∂θ)×((σn+1Trial-αn+1Trial)×nn+1-σy,n)(E+Hiso+Hkin)2---(28);]]>

其中

∂σn+1∂θ=E(∂ϵn+1∂θ-∂ϵnp∂θ)+∂E∂θ(ϵn+1-ϵnp)---(30);]]>

∂αn+1∂θ=∂αn∂θ+∂Hkin∂θ·(Δλ)·nn+1+Hkin·∂(Δλ)∂θ·nn+1---(31);]]>

∂σy,n+1∂θ=∂σy,n∂θ+∂Hiso∂θ×Δλ+Hiso×(Δλ)∂θ---(32);]]>

以上公式中λ为一致性参数，而是离散连续的参数，中上标Trail表示试算、下标n+1表示处于时间t_n+1时的变量，(…)_y下标y表示处于屈服时的变量，另外α表示背应力，ε^e表示弹性部分的应变，ε^p表示塑性部分的应变，E表示材料的杨氏模量，公式(28)中n_n+1代表(σ_n+1-α_n+1)的正负号，其下标n+1表示在n+1时步；

定义屈服函数如下：

f(σ,α,ϵ‾p)=|σ-α|-(σy0+Hisoϵ‾p)---(33);]]>

若{f≤0}则材料为线弹性，公式(12)到公式(17)即为计算结果，否则进入塑性纠正阶段；

公式(24)中H_kin和H_iso表示等向硬化模型参数，公式(25)中sgn(x)表示x的符号；

步骤S4，建立输电导线响应对于其材料、尺寸、外荷载参数的敏感性计算方法

S4-1，输电导线响应关于材料敏感性分析

从公式(5)有

∂(R(un+1(θ)),θ)∂θ|un+1=Ae=1Nel{∫Ω(e)BT∂(σ(ϵn+1(θ)),θ)∂θ|un+1dΩ(e)};]]>

基于材料本构关系，计算

比如上式表示对从e＝1到e＝Nel时的各项集成，公式(24)和(32)可以计算弹塑性材料的应力敏感性；将应力敏感性代回到公式(5)则有即为输电导线内力响应对材料的敏感性；

S4-2，输电导线响应关于几何尺寸敏感性

当敏感性参数θ为坐标x时，不可忽略，x为输电导线位置的坐标，其值大小为输电导线内力对几何位置和尺寸的敏感性；另外，输入的参数和根据已知能确定的初始条件在导线响应敏感性分析过程中即为边界条件；

S4-3，输电导线响应关于外力的敏感性

公式(4)中项中包括了外力的敏感性，当敏感性参数θ为外力时，计算结果为输电导线关于外力的敏感性。

2.根据权利要求1所述的基于直接微分法的输电导线响应敏感性分析方法，其特征是：提出用扰动法计算输电导线响应敏感性的参数扰动范围指标：

对于响应敏感性有

∂r∂θ≈ΔrΔθ---(34);]]>

公式(34)中r表示宏观响应或者细观响应，响应敏感性用有限差分法近似。