CN109726412A - 一种法兰连接螺栓疲劳强度的校核方法 - Google Patents

一种法兰连接螺栓疲劳强度的校核方法 Download PDF

Info

Publication number
CN109726412A
CN109726412A CN201711042093.4A CN201711042093A CN109726412A CN 109726412 A CN109726412 A CN 109726412A CN 201711042093 A CN201711042093 A CN 201711042093A CN 109726412 A CN109726412 A CN 109726412A
Authority
CN
China
Prior art keywords
bolt
load
flange
stress
suffered
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201711042093.4A
Other languages
English (en)
Other versions
CN109726412B (zh
Inventor
邢春雷
严素欣
李云龙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Datang Eryuan Wind Power Co ltd
Original Assignee
BEIJING WANYUAN INDUSTRY CO LTD
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by BEIJING WANYUAN INDUSTRY CO LTD filed Critical BEIJING WANYUAN INDUSTRY CO LTD
Priority to CN201711042093.4A priority Critical patent/CN109726412B/zh
Publication of CN109726412A publication Critical patent/CN109726412A/zh
Application granted granted Critical
Publication of CN109726412B publication Critical patent/CN109726412B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/72—Wind turbines with rotation axis in wind direction

Landscapes

  • Wind Motors (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明属于风力发电技术领域，具体涉及一种法兰连接螺栓疲劳强度的校核方法，是针对现有技术的不足，提出了一种法兰连接螺栓疲劳强度的校核方法。该方法不仅适用风电塔架法兰的连接螺栓，也适用于任意的法兰螺栓连接结构，计算更精确。其特征在于：它包括生成求解文件、得到每个螺栓所受应力与六个载荷分量间的关系式、计算螺栓疲劳损伤。本发明除了可以对风机塔筒法兰连接螺栓进行计算外，可以对任意结构形式的法兰连接螺栓进行计算，适用范围很广。本发明考虑了所有载荷分量对螺栓疲劳寿命的影响，所取得的螺栓所受应力和载荷之间的关系更准确。本发明理论原理清晰，过程简单，便于实施，计算结果精度更高。

Description

一种法兰连接螺栓疲劳强度的校核方法
技术领域
本发明属于风力发电技术领域，具体涉及一种法兰连接螺栓疲劳强度的校核方法。
背景技术
目前风力发电机组塔架多采用全钢塔架和混凝土塔筒——钢塔组成的混合塔架两种形式，全钢塔架和混合塔架的钢塔部分由钢制塔筒构成，各段钢制塔筒间通过钢制法兰和高强螺栓进行连接。法兰连接螺栓的安全与否直接关系到整个机组的生死存亡。法兰连接螺栓的疲劳强度校核工程上普遍采用基于Petersen方法的Schmidt-Neuper方法进行校核，由于Petersen方法中没有考虑预紧力对螺栓变形的影响、没有考虑弯矩对螺栓的影响、只适用于L型和T型法兰连接，Petersen方法对螺栓的计算采用较大的安全系数来弥补理论缺陷。
Schmidt-Neuper方法基于Petersen方法对螺栓极限强度校核结果基础上，只考虑My和Fz两个载荷分量对螺栓应力的影响来计算螺栓的疲劳损伤，又进一步加大了计算误差。
发明内容
本发明的目的是针对现有技术的不足，提出了一种法兰连接螺栓疲劳强度的校核方法。该方法不仅适用风电塔架法兰的连接螺栓，也适用于任意的法兰螺栓连接结构，计算更精确。
本发明是这样实现的：
一种法兰连接螺栓疲劳强度的校核方法，包括如下步骤：
步骤一：将法兰连接结构的装配体在Hypermesh软件中进行处理，建立有限元建模，其中与法兰连接的塔筒假体且长度不小于法兰外径的两倍，上法兰和下法兰间设置为摩擦接触关系，高强垫圈与上法兰之间和高强垫圈与下法兰之间均设置为摩擦接触关系，高强垫圈与螺栓的螺母之间和高强垫圈与螺栓的螺帽之间采用绑定接触方式进行连接，下法兰一侧的塔筒假体设置固定边界约束条件，上法兰一侧的塔筒假体设置MPC约束，在上、下法兰结合面中心分别施加Fx、Fy、Fz、Mx、My、Mz的单位载荷，生成六个求解文件，所施加的单位载荷的单位与法兰中心的时序载荷的单位一致；
步骤二：将求解文件提交Abaqus进行计算，对螺栓进行编号，并提取每个载荷分量下各个螺栓的应力，根据同一个螺栓在六个载荷分量下各自的应力结果，写出每个螺栓所受应力与六个载荷分量间的关系式；
步骤三：在Bladed软件中将上法兰和下法兰结合面中心的时序载荷代入步骤二中螺栓所受应力与的六个载荷分量间的关系式，得出每个螺栓所受应力与时间之间的关系，然后采用Bladed软件进行雨流计数和损伤累积计算，得到每个螺栓的疲劳损伤。
如上所述的步骤二中，每个螺栓所受应力与六个载荷分量间的关系式为：
σ＝Fx*σFx+Fy*σFy+Fz*σFz+Mx*σMx+My*σMy+Mz*σMz
其中：σ是螺栓所受应力，
Fx为载荷分量，σFx为Fx作用下螺栓所受应力，
Fy为载荷分量，σFy为Fx作用下螺栓所受应力，
Fz为载荷分量，σFz为Fx作用下螺栓所受应力，
Mx为载荷分量，σMx为Fx作用下螺栓所受应力，
My为载荷分量，σMy为Fx作用下螺栓所受应力，
Mz为载荷分量，σMz为Fx作用下螺栓所受应力。
本发明的有益效果是：
本发明除了可以对风机塔筒法兰连接螺栓进行计算外，可以对任意结构形式的法兰连接螺栓进行计算，适用范围很广。本发明考虑了所有载荷分量对螺栓疲劳寿命的影响，所取得的螺栓所受应力和载荷之间的关系更准确。本发明理论原理清晰，过程简单，便于实施，计算结果精度更高。
附图说明
图1是本发明中法兰的连接结构示意图。
图中：1.螺母，2.螺栓，3.高强垫圈，4.上法兰，5.下法兰，6.塔筒假体。
具体实施方式
下面结合附图和实施例对本发明进行进一步描述。
本实施例中利用Hypermesh进行有限元前处理操作，利用Abaqus进行单位载荷求解，利用Bladed软件进行雨流计数以及损伤累积计算。
一种法兰连接螺栓疲劳强度的校核方法，包括如下步骤：
步骤一：将法兰连接结构的装配体在Hypermesh软件中进行处理，建立有限元建模，其中与法兰连接的塔筒假体且长度不小于法兰外径的两倍，上法兰和下法兰间设置为摩擦接触关系，高强垫圈与上法兰之间和高强垫圈与下法兰之间均设置为摩擦接触关系，高强垫圈与螺栓的螺母之间和高强垫圈与螺栓的螺帽之间采用绑定接触方式进行连接，下法兰一侧的塔筒假体设置固定边界约束条件，上法兰一侧的塔筒假体设置MPC约束，在上、下法兰结合面中心分别施加Fx、Fy、Fz、Mx、My、Mz的单位载荷，生成六个求解文件，所施加的单位载荷的单位与法兰中心的时序载荷的单位一致；
步骤二：将求解文件提交Abaqus进行计算，对螺栓进行编号，并提取每个载荷分量下各个螺栓的应力，根据同一个螺栓在六个载荷分量下各自的应力结果，写出每个螺栓所受应力与六个载荷分量间的关系式；
步骤三：在Bladed软件中将上法兰和下法兰结合面中心的时序载荷代入步骤二中螺栓所受应力与的六个载荷分量间的关系式，得出每个螺栓所受应力与时间之间的关系，然后采用Bladed软件进行雨流计数和损伤累积计算，得到每个螺栓的疲劳损伤。
如上所述的步骤二中，每个螺栓所受应力与六个载荷分量间的关系式为：
σ＝Fx*σFx+Fy*σFy+Fz*σFz+Mx*σMx+My*σMy+Mz*σMz
其中：σ是螺栓所受应力，
Fx为载荷分量，σFx为Fx作用下螺栓所受应力，
Fy为载荷分量，σFy为Fx作用下螺栓所受应力，
Fz为载荷分量，σFz为Fx作用下螺栓所受应力，
Mx为载荷分量，σMx为Fx作用下螺栓所受应力，
My为载荷分量，σMy为Fx作用下螺栓所受应力，
Mz为载荷分量，σMz为Fx作用下螺栓所受应力。
上面结合实施例对本发明的实施方法作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。本发明说明书中未作详细描述的内容均可以采用现有技术。

Claims (2)

1.一种法兰连接螺栓疲劳强度的校核方法，包括如下步骤：
步骤一：将法兰连接结构的装配体在Hypermesh软件中进行处理，建立有限元建模，其中与法兰连接的塔筒假体(6)且长度不小于法兰外径的两倍，上法兰(4)和下法兰(5)间设置为摩擦接触关系，高强垫圈(3)与上法兰(4)之间和高强垫圈(3)与下法兰(5)之间均设置为摩擦接触关系，高强垫圈(3)与螺栓(2)的螺母(1)之间和高强垫圈与螺栓(2)的螺帽之间采用绑定接触方式进行连接，下法兰(5)一侧的塔筒假体(6)设置固定边界约束条件，上法兰(4)一侧的塔筒假体(6)设置MPC约束，在上、下法兰结合面中心分别施加Fx、Fy、Fz、Mx、My、Mz的单位载荷，生成六个求解文件，所施加的单位载荷的单位与法兰结合面中心的时序载荷的单位一致；
步骤二：将求解文件提交Abaqus进行计算，对螺栓(2)进行编号，并提取每个载荷分量下各个螺栓(2)的应力，根据同一个螺栓(2)在六个载荷分量下各自的应力结果，写出每个螺栓(2)所受应力与六个载荷分量间的关系式；
步骤三：在Bladed软件中将上法兰(4)和下法兰(5)结合面中心的时序载荷代入步骤二中螺栓(2)所受应力与的六个载荷分量间的关系式，得出每个螺栓(2)所受应力与时间之间的关系，然后采用Bladed软件进行雨流计数和损伤累积计算，得到每个螺栓的疲劳损伤。
2.根据权利要求1所述的法兰连接螺栓疲劳强度的校核方法，其特征在于：所述的步骤二中，每个螺栓(2)所受应力与六个载荷分量间的关系式为：
σ＝Fx*σFx+Fy*σFy+Fz*σFz+Mx*σMx+My*σMy+Mz*σMz
其中：σ是螺栓所受应力，
Fx为载荷分量，σFx为Fx作用下螺栓所受应力，
Fy为载荷分量，σFy为Fx作用下螺栓所受应力，
Fz为载荷分量，σFz为Fx作用下螺栓所受应力，
Mx为载荷分量，σMx为Fx作用下螺栓所受应力，
My为载荷分量，σMy为Fx作用下螺栓所受应力，
Mz为载荷分量，σMz为Fx作用下螺栓所受应力。
CN201711042093.4A 2017-10-31 2017-10-31 一种法兰连接螺栓疲劳强度的校核方法 Active CN109726412B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711042093.4A CN109726412B (zh) 2017-10-31 2017-10-31 一种法兰连接螺栓疲劳强度的校核方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711042093.4A CN109726412B (zh) 2017-10-31 2017-10-31 一种法兰连接螺栓疲劳强度的校核方法

Publications (2)

Publication Number Publication Date
CN109726412A true CN109726412A (zh) 2019-05-07
CN109726412B CN109726412B (zh) 2023-01-13

Family

ID=66292993

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711042093.4A Active CN109726412B (zh) 2017-10-31 2017-10-31 一种法兰连接螺栓疲劳强度的校核方法

Country Status (1)

Country Link
CN (1) CN109726412B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110321599A (zh) * 2019-06-12 2019-10-11 许昌许继风电科技有限公司 一种t型法兰地锚螺栓疲劳强度校核方法、装置及介质
CN110348151A (zh) * 2019-07-17 2019-10-18 国电联合动力技术有限公司 风电机组塔筒法兰连接螺栓设计方法及装置
CN110991111A (zh) * 2019-11-27 2020-04-10 南京安维士传动技术股份有限公司 一种基于摩擦接触的风电齿轮箱行星架的疲劳计算方法
CN112699509A (zh) * 2021-01-08 2021-04-23 北京工业大学 一种法兰连接螺栓极限强度的校核方法
CN113849937A (zh) * 2021-10-14 2021-12-28 哈电发电设备国家工程研究中心有限公司 一种由附加力矩引起的螺栓应力计算方法、设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2541457A1 (fr) * 1982-09-23 1984-08-24 Charbonnages Ste Chimique Appareillage pour le controle de la resistance a la pliure de feuilles thermoplastiques
CN103353907A (zh) * 2013-06-17 2013-10-16 沈阳华创风能有限公司 一种法兰连接螺栓极限强度校核的计算方法
CN105022868A (zh) * 2015-07-02 2015-11-04 许继集团有限公司 风机偏航轴承连接螺栓极限强度和疲劳强度的计算方法
CN105160067A (zh) * 2015-08-05 2015-12-16 许继集团有限公司 风机变桨轴承与轮毂和叶片连接螺栓强度性能的校核方法
CN106570326A (zh) * 2016-10-28 2017-04-19 许继集团有限公司 一种l型法兰连接螺栓疲劳强度校核的计算方法及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2541457A1 (fr) * 1982-09-23 1984-08-24 Charbonnages Ste Chimique Appareillage pour le controle de la resistance a la pliure de feuilles thermoplastiques
CN103353907A (zh) * 2013-06-17 2013-10-16 沈阳华创风能有限公司 一种法兰连接螺栓极限强度校核的计算方法
CN105022868A (zh) * 2015-07-02 2015-11-04 许继集团有限公司 风机偏航轴承连接螺栓极限强度和疲劳强度的计算方法
CN105160067A (zh) * 2015-08-05 2015-12-16 许继集团有限公司 风机变桨轴承与轮毂和叶片连接螺栓强度性能的校核方法
CN106570326A (zh) * 2016-10-28 2017-04-19 许继集团有限公司 一种l型法兰连接螺栓疲劳强度校核的计算方法及装置

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110321599A (zh) * 2019-06-12 2019-10-11 许昌许继风电科技有限公司 一种t型法兰地锚螺栓疲劳强度校核方法、装置及介质
CN110321599B (zh) * 2019-06-12 2023-01-17 许昌许继风电科技有限公司 一种t型法兰地锚螺栓疲劳强度校核方法、装置及介质
CN110348151A (zh) * 2019-07-17 2019-10-18 国电联合动力技术有限公司 风电机组塔筒法兰连接螺栓设计方法及装置
CN110991111A (zh) * 2019-11-27 2020-04-10 南京安维士传动技术股份有限公司 一种基于摩擦接触的风电齿轮箱行星架的疲劳计算方法
CN110991111B (zh) * 2019-11-27 2024-01-30 南京安维士传动技术股份有限公司 一种基于摩擦接触的风电齿轮箱行星架的疲劳计算方法
CN112699509A (zh) * 2021-01-08 2021-04-23 北京工业大学 一种法兰连接螺栓极限强度的校核方法
CN113849937A (zh) * 2021-10-14 2021-12-28 哈电发电设备国家工程研究中心有限公司 一种由附加力矩引起的螺栓应力计算方法、设备及存储介质
CN113849937B (zh) * 2021-10-14 2024-06-21 哈电发电设备国家工程研究中心有限公司 一种由附加力矩引起的螺栓应力计算方法、设备及存储介质

Also Published As

Publication number Publication date
CN109726412B (zh) 2023-01-13

Similar Documents

Publication Publication Date Title
CN109726412A (zh) 一种法兰连接螺栓疲劳强度的校核方法
Xia et al. The finite element method for dynamics of FG porous truncated conical panels reinforced with graphene platelets based on the 3-D elasticity
CN105160067A (zh) 风机变桨轴承与轮毂和叶片连接螺栓强度性能的校核方法
CN110795875B (zh) 一种风机轴-齿轮箱法兰连接校核方法
CN107229765A (zh) 分析计算铁塔每根杆件应力的输电铁塔杆件应力计算方法
CN106570326A (zh) 一种l型法兰连接螺栓疲劳强度校核的计算方法及装置
CN109299532A (zh) 一种风机主机架与后机架连接螺栓强度校核方法及系统
Zheng et al. Time-domain fatigue assessment for blade root bolts of floating offshore wind turbine (FOWT)
Qin et al. Design and nonlinear structural responses of multi-bolted joint composite box-beam for sectional wind turbine blades
CN107563085A (zh) 一种高强螺栓连接结构的有限元建模方法及强度评估方法
Paquette et al. Increased Strength in Wind Turbine Blades through Innovative Structural Design.
CN109726410B (zh) 风力机偏航连接螺栓强度校核的计算方法
Rao Advances in aero structures
CN109241546B (zh) 一种风机塔筒连接法兰疲劳强度的校核方法
Qin et al. Design and structural responses of a 38-meter sectional wind turbine blade under extreme static loads
CN113722950A (zh) 一种风力机t型法兰锚栓疲劳强度校核方法
CN110954330A (zh) 一种航空发动机转子紧固界面切向接触阻尼预测方法
CN110937119A (zh) 一种一体式全复合材料连杆结构
Tian et al. Analysis of the nonlinear characteristics of the bonding surfaces of the bolted structure
He et al. The influence of constraint method on the finite element simulation for bolted joint
CN117034711B (zh) 一种基于风力机叶片有限元模型的分布力加载方法
Inomata et al. Measurement of stress on blade of NEDO's 500 kW prototype wind turbine
CN109932241B (zh) 塑性条件下耦合残余应力和拘束效应的蠕变孕育期预测方法
Weber et al. Fatigue Life calculation using winLIFE
Kong et al. Structural Analysis of Horizontal Tidal Current Turbine Main Shaft

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20221229

Address after: 671000 Luoping Mountain, Eryuan County, Dali Bai Autonomous Prefecture, Yunnan Province

Applicant after: Datang Eryuan Wind Power Co.,Ltd.

Address before: 100176 Aerospace Science Park, No. 6, Jinxiu Street, Beijing Economic and Technological Development Zone, Daxing District, Beijing

Applicant before: BEIJING WANYUAN INDUSTRY Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant