首页 / 专利库 / 物理 / 应力 / 剪切应力 / 一种基于应变能理论的横断型斜裂纹转子变刚度特性计算方法

一种基于应变能理论的横断型斜裂纹转子刚度特性计算方法

阅读:2发布:2021-07-16

专利汇可以提供一种基于应变能理论的横断型斜裂纹转子刚度特性计算方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种基于应变能理论的横断型斜裂纹 转子 变 刚度 特性计算方法,其特征在于,该方法包括如下步骤:1)无裂纹转子轴单元弹性应变能计算;2)无裂纹转子轴单元的位移计算;3)三种类型裂纹转子轴应 力 密度 因子计算;4)裂纹引起的转子轴单元外加弹性应变能计算;5)裂纹引起的转子轴单元外加位移计算;6)裂纹转子轴单元柔度系数计算;7)基于静平衡变换的变刚度特性求解。本发明的有益效果是,方法设计合理,计算精确。,下面是一种基于应变能理论的横断型斜裂纹转子刚度特性计算方法专利的具体信息内容。

1.一种基于应变能理论的横断型斜裂纹转子刚度特性计算方法,其特征在于,该方法包括如下步骤:
1)无裂纹转子轴单元弹性应变能计算;
2)无裂纹转子轴单元的位移计算;
3)三种类型裂纹转子轴应密度因子计算;
4)裂纹引起的转子轴单元外加弹性应变能计算;
5)裂纹引起的转子轴单元外加位移计算;
6)裂纹转子轴单元柔度系数计算;
7)基于静平衡变换的变刚度特性求解。
2.根据权利要求1所述的基于应变能理论的横断型斜裂纹转子变刚度
特性计算方法,其特征在于,所述无裂纹转子轴单元弹性应变能计算式为:
式中M1、M2为弯矩,T为扭矩,F为轴向力。G为刚性
模量,E为杨氏弹性模量。/为裂纹截面惯性矩,/0为裂纹截面极惯矩。
3.根据权利要求1所述的基于应变能理论的横断型斜裂纹转子变刚
度特性计算方法,其特征在于,所述无裂纹转子轴单元的位移计算式为:
4.根据权利要求1所述的基于应变能理论的横断型斜裂纹转子变刚度特性计算方法,其特征在于,所述三种类型裂纹转子轴应力密度因子计算式为:
张开模式S/F:
滑移模式S/F:
剪切模式SIF:
其中
5.根据权利要求1所述的基于应变能理论的横断型斜裂纹转子变刚度特性计算方法,其特征在于,所述裂纹引起的转子轴单元外加弹性应变能计算式为:
式中E′=E/(1-v)且η=1+v。
v为泊松比,KIi|E为裂纹位移张开模式的应力密度因子,KIIi|E为裂纹位移滑移模式的应力密度因子,KIIIi|E为裂纹位移剪切模式的应力密度因子,i=1,...,6。
6.根据权利要求1所述的基于应变能理论的横断型斜裂纹转子变刚
度特性计算方法,其特征在于,裂纹引起的转子轴单元外加位移计算式为:
式中J(A)|E为采用欧拉梁单元建模时根据断裂力
学概念给出的应变能密度函数,由权利要去5所述的式求得。
7.根据权利要求1所述的基于应变能理论的横断型斜裂纹转子变刚度特性计算方法,其特征在于,所述裂纹转子轴单元柔度系数gij|E计算式为:
8.根据权利要求1所述的基于应变能理论的横断型斜裂纹转子变刚度特性计算方法,c T
其特征在于,所述基于静平衡变换的变刚度特性求解式为:[K]|E=[T]G|E[T],式中G|E为柔度矩阵,由[G]|E={gij|E},i,j=1,...,6给出,柔度系数gij|E可参照权利要求7所述内容计算得出,T为逆变矩阵,
c
[K]|E为基于静平衡变换的
变刚度矩阵。

说明书全文

一种基于应变能理论的横断型斜裂纹转子刚度特性计算

方法

技术领域

[0001] 本发明基于应变能理论,提出一种横断型斜裂纹转子变刚度特性计算方法。转子轴分别采用欧拉梁与木辛科梁单元建模,并考虑了纵向、弯曲及扭转六个方向的自由度。根据柔度系数计算裂纹单元的刚度矩阵,而柔度系数依据裂纹单元的应变能理论求得。该计算方法为分析与揭示横断型斜裂纹转子系统振动的非线性特性奠定了理论基础

背景技术

