CN113032892B - Skin panel parameter optimization method based on stability algorithm - Google Patents

Skin panel parameter optimization method based on stability algorithm Download PDF

Info

Publication number
CN113032892B
CN113032892B CN202011612917.9A CN202011612917A CN113032892B CN 113032892 B CN113032892 B CN 113032892B CN 202011612917 A CN202011612917 A CN 202011612917A CN 113032892 B CN113032892 B CN 113032892B
Authority
CN
China
Prior art keywords
thickness
skin
sigma
increasing
unit
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
CN202011612917.9A
Other languages
Chinese (zh)
Other versions
CN113032892A (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 First Aircraft Institute
Original Assignee
AVIC First Aircraft Institute
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 First Aircraft Institute filed Critical AVIC First Aircraft Institute
Priority to CN202011612917.9A priority Critical patent/CN113032892B/en
Publication of CN113032892A publication Critical patent/CN113032892A/en
Application granted granted Critical
Publication of CN113032892B publication Critical patent/CN113032892B/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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/28Fuselage, exterior or interior
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Shaping Metal By Deep-Drawing, Or The Like (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a skin panel parameter optimization method based on a stability algorithm, which comprises the following steps: step 1, calculating equivalent stress and critical buckling stress of a skin unit according to an initial model of a skin panel structure of a metal wing, and calculating a stress ratio of the skin unit; step 2, when the stress ratio of the skin unit is greater than or equal to 1, setting a thickness increasing coefficient for increasing the thickness of the skin unit; and 3, determining to adopt the thickness increasing coefficient or the minimum amount of the skin thickness processing technology to carry out thickness increasing treatment on the skin units according to the comparison of the thickness increasing amount formed by the skin units by adopting the thickness increasing coefficient set in the step 2 and the minimum amount of the skin thickness processing technology. The embodiment of the invention solves the problems that the traditional skin panel parameter optimization method highly depends on the experience of designers and has huge manual parameter adjustment workload because the stability is only checked and does not participate in initiative in the parameterization process of the traditional skin panel.

Description

一种基于稳定性算法的蒙皮壁板参数优化方法A Stability Algorithm-Based Parameter Optimization Method for Skin Panels

技术领域technical field

本发明涉及但不限于飞机结构设计及强度分析技术领域,尤指一种基于稳定性算法的蒙皮壁板参数优化方法。The invention relates to, but is not limited to, the technical field of aircraft structure design and strength analysis, in particular to a method for optimizing skin panel parameters based on a stability algorithm.

背景技术Background technique

传统飞机金属机翼蒙皮壁板参数优化过程中,优化思路为:以蒙皮壁板的厚度参数作为优化变量,在制造工艺限制下以满应力水平为约束条件,再校核壁板稳定性,最终以壁板重量最小为优化目标。In the process of traditional aircraft metal wing skin panel parameter optimization, the optimization idea is as follows: the thickness parameter of the skin panel is used as the optimization variable, and the full stress level is the constraint condition under the constraints of the manufacturing process, and then the stability of the panel is checked. , the final optimization goal is to minimize the weight of the siding.

上述壁板参数优化方式中,壁板稳定性仅作为一种被动的强度校核作用,而不是作为主动设计参考;且现有参数化过程仅能手动修改不满足稳定性的单个单元的厚度参数,改动依据需较强经验性,完全依赖设计人员的主观经验,调参工作量大。In the above method of optimizing the panel parameters, the panel stability is only used as a passive strength check, not as an active design reference; and the existing parameterization process can only manually modify the thickness parameters of a single element that does not meet the stability requirements. , the change basis needs strong experience, completely depends on the subjective experience of the designer, and the workload of parameter adjustment is large.

发明内容SUMMARY OF THE INVENTION

本发明的目的为:本发明实施例提供一种基于稳定性算法的蒙皮壁板参数优化方法,以解决传统蒙皮壁板参数化过程中,由于稳定性仅作校核不参与主动,从而使得传统的蒙皮壁板参数优化方法高度依赖设计人员的经验且手动调参工作量巨大的问题。The purpose of the present invention is as follows: the embodiment of the present invention provides a method for optimizing the parameters of the skin panel based on the stability algorithm, so as to solve the problem that in the process of parameterization of the traditional skin panel, the stability is only checked and does not participate in the initiative. This makes the traditional skin panel parameter optimization method highly dependent on the designer's experience and the huge workload of manual parameter adjustment.

本发明的技术方案为:本发明实施例提供一种基于稳定性算法的蒙皮壁板参数优化方法,包括:The technical solution of the present invention is as follows: the embodiment of the present invention provides a method for optimizing the parameters of a skin panel based on a stability algorithm, including:

步骤1,根据金属机翼蒙皮壁板结构的初始模型,计算蒙皮单元的当量应力σ和临界失稳应力σc,并计算所述蒙皮单元的应力比值σ/σcStep 1, according to the initial model of the metal wing skin panel structure, calculate the equivalent stress σ and the critical instability stress σ c of the skin element, and calculate the stress ratio σ/σ c of the skin element;

步骤2,在所述蒙皮单元的应力比值σ/σc大于或等于1时,设置用于对所述蒙皮单元的厚度进行增大的厚度增大系数;Step 2, when the stress ratio σ/σ c of the skin unit is greater than or equal to 1, set a thickness increase coefficient for increasing the thickness of the skin unit;

步骤3,根据采用步骤2中设置的厚度增大系数对所述蒙皮单元形成的厚度增大量与蒙皮厚度加工工艺的最小量的对比,确定采用所述厚度增大系数或蒙皮厚度加工工艺的最小量对所述蒙皮单元进行厚度增大处理。Step 3, according to the comparison between the thickness increase amount formed by the thickness increase coefficient set in step 2 and the minimum amount of the skin thickness processing technology for the skin unit, determine to use the thickness increase coefficient or skin thickness processing The skin element is thickened with a minimum amount of process.

可选地,如上所述的基于稳定性算法的蒙皮壁板参数优化方法中,所述步骤2中设置厚度增大系数的方式为:Optionally, in the above-mentioned method for optimizing the parameters of the skin panel based on the stability algorithm, the method for setting the thickness increase coefficient in the step 2 is:

设置厚度增大系数为(σ/σc)1/3;其中,对每个独立的蒙皮单元,厚度增大(σ/σc)1/3后,相应蒙皮单元的当量应力σ等于临界失稳应力σcSet the thickness increase factor to (σ/σ c ) 1/3 ; where, for each individual skin element, after the thickness is increased by (σ/σ c ) 1/3 , the equivalent stress σ of the corresponding skin element is equal to Critical buckling stress σ c .

可选地,如上所述的基于稳定性算法的蒙皮壁板参数优化方法中,所述步骤3中获取蒙皮单元的厚度增大量的方式为:Optionally, in the above-mentioned method for optimizing the parameters of the skin panel based on the stability algorithm, the method for obtaining the thickness increase of the skin unit in the step 3 is:

步骤31,计算采用厚度增大系数(σ/σc)1/3对所述蒙皮单元进行厚度增大形成的厚度增大量((σ/σc)1/3-1)t;其中,t为所述蒙皮单元的壁板厚度。Step 31: Calculate the thickness increase amount ((σ/σ c ) 1/3 -1)t formed by increasing the thickness of the skin unit by using the thickness increase coefficient (σ/σ c ) 1/3 ; wherein, t is the wall thickness of the skin unit.

可选地,如上所述的基于稳定性算法的蒙皮壁板参数优化方法中,所述步骤3中确定对所述蒙皮单元进行厚度增大的方式,包括:Optionally, in the above-mentioned method for optimizing the parameters of the skin panel based on the stability algorithm, in the step 3, the method of increasing the thickness of the skin unit is determined, including:

步骤32,判断厚度增大量是否大于或等于0.2mm,0.2mm为蒙皮厚度加工工艺的最小量;Step 32, determine whether the thickness increase is greater than or equal to 0.2mm, and 0.2mm is the minimum amount of the skin thickness processing technology;

步骤33,在所述厚度增大量大于或等于0.2mm时,采用所述厚度增大系数对所述蒙皮单元进行厚度增大;Step 33, when the thickness increase is greater than or equal to 0.2 mm, use the thickness increase coefficient to increase the thickness of the skin unit;

步骤34,在所述厚度增大量小于0.2mm时,以0.2mm为厚度增大量对所述蒙皮单元进行厚度增大。Step 34 , when the thickness increase amount is less than 0.2 mm, increase the thickness of the skin unit by taking 0.2 mm as the thickness increase amount.

可选地,如上所述的基于稳定性算法的蒙皮壁板参数优化方法中,所述步骤33中采用所述厚度增大系数对所述蒙皮单元进行厚度增大的方式,包括:Optionally, in the above-mentioned method for optimizing the parameters of the skin panel based on the stability algorithm, the method of increasing the thickness of the skin unit by using the thickness increase coefficient in the step 33 includes:

对所述厚度增大系数设置加速因子k,通过加速因子k将所述厚度增大系数(σ/σc)1/3放大k倍,并通过k(σ/σc)1/3对所述蒙皮单元进行厚度增大。An acceleration factor k is set for the thickness increase coefficient, the thickness increase coefficient (σ / σ c ) 1/3 is amplified by k times by the acceleration factor k, and the The thickness of the skin element is increased.

可选地,如上所述的基于稳定性算法的蒙皮壁板参数优化方法中,还包括:Optionally, in the above-mentioned method for optimizing the parameters of the skin panel based on the stability algorithm, the method further includes:

步骤4,通过重复执行步骤2到步骤3对所述蒙皮单元进行厚度增大,指定所述蒙皮单元的应力比值σ/σc小于1。Step 4, by repeating steps 2 to 3 to increase the thickness of the skin element, specifying that the stress ratio σ/σ c of the skin element is less than 1.

可选地,如上所述的基于稳定性算法的蒙皮壁板参数优化方法中,所述加速因子k的取值为1到2之间;Optionally, in the above-mentioned method for optimizing the parameters of the skin panel based on the stability algorithm, the value of the acceleration factor k is between 1 and 2;

所述蒙皮壁板参数优化方法在每次执行步骤2到步骤3的过程中,通过设置加速因子k的取值,调整所述蒙皮单元的厚度增大量,从而调整蒙皮壁板参数优化的速度。In the process of performing step 2 to step 3 each time in the method for optimizing the parameters of the skin panel, by setting the value of the acceleration factor k, the thickness increase of the skin unit is adjusted, so as to adjust the optimization of the parameters of the skin panel. speed.

可选地,如上所述的基于稳定性算法的蒙皮壁板参数优化方法中,还包括:Optionally, in the above-mentioned method for optimizing the parameters of the skin panel based on the stability algorithm, the method further includes:

步骤2a,在所述蒙皮单元的应力比值σ/σc小于1时,相应的蒙皮单元符合稳定性控制,优化结束。Step 2a, when the stress ratio σ/ σc of the skin element is less than 1, the corresponding skin element conforms to the stability control, and the optimization ends.

本发明的有益效果为:与传统的金属机翼蒙皮壁板参数优化方法相比,本发明实施例提供的一种基于稳定性算法的蒙皮壁板参数优化方法,一方面,结合理论和工程经验提出蒙皮壁板参数的加速优化方法,创造性的将稳定性校核控制推导为稳定性主动设计,所推导出的增大系数为壁板参数优化稳定性控制指明了优化方向;另一方面,加速因子的设置极大的减少了优化迭代次数,将传统优化方法的手动调节几个小时的工作缩短至几分钟,极大地提高了优化效率。本发明实施例提供的基于稳定性算法的蒙皮壁板参数优化方法,将稳定性设计作为一种主动设计因素引入优化流程,设计加速迭代系数,实现参数快速优化。。The beneficial effects of the present invention are: compared with the traditional method for optimizing the parameters of the skin panel of a metal wing, the embodiment of the present invention provides a method for optimizing the parameters of the skin panel based on the stability algorithm. Engineering experience proposes an accelerated optimization method for skin panel parameters, creatively deduces stability check control as stability active design, and the derived increase coefficient points out the optimization direction for panel parameter optimization stability control; another On the one hand, the setting of the acceleration factor greatly reduces the number of optimization iterations, shortens the manual adjustment work of traditional optimization methods from several hours to several minutes, and greatly improves the optimization efficiency. The skin panel parameter optimization method based on the stability algorithm provided by the embodiment of the present invention introduces the stability design as an active design factor into the optimization process, designs the accelerated iteration coefficient, and realizes the rapid optimization of parameters. .

附图说明Description of drawings

附图用来提供对本发明技术方案的进一步理解,并且构成说明书的一部分,与本申请的实施例一起用于解释本发明的技术方案,并不构成对本发明技术方案的限制。The accompanying drawings are used to provide a further understanding of the technical solutions of the present invention, and constitute a part of the specification. They are used to explain the technical solutions of the present invention together with the embodiments of the present application, and do not limit the technical solutions of the present invention.

图1为本发明实施例提供的一种基于稳定性算法的蒙皮壁板参数优化方法的流程图;1 is a flowchart of a method for optimizing skin panel parameters based on a stability algorithm according to an embodiment of the present invention;

图2为本发明实施例提供的另一种基于稳定性算法的蒙皮壁板参数优化方法的流程图。FIG. 2 is a flowchart of another method for optimizing parameters of a skin panel based on a stability algorithm according to an embodiment of the present invention.

具体实施例specific embodiment

为使本发明的目的、技术方案和优点更加清楚明白,下文中将结合附图对本发明的实施例进行详细说明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。In order to make the objectives, technical solutions and advantages of the present invention clearer, the embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be noted that, the embodiments in the present application and the features in the embodiments may be arbitrarily combined with each other if there is no conflict.

上述背景技术中已经说明,壁板稳定性仅作为一种被动的强度校核作用,而不是作为主动设计参考;且现有参数化过程仅能手动修改不满足稳定性的单个单元的厚度参数,改动依据需较强经验性,完全依赖设计人员的主观经验,调参工作量大。It has been explained in the above background art that the stability of the panel is only used as a passive strength check function, rather than as an active design reference; and the existing parameterization process can only manually modify the thickness parameters of a single unit that does not meet the stability. Changes are based on strong experience, completely relying on the subjective experience of designers, and the workload of parameter adjustment is large.

基于上述问题,本发明实施例提供一种基于稳定性算法的蒙皮壁板参数优化方法,将稳定性设计作为一种主动设计因素引入优化流程,设计加速迭代系数,实现参数快速优化。Based on the above problems, embodiments of the present invention provide a skin panel parameter optimization method based on a stability algorithm, which incorporates stability design as an active design factor into the optimization process, designs accelerated iteration coefficients, and realizes rapid parameter optimization.

本发明提供以下几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例不再赘述。The present invention provides the following specific embodiments that can be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments.

图1为本发明实施例提供的一种基于稳定性算法的蒙皮壁板参数优化方法的流程图。本发明实施例提供的基于稳定性算法的蒙皮壁板参数优化方法,可以包括如下步骤:FIG. 1 is a flowchart of a method for optimizing parameters of a skin panel based on a stability algorithm according to an embodiment of the present invention. The skin panel parameter optimization method based on the stability algorithm provided by the embodiment of the present invention may include the following steps:

步骤1,根据金属机翼蒙皮壁板结构的初始模型,计算蒙皮单元的当量应力σ和临界失稳应力σc,并计算蒙皮单元的应力比值σ/σcStep 1, according to the initial model of the metal wing skin panel structure, calculate the equivalent stress σ and critical instability stress σ c of the skin element, and calculate the stress ratio σ/σ c of the skin element;

步骤2,在蒙皮单元的应力比值σ/σc大于或等于1时,设置用于对蒙皮单元的厚度进行增大的厚度增大系数;Step 2, when the stress ratio σ/σ c of the skin element is greater than or equal to 1, set a thickness increase coefficient for increasing the thickness of the skin element;

步骤3,根据采用步骤2中设置的厚度增大系数对蒙皮单元形成的厚度增大量与蒙皮厚度加工工艺的最小量的对比,确定采用厚度增大系数或蒙皮厚度加工工艺的最小量对蒙皮单元进行厚度增大处理。Step 3: According to the comparison between the thickness increase amount formed by the skin unit using the thickness increase coefficient set in step 2 and the minimum amount of the skin thickness processing technology, determine the thickness increase coefficient or the minimum amount of the skin thickness processing technology. Increase the thickness of the skin elements.

本发明实施例的步骤2中,设置厚度增大系数的实现方式为:In step 2 of the embodiment of the present invention, the implementation manner of setting the thickness increase coefficient is as follows:

设置厚度增大系数为(σ/σc)1/3;其中,对每个独立的蒙皮单元,厚度增大(σ/σc)1/3后,相应蒙皮单元的当量应力σ等于临界失稳应力σcSet the thickness increase factor to (σ/σ c ) 1/3 ; where, for each individual skin element, after the thickness is increased by (σ/σ c ) 1/3 , the equivalent stress σ of the corresponding skin element is equal to Critical buckling stress σ c .

本发明实施例的步骤3中获取蒙皮单元的厚度增大量的实施方式包括:The implementation manner of obtaining the thickness increase of the skin unit in step 3 of the embodiment of the present invention includes:

步骤31,计算采用厚度增大系数(σ/σc)1/3对蒙皮单元进行厚度增大形成的厚度增大量((σ/σc)1/3-1)t;其中,t为蒙皮单元的壁板厚度。Step 31: Calculate the thickness increase amount ((σ/σ c ) 1/3 -1)t formed by increasing the thickness of the skin element by using the thickness increase coefficient (σ/σ c ) 1/3 ; where t is The thickness of the siding of the skinned element.

本发明实施例的步骤3中确定对蒙皮单元进行厚度增大的实施方式包括:In the step 3 of the embodiment of the present invention, the implementation manner of determining to increase the thickness of the skin unit includes:

步骤32,判断厚度增大量是否大于或等于0.2mm,0.2mm为蒙皮厚度加工工艺的最小量;Step 32, determine whether the thickness increase is greater than or equal to 0.2mm, and 0.2mm is the minimum amount of the skin thickness processing technology;

步骤33,在厚度增大量大于或等于0.2mm时,采用厚度增大系数对蒙皮单元进行厚度增大;Step 33, when the thickness increase is greater than or equal to 0.2mm, use the thickness increase coefficient to increase the thickness of the skin unit;

步骤34,在厚度增大量小于0.2mm时,以0.2mm为厚度增大量对蒙皮单元进行厚度增大。Step 34 , when the thickness increase is less than 0.2 mm, the thickness of the skin unit is increased by taking 0.2 mm as the thickness increase.

本发明实施例的步骤33中采用厚度增大系数对蒙皮单元进行厚度增大的方式,包括:In step 33 of the embodiment of the present invention, the method of increasing the thickness of the skin unit by the thickness increase coefficient includes:

对厚度增大系数设置加速因子k,通过加速因子k将厚度增大系数 (σ/σc)1/3放大k倍,并通过k(σ/σc)1/3对蒙皮单元进行厚度增大;本发明实施例中通过设置加速因子k的方式,以减少蒙皮壁板参数优化的迭代次数。The acceleration factor k is set for the thickness increase factor, the thickness increase factor (σ/σ c ) 1/3 is amplified by k times through the acceleration factor k, and the thickness of the skin element is increased by k(σ/σ c ) 1/3 . Increase; in the embodiment of the present invention, the acceleration factor k is set to reduce the number of iterations for parameter optimization of the skin panel.

本发明实施例提供的蒙皮壁板参数优化方法还可以包括如下步骤:The skin panel parameter optimization method provided by the embodiment of the present invention may further include the following steps:

步骤4,通过重复执行步骤2到步骤3对蒙皮单元进行厚度增大,指定蒙皮单元的应力比值σ/σc小于1。Step 4, by repeating steps 2 to 3, the thickness of the skin element is increased, and the stress ratio σ/σ c of the specified skin element is less than 1.

需要说明的是,本发明实施例中的加速因子k的取值为1到2之间;另外,蒙皮壁板参数优化方法在每次执行步骤2到步骤3的过程中,通过设置加速因子k的取值,调整蒙皮单元的厚度增大量,从而调整蒙皮壁板参数优化的速度。It should be noted that the value of the acceleration factor k in the embodiment of the present invention is between 1 and 2; in addition, the parameter optimization method of the skin panel sets the acceleration factor during each execution of steps 2 to 3. The value of k adjusts the thickness increase of the skin unit, thereby adjusting the speed of the optimization of the parameters of the skin panel.

本发明实施例提供的蒙皮壁板参数优化方法还可以包括如下步骤:The skin panel parameter optimization method provided by the embodiment of the present invention may further include the following steps:

步骤2a,在蒙皮单元的应力比值σ/σc小于1时,相应的蒙皮单元符合稳定性控制,优化结束。Step 2a, when the stress ratio σ/ σc of the skin element is less than 1, the corresponding skin element conforms to the stability control, and the optimization ends.

与传统的金属机翼蒙皮壁板参数优化方法相比,本发明实施例提供的基于稳定性算法的蒙皮壁板参数优化方法,一方面,结合理论和工程经验提出蒙皮壁板参数的加速优化方法,创造性的将稳定性校核控制推导为稳定性主动设计,所推导出的增大系数为壁板参数优化稳定性控制指明了优化方向;另一方面,加速因子的设置极大的减少了优化迭代次数,将传统优化方法的手动调节几个小时的工作缩短至几分钟,极大地提高了优化效率。本发明实施例提供的基于稳定性算法的蒙皮壁板参数优化方法,将稳定性设计作为一种主动设计因素引入优化流程,设计加速迭代系数,实现参数快速优化。Compared with the traditional metal wing skin panel parameter optimization method, the skin panel parameter optimization method based on the stability algorithm provided by the embodiment of the present invention, on the one hand, combines theory and engineering experience to propose a skin panel parameter optimization method. The acceleration optimization method creatively deduces the stability check control as the stability active design, and the derived increase coefficient points out the optimization direction for the optimization stability control of the panel parameters; on the other hand, the setting of the acceleration factor greatly affects The number of optimization iterations is reduced, and the manual adjustment work of traditional optimization methods is shortened from several hours to several minutes, which greatly improves the optimization efficiency. The skin panel parameter optimization method based on the stability algorithm provided by the embodiment of the present invention introduces the stability design as an active design factor into the optimization process, designs the accelerated iteration coefficient, and realizes the rapid optimization of parameters.

以下对本发明实施例中厚度增大系数(σ/σc)1/3的推导方式进行说明。The derivation method of the thickness increase coefficient (σ/σ c ) 1/3 in the embodiment of the present invention will be described below.

采用本发明实施例提供的方法进行优化前,设当量应力为σ0,临界失稳应力σc0,且当量应力σ0满足以下关系式:Before using the method provided by the embodiment of the present invention for optimization, let the equivalent stress be σ 0 , the critical buckling stress σ c0 , and the equivalent stress σ 0 satisfies the following relationship:

Figure RE-GDA0003074017650000061
Figure RE-GDA0003074017650000061

其中,F为蒙皮单元压缩载荷,b为蒙皮单元加载边宽度;当量应力σ0和临界失稳应力σc0满足以下关系:Among them, F is the compressive load of the skin element, b is the width of the loaded edge of the skin element; the equivalent stress σ 0 and the critical buckling stress σ c0 satisfy the following relationship:

Figure RE-GDA0003074017650000062
Figure RE-GDA0003074017650000062

上述公式中,Kc为压缩临界应力系数,E为材料的弹性模型,μ为材料的弹性泊松比。In the above formula, K c is the compressive critical stress coefficient, E is the elastic model of the material, and μ is the elastic Poisson's ratio of the material.

假设蒙皮单元厚度增大X倍后,理论上蒙皮单元的当量应力等于临界失稳应力,设优化后的当量应力为σ1,临界失稳应力为σc1,根据公式可知:Assuming that the thickness of the skin element is increased by X times, the theoretical equivalent stress of the skin element is equal to the critical buckling stress, and the optimized equivalent stress is σ 1 , and the critical buckling stress is σ c1 , according to the formula:

Figure RE-GDA0003074017650000063
Figure RE-GDA0003074017650000063

σc1=σc0*X2σ c1c0 *X 2 ;

二者相等即可推导出:It can be deduced that the two are equal:

X=(σ0c0)1/3X=(σ 0c0 ) 1/3 .

以下通过一个实施示例对本发明实施例提供的基于稳定性算法的蒙皮壁板参数优化方法具体实施方式进行详细说明。The specific implementation of the method for optimizing the parameters of the skin panel based on the stability algorithm provided by the embodiment of the present invention will be described in detail below through an implementation example.

如图1所示,为本发明实施例提供的基于稳定性算法的蒙皮壁板参数优化方法,图1具体是从CAD模型到CAE模型的快速属性映射方法流程图,该实施示例中蒙皮壁板参数优化方法具体实现原理如下:As shown in FIG. 1, it is a method for optimizing the parameters of a skin panel based on a stability algorithm provided by an embodiment of the present invention. FIG. 1 is a flowchart of a method for fast attribute mapping from a CAD model to a CAE model. In this implementation example, the skin The specific realization principle of the panel parameter optimization method is as follows:

1,对金属机翼蒙皮壁板结构的初始模型,在参数优化过程中,选取蒙皮厚度t作为优化变量,计算蒙皮单元格的当量应力σ,按满应力和稳定性约束,本发明重点为稳定性控制及稳定性反馈设计,计算该蒙皮单元格的临界失稳应力σc,判断比值σ/σc1. For the initial model of the metal wing skin panel structure, in the parameter optimization process, the skin thickness t is selected as the optimization variable, and the equivalent stress σ of the skin cell is calculated. According to the full stress and stability constraints, the present invention The emphasis is on stability control and stability feedback design, calculating the critical buckling stress σ c of the skin cell, and judging the ratio σ/σ c ;

2,若上述比值σ/σc小于1,说明该蒙皮单元符合稳定性控制,满足设计要求,优化结束;2. If the above ratio σ/σ c is less than 1, it means that the skin unit meets the stability control, meets the design requirements, and the optimization ends;

3,若上述比值σ/σc大于或等1,说明该蒙皮单元将失稳,需加厚该蒙皮单元,设计厚度增大系数(σ/σc)1/33. If the above ratio σ/σ c is greater than or equal to 1, it means that the skin element will be unstable, the skin element needs to be thickened, and the design thickness increase factor (σ/σ c ) 1/3 ;

需要说明的是:该增大系数的设计原理为,当量应力σ∝t-1,临界失稳应力σc∝t2(已知公式),理论上对每个单独的蒙皮单元,厚度增大(σ/σc)1/3后,该蒙皮单元的当量应力σ等于临界失稳应力σcIt should be noted that the design principle of the increase factor is that the equivalent stress σ∝t -1 , the critical instability stress σ c ∝t 2 (known formula), theoretically, for each individual skin element, the thickness increases After increasing (σ/σ c ) 1/3 , the equivalent stress σ of the skin element is equal to the critical buckling stress σ c ;

4,考虑制造工艺,判断厚度增大量((σ/σc)1/3-1)t是否大于或等于 0.2mm,若大于或等于0.2mm,则可按步骤3执行增大,若小于0.2mm时则按0.2mm增大;该步骤的制造工艺中蒙皮厚度加工工艺的最小量;4. Considering the manufacturing process, determine whether the thickness increase ((σ/σ c ) 1/3 -1)t is greater than or equal to 0.2mm, if it is greater than or equal to 0.2mm, it can be increased according to step 3, if it is less than 0.2 When it is mm, it is increased by 0.2mm; the minimum amount of skin thickness processing technology in the manufacturing process of this step;

5,实际优化过程中,为加快蒙皮壁板参数优化过程,可设计加速因子k,将(σ/σc)1/3放大k倍,通过k(σ/σc)1/3对蒙皮单元进行厚度增大,从而大大减少迭代次数。5. In the actual optimization process, in order to speed up the optimization process of the parameters of the skin panel, the acceleration factor k can be designed, and the (σ / σ c ) 1/3 is enlarged by k times, and the The thickness of the pico elements is increased, which greatly reduces the number of iterations.

以下对本发明实施例中使用加速因子进行优化的实施方式进行说明,如图 2所示,为本发明实施例提供的另一种基于稳定性算法的蒙皮壁板参数优化方法的流程图。该流程具体说明加速因子k的设置对流程的影响。The following describes the implementation of the optimization using the acceleration factor in the embodiment of the present invention. As shown in FIG. 2, it is a flowchart of another method for optimizing the parameters of the skin panel based on the stability algorithm provided by the embodiment of the present invention. This process specifically illustrates the influence of the setting of the acceleration factor k on the process.

设定初始优化阶段,k为1;Set the initial optimization stage, k is 1;

迭代次数0次:计算初始模块不满足稳定性控制的蒙皮单元格数目m;The number of iterations is 0: calculate the number of skin cells m that the initial module does not meet the stability control;

迭代次数1次:计算首次优化后不满足稳定性控制的蒙皮单元格数目n;The number of iterations is 1: calculate the number n of skin cells that do not meet the stability control after the first optimization;

若n/m大于m-0.125时,认为经过8次迭代后无法完成优化,则设置加速因子k,视优化精度要求在区间(1,2)取值,取值越大,优化速度越快,优化精度越低。If n/m is greater than m -0.125 , it is considered that the optimization cannot be completed after 8 iterations, then the acceleration factor k is set, and the value is in the interval (1, 2) depending on the optimization accuracy requirements. The larger the value, the faster the optimization speed. The lower the optimization precision.

虽然本发明所揭露的实施方式如上,但所述的内容仅为便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式及细节上进行任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。Although the embodiments disclosed in the present invention are as above, the described contents are only the embodiments adopted to facilitate the understanding of the present invention, and are not intended to limit the present invention. Any person skilled in the art to which the present invention belongs, without departing from the spirit and scope disclosed by the present invention, can make any modifications and changes in the form and details of the implementation, but the scope of the patent protection of the present invention still needs to be The scope defined by the appended claims shall prevail.

Claims (3)

1. A skin panel parameter optimization method based on a stability algorithm is characterized by comprising the following steps:
step 1, calculating the equivalent stress sigma and the critical instability stress sigma of a skin unit according to an initial model of a skin panel structure of a metal wing c And calculating the stress ratio sigma/sigma of the skin unit c
Step 2, stress ratio sigma/sigma of the skin unit c When the thickness of the skin unit is greater than or equal to 1, setting a thickness increasing coefficient for increasing the thickness of the skin unit;
step 3, determining that the thickness increasing coefficient or the minimum amount of the skin thickness processing technology is used for performing thickness increasing treatment on the skin unit according to the comparison between the thickness increasing amount formed by the skin unit and the minimum amount of the skin thickness processing technology by using the thickness increasing coefficient set in the step 2;
step 4, the thickness of the skin units is increased by repeatedly executing the steps 2 to 3, and the stress ratio sigma/sigma of the skin units is specified c Less than 1;
wherein, the thickness increasing coefficient is set in the step 2 in a manner that:
setting the thickness increase coefficient to (sigma/sigma) c ) 1/3 (ii) a Wherein the thickness increases (σ/σ) for each individual skin cell c ) 1/3 Then, the equivalent stress σ of the corresponding skin unit is equal to the critical buckling stress σ c
The manner of obtaining the thickness increase amount of the skin unit in the step 3 is as follows:
step 31, calculating the thickness increase coefficient (sigma/sigma) c ) 1/3 An increased thickness amount ((sigma/sigma) formed by increasing the thickness of the skin unit c ) 1/3 -1) t; wherein t is the thickness of the wall plate of the skin unit;
determining a thickness increasing mode of the skin unit in the step 3 includes:
step 32, judging whether the thickness increase is greater than or equal to 0.2mm, wherein 0.2mm is the minimum amount of the skin thickness processing technology;
step 33, when the thickness increase is greater than or equal to 0.2mm, adopting the thickness increase coefficient to increase the thickness of the skin unit;
step 34, when the thickness increment is smaller than 0.2mm, the thickness increment is carried out on the skin unit by taking 0.2mm as the thickness increment;
in the step 33, the increasing the thickness of the skin unit by using the thickness increase coefficient includes:
setting an acceleration factor k to the thickness increasing coefficient by which the thickness increasing coefficient (σ/σ) is increased c ) 1/3 Amplifying by k times and passing k (sigma/sigma) c ) 1/3 And increasing the thickness of the skin unit.
2. The stability algorithm-based skin panel parameter optimization method of claim 1, wherein the acceleration factor k has a value between 1 and 2;
in the process of executing the steps 2 to 3 each time, the skin panel parameter optimization method adjusts the thickness increase amount of the skin unit by setting the value of the acceleration factor k, so that the skin panel parameter optimization speed is adjusted.
3. The stability algorithm-based skin panel parameter optimization method of claim 1, further comprising:
step 2a, stress ratio sigma/sigma of the skin unit c When less than 1, the corresponding maskAnd (5) the skin unit accords with stability control, and optimization is finished.
CN202011612917.9A 2020-12-29 2020-12-29 Skin panel parameter optimization method based on stability algorithm Active CN113032892B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011612917.9A CN113032892B (en) 2020-12-29 2020-12-29 Skin panel parameter optimization method based on stability algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011612917.9A CN113032892B (en) 2020-12-29 2020-12-29 Skin panel parameter optimization method based on stability algorithm

Publications (2)

Publication Number Publication Date
CN113032892A CN113032892A (en) 2021-06-25
CN113032892B true CN113032892B (en) 2022-10-11

Family

ID=76460527

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011612917.9A Active CN113032892B (en) 2020-12-29 2020-12-29 Skin panel parameter optimization method based on stability algorithm

Country Status (1)

Country Link
CN (1) CN113032892B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105184390A (en) * 2015-08-12 2015-12-23 中国运载火箭技术研究院 Integrated optimization method of static strength, rigidity, stability of wallboard structure

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105184390A (en) * 2015-08-12 2015-12-23 中国运载火箭技术研究院 Integrated optimization method of static strength, rigidity, stability of wallboard structure

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于PATRAN/NASTRAN的复合材料机翼蒙皮优化设计;常楠等;《西北工业大学学报》;20060630(第03期);全文 *
复合材料盒段加筋壁板结构屈曲规律研究;谭景磊等;《飞机设计》;20130815(第04期);全文 *
考虑复合材料蒙皮稳定性的飞机翼面结构布局优化设计;刘颖卓等;《航空学报》;20101025(第10期);全文 *

Also Published As

Publication number Publication date
CN113032892A (en) 2021-06-25

Similar Documents

Publication Publication Date Title
CN107391891B (en) Large-aspect-ratio wing optimization design method based on model fusion method
CN107885945A (en) Towards the multiple dimensioned pore space structure lightweight modeling method of 3D printing
CN105183996A (en) Surface element correction and grid beforehand self-adaption calculation method
CN104933250A (en) Wing aerodynamics modeling method
CN105956342B (en) A kind of composite material predeformation cabin door structure optimum design method from locking
CN103488847B (en) A kind of aerodynamic shape optimization method based on Artificial neural network ensemble
CN116628854B (en) Wing section aerodynamic characteristic prediction method, system, electronic equipment and storage medium
CN104573281A (en) Complex space curved surface thin plate forming die face designing method taking springback compensation
CN112528540A (en) Structure optimization method based on force-heat coupling calculation
CN112906331A (en) Standard unit delay model construction method based on logarithmic expansion skewed state distribution
CN113139313B (en) Equivalent lightweight structure analysis method based on finite element theory
CN106372348B (en) Vector fitting model order reducing method based on control errors in a kind of linear system
CN106202693A (en) A kind of Material Stiffened Panel structure anti-vibration fatigue optimization method based on parametric modeling
CN113032892B (en) Skin panel parameter optimization method based on stability algorithm
CN102411652B (en) Method for determining sectional dimension of rectangular beam with lugs
CN111046494B (en) Simplified vehicle body floor design method based on multi-component structural form
CN111539138A (en) Method for solving time domain response sensitivity of structural dynamics peak value based on step function
CN119476016A (en) Mechanical structure reliability analysis method and system based on active learning agent model
CN106971023B (en) Design method of special-shaped hole structure of roulette based on hyperelliptic curve
CN116244894B (en) A power system transient simulation method and system based on large step size
CN106294908B (en) Sound serves as a contrast design method
CN111597724A (en) Structural dynamics topology optimization method and system considering frequency band constraint and application
CN115374545A (en) A Prediction Method and System for Rough Surface Contact Stiffness Considering Elastoplastic State
CN117113589B (en) Multi-objective intelligent optimization method for pre-connection process of weak-rigidity thin-wall structure
CN112580223A (en) Bimetal interface self-adaptive loading simulation method based on selective area stress criterion

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