CN111122346B - Test load processing method for main structure of airfoil - Google Patents
Test load processing method for main structure of airfoil Download PDFInfo
- Publication number
- CN111122346B CN111122346B CN201911346781.9A CN201911346781A CN111122346B CN 111122346 B CN111122346 B CN 111122346B CN 201911346781 A CN201911346781 A CN 201911346781A CN 111122346 B CN111122346 B CN 111122346B
- Authority
- CN
- China
- Prior art keywords
- load
- test
- theoretical
- test load
- state
- 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
Images
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/20—Investigating strength properties of solid materials by application of mechanical stress by applying steady bending forces
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
-
- 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/24—Investigating strength properties of solid materials by application of mechanical stress by applying steady shearing forces
-
- 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/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0216—Finite elements
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a test load processing method for a main structure of an airfoil, provides a set of test load processing method and flow, provides an evaluation method, can objectively evaluate the feasibility of test loading load, ensures the loading accuracy and improves the load processing efficiency.
Description
Technical Field
The invention belongs to the technical field of aviation strength, and particularly relates to a load processing method for an airfoil main structure test.
Background
When the ground static test is carried out on the airplane structure, the theoretical load is often processed to obtain the test load convenient to implement in consideration of the enforceability of the test loading scheme, but the difference between the test load and the theoretical load brings different stress of the airplane body structure, and in order to guarantee the accuracy of force transmission and the safety of the structure, the structural stress distribution under the test load state needs to be checked, and the safety of the local structure of a loading point needs to be checked.
The above process is often a repeated iteration process, and a reasonable scheme can be found between the theoretical precision and the test implementability only by adjusting the test load for multiple times, which is time-consuming and labor-consuming, and the determination of the reasonable scheme has no uniform standard and lacks objectivity.
Disclosure of Invention
The purpose of the invention is as follows: the method comprises the steps of processing theoretical load into test load through an iteration method under the condition of limiting single-point test load for the main structure of the airfoil, checking a load bending-shearing torsion curve, displacement and stress strain to determine whether the test load is feasible, and adjusting the test load if the test load is not feasible until all requirements are met.
The technical scheme of the invention is as follows:
a test load processing method for an airfoil main structure comprises the following steps:
step 1: calculating the test load P of each test loading pointj;
Step 2: checking the test load P of each test load pointjWhether or not the allowable load P is exceededj]If the test load P of the test load point isjExceeding allowable load [ P ]j]I.e. Pj≥[Pj]The excess load amount Δ PjProcessing until the allowable loading load is not exceeded Pj]At the test load point of (a);
and step 3: repeating the step 2 until the test load P on each test loading pointjAll do not exceed allowable load [ Pj]I.e. Pj≤[Pj];
And 4, step 4: calculating the difference delta F on each rib station profile of the bending shear in the theoretical state and the experimental statek、ΔMk、ΔJk;
And 5: calculating the difference delta u of the deflection of the front beam and the rear beam of the wing box on each rib station sectionf,kAnd Δ ur,kAnd the difference in stress Δ σ at the primary force-transmitting locationp;
Step 6: simultaneous determination of Δ Fk、ΔMkAnd Δ JkWhether or not smaller than the error allowable value epsilon respectivelyF,k、εM,k、εJ,k;Δuf,kAnd Δ ur,kWhether or not less than the error allowable value epsilonu,f,kAnd εu,r,k;ΔσpWhether or not less than the error allowable value epsilonσ,pIf not, locally adjusting the load by observing the distribution characteristics of the difference;
and 7: and (5) repeating the steps 4 to 7 until all error requirements are met, and obtaining the final test loading load.
Step 1, calculating the test load P of each test loading pointjThe method specifically comprises the following steps: for a givenLoad case, loading nodal point on airfoil FiProcessing the load to a test loading point according to an energy method, and calculating to obtain a test load Pj。
The method is characterized in that: test load PjThe calculation formula is as follows:
wherein: m is the number of theoretical nodes, Pi,jThe load distributed to the jth test loading point for the ith theoretical node load is calculated according to the following formula:
wherein: j is the number of test load points, coefficient lambdai、λi,xAnd λi,zObtained by solving the following system of equations:
wherein:xjand yjIs the coordinate value of the jth test load point, xiAnd yiAnd the coordinate value of the ith theoretical load point.
Excess load Δ P in step 2j=Pj-[Pj]。
The excess load amount delta P in step 2jProcessing until the allowable loading load is not exceeded Pj]The treatment method used at the test load point of (1) is an energy method.
Step 4, difference delta F on each rib station section of the bending shear torsion curves in the theoretical state and the test statek、ΔMk、ΔJkThe calculation formula is as follows:
wherein m iskThe number of theoretical nodes outside the k section, nkNumber of test load points outside the k section, dk,iIs the distance from the ith theoretical node to the profile k, dk,jDistance from jth test load point to section k, diDistance from the ith theoretical node to the rigid axis of the airfoil, djThe distance from the jth test loading point to the airfoil rigid axis.
Step 5, the difference delta u of the deflection of the front beam and the deflection of the rear beam of the wing box on the section of each rib stationf,kAnd Δ ur,kAnd the difference in stress Δ σ at the primary force-transmitting locationp,Δuf,kAnd Δ ur,kAnd Δ σpCalculated as follows: respectively applying theoretical node load and test loading load by establishing a finite element model of the airfoil structure to obtain the displacement and stress analysis results of two states, whereinAndis the deflection of the front and rear beams at the section k of the theoretical statef,kAnd ur,kThe deflection of the front beam and the back beam at the section k under the test state,and σpStress at the theoretical state and the test state, respectively, at the site p
Test condition of wing box front beam at section kThe difference in deflection from the theoretical state is
The difference of the deflection of the wing box back beam at the section k between the test state and the theoretical state is
The stress difference between the test state and the theoretical state at any main force transmission part p is
k is 1, …, s, s is the number of sections, p is 1, …, t, t is the number of main force transfer parts.
And 5, the main force transmission part is a wall plate.
And 5, the main force transmission part is a wing beam.
The invention has the beneficial effects that: compared with the prior art: the invention provides a set of test load processing method and flow and an evaluation method, which can objectively evaluate the feasibility of test loading load, ensure the loading accuracy and improve the load processing efficiency.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The invention is further described with reference to the accompanying drawings, and the test load processing method of the airfoil main structure comprises the following steps:
step 1: calculating the test load P of each test loading pointjThe method specifically comprises the following steps: for a given load case, loading the nodal point on the airfoil FiProcessing the load to a test loading point according to an energy method, and calculating to obtain a test load PjTest load PjThe calculation formula is as follows:
wherein: m is the number of theoretical nodes, Pi,jThe load distributed to the jth test loading point for the ith theoretical node load is calculated according to the following formula:
wherein: j is the number of test load points, coefficient lambdai、λi,xAnd λi,zObtained by solving the following system of equations:
wherein:xjand yjIs the coordinate value of the jth test load point, xiAnd yiAnd the coordinate value of the ith theoretical load point.
Step 2: checking the test load P of each test load pointjWhether or not the allowable load P is exceededj]If the test load P of the test load point isjExceeding allowable load [ P ]j]I.e. Pj≥[Pj]Then the excess load Δ P is measured by an energy methodjProcessing until the allowable loading load is not exceeded Pj]At the test load point of (1), the excess load amount Δ Pj=Pj-[Pj]。
And step 3: repeating the step 2 until the test load P on each test loading pointjAll do not exceed allowable load [ Pj]I.e. Pj≤[Pj];
And 4, step 4: calculating the difference delta F on each rib station profile of the bending shear in the theoretical state and the experimental statek、ΔMk、ΔJk;
The calculation formula is as follows:
wherein m iskThe number of theoretical nodes outside the k section, nkNumber of test load points outside the k section, dk,iIs the distance from the ith theoretical node to the profile k, dk,jDistance from jth test load point to section k, diDistance from the ith theoretical node to the rigid axis of the airfoil, djThe distance from the jth test loading point to the airfoil rigid axis.
And 5: calculating the difference delta u of the deflection of the front beam and the rear beam of the wing box on each rib station sectionf,kAnd Δ ur,kAnd the difference in stress Δ σ at the primary force-transmitting locationp;
Δσp,Δuf,kAnd Δ ur,kAnd Δ σpCalculated as follows:
respectively applying theoretical node load and test loading load by establishing a finite element model of the airfoil structure to obtain the displacement and stress analysis results of two states, whereinAndis the deflection of the front and rear beams at the section k of the theoretical statef,kAnd ur,kThe deflection of the front beam and the back beam at the section k under the test state,and σpAt the theoretical and experimental states, respectively, at the site pStress of
The difference of the deflection of the wing box front beam at the section k between the test state and the theoretical state is
The difference of the deflection of the wing box back beam at the section k between the test state and the theoretical state is
The stress difference between the test state and the theoretical state at any main force transmission part p is
k is 1, …, s, s is the number of sections, p is 1, …, t, t is the number of main force transfer parts.
The main force transmission part is a wall plate and a wing beam.
Step 6: simultaneous determination of Δ Fk、ΔMkAnd Δ JkWhether or not smaller than the error allowable value epsilon respectivelyF,k、εM,k、εJ,k;Δuf,kAnd Δ ur,kWhether or not less than the error allowable value epsilonu,f,kAnd εu,r,k;ΔσpWhether or not less than the error allowable value epsilonσ,pIf not, locally adjusting the load by observing the distribution characteristics of the difference;
and 7: and (5) repeating the steps 4 to 7 until all error requirements are met, and obtaining the final test loading load. The invention provides a set of test load processing method and flow and an evaluation method, which can objectively evaluate the feasibility of test loading load, ensure the loading accuracy and improve the load processing efficiency.
Claims (8)
1. A test load processing method for an airfoil main structure is characterized by comprising the following steps: the method comprises the following steps:
step 1: calculating the test load P of each test loading pointj;
Step 2: checking the test load P of each test load pointjWhether or not the allowable load P is exceededj]If the test load P of the test load point isjExceeding allowable load [ P ]j]I.e. Pj≥[Pj]The excess load amount Δ PjProcessing until the allowable loading load is not exceeded Pj]At the test load point of (a);
and step 3: repeating the step 2 until the test load P on each test loading pointjAll do not exceed allowable load [ Pj]I.e. Pj≤[Pj];
And 4, step 4: calculating the difference delta F on each rib station profile of the bending shear in the theoretical state and the experimental statek、ΔMk、ΔJk;
And 5: calculating the difference delta u of the deflection of the front beam and the rear beam of the wing box on each rib station sectionf,kAnd Δ ur,kAnd the difference in stress Δ σ at the primary force-transmitting locationp(ii) a Calculated as follows:
respectively applying theoretical node load and test loading load by establishing a finite element model of the airfoil structure to obtain the displacement and stress analysis results of two states, whereinAndis the deflection of the front and rear beams at the section k of the theoretical statef,kAnd ur,kThe deflection of the front beam and the back beam at the section k under the test state,and σpStress at the theoretical state and the test state, respectively, at the site p
The difference of the deflection of the wing box front beam at the section k between the test state and the theoretical state is
The difference of the deflection of the wing box back beam at the section k between the test state and the theoretical state is
The stress difference between the test state and the theoretical state at any main force transmission part p is
k is 1, …, s, s is the number of sections, p is 1, …, t, t is the number of main force transmission parts;
step 6: simultaneous determination of Δ Fk、ΔMkAnd Δ JkWhether or not smaller than the error allowable value epsilon respectivelyF,k、εM,k、εJ,k;Δuf,kAnd Δ ur,kWhether or not less than the error allowable value epsilonu,f,kAnd εu,r,k;ΔσpWhether or not less than the error allowable value epsilonσ,pIf not, locally adjusting the load by observing the distribution characteristics of the difference;
and 7: and (5) repeating the steps 4 to 6 until all error requirements are met, and obtaining the final test loading load.
2. The test load handling method for the main structure of the airfoil according to claim 1, wherein: step 1, calculating the test load P of each test loading pointjThe method specifically comprises the following steps: for a given load case, loading the nodal point on the airfoil FiProcessing the load to a test loading point according to an energy method, and calculating to obtain a test load Pj。
3. The test load handling method for the main structure of the airfoil according to claim 2, wherein: test load PjThe calculation formula is as follows:
wherein: m is the number of theoretical nodes, Pi,jThe load distributed to the jth test loading point for the ith theoretical node load is calculated according to the following formula:
wherein: j is the number of test load points, coefficient lambdai、λi,xAnd λi,zObtained by solving the following system of equations:
xjand yjIs the coordinate value of the jth test load point, xiAnd yiAnd the coordinate value of the ith theoretical load point.
4. The test load handling method for the main structure of the airfoil according to claim 1, wherein: excess load Δ P in step 2j=Pj-[Pj]。
5. The test load handling method for the main structure of the airfoil according to claim 1, wherein: the excess load amount delta P in step 2jProcessing until the allowable loading load is not exceeded Pj]The treatment method used at the test load point of (1) is an energy method.
6. The test load handling method for the main structure of the airfoil according to claim 1, wherein: step 4, the difference delta F of the sections of the theoretical state and the experimental state of the bending shear on the station positions of each ribk、ΔMk、ΔJkThe calculation formula is as follows:
wherein m iskThe number of theoretical nodes outside the k section, nkNumber of test load points outside the k section, dk,iIs the distance from the ith theoretical node to the profile k, dk,jDistance from jth test load point to section k, diDistance from the ith theoretical node to the rigid axis of the airfoil, djDistance from jth test load point to airfoil rigid axis, FiIs the nodal load on the airfoil.
7. The test load handling method for the main structure of the airfoil according to claim 1, wherein: and 5, the main force transmission part is a wall plate.
8. The test load handling method for the main structure of the airfoil according to claim 1, wherein: and 5, the main force transmission part is a wing beam.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911346781.9A CN111122346B (en) | 2019-12-24 | 2019-12-24 | Test load processing method for main structure of airfoil |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911346781.9A CN111122346B (en) | 2019-12-24 | 2019-12-24 | Test load processing method for main structure of airfoil |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111122346A CN111122346A (en) | 2020-05-08 |
CN111122346B true CN111122346B (en) | 2022-04-19 |
Family
ID=70501908
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911346781.9A Active CN111122346B (en) | 2019-12-24 | 2019-12-24 | Test load processing method for main structure of airfoil |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111122346B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5511430A (en) * | 1993-12-06 | 1996-04-30 | Aerospatiale Societe Nationale Industrielle | Method and device for detecting that the design loads of an aircraft have been exceeded |
CN1514213A (en) * | 2002-12-31 | 2004-07-21 | 中国农业机械化科学研究院 | Full machine ground load on site calibration test method and its device |
CN104058101A (en) * | 2014-06-24 | 2014-09-24 | 中国飞机强度研究所 | Method for applying normal load under large deformation condition of wings |
CN108528758A (en) * | 2018-05-04 | 2018-09-14 | 中电科芜湖钻石飞机制造有限公司 | General-purpose aircraft mechanical property testing system |
CN109490116A (en) * | 2018-12-12 | 2019-03-19 | 中国航空工业集团公司西安飞机设计研究所 | A kind of wing Plumb load method of full scale fatigue test |
CN109710982A (en) * | 2018-12-04 | 2019-05-03 | 中国特种飞行器研究所 | Component-level simulation test piece design and verification method based on FEM Simulation |
CN109733638A (en) * | 2018-09-07 | 2019-05-10 | 中国飞机强度研究所 | Load applying method in the case of a kind of long struts undercarriage large deformation |
CN109977556A (en) * | 2019-03-28 | 2019-07-05 | 中国飞机强度研究所 | Load-optimised method based on least square method |
-
2019
- 2019-12-24 CN CN201911346781.9A patent/CN111122346B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5511430A (en) * | 1993-12-06 | 1996-04-30 | Aerospatiale Societe Nationale Industrielle | Method and device for detecting that the design loads of an aircraft have been exceeded |
CN1514213A (en) * | 2002-12-31 | 2004-07-21 | 中国农业机械化科学研究院 | Full machine ground load on site calibration test method and its device |
CN104058101A (en) * | 2014-06-24 | 2014-09-24 | 中国飞机强度研究所 | Method for applying normal load under large deformation condition of wings |
CN108528758A (en) * | 2018-05-04 | 2018-09-14 | 中电科芜湖钻石飞机制造有限公司 | General-purpose aircraft mechanical property testing system |
CN109733638A (en) * | 2018-09-07 | 2019-05-10 | 中国飞机强度研究所 | Load applying method in the case of a kind of long struts undercarriage large deformation |
CN109710982A (en) * | 2018-12-04 | 2019-05-03 | 中国特种飞行器研究所 | Component-level simulation test piece design and verification method based on FEM Simulation |
CN109490116A (en) * | 2018-12-12 | 2019-03-19 | 中国航空工业集团公司西安飞机设计研究所 | A kind of wing Plumb load method of full scale fatigue test |
CN109977556A (en) * | 2019-03-28 | 2019-07-05 | 中国飞机强度研究所 | Load-optimised method based on least square method |
Non-Patent Citations (2)
Title |
---|
基于全机平衡的试验载荷谱优化方法;宋海平 等;《航空安全与装备维修技术——航空安全与装备维修技术学术研讨会论文集》;20140630;第275-279页 * |
复合材料风力机叶片尾缘结构优化设计;郑玉巧 等;《现代制造工程》;20190831(第8期);第56-61页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111122346A (en) | 2020-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107273590B (en) | A kind of Multipurpose Optimal Method determined for composite material constitutive equation parameter | |
WO2019047529A1 (en) | Construction method for dynamic shearing constitutive model of fiber-reinforced composite material | |
Barry et al. | A three‐dimensional element‐free Galerkin elastic and elastoplastic formulation | |
CN105930571B (en) | Big across steel bridge correction method for finite element model based on unit temperature response monitor value | |
Krishnamurthy | Modelling and prediction of steel bolted connection behavior | |
CN108984834B (en) | Wing reliability evaluation system and method based on response surface method | |
CN107766612B (en) | Method for measuring wing load in connecting wing structure form | |
CN111122346B (en) | Test load processing method for main structure of airfoil | |
CN113343559B (en) | Reliability analysis method for response surface of iterative reweighted least square method extreme learning machine | |
CN111159934A (en) | Aluminum matrix composite material dynamics simulation method | |
CN103272982B (en) | Method for determining upsetting direction of rivet for riveting assembly of metal thin-wall part | |
Huang et al. | A Kriging-based non-probability interval optimization of loading path in T-shape tube hydroforming | |
CN109214131B (en) | Error-optimized static test load design method and system | |
Davids et al. | Nonlinear elastic response of locally buckled thin-walled beam-columns | |
Chepurnenko et al. | Features of compressed rods calculations with account of initial imperfections and creep effects | |
CN114993604B (en) | Wind tunnel balance static calibration and measurement method based on deep learning | |
CN111310251B (en) | High-precision structure reliability analysis method based on three-weighted response surface | |
CN112763345B (en) | Airplane airfoil load calibration area loading test method | |
Gonzales et al. | Uncertainty quantification in aircraft load calibration | |
CN112651153B (en) | Method for determining material parameters of crystal plasticity finite element model | |
Su et al. | Mixed mode cracks in Reissner plates | |
CN108169013A (en) | A kind of fatigue life calibration method suitable for Multi-fasteners joint test specimen | |
Sadovský et al. | Imperfection measures of eigen-and periodic modes of axially loaded stringer-stiffened cylindrical shell | |
Franz et al. | Estimation of composite laminate ply angles using an inverse Bayesian approach based on surrogate models | |
CN117235883B (en) | Combined prediction method for aircraft panel assembly deviation |
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 |