CN112711804B - Method for analyzing wave isolation intensity of wall plate of high-lift device - Google Patents
Method for analyzing wave isolation intensity of wall plate of high-lift device Download PDFInfo
- Publication number
- CN112711804B CN112711804B CN202011611555.1A CN202011611555A CN112711804B CN 112711804 B CN112711804 B CN 112711804B CN 202011611555 A CN202011611555 A CN 202011611555A CN 112711804 B CN112711804 B CN 112711804B
- Authority
- CN
- China
- Prior art keywords
- wave
- wall plate
- lift device
- buckling
- load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/28—Fuselage, exterior or interior
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
The invention relates to the technical field of structural strength analysis, in particular to a method for analyzing the wave isolation strength of a wall plate of a high-lift device. The method comprises the following steps: establishing a static finite element model for the wall plate of the high-lift device to obtain a static stress analysis result of the wall plate of the high-lift device, and determining the ratio of the spanwise working load to the chordwise working load according to the stress analysis result; carrying out engineering stress analysis on the known hyperplasia device structure of the high-lift device wall plate, and determining the critical buckling load of the high-lift device wall plate according to the engineering stress analysis result and the ratio; according to the ratio, carrying out structural buckling mode analysis to obtain a structural buckling mode of the high lift device wall plate; determining initial arrangement parameters of the wave isolation member arranged on the wall plate of the high lift device according to the structural buckling mode and the critical buckling load; and (4) carrying out position optimization, configuration optimization and layering parameter optimization on the wave isolation component by adopting a finite element secondary modal analysis method according to the initial arrangement parameters.
Description
Technical Field
The invention relates to the technical field of structural strength analysis, in particular to a method for analyzing the wave isolation strength of a wall plate of a high-lift device.
Background
The wallboard has large structure size, complex stress form, high stress level and large occupied structure weight proportion. The large structure size enables local instability of the structure to occur under low working load, the existing wall plate wave isolation calculation method is mainly determined by means of engineering experience, and the stability design requirement of the existing structure cannot be met.
Disclosure of Invention
The purpose of the invention is as follows: the position optimization, the configuration and the size optimization of the wave-isolating component under the instability condition are carried out on the wall plate under the conditions of large size, large load and complex boundary, and finally the instability control design method of the wall plate under the complex stress state is obtained, so that the structural weight is effectively reduced.
The technical scheme is as follows:
a method for analyzing the wave isolation intensity of a wall plate of a high-lift device comprises the following steps:
establishing a static finite element model for the wall plate of the high-lift device to obtain a static stress analysis result of the wall plate of the high-lift device, and determining the ratio of the spanwise working load to the chordwise working load according to the stress analysis result;
carrying out engineering stress analysis on the known hyperplasia device structure of the high-lift device wall plate, and determining the critical buckling load of the high-lift device wall plate according to the engineering stress analysis result and the ratio;
according to the ratio, carrying out structural buckling mode analysis to obtain a structural buckling mode of the high lift device wall plate;
determining initial arrangement parameters of the wave isolation members arranged on the wall plate of the high lift device according to the structural buckling mode and the critical buckling load;
and (4) carrying out position optimization, configuration optimization and layering parameter optimization on the wave isolation component by adopting a finite element secondary modal analysis method according to the primary arrangement parameters.
And (3) carrying out engineering stress analysis on the known hyperplasia device structure of the high-lift device wall plate, and determining the critical buckling load of the high-lift device wall plate according to the engineering stress analysis result and the ratio, wherein the method comprises the following steps:
because the wall plate of the high lift device bears the biaxial pressure load of the spanwise bending load and the chordwise bending load, the wall plate engineering of the high lift device is simplified into a biaxial pressure bearing rectangular flat plate under the condition of simple boundary;
calculating a bending stiffness coefficient matrix of the composite material according to the layering parameters of the composite material;
under the condition that the ratio is constant, calculating the double-shaft pressure stress rectangular flat plate shaft pressure buckling load, namely the critical buckling load;
wherein D is 11 、D 22 、D 66 The bending rigidity coefficient, m, n, plate length and width direction of the laminated plateA, b refer to the length and width of the plate, nx, ny refer to the buckling load in two directions.
The wave-insulating member must satisfy the buckling condition:
(a) The sectional area of the wave-isolation component is larger than the area of the wall plate by 40 percent;
A min =0.4×A;
A、A min the sectional area of the wave-isolating member and the sectional area of the wall plate member.
(b) The wave isolation component meets the requirement of minimum moment of inertia and is later than the instability of the wallboard structure;
I u 、I min the moment of inertia and the minimum moment of inertia of the wave isolation component are pointed; d. t, h e 、k s The distance between the struts, the thickness of the web plate, the distance between the strip centers of the upper and lower edges and the shear buckling coefficient;
(c) The instability stress of the wave isolation component is as follows:
D 11 、D 66 refers to the flexural rigidity coefficient, b, L, delta, sigma, of the laminate W The width, length and thickness of the wave-isolating component and the working stress of the wave-isolating component are indicated.
(d) The working stress of the wave isolation component is smaller than the allowable strain;
(e) The wave-isolating component meets the requirement of the connection strength of the wall plate.
Adopting a finite element quadratic modal analysis method according to the primary arrangement parameters to carry out position optimization, configuration optimization and layering parameter optimization on the wave isolation component, and the method comprises the following steps:
carrying out integral static stress analysis on the high-lift device wall plate which is preliminarily provided with the wave isolation member;
taking out the parts which are possibly damaged by instability from the whole static stress analysis model, and establishing a local analysis model;
simulating a forced displacement elastic support boundary of a local 'analysis model' according to the analysis result of the overall static stress;
and carrying out buckling modal analysis on the local analysis model, and optimizing the position, the configuration form and the structural parameters of the wave isolation member according to the buckling analysis result.
Optimizing the position, the configuration form and the structural parameters of the wave-isolating component according to the buckling analysis result, wherein the method comprises the following steps:
defining the configuration and the preliminary arrangement parameters of the wave-isolating component, and carrying out position optimization on the wave-isolating component on the local 'analysis model', namely selecting the position of the wave-isolating component with the maximum critical buckling load in different positions as a better position;
limiting the wave isolation member to be at a better position and preliminarily arranging parameters, and optimizing T-shaped and I-shaped configuration forms to obtain a better configuration form;
and limiting the wave isolation member to be at a better position, and carrying out the parameter optimization of the wave isolation member in a better configuration form to obtain the paving parameters meeting the instability condition.
Composite material bending rigidity coefficient matrix D ij ;
In the formula, Q ij 、θ、z k 、z k-1 The positive axis modulus of the main direction of the ply material, the off-axis modulus of a certain ply, the angle of the ply of the composite material, the z coordinate of the K-th layer and the z coordinate of the K-1-th layer are shown.
The preliminary arrangement parameters comprise a preliminary layering parameter, a height parameter and a thickness parameter.
A computer readable storage medium having stored thereon computer instructions which, when executed by a processor, implement the method of any of the above.
Has the advantages that: carrying out wall plate wave isolation instability control on a wall plate structure under the conditions of large size, large load and complex boundary; the position, the configuration and the size of the wave isolation component are optimized, and the weight of the wallboard structure is effectively reduced.
Drawings
FIG. 1 is a schematic view of a large size, high load, complex boundary wall panel construction;
FIG. 2 is a diagram of a rectangular flat plate under biaxial compression;
FIG. 3 is a static finite element model of a wall panel structure;
FIG. 4 is a partial buckling "analytical model" finite element diagram;
FIG. 5 is a comparison graph of the optimization results of the positions of the wave-isolating members;
FIG. 6 is a comparison graph of the optimization results of the wave-insulating member configuration.
Detailed Description
A proliferation device structure is known, which comprises a main bearing part composite material laminated wall plate 1, a wall plate longitudinal supporting part 2 and a wall plate transverse supporting part 3, wherein the supporting parts are combined to form a closed structure form, as shown in figure 1.
(1) Carrying out engineering stress analysis on the structure shown in the figure 1, wherein the composite material laminated wall plate 1 bears biaxial pressure load under the action of spanwise bending and chordwise bending loads, and the engineering is simplified into a biaxial pressure stressed rectangular flat plate under the condition of simple boundary, and the stress form is shown in figure 2;
(2) The bending mode of the laminate is mainly related to the bending stiffness matrix of the composite material layer, and the bending stiffness coefficient Dij of the composite material is calculated;
wherein, qij,θ、z k 、z k-1 The positive axis modulus of the main direction of the ply material, the off-axis modulus of a certain ply, the angle of the ply of the composite material, the z coordinate of the K-th layer and the z coordinate of the K-1 th layer are shown.
(3) Under the condition that the Nx/Ny ratio is certain, calculating the axial compression buckling load of the double-axial compression stressed rectangular flat plate
Wherein D is 11 、D 22 、D 66 The bending stiffness coefficient of the laminated plate is shown, m, n, half wave number in the length direction and the width direction of the plate, a and b are shown as the length and the width of the plate, and Nx and Ny are buckling loads in two directions.
(4) Establishing a static finite element model for the composite laminated wallboard (namely the high-lift device wallboard) to obtain a stress analysis result of the wallboard 1, and obtaining a structural buckling mode according to ratios of spanwise working load/spanwise buckling load and chordwise working load/chordwise buckling load;
(5) Arranging wave isolation members according to the structural buckling mode 4 and the buckling load size 3; the wave-insulating member must satisfy the buckling condition:
a) The sectional area of the wave-insulating component is 40% larger than that of the wall plate (between the two wave-insulating components);
A min =0.4×A
A、A min the sectional area of the wave-isolating member and the sectional area of the wall plate member.
b) The wave isolation component meets the requirement of minimum moment of inertia and is later than the instability of the wallboard structure;
I u 、I min the moment of inertia and the minimum moment of inertia of the wave isolation component are pointed; d. t, h e 、k s The distance between the struts, the thickness of the web plate, the distance between the strip centers of the upper and lower edges and the shear buckling coefficient;
the instability stress of the wave isolation component is as follows:
D 11 、D 66 refers to the flexural rigidity coefficient, b, L, delta, sigma, of the laminate W The width, the length, the thickness and the working stress of the wave-isolating component are indicated.
c) The working stress of the wave isolation component is smaller than the allowable strain;
d) The wave isolation component meets the requirement of the connection strength of the wall plate;
(6) And carrying out strong-adaptability finite element secondary analysis on the simplified engineering analysis wave isolation component to obtain better wave isolation component position, configuration form and structure layering parameters.
a) The overall stress analysis of the overall structure is carried out, see fig. 3;
b) Taking out the parts which are possible to be destabilized and damaged, and establishing a local analysis model (see figure 4);
c) Simulating the forced displacement elastic support boundary of the local 'analysis model' according to the overall stress result, and showing in figure 4;
d) And carrying out local buckling analysis on the local 'analysis model', and optimizing the position, the configuration form and the structural parameters of the wave isolation component according to buckling analysis results.
Limiting the form of the wave-isolating component, optimizing the position, and selecting different positions by using the finite element model shown in FIG. 4 to obtain a preferred position which is the wave-isolating component position 1 for improving the critical buckling load shown in FIG. 5;
defining the position of the wave isolation member, and carrying out configuration form optimization: typical T-type, i-type wave-blocking building block results are shown in figure 6.
And limiting the position and the configuration form of the wave isolation member, and optimizing the parameters of the I-shaped wave isolation member to obtain the paving parameters meeting the instability condition.
A flap hyperplasia device of a certain type of airplane adopts a typical dense-rib composite laminated board structure, and the stability of the lower board is a prominent problem in large size, large load and complex boundary.
According to the engineering analysis result, the structure generates instability at about 40% of load, according to the engineering analysis result, the wave isolation member is arranged as a condition for controlling the instability of the wallboard, and through the position optimization, the structure configuration optimization and the structure parameter optimization of the wave isolation member, the weight of the wallboard structure is finally reduced by 19.5%, and the total weight of the structure is reduced by 12% as shown in table 1.
TABLE 1 results of wave-damping destabilizing element placement
Claims (5)
1. A method for analyzing the wave isolation strength of a wall plate of a high-lift device is characterized by comprising the following steps:
establishing a static finite element model for the wall plate of the high-lift device to obtain a static stress analysis result of the wall plate of the high-lift device, and determining the ratio of the spanwise working load to the chordwise working load according to the stress analysis result;
carrying out engineering stress analysis on the known hyperplasia device structure of the high-lift device wall plate, and determining the critical buckling load of the high-lift device wall plate according to the engineering stress analysis result and the ratio;
according to the ratio, carrying out structural buckling mode analysis to obtain a structural buckling mode of the high lift device wall plate;
determining initial arrangement parameters of the wave isolation member arranged on the wall plate of the high lift device according to the structural buckling mode and the critical buckling load;
performing position optimization, configuration optimization and layering parameter optimization on the wave isolation component by adopting a finite element secondary modal analysis method according to the primary arrangement parameters;
carrying out engineering stress analysis on the known hyperplasia device structure of the high-lift device wallboard, and determining the critical buckling load of the high-lift device wallboard according to the engineering stress analysis result and the ratio, wherein the method comprises the following steps:
because the wall plate of the high lift device bears the biaxial pressure load of the spanwise bending load and the chordwise bending load, the wall plate engineering of the high lift device is simplified into a biaxial pressure bearing rectangular flat plate under the condition of simple boundary;
calculating a bending stiffness coefficient matrix of the composite material according to the layering parameters of the composite material;
under the condition that the ratio is constant, calculating the double-shaft pressure stress rectangular flat plate shaft pressure buckling load, namely the critical buckling load;
wherein D is 11 、D 22 、D 66 The bending rigidity coefficient of the laminated plate is defined, m, n, the length of the plate, the half wave number in the width direction, a, b are the length and the width of the plate, and Nx, ny are the buckling loads in two directions;
the wave-insulating member must satisfy the buckling condition:
(a) The sectional area of the wave-isolation component is larger than the area of the wall plate by 40 percent;
A min =0.4×A;
A、A min is a sectional area of the wave-isolating member, a wall structureA cross-sectional area of the part;
(b) The wave isolation component meets the requirement of minimum moment of inertia and is later than the instability of the wallboard structure;
I u 、I min the moment of inertia and the minimum moment of inertia of the wave isolation component are pointed; d. t, h e 、k s The distance between the struts, the thickness of the web plate, the distance between the strip centers of the upper and lower edges and the shear buckling coefficient;
(c) The instability stress of the wave isolation component is as follows:
D 11 、D 66 refers to the flexural rigidity coefficient, b, L, delta, sigma, of the laminate W The width, the length, the thickness and the working stress of the wave-isolating component are measured;
(d) The working stress of the wave isolation component is smaller than the allowable strain;
(e) The wave isolation component meets the requirement of the connection strength of the wall plate;
composite material bending rigidity coefficient matrix D ij ;
2. The method of claim 1, wherein performing position optimization, configuration optimization and ply parameter optimization of the wave-separating member using finite element quadratic modal analysis based on the preliminary layout parameters comprises:
carrying out integral static stress analysis on the high-lift device wall plate which is preliminarily provided with the wave isolation member;
taking out the parts which are possibly damaged by instability from the whole static stress analysis model, and establishing a local analysis model;
simulating a forced displacement elastic support boundary of a local 'analysis model' according to the analysis result of the overall static stress;
and carrying out buckling modal analysis on the local analysis model, and optimizing the position, the configuration form and the structural parameters of the wave isolation member according to the buckling analysis result.
3. The method of claim 2, wherein the wave-isolating member position, configuration form and structural parameters are optimized according to the buckling analysis result, and the method comprises the following steps:
the configuration and the preliminary arrangement parameters of the wave-isolating component are limited, the position of the wave-isolating component on the local 'analysis model' is optimized, namely the position of the wave-isolating component with the maximum critical buckling load in different positions is selected as a better position;
limiting the wave isolation member to be at a better position and preliminarily arranging parameters, and optimizing T-shaped and I-shaped configuration forms to obtain a better configuration form;
and limiting the wave isolation member to be at a better position, and carrying out the parameter optimization of the wave isolation member in a better configuration form to obtain the paving parameters meeting the instability condition.
4. The method of claim 1, wherein the preliminary placement parameters include a preliminary layup parameter, a height parameter, and a thickness parameter.
5. A computer-readable storage medium having computer instructions stored thereon, wherein the instructions, when executed by a processor, implement the method of any of claims 1-4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011611555.1A CN112711804B (en) | 2020-12-29 | 2020-12-29 | Method for analyzing wave isolation intensity of wall plate of high-lift device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011611555.1A CN112711804B (en) | 2020-12-29 | 2020-12-29 | Method for analyzing wave isolation intensity of wall plate of high-lift device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112711804A CN112711804A (en) | 2021-04-27 |
CN112711804B true CN112711804B (en) | 2022-11-22 |
Family
ID=75547300
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011611555.1A Active CN112711804B (en) | 2020-12-29 | 2020-12-29 | Method for analyzing wave isolation intensity of wall plate of high-lift device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112711804B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113836638A (en) * | 2021-09-18 | 2021-12-24 | 中航西安飞机工业集团股份有限公司 | Layout method for determining flat plate structure under pure bending load |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016141725A1 (en) * | 2015-03-12 | 2016-09-15 | 中车青岛四方机车车辆股份有限公司 | Local structure instability finite element simulation apparatus for high-speed train body and method therefor |
-
2020
- 2020-12-29 CN CN202011611555.1A patent/CN112711804B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016141725A1 (en) * | 2015-03-12 | 2016-09-15 | 中车青岛四方机车车辆股份有限公司 | Local structure instability finite element simulation apparatus for high-speed train body and method therefor |
Non-Patent Citations (3)
Title |
---|
剪切载荷作用下复合材料加筋壁板蒙皮屈曲;李爱环等;《航空计算技术》;20170125(第01期);全文 * |
斜削型筋条复合材料加筋壁板轴压屈曲分析;赵铭等;《科学技术与工程》;20170228(第06期);全文 * |
铝合金加筋板轴压屈曲稳定性的有限元分析;叶广宁等;《机械工程材料》;20130320(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112711804A (en) | 2021-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108153981B (en) | Finite element analysis-based fuselage stiffened wall panel post-buckling analysis method | |
CN108009323A (en) | A kind of composite structure uncertain optimization method based on multi-state failure load | |
JP6446455B2 (en) | Bonded and adjustable composite assembly | |
CN106844874A (en) | A kind of all-aluminium body light-weight design method based on sensitivity Yu CAE analysis | |
CN105184390A (en) | Integrated optimization method of static strength, rigidity, stability of wallboard structure | |
CN112711804B (en) | Method for analyzing wave isolation intensity of wall plate of high-lift device | |
CN203306222U (en) | Wing box formed by composite materials and wing for plane | |
CN105447269A (en) | Non-structural mass trimming method for airplane in steady overload state | |
Nali et al. | Assessments of refined theories for buckling analysis of laminated plates | |
US20140054420A1 (en) | Passive load alleviation for a fiber reinforced wing box of an aircraft with a stiffened shell structure | |
US9026406B2 (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 | |
Moors et al. | Weight trades in the design of a composite wing box: effect of various design choices | |
Zhang et al. | Optimization and analysis of composite sandwich box beam for solar drones | |
CN117352108A (en) | Composite material laminated plate design method based on D-D layering | |
CN108563916B (en) | Initial size optimization design method for thin-wall structure of aircraft wing fuselage | |
Liu et al. | Optimization of composite stiffened panels subject to compression and lateral pressure using a bi-level approach | |
CN110334364B (en) | Composite structural strength checking system | |
Hwu et al. | Aeroelastic divergence of stiffened composite multicell wing structures | |
CN111824356B (en) | Design method of I-shaped sandwich structure under action of wheel print load | |
CN112784369B (en) | Method for checking bearing capacity of composite material reinforced wall plate | |
Liu et al. | Optimum Buckling Design for Composite Wing Cover Panels with Manufacturing Constraints | |
CN113868761A (en) | Optimization design method for composite material airfoil surface skin | |
Shrivastava et al. | Design and optimization of a composite canard control surface of an advanced fighter aircraft under static loading | |
CN116305532A (en) | Variable gradient composite material wallboard strength analysis method | |
CN112560182A (en) | Civil aircraft composite material reinforced wallboard stringer rigidity determination method based on stability characteristics |
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 |