CN113312815B - Method for calculating overall stability of stringers in fuselage structure - Google Patents

Method for calculating overall stability of stringers in fuselage structure Download PDF

Info

Publication number
CN113312815B
CN113312815B CN202110547079.XA CN202110547079A CN113312815B CN 113312815 B CN113312815 B CN 113312815B CN 202110547079 A CN202110547079 A CN 202110547079A CN 113312815 B CN113312815 B CN 113312815B
Authority
CN
China
Prior art keywords
stringer
compressive stress
stringers
stress sigma
fuselage
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
Application number
CN202110547079.XA
Other languages
Chinese (zh)
Other versions
CN113312815A (en
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.)
AVIC Xian Aircraft Industry Group Co Ltd
Original Assignee
AVIC Xian Aircraft Industry Group 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 AVIC Xian Aircraft Industry Group Co Ltd filed Critical AVIC Xian Aircraft Industry Group Co Ltd
Priority to CN202110547079.XA priority Critical patent/CN113312815B/en
Publication of CN113312815A publication Critical patent/CN113312815A/en
Application granted granted Critical
Publication of CN113312815B publication Critical patent/CN113312815B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

According to the invention, the minimum compressive stress of the fuselage stringer is calculated according to the stringer compressive stress calculation formula, and finally the total critical unstability stress of the stringer corresponding to the minimum compressive stress and the ratio of the calculated minimum compressive stress are used to obtain the total stability residual strength of the fuselage stringer.

Description

Method for calculating overall stability of stringers in fuselage structure
Technical Field
The invention belongs to the field of aircraft structural strength design, and particularly relates to a method for calculating overall stability of stringers in a fuselage structure.
Background
The invention mainly researches a method for calculating the overall stability of a stringer when a fuselage skin is mainly subjected to pure shearing.
In the current strength check, a checking method is generally adopted for the fuselage stringer, a stringer unit and left and right skin units under severe working conditions are taken as research objects, axial compressive stress of the stringer rod unit is directly extracted from an overall stress analysis model, maximum shearing stress of the skin unit is calculated, the compressive stress and bending stress of the stringer caused by a tension field after shearing instability of the skin unit are summed with the compressive stress of a stringer body, the overall stability of the fuselage stringer is calculated, and the physical geometric plastic coefficient of the stringer during bending and the influence of the compressive stress and the bending stress of the stringer on the physical geometric plastic coefficient are not considered during calculation.
For the fuselage stringers, the overall stability of the fuselage stringers is checked directly by adopting the stress calculation result of the overall stress analysis model, and the accuracy of the obtained result is lower. Therefore, a new calculation method is needed to be explored, and the overall stability of the stringer of the fuselage structure is calculated conveniently, rapidly and accurately.
Disclosure of Invention
The invention aims to provide a method for calculating the overall stability of a stringer in a fuselage structure, which is used for calculating static strength based on a finite element analysis method and aims to improve the calculation accuracy of the overall stability of the stringer when a skin is mainly sheared purely.
A method for calculating the overall stability of a stringer in a fuselage structure, the fuselage structure comprising a stringer, a skin and a fuselage frame, wherein the design digital-analog and load conditions of the fuselage structure are known, is characterized by comprising the following steps: 1) According to the design counting model of the fuselage structure, taking intersection points of the stringers and fuselage frames as references, establishing a finite element model of the fuselage structure, wherein the skin is simplified into plate units, the stringers are simplified into rod units, the fuselage frames are simplified into beam units, a calculation unit is formed between adjacent beam units by the rod units and the plate units on two sides of the rod units, and each calculation unit comprises the stringers in the middle position and the skin on two sides of the stringers; 2) In the finite element model, the frame beam unit of the machine body close to the machine head is simply supported and restrained to be used as a boundary condition of the overall stress analysis model of the machine body structure; 3) According to the boundary conditions and the load working conditions, stress calculation results of all calculation units of the fuselage structure are obtained; 4) According to the stress calculation results of the calculation units, extracting the calculation units with the skins subjected to pure shearing as analysis units; 5) Calculating each analysis unitThe compressive stress of the stringer is compared and analyzed, and the minimum compressive stress sigma of the stringer is extracted c The method comprises the steps of carrying out a first treatment on the surface of the 6) Calculating the critical instability stress sigma of the stringer by utilizing the stringer geometric parameter corresponding to the minimum compressive stress cf The method comprises the steps of carrying out a first treatment on the surface of the 7) Critical buckling stress sigma of stringers corresponding to minimum compressive stress cf And ratio sigma of minimum compressive stress of stringers c As the overall stability residual strength of the fuselage stringers.
In step 5) above, the stringer compressive stress σ of each analysis unit c Comprises three parts: the first part is the stringer compressive stress sigma calculated in the fuselage structure overall finite element model ch The method comprises the steps of carrying out a first treatment on the surface of the The second part is stringer average compressive stress sigma caused by diagonal stretching after skin shear instability st The method comprises the steps of carrying out a first treatment on the surface of the The third part is the bending stress sigma of the stringer at two fulcrums under inward bending load M
According to the method for calculating the overall stability of the stringer in the fuselage structure, a finite element model of the fuselage structure is established according to a design digital model of the fuselage structure, calculation units are formed by a rod unit and plate units on two sides of the rod unit, stress calculation results of all calculation units of the fuselage structure are obtained according to boundary conditions and load working conditions, geometric parameters of the calculation units corresponding to the maximum shear stress plate units and the maximum compression stress rod units are utilized, the minimum compression stress of the stringer of the fuselage is calculated according to the stringer compression stress calculation formula, and finally the overall stability residual strength of the stringer of the fuselage is obtained according to the ratio of the overall critical unstability stress of the stringer corresponding to the minimum compression stress and the calculated minimum compression stress.
The beneficial effects of this application lie in: 1) According to the invention, by utilizing a finite element analysis method, through accurate simplification and constraint of the fuselage skin, the fuselage frame and the stringers, the overall stability of the stringers in the fuselage structure is calculated, and a calculation unit of the skin subjected to pure shearing is extracted as an analysis unit, so that the method has the advantages of simple algorithm and accurate result; 2) In the calculation process, the plasticity of the material of the stringers in the analysis unit during bending and the stress redistribution condition after the stringers enter the plasticity are fully considered; 3) According to the actual loading condition of the structure, the influence of the compressive stress and the bending stress of the stringer caused by a tension field on the physical geometric plastic coefficient is fully considered, and a correction method of the physical geometric plastic coefficient of the stringer during bending is provided; 4) The overall stability calculation precision of the stringers when the fuselage structure skin is mainly subjected to pure shearing is improved, and the structural strength can be calculated more accurately and the weight reduction design of the structure can be carried out.
The present application is described in further detail below with reference to the examples and the accompanying drawings.
Drawings
FIG. 1 is a schematic illustration of a typical airframe structure;
FIG. 2 is a schematic diagram of a finite element analysis model of a typical fuselage structure;
fig. 3 is a schematic diagram of a computing unit.
The numbering in the figures illustrates: 1 skin, 2 stringers, 3 fuselage frames, 4 calculation units.
Detailed Description
Referring to the drawings, a typical fuselage structure contains a monolithic skin 1, a plurality of longitudinal stringers 2 and a plurality of transverse fuselage frames 3. The method for calculating the overall stability of the stringers in the fuselage structure comprises the following specific steps:
the static strength calculation is carried out based on a finite element analysis method, and the purpose is to improve the calculation accuracy of the overall stability of the fuselage stringer when the skin is mainly subjected to pure shearing.
According to the design digital model of the fuselage structure, taking the intersection point of the stringer 2 and the fuselage frame 3 as a reference, a finite element model of the fuselage structure is established, wherein the skin 1 is simplified into a plate unit, the stringer 2 is simplified into a rod unit, the fuselage frame 3 is simplified into a beam unit, the rod unit and the adjacent plate units form calculation units 4, and each calculation unit comprises the stringer 2 at the middle position and the skins 1 at two sides of the stringer 2, as shown in figure 3.
During analysis, the frame 3 at the end of the fuselage is simply supported and constrained to serve as a boundary condition of the overall stress analysis model of the fuselage structure.
And obtaining the stress calculation result of each calculation unit of the fuselage structure according to the boundary condition and the load working condition.
Extracting different load conditionsThe lower skin 1 is mainly subjected to pure shearing by a calculating unit 4, the shearing stress of the skin 1 and the calculating stress of the stringers 2 are extracted, and the compressive stress of the stringers calculated by the overall finite element model is recorded as sigma ch The method comprises the steps of carrying out a first treatment on the surface of the The average compressive stress of the stringers caused by diagonal tension after shear instability of the skin is noted as sigma st The method comprises the steps of carrying out a first treatment on the surface of the The third part is that the bending stress of the stringer under inward bending load at two fulcrums is recorded as sigma M Ultimate load bearing stress of the stringer is sigma b
According to the geometric plastic coefficient of various section bending of chapter 17-22 of book 9 of the aircraft design manual, the geometric plastic coefficient alpha of the I-beam during bending is obtained, and according to the physical geometric plastic coefficient (K) of chapter 17 of book 9 of the aircraft design manual, chapter 17 of the drawing 17-19 P - α curve) "to derive the physical geometric plasticity coefficient K of the stringer upon bending P
Considering stringer compressive stress sigma ch And diagonal tension induced stringer average compressive stress sigma st The maximum bending stress that the stringer can withstand when bending is:
σ max =K P σ bstch
considering stringer compressive stress sigma ch And diagonal tension induced stringer average compressive stress sigma st After the influence of (a), the physical geometric plastic correction coefficient of the stringer when bending is:
calculating a minimum calculated stress sigma for the stringers according to the following formula c
According to the general stability formula of section bar of chapter 21.4 section bar of book 9 of aircraft design manual, the geometric parameters of stringer 2 are used to calculate the general critical instability stress sigma of the stringer cf
Finally, the minimum stress corresponds to the total critical buckling stress sigma of the stringer 2 cf And minimum calculated stress sigma c To yield the overall stability residual strength of the fuselage stringers.

Claims (1)

1. A method for calculating the overall stability of a stringer in a fuselage structure, the fuselage structure comprising a stringer, a skin and a fuselage frame, wherein the design digital-analog and load conditions of the fuselage structure are known, is characterized by comprising the following steps: 1) According to the design counting model of the fuselage structure, taking intersection points of the stringers and fuselage frames as references, establishing a finite element model of the fuselage structure, wherein the skin is simplified into plate units, the stringers are simplified into rod units, the fuselage frames are simplified into beam units, a calculation unit is formed between adjacent beam units by the rod units and the plate units on two sides of the rod units, and each calculation unit comprises the stringers in the middle position and the skin on two sides of the stringers; 2) In the finite element model, the frame beam unit of the machine body close to the machine head is simply supported and restrained to be used as a boundary condition of the overall stress analysis model of the machine body structure; 3) According to the boundary conditions and the load working conditions, stress calculation results of all calculation units of the fuselage structure are obtained; 4) According to the stress calculation results of the calculation units, extracting the calculation units with the skins subjected to pure shearing as analysis units; 5) Calculating the compressive stress of the stringers of each analysis unit, comparing and analyzing the calculation results of the compressive stress of the stringers, and extracting the minimum compressive stress sigma of the stringers c Stringer compressive stress sigma for each analysis unit c Comprises three parts: the first part is the stringer compressive stress sigma calculated in the fuselage structure overall finite element model ch The method comprises the steps of carrying out a first treatment on the surface of the The second part is stringer average compressive stress sigma caused by diagonal stretching after skin shear instability st The method comprises the steps of carrying out a first treatment on the surface of the The third part is the bending stress sigma of the stringer at two fulcrums under inward bending load M The method comprises the steps of carrying out a first treatment on the surface of the 6) Calculating the critical instability stress sigma of the stringer by utilizing the stringer geometric parameter corresponding to the minimum compressive stress cf The method comprises the steps of carrying out a first treatment on the surface of the 7) Critical buckling stress sigma of stringers corresponding to minimum compressive stress cf And minimum compressive stress sigma of stringers c As the overall stability residual strength of the fuselage stringer, the stringer compressive stress sigma is calculated c Consideration should be given to the timePhysical geometric plasticity coefficient K of stringer upon bending P For bending stress sigma M And stringer compressive stress sigma calculated in the fuselage structure overall finite element model ch And stringer mean compressive stress sigma due to diagonal tensile load after skin shear failure st Coefficient of physical geometric plasticity K for stringers when curved P The ultimate load bearing stress of the stringer is sigma b Stringer compressive stress sigma c The specific calculation method of (2) is as follows:
first consider the stringer compressive stress sigma ch And diagonal tension induced stringer average compressive stress sigma st The maximum bending stress sigma that can be sustained when the stringer is bent is calculated according to the following formula max
σ max =K P σ bstch
Considering the stringer compressive stress sigma ch And diagonal tension induced stringer average compressive stress sigma st The physical geometric plastic correction coefficient K of the stringer during bending is calculated according to the following formula P1
Calculating a minimum calculated stress sigma for the stringers according to the following formula c
CN202110547079.XA 2021-05-19 2021-05-19 Method for calculating overall stability of stringers in fuselage structure Active CN113312815B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110547079.XA CN113312815B (en) 2021-05-19 2021-05-19 Method for calculating overall stability of stringers in fuselage structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110547079.XA CN113312815B (en) 2021-05-19 2021-05-19 Method for calculating overall stability of stringers in fuselage structure

Publications (2)

Publication Number Publication Date
CN113312815A CN113312815A (en) 2021-08-27
CN113312815B true CN113312815B (en) 2024-04-02

Family

ID=77373627

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110547079.XA Active CN113312815B (en) 2021-05-19 2021-05-19 Method for calculating overall stability of stringers in fuselage structure

Country Status (1)

Country Link
CN (1) CN113312815B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011030079A2 (en) * 2009-09-14 2011-03-17 Airbus Operation (S.A.S) Method for the structural analysis of panels consisting of an isotropic material, and stiffened by triangular pockets
WO2018028284A1 (en) * 2016-08-09 2018-02-15 苏州数设科技有限公司 Method and device for creating strength model of aircraft structure
WO2019050176A1 (en) * 2017-09-07 2019-03-14 주식회사 엘지화학 Structural analysis tool for mono frame and method for designing mono frame
CN111914351A (en) * 2020-07-06 2020-11-10 西安飞机工业(集团)有限责任公司 Method for calculating overall stability of reinforced wall plate of fuselage structure

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011030079A2 (en) * 2009-09-14 2011-03-17 Airbus Operation (S.A.S) Method for the structural analysis of panels consisting of an isotropic material, and stiffened by triangular pockets
WO2018028284A1 (en) * 2016-08-09 2018-02-15 苏州数设科技有限公司 Method and device for creating strength model of aircraft structure
WO2019050176A1 (en) * 2017-09-07 2019-03-14 주식회사 엘지화학 Structural analysis tool for mono frame and method for designing mono frame
CN111914351A (en) * 2020-07-06 2020-11-10 西安飞机工业(集团)有限责任公司 Method for calculating overall stability of reinforced wall plate of fuselage structure

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Cause analysis and preventive measures for longitudinal cracks in bottom slab of trough girder bridges;Wu Xun;IEEE;20110527;全文 *
基于有限元方法的飞机壁板结构优化设计;姚福印;;科技经济导刊;20170905(第25期);全文 *
民用飞机壁板蒙皮及长桁布置结构优化设计;方阳等;《民用飞机设计与研究》;20170630(第02期);全文 *

Also Published As

Publication number Publication date
CN113312815A (en) 2021-08-27

Similar Documents

Publication Publication Date Title
CN111914351B (en) Method for calculating overall stability of reinforced wallboard of fuselage structure
CN109344437B (en) Reinforced concrete complex stress component reinforcement design method based on force transmission path
CN104330253B (en) A kind of Material Stiffened Panel damage tolerance characteristic analysis method
CN109033526B (en) Calculation method for connecting load of wing rib and skin rivet
CN107506529A (en) A kind of Composite Material Stiffened Panel Axial Compression Stability computational methods
CN112784359B (en) Iterative calculation method for constraint torsion limit bearing capacity of thin-wall beam
CN112989659A (en) Method for establishing surface crack strength factor database based on point weight function method
CN113312815B (en) Method for calculating overall stability of stringers in fuselage structure
Pevzner et al. Calculation of the collapse load of an axially compressed laminated composite stringer-stiffened curved panel–An engineering approach
Kumar et al. Static & dynamic analysis of a typical aircraft wing structure using Msc Nastran
US20130160295A1 (en) Method, apparatus and computer program product for determining the strain induced at a selected point in a stiffened panel structure in response to a load, taking into account one or more out of plane (oop) effects
CN108334690A (en) The reaction beam construction design method of more anchor pole reaction beam load tests and more anchor pole reaction beam load test design methods
CN112504589B (en) Helicopter composite material main blade airfoil section static strength test system and method
CN106326551B (en) Method for calculating effective width of skin in stiffened wall plate structure
CN115482888A (en) Method for predicting crack propagation life under action of pressure-pressure cyclic load
CN113051657B (en) Method for calculating bearing capacity of closed frame beam type fuselage
CN104648689B (en) A kind of static(al)/fatigue test method being applicable to reinforced frame structure
Mert et al. Post-Buckling Load Redistribution of Stiffened Panels in Aircraft Wingbox Structures
CN112699471B (en) Method and device for calculating effective width of skin under axial compression load of fuselage wallboard
CN106934181A (en) A kind of computational methods for wearing the huge component computational length coefficient of layer
He et al. FEA of in-plane compression of aluminum alloy honeycomb panels
Andreyachshenko et al. Simulation of fullering technology as a plastic deformation method for high quality forgings production
CN113049360B (en) Method for determining compression allowable strain value of aircraft composite material reinforced wallboard
RU2243525C1 (en) Method of modeling stressed-deformed state in aviation panel
Shankar et al. Fatigue analysis on wing structure

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