CN112364546B - Fiber reinforced composite material structure optimization design method based on bilateral filtering - Google Patents
Fiber reinforced composite material structure optimization design method based on bilateral filtering Download PDFInfo
- Publication number
- CN112364546B CN112364546B CN202011403764.7A CN202011403764A CN112364546B CN 112364546 B CN112364546 B CN 112364546B CN 202011403764 A CN202011403764 A CN 202011403764A CN 112364546 B CN112364546 B CN 112364546B
- Authority
- CN
- China
- Prior art keywords
- design
- unit
- composite material
- reinforced composite
- fiber reinforced
- 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
- 238000013461 design Methods 0.000 title claims abstract description 77
- 238000005457 optimization Methods 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000001914 filtration Methods 0.000 title claims abstract description 34
- 239000003733 fiber-reinforced composite Substances 0.000 title claims abstract description 30
- 230000002146 bilateral effect Effects 0.000 title claims abstract description 28
- 239000000463 material Substances 0.000 title claims abstract description 26
- 239000000835 fiber Substances 0.000 claims abstract description 48
- 230000035945 sensitivity Effects 0.000 claims abstract description 39
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 26
- 238000006073 displacement reaction Methods 0.000 claims abstract description 18
- 239000002131 composite material Substances 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 238000012545 processing Methods 0.000 claims abstract description 4
- 238000004519 manufacturing process Methods 0.000 description 13
- 230000007547 defect Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
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]
-
- 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
Abstract
The invention belongs to the technical field related to the structural optimization design of fiber reinforced composite materials, and discloses a structural optimization design method of a fiber reinforced composite material based on bilateral filtering, which comprises the following steps: (1) Dividing a design domain of the fiber reinforced composite material structure into N units, and taking a fiber angle theta e at the center point of each unit as a design variable; (2) Establishing a unit stiffness matrix K e, and obtaining an overall displacement vector U and a unit displacement vector U e according to finite element analysis and calculation; (3) Calculating an objective function value according to an objective function formula c=f T U, and calculating a sensitivity value of the objective function with respect to the design variable θ e (4) Performing bilateral filtering processing on the sensitivity value at the center point of the unit e in a circular area N e; (5) Updating the design variable theta e by using a sensitivity-based optimization algorithm; (6) And (3) repeating the steps (2) to (5) until the optimization termination condition is met, thereby completing the optimization design conforming to the composite material structure. The invention improves the optimization quality.
Description
Technical Field
The invention belongs to the technical field related to the structural optimization design of fiber reinforced composite materials, and particularly relates to a structural optimization design method of a fiber reinforced composite material based on bilateral filtering.
Background
Compared with the conventional metal material, the fiber reinforced composite material has the advantages of high specific strength, large specific modulus, good fatigue resistance, corrosion resistance, good shock absorption performance and the like, and is widely applied to the fields of aerospace, automobiles, buildings and the like. With the continuous development of automatic tape laying/wire laying technology, the manufacture of the variable-rigidity composite material is possible, and fibers can be laid along a designed curve path. Compared with the linear fiber reinforced constant-stiffness structural design, the variable-stiffness composite structural design has better mechanical properties, and a designer can obtain a structure with better properties by optimally designing the fiber laying angle or path.
The fiber laying angle or path is reasonable, the structural mechanical property is met, and the fiber laying angle or path is an important content of the structural optimization design of the composite material. Existing optimization techniques have proven effective in providing design solutions with significantly enhanced structural performance, but inherent problems are that the fiber paths are not parallel, i.e. spatial continuity of fiber angles cannot be guaranteed, and adjacent fiber angles may be abrupt (vary widely), thus becoming one of the main sources of gaps or overlaps. In the paths set in these design methods, the amount of gaps, overlaps, can affect the structural response, manufacturing time, surface quality of the finished product, and can also present difficulties in modeling in numerical analysis. It is particularly important to consider the manufacturability of the structure during the optimization design stage, so as to minimize manufacturing defects such as gaps or overlapping caused by non-parallel fiber paths, and to ensure the manufacturing quality of the composite material.
Disclosure of Invention
Aiming at the defects or improvement demands of the prior art, the invention provides a fiber reinforced composite material structure optimization design method based on bilateral filtering, which divides a design domain of a fiber reinforced composite material structure into a limited number of units, sets a specific circular area for each unit, takes a unit center point as a circle center, has a radius not smaller than 10 times of the unit size, carries out bilateral filtering treatment on a sensitivity value at a unit e center point in the area to obtain a filtered sensitivity value, takes the filtered sensitivity value as the sensitivity value at the unit e center point, and then utilizes an optimization algorithm based on sensitivity information to update design variables until an optimization termination condition is met, thereby obtaining the fiber reinforced composite material structure design meeting manufacturing requirements, and greatly improving the manufacturing quality of the composite material.
In order to achieve the above object, according to one aspect of the present invention, there is provided a method for optimizing a fiber reinforced composite structure based on bilateral filtering, the method comprising the steps of:
(1) Taking a fiber reinforced composite material structure as a design domain, dividing the design domain into N units, taking a fiber angle theta e at the center point of each unit as a design variable, and setting an initial value of the fiber angle at the center point of each unit, wherein e=1, 2, … and N;
(2) Establishing a unit stiffness matrix K e, and obtaining an overall displacement vector U and a unit displacement vector U e according to finite element analysis and calculation, wherein the unit stiffness matrix depends on a fiber angle theta e;
(3) Calculating an objective function value c according to an objective function formula c=f T U, and calculating a sensitivity value of the objective function with respect to the design variable θ e
(4) Setting a circular area N e with the center point of the unit e as the center, and carrying out bilateral filtering processing on the sensitivity value at the center point of the unit e in the circular area N e to obtain a filtered sensitivity value
(5) Updating the design variable theta e by using a sensitivity-based optimization algorithm;
(6) And (3) repeating the steps (2) to (5) until the optimization termination condition is met, thereby completing the optimization design conforming to the composite material structure.
Further, by the calculation formulaA cell stiffness matrix K e is established that depends on the fiber angle θ e, where dΩ is the area infinitesimal, Ω e is the area of the cell e, B is the displacement strain matrix, and D (θ e) is the cell elastic matrix that depends on the fiber angle θ e.
Further, the calculation formula of D (θ e) is: d (θ e)=T(θe)D0T(θe)T, where D 0 is the original elastic matrix when the fiber is not rotating, and T (θ e) is the rotation matrix.
Further, in step (3), the sensitivity value of the objective function with respect to the design variable is derived using the relationship between the unit stiffness matrix K e and θ e The calculation formula is/>Where U e is a unit displacement vector corresponding to unit e, and the number of unit e is n e, U e is a column vector composed of elements in the 2n e th and 2n e +1 th rows of the overall displacement vector U.
Further, the filtered sensitivity valueThe calculation formula is as follows:
Wherein θ p is the fiber angle at any cell center point within circular region N e; h (e) is a weight function, and the calculation formula is W c (p-e) is a domain filter function, and p-e represents the Euclidean distance between two unit center points; w s(|θp-θe |) is a range filter function, |θ p-θe | represents the absolute value of the difference in fiber angle at the center point of the two units; the domain filter function W c (x) is calculated/>, using the following formulaWherein σd is the domain filter coefficient; the range filter function W s (x) is calculated/>, using the following formulaWhere σr is the range filter coefficient.
Further, the sensitivity value of the processed objective function with respect to the design variable is used to determine the sensitivity value of the design variable using the objective function value cAnd upper and lower bounds θ max and θ min of the design variables, and updating the design variables in combination with a moving asymptote algorithm.
Further, θ max and θ min are respectively:
θ max=90°-ε,θmin = -90 °, where ε is a very small positive value.
Further, the mathematical expression of the optimal design model adopted by the method is as follows:
findθe(e=1,2,...,1600)
min c=FTU
s.t. KU=F
θmin≤θe≤θmax
Where the fiber angle value θ e at the cell center point is a design variable, the objective function is compliance c, the design objective minimizes the compliance c of the structure, and the constraints include the equilibrium equation ku=f, the upper and lower bounds θ max and θ min of θ e.
In general, compared with the prior art, the fiber reinforced composite material structure optimization design method based on bilateral filtering mainly has the following beneficial effects:
1. Dividing the design domain of the fiber reinforced composite material structure into a limited number of units, taking the center point of each unit as the center of a circle, setting a specific circular area, carrying out bilateral filtering treatment on the sensitivity value at the center point of the unit e in the corresponding area to obtain a filtered sensitivity value as the sensitivity value at the center point of the unit e, and updating design variables by using an optimization algorithm based on sensitivity information until the optimization termination condition is met, so that the fiber reinforced composite material structure design meeting the manufacturing requirement is obtained, and the manufacturing quality of the composite material is greatly improved.
2. The fiber reinforced composite material obtained by adopting the method of the invention is optimally designed, so that the manufacturing defects of gaps and overlapping caused by unparallel fiber ribbons can be avoided as much as possible, and the manufacturing quality of the composite material is improved to a great extent.
3. The method provided by the invention is easy to implement, has strong applicability and is beneficial to popularization and application.
Drawings
FIG. 1 is a schematic flow chart of a fiber reinforced composite structure optimization design method based on bilateral filtering;
FIG. 2 is a schematic diagram of an example of an optimized design of a planar simply supported beam structure provided by an embodiment of the present invention;
FIG. 3 is a schematic view of an arrangement with respect to a circular region N e;
FIG. 4 is a schematic view of the fiber angle layout at the cell center point obtained using the method provided by the present invention for the example of FIG. 2;
FIG. 5 is a schematic view of the fiber angle layout at the center point of the cell calculated without the method provided by the present invention for the example of FIG. 2.
Detailed Description
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention. In addition, the technical features of the embodiments of the present invention described below may be combined with each other as long as they do not collide with each other.
Referring to fig. 1, the method for optimizing the design of the fiber reinforced composite material structure based on bilateral filtering provided by the invention mainly comprises the following steps:
Step one, taking a fiber reinforced composite material structure as a design domain, dividing the design domain into N units, taking a fiber angle theta e at the center point of each unit as a design variable, and setting an initial value of the fiber angle at the center point of each unit, wherein e=1, 2, … and N.
Step two, a unit stiffness matrix K e is established, and then an overall displacement vector U and a unit displacement vector U e are obtained according to finite element analysis and calculation, wherein the unit stiffness matrix depends on a fiber angle theta e.
Step three, calculating an objective function value c according to an objective function formula c=f T U, and calculating a sensitivity value of the objective function with respect to the design variable θ e
Wherein the sensitivity valueThe calculation formula of (2) is as follows:
Setting a circular area N e with the center point of the unit e as the center, and carrying out bilateral filtering processing on the sensitivity value at the center point of the unit e in the circular area N e to obtain a filtered sensitivity value
Wherein the radius of the circular region N e is not less than 10 times the cell size; filtered sensitivity valuesThe following calculation formula is adopted:
Wherein θ p is the fiber angle at any cell center point within circular region N e; h (e) is a weight function, H -1 (e) is the reciprocal thereof to ensure that the output of the bilateral filtering is the average value of sensitivity near similar values, and the calculation formula is W c (p-e) is a domain filter function, and p-e represents the Euclidean distance between two unit center points; w s(|θp-θe |) is a range filter function, |θ p-θe | represents the absolute value of the difference in fiber angle at the center point of the two units.
The domain filter function W c (x) is calculated using the following formula:
Wherein σd is a domain filter coefficient, and the selection of this parameter directly affects the performance of bilateral filtering, and needs to be adjusted according to different calculation examples.
The range filter function W s (x) is calculated using the following formula:
Wherein σr is a range filter coefficient, and the selection of the parameter directly affects the performance of bilateral filtering and needs to be adjusted according to different calculation examples.
And fifthly, updating the design variable theta e by using a sensitivity-based optimization algorithm.
And step six, repeating the steps 2 to 5 until the optimization termination condition is met, thereby completing the optimization design conforming to the material structure.
The present invention will be described in further detail with reference to the following examples.
Examples
Referring to fig. 2 to 5, the present embodiment illustrates the present invention by taking an optimization problem of minimizing the flexibility of the simple beam structure with in-plane load as an example, and given a rectangular design domain with a size of 1m×4m, the lower left corner and the lower right corner are fixed supporting points, and the lower center receives a downward concentrated load F with a size of 1N.
And (3) performing fiber angle layout optimization on the fiber reinforced composite material simply supported beam structure to minimize the flexibility, wherein the method comprises the following specific steps:
Step one, dividing the design domain of the fiber reinforced composite structure into 20×80 square units, wherein the unit side length is 0.05m, and setting the initial value of the fiber angle theta e at the center point of each unit to be 0 °, namely the fiber level, wherein e=1, 2, …,1600.
Step two, through calculation formulaEstablishing a cell stiffness matrix K e depending on a fiber angle theta e, wherein dΩ is an area infinitesimal, Ω e is an area of a cell e, B is a displacement strain matrix, D (theta e) is a cell elastic matrix depending on the fiber angle theta e, a calculation formula is D (theta e)=T(θe)D0T(θe)T,D0 is an original elastic matrix when the fiber is not rotated, and T (theta e) is a rotation matrix, and the calculation formula is respectively adopted as follows:
Where E x and E y are Young's modulus, G xy is shear modulus, v xy and v yx are Poisson's ratio, satisfying v xyEy=νyxEx. The overall rigidity matrix K is obtained through assembling K e, belongs to basic operation in a finite element analysis method, and then the overall displacement field U is obtained through solving according to a formula KU=F, wherein F is an external vector.
Step three, calculating an objective function value c according to c=f T U, and simultaneously deducing a sensitivity value of the objective function on the design variable by utilizing the relation between the unit stiffness matrix K e and the theta e The calculation formula is/>In the formula, U e is a unit displacement vector corresponding to a unit e, and the number of the unit e is n e, and U e is a column vector formed by elements of the 2n e th row and the 2n e +1 th row of the overall displacement vector U.
Fourth, referring to fig. 3, a circular area N e is set with a radius of 0.5m (i.e. 10 times of the unit side length) around the center point of the unit e, and the sensitivity value at the center point of the unit e is filtered in the area to obtain a filtered sensitivity valueThe calculation formula is as follows:
Wherein θ p is the fiber angle at any cell center point within circular region N e; h (e) is a weight function, H -1 (e) is the reciprocal thereof to ensure that the output of the bilateral filtering is the average value of sensitivity near similar values, and the calculation formula is W c (p-e) is a domain filter function, and p-e represents the Euclidean distance between two unit center points; w s(|θp-θe |) is a range filter function, |θ p-θe | represents the absolute value of the difference in fiber angle at the center point of the two units; the domain filter function W c (x) is calculated using the following formula: /(I)Wherein σd is a domain filtering coefficient, and the selection of the parameter directly affects the performance of bilateral filtering and needs to be adjusted according to different calculation examples; the range filter function W s (x) is calculated using the following formula: /(I)Where σr is a range filter coefficient, and the selection of this parameter directly affects the performance of bilateral filtering, and needs to be adjusted according to different calculation examples.
Step five, using the objective function value c, the sensitivity value of the processed objective function with respect to the design variableAnd upper and lower bounds θ max and θ min of the design variables, θ max=90°-ε,θmin = -90 ° where ε is a very small positive value, here 1×10 -8, and updating the design variables in combination with a moving asymptote algorithm (Method of Moving Asymptotes, abbreviated as MMA).
Step six, repeating the steps two to five, wherein each repetition is called an iteration process, until the optimization termination condition is met, completing the whole optimization design, and the optimization termination condition is simply set to be that the iteration times reach 50 times.
In summary, the fiber reinforced composite structure optimization design model based on bilateral filtering can be summarized as follows:
findθe(e=1,2,...,1600)
min c=FTU
s.t. KU=F
θmin≤θe≤θmax
Wherein the fiber angle value θ e at the cell center point is a design variable, the objective function is compliance c, the design objective minimizes the compliance c of the structure, the constraint includes the equilibrium equation ku=f, and the upper and lower bounds of θ e are θ max and θ min, respectively.
The optimization result of the invention is as follows: the fiber angle layout obtained after filtering by the method is shown in figure 4, and the flexibility value of the objective function is 99.3; the fiber angle layout obtained without any filtering method is shown in fig. 5, the flexibility value of the objective function is 36.7, and the structural rigidity obtained by filtering by the method is smaller than that obtained without any filtering method, but the structural fiber angle layout obtained by filtering by the method is reasonable, is more convenient to manufacture, can meet the manufacturing requirement, and greatly improves the manufacturing quality of the composite material.
It will be readily appreciated by those skilled in the art that the foregoing description is merely a preferred embodiment of the invention and is not intended to limit the invention, but any modifications, equivalents, improvements or alternatives falling within the spirit and principles of the invention are intended to be included within the scope of the invention.
Claims (6)
1. A fiber reinforced composite material structure optimization design method based on bilateral filtering is characterized in that:
(1) Taking a fiber reinforced composite material structure as a design domain, dividing the design domain into N units, taking a fiber angle theta e at the center point of each unit as a design variable, and setting an initial value of the fiber angle at the center point of each unit, wherein e=1, 2, … and N;
(2) Establishing a unit stiffness matrix K e, and obtaining an overall displacement vector U and a unit displacement vector U e according to finite element analysis and calculation, wherein the unit stiffness matrix depends on a fiber angle theta e;
(3) Calculating an objective function value c according to an objective function formula c=f T U, and calculating a sensitivity value of the objective function with respect to the design variable θ e
(4) Setting a circular area N e with the center point of the unit e as the center, and carrying out bilateral filtering processing on the sensitivity value at the center point of the unit e in the circular area N e to obtain a filtered sensitivity value
(5) Updating the design variable theta e by using a sensitivity-based optimization algorithm;
(6) Repeating the steps (2) to (5) until the optimization termination condition is met, thereby completing the optimization design conforming to the composite material structure;
balance equation ku=f; filtered sensitivity values The calculation formula is as follows:
Wherein θ p is the fiber angle at any cell center point within circular region N e; h (e) is a weight function, and the calculation formula is W c (p-e) is a domain filter function, and p-e represents the Euclidean distance between two unit center points; w s(|θp-θe |) is a range filter function, |θ p-θe | represents the absolute value of the difference in fiber angle at the center point of the two units; the domain filter function W c (x) is calculated/>, using the following formulaWherein σd is the domain filter coefficient; the range filter function W s (x) is calculated/>, using the following formulaWhere σr is the range filter coefficient;
Sensitivity values of the processed objective function with respect to the design variables using the objective function value c And upper and lower bounds θ max and θ min of the design variables, and updating the design variables in combination with a moving asymptote algorithm.
2. The fiber reinforced composite material structure optimization design method based on bilateral filtering as claimed in claim 1, wherein the method comprises the following steps: by calculationA cell stiffness matrix K e is established that depends on the fiber angle θ e, where dΩ is the area infinitesimal, Ω e is the area of the cell e, B is the displacement strain matrix, and D (θ e) is the cell elastic matrix that depends on the fiber angle θ e.
3. The fiber reinforced composite material structure optimization design method based on bilateral filtering as claimed in claim 2, wherein the method is characterized in that: d (θ e) has the formula: d (θ e)=T(θe)D0T(θe)T, where D 0 is the original elastic matrix when the fiber is not rotating, and T (θ e) is the rotation matrix.
4. The fiber reinforced composite material structure optimization design method based on bilateral filtering as claimed in claim 1, wherein the method comprises the following steps: in step (3), the sensitivity value of the objective function with respect to the design variable is derived using the relationship between the unit stiffness matrix K e and θ e The calculation formula is/>Where U e is a unit displacement vector corresponding to unit e, and the number of unit e is n e, U e is a column vector composed of elements in the 2n e th and 2n e +1 th rows of the overall displacement vector U.
5. The fiber reinforced composite material structure optimization design method based on bilateral filtering as claimed in claim 1, wherein the method comprises the following steps: θ max and θ min are respectively:
θ max=90°-ε,θmin = -90 °, where ε is a very small positive value.
6. The fiber reinforced composite structure optimization design method based on bilateral filtering as in any one of claims 1-4, wherein: the mathematical expression of the optimal design model adopted by the method is as follows:
findθe,e=1,2,...,1600
min c=FTU
s.t.KU=F
θmin≤θe≤θmax
Where the fiber angle value θ e at the cell center point is a design variable, the objective function is compliance c, and the design objective minimizes the compliance c of the structure, the upper and lower bounds θ max and θ min of θ e.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011403764.7A CN112364546B (en) | 2020-12-03 | 2020-12-03 | Fiber reinforced composite material structure optimization design method based on bilateral filtering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011403764.7A CN112364546B (en) | 2020-12-03 | 2020-12-03 | Fiber reinforced composite material structure optimization design method based on bilateral filtering |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112364546A CN112364546A (en) | 2021-02-12 |
CN112364546B true CN112364546B (en) | 2024-04-19 |
Family
ID=74535938
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011403764.7A Active CN112364546B (en) | 2020-12-03 | 2020-12-03 | Fiber reinforced composite material structure optimization design method based on bilateral filtering |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112364546B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113515850B (en) * | 2021-05-20 | 2024-02-06 | 大连理工大学 | Fiber-reinforced composite material structural layout optimization design method considering fiber continuity |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590325A (en) * | 2017-08-29 | 2018-01-16 | 华中科技大学 | A kind of fiber-reinforced composite materials structures optimization method based on Shepard interpolation |
CN109460577A (en) * | 2018-10-10 | 2019-03-12 | 华中科技大学 | A kind of variation rigidity composite structural optimization design method guaranteeing manufacturing quality |
CN111444579A (en) * | 2020-03-11 | 2020-07-24 | 华中科技大学 | Composite material structure optimization design method considering manufacturability |
WO2020211012A1 (en) * | 2019-04-17 | 2020-10-22 | 大连理工大学 | Quick collaborative optimization method for hybrid fiber composite plate-shell structure |
-
2020
- 2020-12-03 CN CN202011403764.7A patent/CN112364546B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107590325A (en) * | 2017-08-29 | 2018-01-16 | 华中科技大学 | A kind of fiber-reinforced composite materials structures optimization method based on Shepard interpolation |
CN109460577A (en) * | 2018-10-10 | 2019-03-12 | 华中科技大学 | A kind of variation rigidity composite structural optimization design method guaranteeing manufacturing quality |
WO2020211012A1 (en) * | 2019-04-17 | 2020-10-22 | 大连理工大学 | Quick collaborative optimization method for hybrid fiber composite plate-shell structure |
CN111444579A (en) * | 2020-03-11 | 2020-07-24 | 华中科技大学 | Composite material structure optimization design method considering manufacturability |
Non-Patent Citations (2)
Title |
---|
基于灵敏度分析的三维结构等几何形状优化方法;陈龙 等;中国机械工程;20170724(第14期);全文 * |
复合材料框架结构刚度几何双尺度优化设计;段尊义;阎军;张瑞;赵国忠;;计算力学学报;20160215(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112364546A (en) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110955941B (en) | Vector field-based composite material structure optimization design method and device | |
CN108629140B (en) | Porous composite material structure design optimization method based on geodesic distance | |
CN108829914B (en) | Structure and process integrated design method of FRP structural member | |
CN107729648A (en) | A kind of wavy fiber composite structural design Waterfall type multilevel optimization method based on Shepard interpolation | |
CN107590325A (en) | A kind of fiber-reinforced composite materials structures optimization method based on Shepard interpolation | |
CN108763658B (en) | Combined thin-wall structure natural frequency design method based on isogeometric method | |
CN109460577B (en) | Variable-rigidity composite material structure optimization design method for guaranteeing manufacturing quality | |
CN112446163B (en) | Energy finite element topological optimization method based on parameterized level set | |
CN111723457B (en) | Level set method for optimization design of fiber curve laying variable-stiffness structure | |
CN112364546B (en) | Fiber reinforced composite material structure optimization design method based on bilateral filtering | |
CN108416083B (en) | Two-dimensional dynamic model analysis method and system for towering television tower structure | |
CN112836411A (en) | Method and device for optimizing structure of stiffened plate shell, computer equipment and storage medium | |
CN108180918A (en) | A kind of cloud geodetic path forward trace generation method and device | |
CN111444579B (en) | Composite material structure optimization design method considering manufacturability | |
CN109255141A (en) | A kind of body of a motor car forward direction conceptual design cross sectional shape optimization method | |
CN113191048A (en) | Optimal design method of equidistant curve fiber reinforced structure | |
CN113505405A (en) | Equivalent load obtaining method, and topology optimization method and system based on equivalent load | |
CN113742970B (en) | Density changing method based on level set for lattice structure optimization design | |
CN113779825A (en) | Structure optimization design method based on unit decomposition and mobile control point | |
CN113515850B (en) | Fiber-reinforced composite material structural layout optimization design method considering fiber continuity | |
CN112818576B (en) | Multi-level optimization method for curve fiber composite structure design | |
Belblidia et al. | Integrated design optimization of stiffened plate structures | |
CN112818583A (en) | Equivalent dead load obtaining method, topology optimization method and system | |
CN116611213A (en) | Topological optimization method and system for linear elastic structure | |
CN115563838B (en) | COMSOL-based seismic metamaterial energy band structure and frequency domain analysis method |
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 |