[0002] 裂纹是一种重要的转子系统故障,如果不及时加以监测将导致灾难性后果。现代的大型旋转机械设备中,转子系统通常工作于严酷的热应与机械应力工况下,裂纹往往频繁发生,因此有关裂纹建模与故障诊断问题得到广泛重视和研究。已有的裂纹模型可分[1] [2]为两大类,即常开裂纹 与呼吸裂纹 。常开裂纹将造成转子轴刚度的一个局部非对称常值削弱,适用于拉应力始终作用于裂纹面的情形。在呼吸裂纹情况下,拉应力与压应力随时[3]
间变化交替作用于裂纹面上 ,可以观察到裂纹呈现出不断张开、闭合的呼吸形态,从而引起转子轴刚度的周期变化。目前,有关裂纹建模与故障诊断研究大都针对最简单的横向表[4-8]
面裂纹 ,这限制了研究成果的实际应用。由于问题的复杂性,已知的斜裂纹研究比较少。
Darpe建立了包含斜裂纹的转子有限元模型,并基于断裂力学推导出一个新的柔度系数矩阵,该矩阵中包含由裂纹方向引入的外加应力密度因子。进而,对比了斜裂纹与横向裂纹转[9]
子不同的刚度参数以及耦合振动响应特性 。Sekhar等对比了斜裂纹与横向裂纹转子的[10]
模型特性,并对不同深度、不同位置的裂纹进行了诊断研究 。
[0003] 裂纹源有多种,比如穿晶裂纹、机械疲劳裂纹、汇流裂纹、晶界裂纹、热疲劳裂纹、淬火裂纹、龟裂等,产生原因十分复杂。根据失效分析理论,不同的裂纹源以及不同的材料和形状其裂纹生成与扩展方向均不相同。Darpe与Sekhar研究的是垂直型斜裂纹,即裂纹面垂直于xoy平面,如图2所示。
[0004] 实际上,另一种类型的斜裂纹对于失效分析也是十分重要的,即横断型斜裂纹,裂纹表面横断线垂直于转子轴轴线但与xoy面不垂直,如图3所示。

发明内容

[0005] 本发明的目的是为了解决上述问题,设计了一种基于应变能理论的横断型斜裂纹转子变刚度特性计算方法
[0006] 实现上述目的本发明的技术方案为,一种基于应变能理论的横断型斜裂纹转子变刚度特性计算方法,其特征在于,该方法包括如下步骤:
[0007] 1)无裂纹转子轴单元弹性应变能计算;
[0008] 2)无裂纹转子轴单元的位移计算;
[0009] 3)三种类型裂纹转子轴应力密度因子计算;
[0010] 4)裂纹引起的转子轴单元外加弹性应变能计算;
[0011] 5)裂纹引起的转子轴单元外加位移计算;
[0012] 6)裂纹转子轴单元柔度系数计算;
[0013] 7)基于静平衡变换的变刚度特性求解。
[0014] 所述无裂纹转子轴单元弹性应变能计算式为:
[0015] 式中M1、M2为弯矩,T为扭矩,F为轴向力。G为刚性模量,E为杨氏弹性模量。I为裂纹截面惯性矩,I0为裂纹截面极惯矩。
[0016] 所述无裂纹转子轴单元的位移计算式为:
[0017]
[0018] 所述三种类型裂纹转子轴应力密度因子计算式为:
[0019] 张开模式S/F:
[0020]
[0021] 滑移模式S/F:
[0022]
[0023]
[0024] 剪切模式SIF:
[0025]
[0026] 其中
[0027]
[0028] 所述裂纹引起的转子轴单元外加弹性应变能计算式为:
[0029] 式中E′=E/(1-v)且η=1+v。v为泊松比,KIi|E为裂纹位移张开模式的应力密度因子,KIIi|E为裂纹位移滑移模式的应力密度因子,KIIIi|E为裂纹位移剪切模式的应力密度因子,i=1,...,6。
[0030] 裂纹引起的转子轴单元外加位移计算式为:
[0031] 式中J(A)|E为采用欧拉梁单元建模时根据断裂力学概念给出的应变能密度函数,由权利要去5所述的式求得。
[0032] 所述裂纹转子轴单元柔度系数gij|E计算式为:
[0033]
[0034] 所述基于静平衡变换的变刚度特性求解式为:[K]c|E=[T]G|E[T]T,式中G|E为柔度矩阵,由[G]|E={gij|E},i,j=1,…,6给出,柔度系数gij|E可参照权利要求7所述内容计算得出,T为逆变矩阵,
[0035] [K]c|E为基于静平衡变换的变刚度矩阵。
[0036] 利用本发明的技术方案制作的基于应变能理论的横断型斜裂纹转子变刚度特性计算方法,设计合理,计算结果准确度高,可很好的计算出基于应变能理论的横断型斜裂纹转子变刚度特性。附图说明
[0037] 图1是本发明所述基于应变能理论的横断型斜裂纹转子变刚度特性计算方法的流程示意图:
[0038] 图2垂直型斜裂纹转子受力坐标系统及裂纹面;
[0039] 图3横向型斜裂纹转子受力坐标系统及裂纹面。

具体实施方式

