CN112710539B - 含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法 - Google Patents
含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法 Download PDFInfo
- Publication number
- CN112710539B CN112710539B CN202011527291.1A CN202011527291A CN112710539B CN 112710539 B CN112710539 B CN 112710539B CN 202011527291 A CN202011527291 A CN 202011527291A CN 112710539 B CN112710539 B CN 112710539B
- Authority
- CN
- China
- Prior art keywords
- fatigue
- test piece
- test
- main beam
- normalized
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000007547 defect Effects 0.000 title claims abstract description 40
- 238000012360 testing method Methods 0.000 claims abstract description 128
- 238000009661 fatigue test Methods 0.000 claims abstract description 31
- 230000008569 process Effects 0.000 claims abstract description 18
- 238000010586 diagram Methods 0.000 claims abstract description 7
- 239000000463 material Substances 0.000 claims description 28
- 230000037303 wrinkles Effects 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 10
- 239000003365 glass fiber Substances 0.000 claims description 7
- 230000006641 stabilisation Effects 0.000 claims description 7
- 238000011105 stabilization Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 5
- 238000009825 accumulation Methods 0.000 claims description 4
- 239000003822 epoxy resin Substances 0.000 claims description 4
- 229920000647 polyepoxide Polymers 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 102100040287 GTP cyclohydrolase 1 feedback regulatory protein Human genes 0.000 claims 1
- 101710185324 GTP cyclohydrolase 1 feedback regulatory protein Proteins 0.000 claims 1
- 239000002131 composite material Substances 0.000 description 7
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 239000000805 composite resin Substances 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- LNEPOXFFQSENCJ-UHFFFAOYSA-N haloperidol Chemical compound C1CC(O)(C=2C=CC(Cl)=CC=2)CCN1CCCC(=O)C1=CC=C(F)C=C1 LNEPOXFFQSENCJ-UHFFFAOYSA-N 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 239000012779 reinforcing material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 238000001931 thermography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/02—Details
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0016—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings of aircraft wings or blades
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/02—Details
- G01N3/04—Chucks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/32—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2240/00—Components
- F05B2240/20—Rotors
- F05B2240/30—Characteristics of rotor blades, i.e. of any element transforming dynamic fluid energy to or from rotational energy and being attached to a rotor
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2260/00—Function
- F05B2260/82—Forecasts
- F05B2260/821—Parameter estimation or prediction
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2260/00—Function
- F05B2260/83—Testing, e.g. methods, components or tools therefor
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2260/00—Function
- F05B2260/84—Modelling or simulation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/303—Temperature
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
- F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
- F05B2270/00—Control
- F05B2270/30—Control parameters, e.g. input parameters
- F05B2270/332—Maximum loads or fatigue criteria
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0001—Type of application of the stress
- G01N2203/0005—Repeated or cyclic
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0014—Type of force applied
- G01N2203/0016—Tensile or compressive
- G01N2203/0017—Tensile
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0073—Fatigue
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/0641—Indicating or recording means; Sensing means using optical, X-ray, ultraviolet, infrared or similar detectors
-
- 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
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Health & Medical Sciences (AREA)
- Sustainable Energy (AREA)
- Sustainable Development (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Wind Motors (AREA)
- Fluid Mechanics (AREA)
Abstract
本发明涉及风力发电机叶片的疲劳寿命预测技术领域,尤其涉及一种含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法。该预测方法更接近试验和实际情况,预测结果更加科学而准确。包括步骤:S1、对待测的含褶皱缺陷的主梁试件进行拉伸性能测试;S2、根据S1获得的试件表面温度数据,计算主梁试件在不同加载应力下固有耗散能量;S3、绘制试件的固有耗散与相应的UTS水平的关系图;S4:基于疲劳过程中的固有耗散变化建立含待确定参数的归一化剩余刚度模型,将疲劳试验数据代入模型;S5:对于含褶皱缺陷的主梁试件,根据确定参数后的归一化剩余刚度模型推导出疲劳寿命预测模型;S6:求出归一化失效刚度。
Description
技术领域
本发明涉及风力发电机叶片的疲劳寿命预测技术领域,同时也适用于其他树脂基复合材料的疲劳寿命预测的方法;尤其涉及一种含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法。
背景技术
风能是一种绿色、清洁、开发和利用成本相对较低的可再生能源之一,已成为我国最具广阔商业开发潜力和发展前景的新能源。风能的主要利用方式是通过风力发电机将风能转换为电能。风力发电机是典型的机电一体化系统,叶片是风力机接收风能最重要的组件,其成本约占风力机总成本的15-20%。然而,叶片长期承受着恶劣的工作环境和复杂的交变载荷,在运行工况下疲劳损伤逐渐累积,对风电场的安全运行造成极大隐患与威胁。为了确保风力机可以满足20年的使用寿命,在设计和生产制造之前对叶片的疲劳性能进行有效的评估对于提高叶片的可靠性、经济性和避免叶片的灾难性事故具有重要意义。
目前,风力机叶片制造普遍采用的是以玻璃纤维为增强材料、以树脂为基体的玻璃纤维增强复合材料(GFRP)。主梁是叶片的主要承力结构,是引发叶片断裂事故的主要部位。大多数主梁的主要结构采用单向GFRP层合板,以00的玻璃纤维作为铺层。褶皱缺陷已被风能行业视为一种在主梁结构中经常出现,但在制造过程中难以预防、难以检测且可能引起意外故障的缺陷。具有缺陷的主梁材料在交变载荷的作用下,疲劳损伤会加速累积,导致使用寿命缩短。目前,对不含缺陷的复合材料进行疲劳寿命预测的相关研究较多,然而这类模型不适用于预测含有缺陷的复合材料的疲劳寿命。因此,一个能够准确预测含褶皱缺陷的主梁材料的疲劳寿命模型是十分必要的。如果检测到有不满足其使用寿命的,可提前有规划地进行叶片的更换和维护,有效减少叶片的运营和维护成本。
目前,常用的疲劳试验研究方法是通过在不同的载荷水平下测试大量的试件,记录和统计试件在不同载荷下的疲劳寿命,进而评估其疲劳性能。但是传统的试验方法存在着试验周期长、试件消耗大、数据分散等一系列不足。由于材料的疲劳损伤是伴随温度变化的能量耗散过程,近年来,一系列基于试件在疲劳载荷下能量耗散的复合材料疲劳性能分析检测方法被提出,被称为能量方法。红外热像技术作为一种常用的温度无损检测手段,已经被用于从能量的角度来建立疲劳寿命预测模型。与传统试验方法相比,能量方法在很大程度上缩短了试验周期,减少了试验材料的数量,是一种快速、经济、准确的疲劳寿命预测技术。
以能量耗散因素作为疲劳损伤指标的模型,本质上反映了复合材料疲劳失效过程的外在表现形式,但是没有考虑疲劳的物理性能变化。复合材料疲劳过程中在宏观上除了会导致能量耗散的变化外,还会引起材料的强度、刚度、泊松比和阻尼系数等物理性能的变化。以复合材料疲劳过程中力学性能影响为基础的寿命预测模型的研究大多是基于剩余刚度因素的损伤变量。然而,现有的寿命预测模型没有能够建立材料的疲劳性能退化与其宏观可测力学性能变化之间的关系,都仅仅考虑了一个因素的损伤指标。如果可以综合考虑两个以上因素各自的优势,就可以建立一个全面准确的含褶皱缺陷的主梁试件的疲劳寿命预测模型,具有重大的工程应用价值。
发明内容
本发明就是针对现有技术存在的缺陷,提供一种含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法,其在寿命预测过程中,分析了主梁内部缺陷对疲劳寿命的影响,充分考虑了主梁试件疲劳过程中外部固有能量耗散与剩余刚度变化的复杂性,可以在短时间内精确快速地预测疲劳寿命。本发明的预测方法更接近试验和实际情况,预测结果更加科学而准确。
为实现上述目的,本发明采用如下技术方案,包括以下步骤:
S1、对待测的含褶皱缺陷的主梁试件在万能疲劳试验机上进行拉伸性能测试,测出其最大载荷应力UTS;疲劳测试时使用红外热像仪同步观察和记录试件表面的温度变化,试验结束后统计并记录试件在不同循环次数下的表面温度、应力、应变和刚度数据;
S2、根据S1获得的试件表面温度数据,计算主梁试件在不同加载应力下温度稳定阶段时的固有耗散能量;
S3、绘制试件的固有耗散dista与相应的UTS水平的关系图,利用双直线方法对关系图中具有两种不同斜率的数据趋势线进行插值,由两条直线交点的σmax坐标值确定出试件的疲劳极限;把疲劳极限以上载荷认定为过载区域,即试件发生失效的区域;而低于疲劳极限以下载荷的区域在试验过程中没有失效试件;
步骤S4:基于疲劳过程中的固有耗散变化建立含待确定参数的归一化剩余刚度模型,将步骤S1、S2记录的疲劳试验数据代入模型,确定模型中未知参数a,b,p,q的值;a和b是与材料内部褶皱缺陷有关的参数;p和q是独立于固有耗散能量和加载周期的两个材料参数;
步骤S5:对于含褶皱缺陷的主梁试件,根据确定参数后的归一化剩余刚度模型推导出疲劳寿命预测模型;
步骤S6:定义试件可以承受最大疲劳载荷的循环次数,以此求出归一化失效刚度;再将确定归一化失效刚度后的疲劳寿命预测模型用于待测试件在过载区域内的疲劳寿命预测,以获得含褶皱缺陷的主梁试件的预测S-N曲线。
进一步地,S1步骤中,疲劳测试的最大加载应力在UTS的20%到90%之间以5%为步长分别对试件施加;每次仅在8000个加载周期内对试件进行测试。
进一步地,所述步骤S1中的疲劳试验机为电动液压伺服万能疲劳试验机(MTS810);进行拉伸性能试验,拉伸速率为2毫米/分钟;进行疲劳试验,试验机的疲劳载荷为拉-拉载荷,应力比为0.1,振动频率为10Hz。
进一步地,所述步骤S2中计算主梁试件在温度稳定阶段时的固有耗散能量的步骤包括:
主梁承受均匀的单轴循环载荷时,内变量和温度之间的内部耦合源可以忽略;在疲劳过程中温度变化的稳定阶段,当试件的每一周循环结束时,热损失与固有耗散能量之间存在平衡,该阶段稳定的固有耗散计算模型可表达为:
distab=-div(kgradTstab) (1)
其中,Tstab是试件在表面温度稳定阶段下的温度,其单位为℃,k是导热系数,其单位是W·m-1·K-1;div(kgradTstab)表示由传导引起的热损失率。
进一步地,所述步骤S4中归一化剩余刚度模型的建立过程如下:
步骤S41:在疲劳试验过程中,疲劳损伤指数D(n)通过疲劳载荷下的剩余刚度表示如下:
其中,n代表循环次数;E(n)代表第n个循环次数相对应的刚度;E0是沿x轴方向的初始有效刚度;E’(n)是第n个循环次数相应的归一化刚度,定义为E(n)/E0;Ef代表失效刚度;E′f是最终循环Nf时的归一化失效刚度,定义为Ef/E0;
步骤S42:由于稳定的固有耗散能量值distab与损伤累积之间存在一定的关系,通过引入褶皱缺陷中高宽比的影响,将疲劳加载x方向下的归一化剩余刚度E′(n)与循环次数的关系式建立如下:
其中,A和L分别是平面外褶皱缺陷的高度和宽度;a和b是与材料内部褶皱缺陷有关的参数;p和q是独立于固有耗散能量和加载周期的两个材料参数;归一化刚度E’(n)是无量纲的,为了保持量纲的统一,p的单位是(J·m-3·s-1)-1,a,b和q被定义为无量纲参数;参数a,b的作用是反映褶皱几何形状的影响;参数q的作用是控制函数的形状;参数p是用来调节distab的影响,因为疲劳过程中的distab取决于试验材料。
进一步地,所述步骤S5中疲劳寿命预测模型的建立过程如下:
步骤S51:结合公式(2)和公式(3),疲劳损伤指数可以表示为:
步骤S52:达到最终失效循环数Nf时的损伤累积指数必须等于1,即D(Nf)=1,经过推导后的疲劳寿命预测模型可以表达如下:
进一步地,选用的待测风力机叶片主梁是由单向玻璃纤维增强环氧树脂GFRP材料制成的。
进一步地,疲劳寿命预测结果可以通过传统的疲劳试验得到,也可以通过提出的寿命预测模型进行估算。
与现有技术相比本发明有益效果。
本发明通过引入未知参数,将疲劳过程的损伤变量与宏观的能量耗散和剩余刚度建立了直观联系,提出了一种含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法。与现有预测方法相比,本发明具有测试时间短、节省测试成本、试验结果准确等优点,是一种可以预测含缺陷的主梁材料疲劳性能的新技术。同时,本发明的物理意义明确,且技术简单,方便工程技术人员的使用并为主梁材料的制作、检测和修复提供参考。
附图说明
下面结合附图和具体实施方式对本发明做进一步说明。本发明保护范围不仅局限于以下内容的表述。
图说明
图1为本发明的疲劳寿命预测方法的流程示意图。
图2为本发明实施例的表面温度变化曲线图。
图3为本发明实施例的最大加载应力σmax与固有耗散distab的关系图。
图4为本发明实施例的循环次数n、归一化剩余刚度E′(n)与固有耗散distab三维数据图。
图5为本发明实施例的疲劳试验预测S-N曲线结果与传统试验疲劳寿命结果图。
具体实施方式
为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及具体实施方式,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
如图1所示,提供了本发明的疲劳寿命预测方法的流程示意图,包括:
步骤一:对待测的含褶皱缺陷的主梁试件在万能疲劳试验机上首先进行拉伸性能测试,测出其最大载荷应力(UTS)。疲劳测试的最大加载应力在UTS的20%到90%之间以5%为步长分别对试件施加。每次仅在8000个加载周期内对试件进行测试。疲劳测试时使用红外热像仪同步观察和记录试件表面的温度变化,获取与疲劳寿命相关的试验数据。
步骤二:根据步骤S1获得的试件表面温度数据,计算主梁试件在不同加载应力下温度稳定阶段时的固有耗散能量。
步骤三:绘制试件的固有耗散dista与相应的UTS水平的关系图,利用双直线方法对图中具有两种不同斜率的数据趋势线进行插值,由两条直线交点的σmax坐标值确定出试件的疲劳极限。把疲劳极限以上载荷认定为过载区域,即试件发生失效的区域;而低于疲劳极限以下载荷的区域在试验过程中没有失效试件。
步骤四:基于疲劳过程中的固有耗散变化建立含待确定参数的归一化剩余刚度模型。将步骤S1、S2记录的疲劳试验数据代入模型,确定模型中的相关参数值。
步骤五:对于含褶皱缺陷的主梁试件,根据确定参数后的归一化剩余刚度模型推导出疲劳寿命预测模型。
步骤六:在本申请中,试件可以承受最大疲劳载荷的循环次数定义为106次,以此求出归一化失效刚度。再将确定归一化失效刚度后的疲劳寿命预测模型用于待测试件在过载区域内的疲劳寿命预测,以获得含褶皱缺陷的主梁试件的预测S-N曲线。
本实施例中,试件所采用的材料选取自实际大型叶片生产公司,包括PPG-Devold公司生产的1250gsm型号的玻璃纤维材料和Hexion RIM公司生产的135型号的环氧树脂材料。本实施例制作同实际主梁结构完全相同的单向铺层为10层的层合板试件,试件尺寸以255mm×30mm×10mm作为层合板的长度、宽度和厚度。层合板试件制作时,在抽真空之前将底层纤维布下面的中心处插入一根圆柱形塑料棒以引入褶皱缺陷。本实施例的层合板制作为中心位置具有高宽比为0.4的含褶皱缺陷的试件。玻璃纤维和环氧树脂的性能参数已从供应商数据表中获得,表1显示了这些材料的弹性参数。试件的几何特性如表2所示。
表1.材料的弹性力学性能
表2.褶皱缺陷试件的几何特性
疲劳试验设备为电动液压伺服万能疲劳试验机(MTS 810),首先对试件进行拉伸性能试验,其中,拉伸速率为2毫米/分钟。试验至试件最终断裂,确定了试件的极限拉伸应力(UTS)和极限拉伸应变分别为0.92GPa和0.1625mm。然后对试件进行疲劳试验,疲劳载荷为拉-拉载荷,应力比为0.1,振动频率为10Hz。本实施例中,施加的最大载荷应力水平以5%的间隔在UTS的20%到90%之间变化。对于每个最大加载应力,仅在8000个加载周期内对试样进行测试,不必做疲劳试验直到试件失效,足以获得最大载荷应力σmax与稳定的固有耗散distab之间的关系曲线。试验过程中用红外热像仪对试件表面的温度进行监测记录,红外热像仪与试件的表面相距400 mm。同时,将应变片粘贴在试件的两端,用应变仪来测量试件加载方向的应变以获得剩余刚度。
在不同的最大应力水平加载下,由红外热像仪记录的试件表面平均温度变化如图2所示。将试件在每个应力水平下第II阶段的表面稳定温度数据T代入如下的固有耗散计算表达式:
distab=-div(kgradT) (1)
其中,k是导热系数,k=1.5W/m·k。公式(1)计算后得到的最大加载应力σmax与固有耗散distab的关系如图3所示。
图3显示出数据具有两种不同斜率的线性趋势,通过最小二乘法拟合直线确定出低于和高于疲劳极限的两条直线。两条线性回归直线中数据交点的横坐标对应为试件的疲劳极限。从图3可以将试件的疲劳极限确定为80.8%UTS。
本发明将疲劳极限以上载荷认定为过载区域,即发生失效的区域,低于疲劳极限的载荷下没有失效试件,故只需对失效区域下的疲劳寿命进行预测。本发明的疲劳极限确定方法只需要采集试件在不同加载应力下8000次循环内的平均稳定温升,克服了传统疲劳方法需要大量试件和较长试验时间的局限性。
根据疲劳加载方向下的归一化剩余刚度与循环次数建立的关系式:
将试件的几何参数A/L=0.4以及图4中不同循环次数n下的distab和E′(n)数据代入公式(2),求出a,b,p和q的参数值。可以求出四个未知参数值:a=-5.31034,b=5.09028,p=0.00336,q=-2.18267E17。参数值确定后代入公式(2),由于对应于疲劳极限的失效循环次数取n=106,此时distab=41.3kJ·m-3,再次根据公式(2)求出试件失效时的归一化刚度
根据疲劳寿命预测模型表达式:
在已知a,b,p,q和后,代入不同应力水平下的distab,可以根据公式(3)预测出试件的整个S-N曲线。本发明实施例得到的疲劳试验预测S-N曲线结果如图5所示。
试件的S-N曲线既可以通过提出的寿命预测模型预测出,也可以通过传统的疲劳试验得到。本实施例还对试件进行了传统的疲劳试验,分别对试件施加五种最大应力水平:极限应力(UTS)的95%,90%,85%,80%和75%,在载荷控制下加载直至每个试件失效或循环至106次。图5显示了试件的传统疲劳试验结果。预测的S-N曲线与传统的试验寿命结果拟合度非常高,且预测的寿命结果都在传统结果95%的置信区间,符合风力发电机叶片的通用认证规则。对应于本发明失效循环次数取106下预测的S-N曲线没有被高估,预测结果基本都是小于传统试验结果的,是相对保守的,符合实际工程的安全考虑。
可以理解的是,以上关于本发明的具体描述,仅用于说明本发明而并非受限于本发明实施例所描述的技术方案,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换,以达到相同的技术效果;只要满足使用需要,都在本发明的保护范围之内。
Claims (2)
1.含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法,其特征在于,包括以下步骤:S1、对待测的含褶皱缺陷的主梁试件在万能疲劳试验机上进行拉伸性能测试,测出其最大载荷应力UTS;疲劳测试时使用红外热像仪同步观察和记录试件表面的温度变化,试验结束后统计并记录试件在不同循环次数下的表面温度、应力、应变和刚度数据;
S2、根据S1获得的试件表面温度数据,计算主梁试件在不同加载应力下温度稳定阶段时的固有耗散能量;
S3、绘制试件的固有耗散dista与相应的UTS水平的关系图,利用双直线方法对关系图中具有两种不同斜率的数据趋势线进行插值,由两条直线交点的σmax坐标值确定出试件的疲劳极限;把疲劳极限以上载荷认定为过载区域,即试件发生失效的区域;而低于疲劳极限以下载荷的区域在试验过程中没有失效试件;
步骤S4:基于疲劳过程中的固有耗散变化建立含待确定参数的归一化剩余刚度模型,将步骤S1、S2记录的疲劳试验数据代入模型,确定模型中未知参数a,b,p,q的值;a和b是与材料内部褶皱缺陷有关的参数;p和q是独立于固有耗散能量和加载周期的两个材料参数;
步骤S5:对于含褶皱缺陷的主梁试件,根据确定参数后的归一化剩余刚度模型推导出疲劳寿命预测模型;
步骤S6:定义试件可以承受最大疲劳载荷的循环次数,以此求出归一化失效刚度;再将确定归一化失效刚度后的疲劳寿命预测模型用于待测试件在过载区域内的疲劳寿命预测,以获得含褶皱缺陷的主梁试件的预测S-N曲线;
S1步骤中,疲劳测试的最大加载应力在UTS的20%到90%之间以5%为步长分别对试件施加;每次仅在8000个加载周期内对试件进行测试;
所述步骤S1中的疲劳试验机为电动液压伺服万能疲劳试验机;进行拉伸性能试验,拉伸速率为2毫米/分钟;进行疲劳试验,试验机的疲劳载荷为拉-拉载荷,应力比为0.1,振动频率为10Hz;
所述步骤S2中计算主梁试件在温度稳定阶段时的固有耗散能量的步骤包括:
主梁承受均匀的单轴循环载荷时,内变量和温度之间的内部耦合源可以忽略;在疲劳过程中温度变化的稳定阶段,当试件的每一周循环结束时,热损失与固有耗散能量之间存在平衡,该阶段稳定的固有耗散计算模型可表达为:
distab=-div(kgradTstab) (1)
其中,Tstab是试件在表面温度稳定阶段下的温度,其单位为℃,k是导热系数,其单位是W·m-1·K-1;div(kgradTstab)表示由传导引起的热损失率;
所述步骤S4中归一化剩余刚度模型的建立过程如下:
步骤S41:在疲劳试验过程中,疲劳损伤指数D(n)通过疲劳载荷下的剩余刚度表示如下:
其中,n代表循环次数;E(n)代表第n个循环次数相对应的刚度;E0是沿x轴方向的初始有效刚度;E’(n)是第n个循环次数相应的归一化刚度,定义为E(n)/E0;Ef代表失效刚度;E′f是最终循环Nf时的归一化失效刚度,定义为Ef/E0;
步骤S42:由于稳定的固有耗散能量值distab与损伤累积之间存在一定的关系,通过引入褶皱缺陷中高宽比的影响,将疲劳加载x方向下的归一化剩余刚度E′(n)与循环次数的关系式建立如下:
其中,A和L分别是平面外褶皱缺陷的高度和宽度;a和b是与材料内部褶皱缺陷有关的参数;p和q是独立于固有耗散能量和加载周期的两个材料参数;归一化刚度E’(n)是无量纲的,为了保持量纲的统一,p的单位是(J·m-3·s-1)-1,a,b和q被定义为无量纲参数;参数a,b的作用是反映褶皱几何形状的影响;参数q的作用是控制函数的形状;参数p是用来调节distab的影响,因为疲劳过程中的distab取决于试验材料;
所述步骤S5中疲劳寿命预测模型的建立过程如下:
步骤S51:结合公式(2)和公式(3),疲劳损伤指数可以表示为:
步骤S52:达到最终失效循环数Nf时的损伤累积指数必须等于1,即D(Nf)=1,经过推导后的疲劳寿命预测模型可以表达如下:
2.根据权利要求1所述的含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法,其特征在于:选用的待测风力机叶片主梁是由单向玻璃纤维增强环氧树脂GFRP材料制成的。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011527291.1A CN112710539B (zh) | 2020-12-22 | 2020-12-22 | 含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法 |
US17/544,607 US11566606B2 (en) | 2020-12-22 | 2021-12-07 | Method for quickly predicting fatigue life of wrinkle defect-containing main spar in wind turbine blade |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011527291.1A CN112710539B (zh) | 2020-12-22 | 2020-12-22 | 含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112710539A CN112710539A (zh) | 2021-04-27 |
CN112710539B true CN112710539B (zh) | 2024-05-17 |
Family
ID=75545109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011527291.1A Active CN112710539B (zh) | 2020-12-22 | 2020-12-22 | 含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11566606B2 (zh) |
CN (1) | CN112710539B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113378327B (zh) * | 2021-07-02 | 2022-09-02 | 吉林重通成飞新材料股份公司 | 风电叶片后缘抗开裂设计方法、装置、设备及存储介质 |
CN114065556B (zh) * | 2022-01-14 | 2022-05-03 | 湖南大学 | 一种基于耗散能的变幅疲劳寿命预测方法 |
JP7487248B2 (ja) | 2022-03-28 | 2024-05-20 | 株式会社豊田中央研究所 | 疲労損傷度特定装置及び疲労損傷度の特定方法 |
CN115408893B (zh) * | 2022-07-01 | 2023-10-13 | 重庆大学 | 一种基于疲劳寿命预测的电池包设计方法 |
CN115438556B (zh) * | 2022-11-10 | 2022-12-30 | 北京理工大学 | 柔性充气飞行器结构刚度退化率预测方法、装置及设备 |
CN116579193B (zh) * | 2022-12-20 | 2024-03-19 | 湘潭大学 | 挠度应变能等效的梁结构局部等效刚度计算方法 |
CN117150685B (zh) * | 2023-10-26 | 2024-02-06 | 国合通用(青岛)测试评价有限公司 | 一种紧固螺栓的疲劳寿命评估系统及方法 |
CN117436344A (zh) * | 2023-11-10 | 2024-01-23 | 沈阳工业大学 | 一种基于参数化描述的风力机叶片结构优化设计方法 |
CN117454719B (zh) * | 2023-12-21 | 2024-02-23 | 无锡阳光精机股份有限公司 | 一种机床主轴疲劳强度分析方法、装置、设备及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007107955A (ja) * | 2005-10-12 | 2007-04-26 | Ntn Corp | 捩り疲労試験機 |
CN103344514A (zh) * | 2013-07-05 | 2013-10-09 | 北京航空航天大学 | 一种基于名义应力法的高周疲劳和低强度冲击耦合的损伤计算方法 |
CN105447316A (zh) * | 2015-12-01 | 2016-03-30 | 北京航空航天大学 | 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法 |
WO2018107730A1 (zh) * | 2016-12-16 | 2018-06-21 | 海口未来技术研究院 | 复合材料的疲劳寿命预测方法及预测系统 |
CN111323316A (zh) * | 2020-01-06 | 2020-06-23 | 湖南大学 | 多轴疲劳寿命预测方法以及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9476815B2 (en) * | 2011-02-25 | 2016-10-25 | Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College | Fatigue monitoring for composite materials |
US10744727B2 (en) * | 2017-03-21 | 2020-08-18 | Textron Innovations Inc. | Methods of making a specimen with a predetermined wrinkle defect |
-
2020
- 2020-12-22 CN CN202011527291.1A patent/CN112710539B/zh active Active
-
2021
- 2021-12-07 US US17/544,607 patent/US11566606B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007107955A (ja) * | 2005-10-12 | 2007-04-26 | Ntn Corp | 捩り疲労試験機 |
CN103344514A (zh) * | 2013-07-05 | 2013-10-09 | 北京航空航天大学 | 一种基于名义应力法的高周疲劳和低强度冲击耦合的损伤计算方法 |
CN105447316A (zh) * | 2015-12-01 | 2016-03-30 | 北京航空航天大学 | 一种预测含初始缺陷复合材料桨叶剩余寿命的渐进损伤方法 |
WO2018107730A1 (zh) * | 2016-12-16 | 2018-06-21 | 海口未来技术研究院 | 复合材料的疲劳寿命预测方法及预测系统 |
CN111323316A (zh) * | 2020-01-06 | 2020-06-23 | 湖南大学 | 多轴疲劳寿命预测方法以及装置 |
Non-Patent Citations (2)
Title |
---|
基于能量耗散的不锈钢Fv520B疲劳性能分析;蔺宏岩;王振芬;;钢结构;20180822(第08期);49-52 * |
微缺陷风力机叶片复合材料的疲劳损伤及疲劳寿命研究;王琳琳;陈长征;周勃;康爽;;玻璃钢/复合材料;20191228(第12期);20-25 * |
Also Published As
Publication number | Publication date |
---|---|
US20220195991A1 (en) | 2022-06-23 |
CN112710539A (zh) | 2021-04-27 |
US11566606B2 (en) | 2023-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112710539B (zh) | 含褶皱缺陷的风力机叶片主梁疲劳寿命的快速预测方法 | |
Vassilopoulos | The history of fiber-reinforced polymer composite laminate fatigue | |
Nijssen | Fatigue life prediction and strength degradation of wind turbine rotor blade composites | |
McGugan et al. | Damage tolerance and structural monitoring for wind turbine blades | |
EP2390074B1 (en) | Apparatus and closed-loop method for the manufacture of prepregs and/or laminates comprising the same | |
Nachtane et al. | Damage prediction of horizontal axis marine current turbines under hydrodynamic, hydrostatic and impacts loads | |
EP2795309B1 (en) | A method for calibrating an on-line control of a manufacturing process for a multicomponent sheet material | |
Huang et al. | Statistical methodology for assessing manufacturing quality related to transverse cracking in cross ply laminates | |
Krueger et al. | Fatigue life methodology for bonded composite skin/stringer configurations | |
CN109855991B (zh) | 一种预测复合材料疲劳寿命的方法 | |
CN108228979B (zh) | 一种用于发电机结构的低周疲劳寿命分析方法 | |
Sala et al. | Fibre optics health monitoring for aeronautical applications | |
Khan et al. | A review of physics-based models in prognostics and health management of laminated composite structures | |
Cairns et al. | Wind turbine composite blade manufacturing: the need for understanding defect origins, prevalence, implications and reliability. | |
Sayer et al. | ReliaBlade project: A material’s perspective towards the digitalization of wind turbine rotor blades | |
Vassilopoulos | Fatigue life prediction of wind turbine blade composite materials | |
ITTO20080908A1 (it) | Procedimento per la prognostica di una struttura sottoposta a carichi. | |
Van Paepegem | Fatigue testing methods for polymer matrix composites | |
Pawar et al. | Fuzzy-logic-based health monitoring and residual-life prediction for composite helicopter rotor | |
Philippidis et al. | THIN/PROBEAM: a software for fatigue design and analysis of composite rotor blades | |
Yu et al. | Fatigue life prediction for the main spar with wrinkle defects of a wind turbine blade | |
Hugaas et al. | Estimating SN curves for local fiber dominated fatigue failure in ring specimens representing filament wound pressure vessels with damage | |
Bacharoudis et al. | Probabilistic analysis of wind turbine blades considering stiffness, strength and stability under extreme and fatigue loading | |
Voronkov et al. | Experimental study of the stress-strain state features of outlet guide vane made from polymer composite material using fiber optic sensors | |
Koch et al. | Analysis of the micro-cracking behaviour of carbon fibre reinforced flywheel rotors considering residual stresses |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |