CN115146377A - Coordinated design method for structural rigidity of aircraft with layout of connecting wings - Google Patents

Coordinated design method for structural rigidity of aircraft with layout of connecting wings Download PDF

Info

Publication number
CN115146377A
CN115146377A CN202210726425.5A CN202210726425A CN115146377A CN 115146377 A CN115146377 A CN 115146377A CN 202210726425 A CN202210726425 A CN 202210726425A CN 115146377 A CN115146377 A CN 115146377A
Authority
CN
China
Prior art keywords
wing
rigidity
structural
design
aircraft
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
CN202210726425.5A
Other languages
Chinese (zh)
Other versions
CN115146377B (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202210726425.5A priority Critical patent/CN115146377B/en
Publication of CN115146377A publication Critical patent/CN115146377A/en
Application granted granted Critical
Publication of CN115146377B publication Critical patent/CN115146377B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64FGROUND 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/00Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
    • 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/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • 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)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Manufacturing & Machinery (AREA)
  • Transportation (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a structural rigidity coordination design method of a connecting wing layout airplane, which adopts a piecewise quadratic function to describe bending/torsion rigidity curves of front and rear wings, utilizes a gradient method to quickly optimize a rigidity curve representation parameter, and further determines the structural configuration and the size based on an evolutionary algorithm. The optimization method has the advantages that the first-stage design variables are few, the second stage is pure numerical calculation, the optimization period is short, and the optimal initial structure scheme of the tie wing layout aircraft with the rigidity coordination characteristic can be quickly obtained.

Description

Coordinated design method for structural rigidity of aircraft with layout of connecting wings
Technical Field
The invention belongs to the field of aircraft structure design methods, and particularly relates to a coordinated design method for structural rigidity of a connection wing layout aircraft.
Background
The structural design of the high-aspect-ratio wing generally adopts the thin-wall engineering beam theory, namely the internal force (bending moment, shearing force and torque) of the airfoil section can be calculated through the aerodynamic load, the inertial load and the concentrated force, and the size of the section element can be further determined according to a normal stress and shear flow calculation formula. However, for the joined wing layout aircraft, the front wing and the rear wing integrally form a statically indeterminate double-girder structure, the internal force on the front wing and the rear wing can be influenced by the support reaction force at the joint, and the support reaction force at the joint depends on the aerodynamic load distribution, the inertial load distribution and the rigidity distribution of the front wing and the rear wing, so that the internal force of the front wing section and the rear wing section of the joined wing layout aircraft is related to the rigidity distribution of the front wing and the rear wing. At the moment, the internal force cannot be solved, so that the element size cannot be determined, and the traditional theoretical structural design method of the thin-wall engineering beam cannot be applied to the structural design of the connecting wing.
Therefore, the novel optimization design method is crucial to the structural design of the aircraft with the joined wing layout, and the corresponding optimization method can be particularly adopted to design the structural size of the front wing and the rear wing, but the size is taken as the optimization design variable, so that the design variable is more. If the traditional gradient optimization algorithm is adopted for optimization, a local optimal solution is inevitably trapped; if the genetic algorithm is adopted for optimization, the global optimal solution is easy to guarantee, however, numerous design variables bring huge challenges to the evolutionary algorithm and the engineering application thereof, the optimization design period is long, and the method is not suitable for the initial scheme design stage. Therefore, an efficient optimization design model and a solving method suitable for structural rigidity coordination design of the aircraft with the joined wing layout are urgently needed to be developed.
Disclosure of Invention
Aiming at the problems that the traditional design method of the theoretical structure of the thin-wall engineering beam cannot be applied to the structural design of the aircraft with the layout of the connecting wings, and the traditional optimization method has many design variables and long design period, the invention provides a coordinated design method of the structural rigidity of the aircraft with the layout of the connecting wings, which adopts two-stage optimization: in the first-stage optimization, the bending stiffness and the torsional stiffness curves of the front wing and the rear wing are used as design variables, the number of the design variables is small, the optimization period is short, in the second-stage optimization, mechanical analysis or finite element calling is not needed, pure numerical calculation is carried out, and the calculated amount is small; the high-efficiency coordinated design of the structural rigidity of the aircraft with the joined wing layout is realized.
The technical scheme of the invention is as follows:
the coordinated design method for the structural rigidity of the aircraft with the joined wing layout comprises the following steps:
step 1: finite element modeling is carried out on the structure of the aircraft with the layout of the connecting wings; the bending stiffness curve and the torsional stiffness curve of the front wing and the rear wing are fitted by adopting a piecewise quadratic function;
and 2, step: applying the set pneumatic load and the set inertial load to the finite element model nodes established in the step 1;
and 3, step 3: optimally designing by taking the rigidity description parameters of the front wing and the rear wing as design variables and taking the full-aircraft displacement and the minimum torsion angle as design targets to obtain an optimal bending rigidity curve and an optimal torsion rigidity curve of the front wing and the rear wing;
and 4, step 4: and 3, taking the optimal bending stiffness curve and the optimal torsional stiffness curve of the front wing and the rear wing obtained in the step 3 as design targets, taking the structural forms and the structural sizes of the main beams of the front wing and the rear wing as design variables, and taking the set structural sizes and the set weight as constraint conditions to carry out optimal design so as to obtain the optimal structural forms and the optimal structural sizes of the main beams of the front wing and the rear wing.
Further, in step 1, the finite element modeling of the structure of the aircraft with the bonded wing layout comprises: the main beams of the front wing and the rear wing adopt I-shaped cantilever beam models, and the wing rib section is modeled by adopting concentrated mass points and the beam models to obtain a full-aircraft roof beam fishbone structure model.
Further, in step 1, the bending stiffness curve and the torsional stiffness curve of the front wing and the rear wing are expressed as:
Figure BDA0003711127130000021
Figure BDA0003711127130000022
Figure BDA0003711127130000023
Figure BDA0003711127130000024
Figure BDA0003711127130000025
Figure BDA0003711127130000031
where EI denotes the bending stiffness of the front or rear wing, GJ denotes the torsional stiffness of the front or rear wing, EI root ,GJ root And x root Respectively representing the bending stiffness, torsional stiffness and spanwise attitude, EI, of the front or rear wing root section tip ,GJ tip And x tip Respectively representing the bending stiffness, torsional stiffness and spanwise attitude, EI, of the front or rear wing tip section extreme And x EI,extreme Respectively representing the bending rigidity and the spanwise station, GJ, of the extreme bending rigidity profile of the front wing or the rear wing extreme And x GJ,extreme Respectively representing the torsional rigidity and the spanwise station position of the extreme torsional rigidity profile of the front wing or the rear wing.
Furthermore, in the step 2, the pneumatic load is applied to the finite element model nodes by adopting a three-point picking method, a multi-point picking method, a spline interpolation method or a proxy model method.
Further, in the step 2, in the process of applying the pneumatic load to the finite element model nodes, a CR column algorithm of three-dimensional beam unit analysis based on an updated Lagrange method is adopted, the minimum potential energy principle is used for solving, the structural deformation is calculated, and then the unit internal force is obtained.
Further, in step 3, the stiffness description parameters of the front wing and the rear wing as design variables comprise component stiffness, section stiffness, strain energy density and structural equivalent stiffness of the component.
Further, in step 3, an optimization algorithm for optimizing based on gradient information is adopted for optimization design, and the optimization algorithm for optimizing based on gradient information comprises a feasible direction method or a sequential quadratic programming method.
Further, in step 4, the structural types of the main beams serving as design variables comprise circular tube beams and box-type beams; when the circular tubular beam is adopted as a main beam structure type, the main beam structure size design variables are the outer diameter and the wall thickness; when box beams are used as the main beam structural style, the main beam structural dimensional design variables are width, height, web thickness, and flange thickness.
Further, in step 4, an evolutionary algorithm is adopted for optimization design, wherein the evolutionary algorithm comprises a genetic algorithm, a particle swarm algorithm or a fish swarm algorithm.
Advantageous effects
The invention has the beneficial effects that: according to the structural rigidity coordination design method of the aircraft with the coupled wing layout, the bending/torsion rigidity curves of the front wing and the rear wing are described by adopting a piecewise quadratic function, the characteristic parameters of the rigidity curves are rapidly optimized by utilizing a gradient method, and further structural configuration and size are determined based on an evolutionary algorithm. The optimization method has the advantages that the first-stage design variables are few, the second stage is pure numerical calculation, the optimization period is short, and the optimal initial structure scheme of the tie wing layout aircraft with the rigidity coordination characteristic can be quickly obtained.
The invention has stronger applicability: (1) In the first-stage optimization, the bending stiffness and torsional stiffness curves of the front wing and the rear wing are used as design variables, the number of the design variables is small, and the optimization period is short; (2) In the second-stage optimization, mechanical analysis or finite element calling is not needed, pure numerical calculation is achieved, and the calculated amount is small. The two-stage optimization concept is not only suitable for the structural design of the aircraft with the layout of the connecting wings, but also suitable for the structural design of the wings of the conventional aircraft.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The above and/or additional aspects and advantages of the present invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
FIG. 1: coupling wing layout aerodynamic layout example of the present invention
FIG. 2: the invention relates to a method for describing a bending and torsion stiffness curve by a piecewise quadratic function
FIG. 3: front wing beam fish bone structure model of optimized embodiment of the invention
FIG. 4 is a schematic view of: rear wing spar type fishbone structure model of optimized embodiment of the invention
FIG. 5: 'multipoint picking' pneumatic load applying method of optimized embodiment of the invention
FIG. 6: schematic diagram of circular tube main beam structure scheme of optimized embodiment of the invention
FIG. 7: schematic diagram of box type main beam structure scheme of optimized embodiment of invention
FIG. 8: genetic algorithm implementation procedure for optimized embodiments of the invention
In the figure: 1. a front wing outer section; 2. a rear wing; 3. hanging the tail; 4. an end plate; 5. a radar antenna; 6. a front wing.
Detailed Description
The following detailed description of embodiments of the invention is intended to be illustrative, and not to be construed as limiting the invention.
Fig. 1 shows a typical joined wing aircraft, which includes a front wing and a rear wing connected by an end plate. Unlike conventional large aspect ratio wing engineering beam force transfer, the coupling wing layout front wing and rear wing force transfer paths are unclear. The bending stiffness and the torsional stiffness of the conventional large-exhibition specific wing are a curve which monotonically decreases, but the optimal stiffness curve mathematical characteristic of the layout of the tie wing is unknown, and the monotonically decreasing function description is not suitable for the stiffness curves of the front wing and the rear wing of the layout of the tie wing. Therefore, in order to better adapt to the rigidity matching requirement of the front wing and the rear wing and widen the design space, when the finite element modeling of the structure of the aircraft with the layout of the connecting wings is carried out, the bending rigidity curve and the torsional rigidity curve of the front wing and the rear wing adopt the piecewise quadratic function fitting and are defined by parameterization, and the considered parameters are the position of the segmentation point of the piecewise quadratic function, the front slope and the rear slope of the segmentation point and the positioning parameters of the quadratic functions in different segmentation areas.
The specific piecewise quadratic function is:
Figure BDA0003711127130000051
Figure BDA0003711127130000052
Figure BDA0003711127130000053
Figure BDA0003711127130000054
Figure BDA0003711127130000055
Figure BDA0003711127130000056
where EI denotes the bending stiffness of the front or rear wing, GJ denotes the torsional stiffness of the front or rear wing, EI root ,GJ root And x root Respectively representing the bending stiffness, torsional stiffness and spanwise attitude, EI, of the front or rear wing root section tip ,GJ tip And x tip Respectively representing the bending stiffness, torsional stiffness and spanwise attitude, EI, of the front or rear wing tip section extreme And x EI,extreme Respectively representing the bending rigidity and the spanwise station, GJ, of the extreme bending rigidity profile of the front wing or the rear wing extreme And x GJ,extreme Respectively representing the torsional rigidity and the spanwise station position of the extreme torsional rigidity profile of the front wing or the rear wing.
It can be seen that the characteristic variables of the bending stiffness curves of the front wing and the rear wing are as follows: EI (electronic instrument) root 、EI tip 、x EI,extreme 、EI extreme The torsional rigidity curve characterization variables of the front wing and the rear wing are respectively as follows: GJ root 、GJ tip 、x GJ,extreme 、GJ extreme A total of 16 characterizing variables.
In the specific finite element modeling of the structure of the aircraft with the layout of the connecting wings, the main beams of the front wing and the rear wing adopt I-shaped cantilever beam models, and the wing rib section is modeled by adopting concentrated mass points and beam models to obtain a full-aircraft Liangshi fishbone structure model. The front wing and the rear wing structural models are shown in fig. 3 and 4.
Then applying the set pneumatic load and the set inertial load to the nodes of the established finite element model; the pneumatic load can be applied to the finite element model nodes by adopting a three-point picking method, a multi-point picking method, a spline interpolation method or a proxy model method. In the process of applying the pneumatic load to the finite element model nodes, a CR column algorithm of three-dimensional beam unit analysis based on an updated Lagrange method is adopted, the minimum potential energy principle is applied to solve, structural deformation is calculated, and then unit internal force is obtained.
In this embodiment, a "multi-point picking" method is used to apply a pneumatic load to a structural node of a finite element model, and the specific process is as follows:
assuming that a virtual beam element exists between all the structural points and the pneumatic points, as shown in fig. 7, the virtual beam element is a cantilever beam with one end of the pneumatic point being a fixed support, and the structural point at the free end of the cantilever beam is distributed to the load P i The deformation energy is
Figure BDA0003711127130000061
Where EJ is the bending stiffness of the virtual beam element. l i Is the length of the ith cantilever beam, u i Is the deformation energy of the ith cantilever beam。
The deformation of the entire system can then be
Figure BDA0003711127130000062
Wherein U is the total deformation energy of the system, and n is the total number of the cantilever beams.
The load distributed to the structural points is such that the deformation energy of the system is minimal and the static equivalence conditions are met, i.e.
Figure BDA0003711127130000063
In the formula, X i ,Z i The x and z coordinates of node i, respectively.
Lagrange multiplier method for establishing Lagrange function
Figure BDA0003711127130000071
Figure BDA0003711127130000072
In the formula, λ xz Is a Lagrange multiplier;
Figure BDA0003711127130000073
making 3EJ =1, namely
Figure BDA0003711127130000074
Substituting the formula (6) into the formula (3) to obtain
Figure BDA0003711127130000075
From formula (7) to give λ, λ xz And (4) then, replacing the structure points with the formula (6) to obtain the loads distributed by the structure points, wherein the loads of each pneumatic force point on the whole pneumatic network point group are distributed to the structure point group according to the method to obtain the node pneumatic force load array.
And then, taking rigidity description parameters of the front wing and the rear wing as design variables, including component rigidity, section rigidity, strain energy density and structural equivalent rigidity of each part, and taking the minimum of full-aircraft displacement and torsion angle as a design target, and performing optimization design by adopting an optimization algorithm for optimizing based on gradient information, such as a feasible direction method or a sequential quadratic programming method, to obtain a bending rigidity curve and a torsion rigidity curve of the front wing and the rear wing, which are matched with the optimal front wing and the optimal rear wing.
And then, performing optimal design by taking the obtained optimal bending stiffness curve and torsional stiffness curve of the front wing and the rear wing as design targets, taking the structural types and the structural sizes of the main beams of the front wing and the rear wing as design variables and taking the set structural sizes and weights as constraint conditions to obtain the optimal structural types and structural sizes of the main beams of the front wing and the rear wing.
The structural types of the main beam as design variables comprise a round tubular beam and a box beam; when the circular tubular beam is adopted as a main beam structure, the design variables of the main beam structure size are the outer diameter D and the wall thickness t; when a box beam is used as the main beam structure, the design variables of the main beam structure are width b, height h, web thickness t1 and flange plate thickness t2. And the evolutionary algorithm can adopt a genetic algorithm, a particle swarm algorithm or a fish swarm algorithm. The evolutionary algorithm adopted in the embodiment is a real number coding genetic algorithm, the crossover operator adopts analog binary crossover, the mutation operator adopts polynomial mutation and keeps consistent with the reference setting in the NSGA-II algorithm, and the implementation steps of the genetic algorithm are shown in FIG. 8.
In the embodiment, the final optimization result is reduced by 60% compared with the initial scheme, the torsion angle is reduced by 80%, the first-stage optimization is performed, the design variables are reduced from hundreds of dimensions to more than ten, the optimization design period is reduced by 80%, and the efficient coordinated design of the structural rigidity of the aircraft with the joined wing layout is realized
Although embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are exemplary and not to be construed as limiting the present invention, and that those skilled in the art may make variations, modifications, substitutions and alterations within the scope of the present invention without departing from the spirit and scope of the present invention.

Claims (9)

1. A coordinated design method for structural rigidity of a connecting wing layout aircraft is characterized by comprising the following steps: the method comprises the following steps:
step 1: finite element modeling of the structure of the aircraft with the joined wing layout is carried out; the bending stiffness curve and the torsional stiffness curve of the front wing and the rear wing are fitted by adopting a piecewise quadratic function;
and 2, step: applying the set pneumatic load and the set inertial load to the finite element model nodes established in the step 1;
and step 3: optimally designing by taking the rigidity description parameters of the front wing and the rear wing as design variables and taking the full-aircraft displacement and the minimum torsion angle as design targets to obtain an optimal bending rigidity curve and an optimal torsion rigidity curve of the front wing and the rear wing;
and 4, step 4: and (4) taking the optimal bending stiffness curve and torsional stiffness curve of the front wing and the rear wing obtained in the step (3) as design targets, taking the structural types and the structural sizes of the main beams of the front wing and the rear wing as design variables, and taking the set structural sizes and weights as constraint conditions to carry out optimal design, so as to obtain the optimal structural types and the structural sizes of the main beams of the front wing and the rear wing.
2. A method for coordinately designing structural rigidity of a tie wing aircraft according to claim 1, wherein: in step 1, the finite element modeling of the structure of the aircraft with the layout of the connecting wing comprises the following steps: the main beams of the front wing and the rear wing adopt I-shaped cantilever beam models, and the wing rib section is modeled by adopting concentrated mass points and the beam models to obtain a full-aircraft roof beam fishbone structure model.
3. A coordinated design method for structural rigidity of a join wing layout aircraft according to claim 1 or 2, characterized in that: in step 1, the bending stiffness curve and the torsional stiffness curve of the front wing and the rear wing are expressed as:
Figure FDA0003711127120000011
Figure FDA0003711127120000012
Figure FDA0003711127120000013
Figure FDA0003711127120000014
Figure FDA0003711127120000021
Figure FDA0003711127120000022
where EI denotes the bending stiffness of the front or rear wing, GJ denotes the torsional stiffness of the front or rear wing, EI root ,GJ root And x root Respectively, the bending rigidity, torsional rigidity and the spanwise position, EI, of the root section of the front wing or the rear wing tip ,GJ tip And x tip Respectively, the bending rigidity, torsional rigidity and the spanwise position, EI, of the front wing or rear wing tip section extreme And x EI,extreme Respectively representing the bending rigidity and the spanwise station, GJ, of the extreme bending rigidity profile of the front wing or the rear wing extreme And x GJ,extreme Respectively representing the torsional rigidity and the spanwise station position of the extreme torsional rigidity profile of the front wing or the rear wing.
4. A method for coordinately designing structural rigidity of a tie wing aircraft according to claim 1, wherein: and step 2, applying the pneumatic load to the finite element model nodes by adopting a three-point picking method, a multi-point picking method, a spline interpolation method or a proxy model method.
5. A coordinated design method for structural rigidity of a join wing layout aircraft according to claim 4 is characterized in that: in the step 2, in the process of applying the pneumatic load to the finite element model nodes, a CR column algorithm of three-dimensional beam unit analysis based on an updated Lagrange method is adopted, the minimum potential energy principle is applied to solve, structural deformation is calculated, and then unit internal force is obtained.
6. A method for coordinately designing structural rigidity of a tie wing aircraft according to claim 1, wherein: in step 3, the rigidity description parameters of the front wing and the rear wing as design variables comprise component rigidity, section rigidity, strain energy density and structural equivalent rigidity of the part.
7. A method for coordinately designing structural rigidity of a tie wing aircraft according to claim 1, wherein: in step 3, an optimization algorithm for optimizing based on gradient information is adopted for optimization design, wherein the optimization algorithm for optimizing based on gradient information comprises a feasible direction method or a sequence quadratic programming method.
8. A method for coordinately designing structural rigidity of a tie wing aircraft according to claim 1, wherein: in step 4, the structural forms of the main beams serving as design variables comprise circular tube beams and box-type beams; when the circular tubular beam is adopted as a main beam structure type, the main beam structure size design variables are the outer diameter and the wall thickness; when box beams are used as the main beam structural style, the main beam structural dimensional design variables are width, height, web thickness, and flange thickness.
9. A coordinated design method for structural rigidity of a tie wing aircraft according to claim 1, characterized in that: and 4, performing optimization design by adopting an evolutionary algorithm, wherein the evolutionary algorithm comprises a genetic algorithm, a particle swarm algorithm or a fish swarm algorithm.
CN202210726425.5A 2022-06-23 2022-06-23 Coordinated design method for structural rigidity of aircraft with layout of connecting wings Active CN115146377B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210726425.5A CN115146377B (en) 2022-06-23 2022-06-23 Coordinated design method for structural rigidity of aircraft with layout of connecting wings

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210726425.5A CN115146377B (en) 2022-06-23 2022-06-23 Coordinated design method for structural rigidity of aircraft with layout of connecting wings

Publications (2)

Publication Number Publication Date
CN115146377A true CN115146377A (en) 2022-10-04
CN115146377B CN115146377B (en) 2023-02-24

Family

ID=83408897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210726425.5A Active CN115146377B (en) 2022-06-23 2022-06-23 Coordinated design method for structural rigidity of aircraft with layout of connecting wings

Country Status (1)

Country Link
CN (1) CN115146377B (en)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103543630A (en) * 2012-07-17 2014-01-29 动力专家有限公司 Method for improving concentricity when using spiral spring mechanical watch and spiral spring
CN105975706A (en) * 2016-05-16 2016-09-28 中国航空工业集团公司西安飞机设计研究所 Scheme stage wing parameter estimation method
CN108491576A (en) * 2018-02-12 2018-09-04 北京航空航天大学 A kind of optimum design method of composite wing Cutout reinforcement
CN109409009A (en) * 2018-11-30 2019-03-01 中国航空工业集团公司沈阳飞机设计研究所 A kind of Torsion Coupling wing structure optimization method
CN110334427A (en) * 2019-06-27 2019-10-15 南京航空航天大学 A kind of TBW airplane wing structure analysis method based on engineering beam theory
CN113044200A (en) * 2021-04-20 2021-06-29 西北工业大学 Novel solar unmanned aerial vehicle layout with connection wings close to space
CN113076677A (en) * 2021-04-15 2021-07-06 朱礼云 Flexible body structure high-order nonlinear finite element numerical simulation method based on quintic Hermite shape function
CN113361017A (en) * 2021-07-03 2021-09-07 西北工业大学 Pneumatic design method for secondary folding wings of pipe-jet unmanned aerial vehicle considering mechanism constraint
CN113591307A (en) * 2021-07-30 2021-11-02 郑州轻工业大学 Double-layer optimization method for integrating appearance and internal mechanism of underwater glider by wing body
US20210363961A1 (en) * 2018-04-28 2021-11-25 The Research Foundation For The State University Of New York Flexible wind turbine blade with actively variable twist distribution

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103543630A (en) * 2012-07-17 2014-01-29 动力专家有限公司 Method for improving concentricity when using spiral spring mechanical watch and spiral spring
CN105975706A (en) * 2016-05-16 2016-09-28 中国航空工业集团公司西安飞机设计研究所 Scheme stage wing parameter estimation method
CN108491576A (en) * 2018-02-12 2018-09-04 北京航空航天大学 A kind of optimum design method of composite wing Cutout reinforcement
US20210363961A1 (en) * 2018-04-28 2021-11-25 The Research Foundation For The State University Of New York Flexible wind turbine blade with actively variable twist distribution
CN109409009A (en) * 2018-11-30 2019-03-01 中国航空工业集团公司沈阳飞机设计研究所 A kind of Torsion Coupling wing structure optimization method
CN110334427A (en) * 2019-06-27 2019-10-15 南京航空航天大学 A kind of TBW airplane wing structure analysis method based on engineering beam theory
CN113076677A (en) * 2021-04-15 2021-07-06 朱礼云 Flexible body structure high-order nonlinear finite element numerical simulation method based on quintic Hermite shape function
CN113044200A (en) * 2021-04-20 2021-06-29 西北工业大学 Novel solar unmanned aerial vehicle layout with connection wings close to space
CN113361017A (en) * 2021-07-03 2021-09-07 西北工业大学 Pneumatic design method for secondary folding wings of pipe-jet unmanned aerial vehicle considering mechanism constraint
CN113591307A (en) * 2021-07-30 2021-11-02 郑州轻工业大学 Double-layer optimization method for integrating appearance and internal mechanism of underwater glider by wing body

Non-Patent Citations (11)

* Cited by examiner, † Cited by third party
Title
YANG ZHANG ET AL.: ""Aerodynamic performance of a low-Reynolds UAV with leading-edge protuberances inspired by humpback whale flippers"", 《CHINESE JOURNAL OF AERONAUTICS》 *
乔巍: ""联结翼复合材料机翼结构优化设计与重量特性分析"", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *
何萌等: ""翼型气动优化设计问题多极值特性研究"", 《航空工程进展》 *
刘波等: ""复合材料无人机机翼布局的二级优化方法"", 《玻璃钢/复合材料》 *
尹钧等: ""联结翼飞机主要布局参数对全机升阻特性影响研究"", 《航空科学技术》 *
张仲桢等: ""飞翼式飞行器结构布局与构件尺寸的两级优化"", 《计算机辅助工程》 *
王一为等: ""飞翼布局飞行器结构拓扑优化设计"", 《北京航空航天大学学报》 *
石运国: ""试验火箭尾翼的气动—结构协同优化设计研究"", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *
石钢等: ""基于工程梁理论的联结翼结构重量特性"", 《南京航空航天大学学报》 *
金朋等: ""基于整数编码并行遗传算法的复合材料螺旋桨结构优化设计"", 《机械强度》 *
陈颂: ""基于梯度的气动外形优化设计方法及应用"", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 *

Also Published As

Publication number Publication date
CN115146377B (en) 2023-02-24

Similar Documents

Publication Publication Date Title
CN109711048B (en) Wing design method considering aerodynamic and structural multi-type design variables
CN112016167B (en) Aircraft aerodynamic shape design method and system based on simulation and optimization coupling
CN106845019B (en) Self-adaptive wing profile design method
Vasista et al. Compliant structures-based wing and wingtip morphing devices
Bottasso et al. Integrated aero-structural optimization of wind turbines
CN115374543B (en) Aerodynamic/structural multidisciplinary design optimization method for Lambda wings
CN115630542B (en) Reinforced layout optimization method of thin-wall reinforced structure
CN111027250B (en) Modeling method for special-shaped curved surface reinforcement shell based on grid deformation technology
CN115146377B (en) Coordinated design method for structural rigidity of aircraft with layout of connecting wings
Wu et al. Aero-structure-stealth coupled optimization for high aspect ratio wing using adaptive metamodeling method
Gandhi et al. Rapid design generation and multifidelity analysis of aircraft structures
CN110096798B (en) Multi-state finite element model correction method
CN114282410B (en) Deformation wing finite element modeling method considering aerodynamic structure coordination updating
CN113361017B (en) Pneumatic design method for secondary folding wings of pipe-jet unmanned aerial vehicle considering mechanism constraint
De Gaspari et al. Final validation and performance assessment of a full-scale morphing droop nose demonstrator
Rumpfkeil et al. Multi-fidelity, Aeroelastic Analysis and Optimization with Control Surface Deflections of an Efficient Supersonic Air Vehicle
Winter et al. NURBS-based shape and parameter optimization of structural components with an adaptive amount of control points
Zhang et al. Combining density-based approach and optimization refinement in the design of morphing airfoil structures
Amrit et al. Aerodynamic design exploration through point-by-point pareto set identification using local surrogate models
De Gaspari et al. A two levels approach for the optimal design of morphing airfoils
Cavalieri et al. Aero-Structural Design Optimization of a Morphing Aileron Considering Actuation Aspects
CN118504144B (en) Variable camber pneumatic optimization design method for aircraft wing with multiple trailing edge control surfaces
CN117669016B (en) Global/local coupled wide-speed-domain waverider-wing layout optimization design method
CN115688528A (en) Carrier rocket aerodynamic-load-structure finite element optimization design method
Barnewitz Flexible wing optimisation based on shapes and structures

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