[0040] 下面结合附图对本发明进行具体描述,如图1是本发明所述基于应变能理论的横断型斜裂纹转子变刚度特性计算方法的流程示意图,如图所示,转子轴分别采用欧拉梁单元与铁木辛科梁单元建模两种,并考虑了纵向、弯曲及扭转六个方向的自由度如图3为一段包含横断型斜裂纹的转子轴,转盘质量为m,转轴直径为D,长度为L。裂纹位于左端x处,裂纹深度为a。单元受力情况如下:剪力P2、P3与P8、P9,弯矩P5、P6与P11、P12,轴向力P1与P7,扭矩P4与P10。
[0041] 本技术方案以欧拉梁单元建模转子轴为例子阐述求解基于静平衡变换的变刚度特性的过程,其基本求解原理为:根据柔度系数计算裂纹单元的刚度矩阵,而柔度系数依据裂纹单元的Castingliano应变能理论求得,即
[0042]
[0043] 式中U0为无裂纹转子轴的应变能,Uc为裂纹导致的附加应变能。ui与Pi分别为沿着第i个坐标方向的位移和力。
[0044] 实施例1
[0045] 基于欧拉梁单元建模的横断型斜裂纹转子变刚度特性计算
[0046] 采用欧拉梁单元建模时,忽略了裂纹面上剪切变形,只考虑轴向力、扭矩及弯矩作用。这种情况下无裂纹转子轴单元的弹性应变能U0|E可表达为
[0047]
[0048] 式中M1、M2为弯矩,T为扭矩,F为轴向力。G为刚性模量,E为杨氏弹性模量。I为裂纹截面惯性矩,I0为裂纹截面极惯矩。
[0049] 由图2,有M1=P2x-P6,M2=P3x+P5,T=P4以及F=P1。从而式(2)变为[0050]
[0051] 无裂纹欧拉梁单元位移 为
[0052]
[0053] 裂纹引起的外加位移为
[0054] Uc|E=∫AJ(A)|EdA.
[0055] 式中J(A)|E为采用欧拉梁单元建模时根据断裂力学概念给出的应变能密度函数,由下式给出
[0056]
[0057] 式中E′=E/(1-v)且η=1+v。v为泊松比,JIi|E为裂纹位移张开模式的应力密度因子,KIIi|E为裂纹位移滑移模式的应力密度因子,KIIIi|E为裂纹位移剪切模式的应力密度因子,i=1,...,6。
[0058] 这些应力密度因子(SIF)由下式给出:
[0059] (1)张开模式SIF:
[0060]
[0061] (2)滑移模式SIF:
[0062]
[0063]
[0064] (3)剪切模式SIF:
[0065]
[0066] 其中
[0067]
[0068] 应用这些应力密度因子(SIF)公式以及式(6)的应变能J(A)|E表达式,求得裂纹c引起的外加位移ui|E,i=1,…,6。进而,整体位移ui|E可以通过加上无裂纹转子单元位移
0 0 c
ui|E求得,可以写成矩阵形式ui|E=ui|E+ui|E.=[G]|EPi,其中G|E为柔度矩阵,由[G]|E={gij|E},i,j=1,…,6给出。柔度系数gij|E为
[0069]
[0070] 考虑有限元静平衡条件,采用下式由柔度矩阵通过变换T求解刚度矩阵:T T
[0071] {q1-12} =[T]{q1-6},
[0072] 式中变换矩阵为
[0073]
[0074] 从而横断式斜裂纹转子欧拉梁单元刚度矩阵为
[0075] [K]c|E=[T]G|E[T]T
[0076] 实施例2
[0077] 基于铁木辛科梁单元建模的横断型斜裂纹转子变刚度特性计算
[0078] 对于铁木辛科梁单元,裂纹面上的剪切变形必须加以考虑,从而无裂纹转子轴单元的弹性应变能U0|T可写为
[0079]
[0080] 式中V1、V2为剪切力,且V1=P2,V2=P3。αS为剪切变形因子。0
[0081] 无裂纹欧拉梁单元位移ui|T为
[0082]
[0083] 当分别采用欧拉梁与铁木辛科梁单元进行建模时,应力密度因子(SIF)有所不同,不同之处由下式给出:
[0084] KII2|E=KIII3|E=0.欧拉梁单元
[0085] 铁木辛科梁单元
[0086] 式中k为剪力分布修正因子。其他的应力密度因子是一样的,可参考前面的公式。
[0087] 最终,横断型斜裂纹下铁木辛科梁单元与欧拉梁单元在柔度系数上的差异由下式给出:
[0088] 欧拉梁单元
[0089] 铁木辛科梁单元
[0090] 其他的柔度系数是一样的,可参考前面的公式。计算出柔度系数gij|T后,即可参c照式实施例1中的公式计算横断型斜裂纹转子铁木辛科梁单元的刚度矩阵[K]|T。
[0091] 上述技术方案仅体现了本发明技术方案的优选技术方案,本技术领域的技术人员对其中某些部分所可能做出的一些变动均体现了本发明的原理,属于本发明的保护范围之内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