CN112182933A - Modeling and optimal design method of variable-rigidity composite material cylindrical shell - Google Patents
Modeling and optimal design method of variable-rigidity composite material cylindrical shell Download PDFInfo
- Publication number
- CN112182933A CN112182933A CN202011038065.7A CN202011038065A CN112182933A CN 112182933 A CN112182933 A CN 112182933A CN 202011038065 A CN202011038065 A CN 202011038065A CN 112182933 A CN112182933 A CN 112182933A
- Authority
- CN
- China
- Prior art keywords
- model
- cylindrical shell
- variable
- analysis
- composite material
- 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.)
- Pending
Links
Images
Classifications
-
- 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/26—Composites
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
Abstract
The invention discloses a modeling and optimal design method of a variable-rigidity composite material cylindrical shell; carrying out unit discretization design on a geometric model of the cylindrical shell; imparting variable stiffness composite properties to discrete units of a geometric model of the cylindrical shell; carrying out variable angle parameter assignment on discrete units in the model and establishing a finite element analysis model; and carrying out a small amount of finite element analysis on the parameterized analysis model of the discrete units in the model, extracting sample points, and establishing an approximate model of buckling analysis and optimization design. The modeling and optimization design method realizes the establishment of a simulation analysis model of the variable-stiffness composite material cylindrical shell through a discrete finite element method and parametric programming, introduces a neural network approximation model in the design optimization process, and performs quick response analysis and optimization design on the mechanical property of the cylindrical shell by taking the fiber placement starting angle and the layering quantity sequence as variables, thereby greatly improving the optimization design efficiency in the application of the underwater pressure-resistant cylindrical shell.
Description
Technical Field
The invention relates to the technical field of composite material revolution body modeling, in particular to a method for modeling and quickly optimizing design of a variable-rigidity composite material cylindrical shell.
Background
As a core component of an underwater vehicle, the cylindrical pressure hull not only needs to have greater strength and greater stability to ensure the safety of its internal instrumentation, but also must have a greater buoyancy/gravity ratio to provide a buoyancy reserve for the underwater vehicle. Compared with the traditional isotropic pressure-resistant shell structure such as high-strength steel, aluminum alloy, titanium alloy and the like, the composite pressure-resistant shell has the characteristics of light weight, high specific stiffness, high specific strength, good corrosion resistance, designability of structure and good low-noise concealment, can effectively improve the load capacity of the underwater vehicle, and has a great application prospect in the aspect of the pressure-resistant structure of the underwater vehicle.
The optimal design of the composite material structure is the key to reducing the design period and the production cost, and the emphasis is on how to fully utilize the directionality of the performance of the composite material and the designability of the structure. The variable stiffness is that the fiber composite material is laid in a curve manner in the laying process, the ply angle and the stiffness are changed along with the change of the fiber angle, the performance of the composite material can be exerted to the maximum extent by optimally designing the path and the ply parameters of the fiber composite material, and compared with a constant stiffness composite material structure, the structural performance of the variable stiffness composite material can be improved by 10-30%.
At present, the optimized design of the pressure shell made of the variable-rigidity composite material is mainly carried out by carrying out preliminary design through finite element analysis and carrying out verification through test. In the aspect of simulation modeling analysis of the pressure shell made of the variable-stiffness composite material, the uncertainties of continuously-changed fiber laying angles and various path laying modes cause that a variable-stiffness structure is difficult to directly depend on the existing commercial software for modeling; on the other hand, the optimal design process of the variable-stiffness composite material structure usually seeks an optimal solution based on repeated finite element calculation and iteration, so the design analysis period is long and the efficiency is not high. Therefore, the finite element analysis model of the variable-rigidity composite material cylindrical shell is accurately and reasonably established, the optimized design parameters of the finite element analysis model are quickly obtained, and the method has important significance for the analysis and design of the variable-rigidity underwater pressure shell.
Disclosure of Invention
In order to avoid the defects in the prior art, the invention provides a modeling and optimal design method of a variable-rigidity composite material cylindrical shell.
The idea of the invention is that the establishment of a simulation analysis model of the variable-stiffness composite material cylindrical shell is realized by a discrete finite element method and parametric programming, a neural network approximation model is introduced in the design optimization process, and the rapid response analysis and optimization design of the mechanical property of the cylindrical shell are carried out by taking the fiber laying starting angle and the layering quantity sequence as variables, so that the problem of overlong analysis period of the variable-stiffness composite material in the application of the underwater pressure-resistant cylindrical shell is solved, and the optimization design efficiency is greatly improved; the method has wide popularization significance in the aspects of improving the structural design and optimizing the efficiency of the variable-rigidity composite material cylindrical shell.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a modeling and optimal design method of a variable-rigidity composite material cylindrical shell is characterized by comprising the following steps:
step 1, carrying out unit discretization design on a geometric model of a cylindrical shell;
discretizing the geometric model according to a mathematical model of linear change of the fiber path of the composite material, wherein the mathematical model of the fiber angle change is expressed as follows:
in the formula, theta is the fiber angle, T0Is the angle of the fiber at the midpoint of the cylinder, T1Is the fiber angle at the cylindrical end point, x is the axial coordinate, and L is the cylindrical shell length; according to the expression, three of the cylindrical shell are arrangedThe dimensional geometric model is evenly divided into a plurality of circular ring structures along the axial coordinate, and each circular ring is further refined into discrete units according to the grid size of the designed finite element analysis;
step 2, endowing the discrete units with variable-rigidity composite material properties;
dispersing continuous fibers into a plurality of independent units, and assigning the variable angles and the same material properties to the discrete units according to the principle that the fiber angles corresponding to the same axial coordinates in a translation method are the same, so as to complete the simulation modeling of the continuous variable-angle fiber laying path;
step 3, carrying out variable angle parameter assignment on the discrete units in the model in the step 2 and establishing a finite element model;
step 4, carrying out a small amount of finite element analysis on the parametric analysis model, extracting sample points, and establishing an approximate model of buckling analysis and optimal design; establishing an approximate analysis model by a Latin hypercube method and a radial basis function neural network agent model method; the approximate model method is based on Gaussian radial basis functions and is realized by self-defined programming;
step 5, carrying out optimization design on the established approximate model; the optimization design is carried out aiming at the approximate model established in the step 4, and a multi-island genetic algorithm is adopted for optimization design.
Advantageous effects
The invention provides a modeling and optimal design method of a variable-rigidity composite material cylindrical shell; firstly, carrying out unit discretization design on a geometric model of a cylindrical shell; then, endowing discrete units of the geometric model of the cylindrical shell with variable-rigidity composite material properties; variable angle parameter assignment is carried out on discrete units in the model, and a finite element analysis model is established; secondly, carrying out a small amount of finite element analysis on the parameterized analysis model of the discrete units in the model, extracting sample points, and establishing an approximate model of buckling analysis and optimal design.
The modeling and optimization design method of the variable-stiffness composite material cylindrical shell realizes the establishment of a simulation analysis model of the variable-stiffness composite material cylindrical shell through a discrete finite element method and parametric programming, then introduces a neural network approximate model in the design optimization process, and takes the sequence of fiber laying starting and stopping angles and the number of layers as variables to perform quick response analysis and optimization design of the mechanical property of the cylindrical shell, thereby overcoming the problem of overlong analysis period of the variable-stiffness composite material in the application of the underwater pressure-resistant cylindrical shell and greatly improving the optimization design efficiency; the method has wide application and popularization significance in the aspects of improving the structural design and optimizing the efficiency of the variable-rigidity composite material cylindrical shell.
Drawings
The following describes the modeling and optimization design method of the variable-stiffness composite material cylindrical shell in detail with reference to the accompanying drawings and embodiments.
FIG. 1 is a two-dimensional schematic of a composite cylindrical shell geometry.
FIG. 2 is a schematic diagram of material distribution of a variable-stiffness cylindrical shell simulation analysis three-dimensional model.
Fig. 3 is a flow chart of modeling and optimization design of the variable stiffness pressure shell.
Detailed Description
The embodiment is a modeling and optimal design method of a variable-rigidity composite material cylindrical shell.
The modeling and optimal design method of the variable-rigidity composite material cylindrical shell comprises the steps of carrying out unit discretization design on a geometric model of the cylindrical shell; then, endowing discrete units of the geometric model of the cylindrical shell with variable-rigidity composite material properties; carrying out variable angle parameter assignment on discrete units in the model and establishing a finite element analysis model; secondly, carrying out a small amount of finite element analysis on the parameterized analysis model of the discrete units in the model, extracting sample points, and establishing an approximate model of buckling analysis and optimal design. The modeling and rapid optimization design method realizes the establishment of a simulation analysis model of the variable-stiffness composite material cylindrical shell through a discrete finite element method and parametric programming, introduces a neural network approximate model in the design optimization process, and performs rapid response analysis and optimization design on the mechanical property of the cylindrical shell by taking the fiber laying starting and stopping angle and the layering quantity sequence as variables.
Referring to fig. 1 to 3, a variable-stiffness composite cylindrical shell model with a certain geometric dimension is taken as an example; firstly, carrying out unit discretization design on the three-dimensional geometric model of the cylindrical shell;
the mathematical model expression for the fiber angle change is as follows:
in the formula, theta is the fiber angle, T0Is the angle of the fiber at the midpoint of the cylinder, T1Is the fiber angle at the cylindrical end point, x is the axial coordinate, and L is the cylindrical shell length; according to the expression, a three-dimensional geometric model of the cylindrical shell is evenly divided into a plurality of circular ring structures along an axial coordinate, namely an x axis, each circular ring is further refined into discrete units according to the expected grid size of finite element analysis, and the discrete units are grouped and numbered by utilizing programming codes.
According to the three-dimensional cylindrical shell structure and the grid serial number after grid discretization, a programming language is utilized and the end point and midpoint angle values in the formula are combined, the attribute of the analysis material to be used and the fiber angle corresponding to the grouping serial number are sequentially given to a simulation analysis model, and simulation modeling of the variable-angle fiber laying path is completed; the continuous fiber is dispersed into a plurality of independent units, the variable angle and the same material property are assigned to the discrete units according to the principle that the fiber angles corresponding to the same axial coordinate in the translation method are the same, and the purpose of simulating continuous variable-angle fiber laying is achieved by the method in the embodiment due to the fact that the variable span of the fiber angles of adjacent units is small.
Then according to the Gaussian radial basis function and the method for establishing an approximate model, the fiber angle T at the midpoint of the cylinder0And the fiber angle T at the end point of the cylinder1And establishing and optimizing an approximate model by using the variable and taking the buckling load response and the buckling maximization optimization design of the structure as targets. This section is implemented in terms of gaussian radial basis functions and programming language methods; the specific mathematical expression of the radial basis function is as follows:
in the formula, cjAndjrespectively representing the central parameter and the width parameter of the jth hidden layer of the neural network. The following equation represents the loss function of the approximate model system; after the sample point and the response value are input, when the loss function is smaller than the given precision, the approximate model building is completed.
In the formula (I), the compound is shown in the specification,
is the sample label, i.e. the response value of the sample point.
And then, according to the radial basis function and the loss function, carrying out a small amount of finite element analysis on the structural model under the working condition of bearing the combined pressure of the axial pressure and the lateral pressure, acquiring sample data of critical buckling response, and constructing and establishing an approximate model through a programming language.
After the approximate model meeting the precision requirement is obtained, the multi-island genetic algorithm is selected to carry out optimization design of the maximum sample response value on the approximate model. The results of the optimization and the comparison with the traditional finite element optimization method are shown in the table below, and the results show that the accuracy of the results is high, and meanwhile, the method greatly reduces the time of optimization design and improves the efficiency.
Modeling and optimization design of variable-rigidity cylindrical shell and efficiency result comparison
Claims (1)
1. A modeling and optimal design method of a variable-rigidity composite material cylindrical shell is characterized by comprising the following steps:
step 1, carrying out unit discretization design on a geometric model of a cylindrical shell;
discretizing the geometric model according to a mathematical model of linear change of the fiber path of the composite material, wherein the mathematical model of the fiber angle change is expressed as follows:
in the formula, theta is the fiber angle, T0Is the angle of the fiber at the midpoint of the cylinder, T1Is the fiber angle at the cylindrical end point, x is the axial coordinate, and L is the cylindrical shell length; according to the expression, the three-dimensional geometric model of the cylindrical shell is evenly divided into a plurality of circular ring structures along the axial coordinate, and each circular ring is further refined into discrete units according to the designed grid size of finite element analysis;
step 2, endowing the discrete units with variable-rigidity composite material properties;
dispersing continuous fibers into a plurality of independent units, and assigning the variable angles and the same material properties to the discrete units according to the principle that the fiber angles corresponding to the same axial coordinates in a translation method are the same, so as to complete the simulation modeling of the continuous variable-angle fiber laying path;
step 3, carrying out variable angle parameter assignment on the discrete units in the model in the step 2 and establishing a finite element model;
step 4, carrying out a small amount of finite element analysis on the parametric analysis model, extracting sample points, and establishing an approximate model of buckling analysis and optimal design; establishing an approximate analysis model by a Latin hypercube method and a radial basis function neural network agent model method; the approximate model method is based on Gaussian radial basis functions and is realized by self-defined programming;
step 5, carrying out optimization design on the established approximate model; the optimization design is carried out aiming at the approximate model established in the step 4, and a multi-island genetic algorithm is adopted for optimization design.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011038065.7A CN112182933A (en) | 2020-09-28 | 2020-09-28 | Modeling and optimal design method of variable-rigidity composite material cylindrical shell |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011038065.7A CN112182933A (en) | 2020-09-28 | 2020-09-28 | Modeling and optimal design method of variable-rigidity composite material cylindrical shell |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112182933A true CN112182933A (en) | 2021-01-05 |
Family
ID=73944328
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011038065.7A Pending CN112182933A (en) | 2020-09-28 | 2020-09-28 | Modeling and optimal design method of variable-rigidity composite material cylindrical shell |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112182933A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112255963A (en) * | 2020-09-24 | 2021-01-22 | 江苏理工学院 | Laying track planning method and device for fiber composite material on cylindrical surface |
| CN112836409A (en) * | 2021-02-03 | 2021-05-25 | 浙江工业大学 | Optimal design method of bistable composite shell |
| CN115495853A (en) * | 2022-11-18 | 2022-12-20 | 北京航空航天大学 | Particle blocking variable stiffness module parameter optimization method and device |
| CN116522420A (en) * | 2023-06-25 | 2023-08-01 | 山东石油化工学院 | Mechanical property simulation design method and system for curved surface column type composite material |
| CN117648846A (en) * | 2024-01-30 | 2024-03-05 | 中国空气动力研究与发展中心计算空气动力研究所 | Image sample generation method for composite material performance prediction modeling |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4938824A (en) * | 1987-01-23 | 1990-07-03 | Thiokol Corporation | Method for making a composite component using a transverse tape |
| US5749985A (en) * | 1989-06-20 | 1998-05-12 | Institut Francais Du Petrole | Process for optimizing multilayered tubes made of composite materials and tubes obtained through the process |
| US20120053907A1 (en) * | 2010-08-25 | 2012-03-01 | Livermore Software Technology Corporation | Efficient Data Management For Shell Finite Elements Representing Layered Composite Materials |
| CN107220404A (en) * | 2017-04-20 | 2017-09-29 | 江苏理工学院 | Composite material automobile accumulator housing design method based on multi-stage optimization |
| CN107526898A (en) * | 2017-09-13 | 2017-12-29 | 大连理工大学 | A kind of variation rigidity composite panel shell structure Accurate Model analysis and reliability-based optimization integrated design method |
| CN107563094A (en) * | 2017-09-21 | 2018-01-09 | 上海交通大学 | Three-dimensional woven carbon fibre composite fender optimization method |
| CN110321613A (en) * | 2019-06-25 | 2019-10-11 | 西北工业大学 | A kind of optimum design of laminate layup method of composite material pressure hull |
-
2020
- 2020-09-28 CN CN202011038065.7A patent/CN112182933A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4938824A (en) * | 1987-01-23 | 1990-07-03 | Thiokol Corporation | Method for making a composite component using a transverse tape |
| US5749985A (en) * | 1989-06-20 | 1998-05-12 | Institut Francais Du Petrole | Process for optimizing multilayered tubes made of composite materials and tubes obtained through the process |
| US20120053907A1 (en) * | 2010-08-25 | 2012-03-01 | Livermore Software Technology Corporation | Efficient Data Management For Shell Finite Elements Representing Layered Composite Materials |
| CN107220404A (en) * | 2017-04-20 | 2017-09-29 | 江苏理工学院 | Composite material automobile accumulator housing design method based on multi-stage optimization |
| CN107526898A (en) * | 2017-09-13 | 2017-12-29 | 大连理工大学 | A kind of variation rigidity composite panel shell structure Accurate Model analysis and reliability-based optimization integrated design method |
| CN107563094A (en) * | 2017-09-21 | 2018-01-09 | 上海交通大学 | Three-dimensional woven carbon fibre composite fender optimization method |
| CN110321613A (en) * | 2019-06-25 | 2019-10-11 | 西北工业大学 | A kind of optimum design of laminate layup method of composite material pressure hull |
Non-Patent Citations (8)
| Title |
|---|
| KECHUN SHEN 等: "Design Optimization of Composite Cylindrical Shell Under Hydrostatic Pressure", 《2018 OCEANS - MTS/IEEE KOBE TECHNO-OCEANS (OTO)》 * |
| WEI, RANFENG 等: "An efficient approach for stacking sequence optimization of symmetrical laminated composite cylindrical shells based on a genetic algorithm", 《THIN-WALLED STRUCTURES》 * |
| 吴双华 等: "变刚度圆柱壳屈曲分析及铺丝路径优化", 《复合材料科学与工程》 * |
| 廖宝华 等: "复合材料格栅加筋圆柱壳结构优化设计", 《复合材料：创新与可持续发展(下册)》 * |
| 张冰 等: "弯扭载荷作用变刚度复合材料圆柱壳的屈曲优化", 《南昌航空大学学报(自然科学版)》 * |
| 潘光 等: "复合材料圆柱壳体水下非线性屈曲数值分析", 《哈尔滨工程大学学报》 * |
| 肖泽平 等: "基于径向基函数神经网络的轨道板运输车优化设计", 《现代制造工程》 * |
| 顾杰斐 等: "考虑制造因素的变刚度层合板的抗屈曲铺层优化设计", 《复合材料学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112255963A (en) * | 2020-09-24 | 2021-01-22 | 江苏理工学院 | Laying track planning method and device for fiber composite material on cylindrical surface |
| CN112836409A (en) * | 2021-02-03 | 2021-05-25 | 浙江工业大学 | Optimal design method of bistable composite shell |
| CN115495853A (en) * | 2022-11-18 | 2022-12-20 | 北京航空航天大学 | Particle blocking variable stiffness module parameter optimization method and device |
| CN116522420A (en) * | 2023-06-25 | 2023-08-01 | 山东石油化工学院 | Mechanical property simulation design method and system for curved surface column type composite material |
| CN116522420B (en) * | 2023-06-25 | 2023-08-29 | 山东石油化工学院 | Mechanical property simulation design method and system for curved surface column type composite material |
| CN117648846A (en) * | 2024-01-30 | 2024-03-05 | 中国空气动力研究与发展中心计算空气动力研究所 | Image sample generation method for composite material performance prediction modeling |
| CN117648846B (en) * | 2024-01-30 | 2024-04-26 | 中国空气动力研究与发展中心计算空气动力研究所 | Image sample generation method for composite material performance prediction modeling |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112182933A (en) | Modeling and optimal design method of variable-rigidity composite material cylindrical shell | |
| CN108984910B (en) | A kind of parametric modeling method of Three-Dimensional Braided Composites | |
| CN104866673A (en) | Opening reinforcement method of shaft pressing reinforced cylindrical shell | |
| CN107169191A (en) | A kind of fan blade modeling method | |
| CN103729694B (en) | The method that improvement GA based on polychromatic sets hierarchical structure solves Flexible workshop scheduling | |
| CN111898295A (en) | Finite element modeling method of variable-rigidity composite laminated plate | |
| CN109145393B (en) | Bionic lightweight design method for structure | |
| CN110705032B (en) | Rapid modeling method for armored vehicle target vulnerability model | |
| CN108319799A (en) | A kind of more fidelity optimum design methods of the shape of Autonomous Underwater Vehicle | |
| CN109117504B (en) | Bidirectional functional gradient curved shell vibration analysis method | |
| CN107038303B (en) | Double-layer experimental design method based on proxy model and used for mechanical reliability analysis and design | |
| CN105005675A (en) | Composite insulator electric field optimization method based on multi-target genetic algorithm | |
| CN104750948A (en) | Optimization method for processing multiple extreme values and multiple restricted problems in flight vehicle design | |
| CN114756974B (en) | Wall surface distance calculation method considering object surface normal information | |
| CN108459993B (en) | Complex high-dimensional system optimization method based on rapid peak-tracking sampling | |
| CN111027250B (en) | Modeling method for special-shaped curved surface reinforcement shell based on grid deformation technology | |
| CN112417738A (en) | Numerical calculation method of spherical pressure-resistant shell containing random pit pitting defects | |
| CN108491612A (en) | The Finite Element Method of scheme of material selection is provided for multiple tube hydraulic bulging process | |
| CN108491654B (en) | Three-dimensional entity structure topology optimization method and system | |
| CN105468826A (en) | Design method of composite material | |
| CN104933261A (en) | High efficient sequential maximin latin hypercube design method | |
| CN110245408A (en) | A kind of steam turbine list circular arc pressure face Blade Design Method | |
| CN114996880A (en) | Composite armor structure optimization method based on ANSYS secondary development | |
| CN111666689B (en) | Feature line tracking method, reactor core neutron physical calculation method and reactor core neutron physical calculation device | |
| CN115019910A (en) | Modeling of special-shaped strand steel wire rope and load calculation method based on Adams and Abaqus combined simulation |
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 |