CN117910323A - Double-Double layer composite material C-type Liang Bianhou DEG optimization method - Google Patents
Double-Double layer composite material C-type Liang Bianhou DEG optimization method Download PDFInfo
- Publication number
- CN117910323A CN117910323A CN202410294681.0A CN202410294681A CN117910323A CN 117910323 A CN117910323 A CN 117910323A CN 202410294681 A CN202410294681 A CN 202410294681A CN 117910323 A CN117910323 A CN 117910323A
- Authority
- CN
- China
- Prior art keywords
- double
- composite material
- shaped beam
- region
- bianhou
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 122
- 238000005457 optimization Methods 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 67
- 238000005452 bending Methods 0.000 claims abstract description 31
- 238000004088 simulation Methods 0.000 claims abstract description 21
- 238000013178 mathematical model Methods 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 5
- 230000008676 import Effects 0.000 claims description 3
- 238000012805 post-processing Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 description 6
- 239000000463 material Substances 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000000265 homogenisation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- 239000011157 advanced composite material Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009830 intercalation Methods 0.000 description 1
- 230000002687 intercalation Effects 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
Abstract
The invention discloses a Double-Double layer composite material C-type Liang Bianhou DEG optimization method, belonging to the technical field of composite material C-type beam layer optimization; the method comprises the following specific steps: equivalently converting the laminated composite material of the C-shaped beam into a Double-Double laminated composite material; dividing the equivalent converted C-shaped beam into areas along the length direction; establishing a finite element simulation model of the bending process of the C-shaped beam; completing parameterized modeling of the bending process of the Double-Double layer composite material C-shaped beam based on a finite element simulation model; establishing a Double-Double layer composite material C-type Liang Bianhou DEG optimization mathematical model by taking the weight G minimum of the C-type beam structure as a target; and according to the mathematical model and the finite element simulation model, performing variable thickness optimization on the Double-Double layering C-shaped beam. The invention solves the problem of discontinuity such as layer interruption, bridging between layers and the like generated in the composite material.
Description
Technical Field
The invention belongs to the technical field of composite material C-shaped beam layering optimization, and particularly relates to a Double-Double layering composite material C-shaped Liang Bianhou-degree optimization method.
Background
Currently, in advanced composite materials used in the aerospace field, carbon fiber reinforced resin matrix composite materials are widely used, and four angles of 0 degree (+/-45 degrees) and 90 degrees are main choices of traditional layering composite materials, however, the fixed fiber angle direction limits the development of the composite materials. In this context, tsai proposes a Double-Double layering method in the form ofWherein/>And/>The number of times the sub-ply set is repeated, r representing the number of times the Double-Double ply composite material has only one optimal ply, i.e./>, at any angleAnd does not require central symmetry.
The method is characterized in that geometric modeling, finite element calculation and layering sequence optimization are serially connected, layering smoothening problems are converted into intercalation position and angle optimization problems in the optimization process, references are provided for composite material layering design, and the method has a certain practical value. Compared with the traditional layering method, the Double-Double layering method has a series of advantages of homogenization, light weight, strong designability, multiple layering types which can be selected and the like, and has great application potential in the field of composite materials. However, related researches on Double-Double layering are still in a starting stage, and researches on the thickening optimization design of the layering method are not reported at present.
Therefore, the invention optimally designs the Double-layer composite material, provides a Double-layer composite material C-type Liang Bianhou DEG optimizing method, and aims at the minimum total weight of the C-type beam on the premise of ensuring rigidity to realize the layer optimization of the Double-layer composite material C-type beam.
Disclosure of Invention
The technical problems to be solved are as follows:
In order to avoid the defects of the prior art, the invention provides the Double-Double composite material C-shaped Liang Bianhou DEG optimizing method, which can realize the ply optimization of the Double-Double composite material C-shaped beam by taking the minimum total weight of the C-shaped beam as the aim on the premise of ensuring the rigidity. The method solves the problem of discontinuous generation of layering interruption, bridging between layers and the like in the composite material, is not limited to layering optimization with a simple structure, can obtain a plurality of groups of feasible solutions and relatively optimal solutions, can realize layering optimization of Double-Double layering composite materials, and has a certain practical value.
The technical scheme of the invention is as follows: a Double-Double layer composite material C-type Liang Bianhou DEG optimization method comprises the following specific steps:
Equivalently converting the laminated composite material of the C-shaped beam into a Double-Double laminated composite material;
dividing the equivalent converted C-shaped beam into regions along the length direction, namely I represents the i-th region of the C-shaped beam, x i represents the number of layering layers of the i-th region of the C-shaped beam, and n represents the total number of regions of the C-shaped beam;
Establishing a finite element simulation model of the bending process of the C-shaped beam;
completing parameterized modeling of the bending process of the Double-Double layer composite material C-shaped beam based on a finite element simulation model;
Taking the total weight Gmin of the C-shaped beam as a target, and establishing a Double-Double layer composite material C-shaped Liang Bianhou-degree optimized mathematical model:
Wherein x i represents the number of layering layers in the ith region of the C-beam; i represents the i-th region of the C-beam, and n represents the total number of regions of the C-beam; g i represents the weight of the ith region of the C-beam as a function of the number of plies x i; g represents the total weight of the C-shaped beam; the constraint condition is that the number of layers of each region is required to be within a feasible region, and the bending rigidity K of the C-shaped beam is not lower than a specified value K 0;
And according to the mathematical model and the finite element simulation model, performing variable-thickness optimization on the Double-Double layer composite material C-shaped beam.
The invention further adopts the technical scheme that: the Double-Double layer composite material of the C-shaped beam is obtained by equivalent conversion of a traditional layer composite material through a tensile rigidity matrix, and the layer angle of the Double-Double layer composite material equivalent to the tensile rigidity of the traditional layer composite material is obtained; the traditional ply composite material comprises a ply composite material with four angles of 0 degree (+/-45 degrees) and 90 degrees.
The invention further adopts the technical scheme that: the calculation formula of the equivalent conversion is as follows:
wherein, Respectively represent the layering angles of Double-Double layering composite materials,/>For calculating the parameters, the formula is as follows:
Wherein a, b and c are the proportions of four angles of 0 degree, plus or minus 45 degrees and 90 degrees in the traditional layered composite material respectively;
The Double-Double layer composite material obtained by equivalent conversion is Where r represents the number of repetitions of the sub-ply and t represents the total laminate.
The invention further adopts the technical scheme that: the parameterized modeling method of the Double-Double laminated composite material C-shaped beam bending process comprises the following steps: firstly, generating script. Py initial script by utilizing script generated by preprocessing on the basis of a finite element simulation model established by Abaqus; secondly, simplifying the initial script, reserving effective operation, and deleting scaling and rotation ineffective operation; setting an import path and a position for parameters to be modified for subsequent direct calling; debugging the initial script after simplifying processing, finally writing the post-processing script, and extracting the needed data into the designated file.
The invention further adopts the technical scheme that: constraint conditions of the Double-Double layered composite material C-type Liang Bianhou-degree optimized mathematical model are as follows:
The number of layers in each region of the C-shaped beam should be a multiple of the number of sub-layers, namely a multiple of 4;
the difference between the layers of adjacent areas of the C-shaped beam is at most 4 layers;
The layer number of the C-shaped beam near the middle area is not lower than that of the edge area;
The ratio of the number of layers in the thickest region to the thinnest region of the C-beam is not greater than 2.
The invention further adopts the technical scheme that: the number of layers of each area of the C-shaped beam is not less than 20.
The invention further adopts the technical scheme that: the feasible region is R, expressed as:
wherein m represents the number of sub-plies of each region of the C-beam, each sub-ply comprising 4 plies; j represents the j-th region of the C-beam near the middle region.
The invention further adopts the technical scheme that: the variable thickness optimization method for the Double-Double laminated composite material C-shaped beam comprises the following steps:
Setting the initial thickness of a C-shaped beam according to the constraint condition of the Double-Double layered composite material C-shaped Liang Bianhou DEG optimization mathematical model so as to leave an optimization space;
Outputting script/py script from the initial thickness, submitting Abaqus to calculate to obtain load deflection curve and node stress-strain data of the C-shaped beam, and judging whether the rigidity of the C-shaped beam meets the requirement; entering the next step after meeting;
calculating strain energy data of each region of the C-shaped beam, and reducing the thickness of the minimum strain energy region according to the constraint condition;
And obtaining all feasible solutions in the feasible region R by adopting an exhaustion method, updating thickness information of the C-shaped beam according to the feasible solutions, changing a script, re-submitting Abaqus calculation until the bending rigidity of the C-shaped beam is lower than an optimization target, and outputting the optimization result of the previous step as a final result.
The invention further adopts the technical scheme that: the strain energy data calculation formula of the ith area of the Double-Double laminated composite material C-shaped beam is as follows:
Wherein: n i and x i represent the number of grid cells and the number of layers in the ith region of the C-beam, Is positively strained,/>Bending stiffness of p-th grid cell of q-th layer of C-shaped beam,/>Thickness of p-th grid cell of q-th layer of C-shaped beam,/>The area of the p-th grid cell of the q-th layer is represented, and T represents the transpose.
The invention further adopts the technical scheme that: the rule of reducing the thickness of the strain energy minimum area according to the constraint condition is to take four layers as a unit and reduce four layers at a time.
Advantageous effects
The invention has the beneficial effects that: in the invention, firstly, the traditional pavement composite material is equivalently converted into the Double-Double pavement composite material based on the tensile stiffness matrix, and the pavement angle of the Double-Double pavement composite material equivalent to the tensile stiffness of the traditional pavement composite material is obtained through the equivalent conversion. Only by converting the composite material of the C beam into Double-Double composite material, the problems of discontinuous layering, bridging between layers and the like can be solved, otherwise, the optimization cannot be performed through a thickness-variable method after the optimization.
Then, after the C-shaped beam is subjected to regional division, a finite element simulation model is established, and parameterized modeling of the Double-Double layer composite material C-shaped beam bending process is realized through Python language; the parametric modeling is adopted, so that the modeling efficiency can be improved, the complexity of the model is simplified, and the improvement and optimization of the finite element simulation model are promoted.
Then taking the total weight Gmin of the C-shaped beam as a target, establishing a Double-Double layer composite material C-shaped Liang Bianhou DEG optimization mathematical model, and giving specific constraint conditions and feasible domain expressions; the step can reduce the weight of the whole structure to the maximum extent under the condition of meeting the bending rigidity of the C-shaped beam, simplify the actual production and manufacturing process and improve the production efficiency.
And finally, performing variable-thickness optimization on the Double-Double layered composite material C-shaped beam according to the established mathematical model and the finite element simulation model. The greater the stored strain energy area, the stronger the load-bearing contribution to the C-beam under flexural loading. Therefore, according to the characteristic that the local thickness of the structure is changed under the condition that the material performance of the Double-Double layer can be kept unchanged, the thickness of the region with large strain energy of the C-shaped beam can be increased, and the thickness of the region with smaller strain energy is reduced, so that the use efficiency of the material performance in the structure is improved, and finally the aim of reducing the weight is fulfilled.
The method can realize the layering optimization of the Double-Double layering composite material C-shaped beam by taking the minimum total weight of the C-shaped beam as a target on the premise of ensuring the rigidity. The method can not cause discontinuous problems such as layer interruption, bridging between layers and the like in the composite material, is not limited to layer optimization with a simple structure, can obtain a plurality of groups of feasible solutions and relatively optimal solutions, can realize layer optimization of Double-Double layer composite material, and has a certain practical value. In addition, when the traditional pavement composite material is optimized by adopting the prior art means, the two sides of the symmetry plane need to be simultaneously deleted, and the DD pavement composite material is repeatedly formed by sub-pavement without deleting the two sides, so that the thickness of the DD pavement composite material cannot be optimized by adopting the prior art means. Meanwhile, the method can be widely applied to composite materials with other shapes and sizes, and has high universality.
Drawings
FIG. 1 is a flow chart of a method for optimizing a Double-Double lay-up composite C-type Liang Bianhou DEG in an embodiment of the invention;
FIG. 2 is a graph of a C-beam region division of a Double-Double lay-up composite C-Liang Bianhou degree optimization method in accordance with an embodiment of the present invention;
FIG. 3 is a parametric modeling flow chart of a Double-Double lay-up composite C-type Liang Bianhou DEG optimization method in an embodiment of the invention;
FIG. 4 is an optimization flow chart of a Double-Double lay-up composite C-type Liang Bianhou DEG optimization method in an embodiment of the invention;
FIG. 5 is a graph of load-deflection of a C-beam before and after optimization of a Double-Double lay-up composite C-Liang Bianhou degree optimization method in an embodiment of the invention.
Detailed Description
The embodiments described below by referring to the drawings are illustrative and intended to explain the present invention and should not be construed as limiting the invention.
Based on the problem that discontinuity such as layer interruption, bridging between layers and the like can be caused in the composite material when the traditional layer composite material is optimized in the prior art, the invention provides a Double-Double layer composite material C-type Liang Bianhou DEG optimizing method, which comprises the following specific steps:
Step 1: equivalently converting the laminated composite material of the C-shaped beam into a Double-Double laminated composite material;
step 2: dividing the equivalent converted C-shaped beam into regions along the length direction, namely I represents the i-th region of the C-shaped beam, x i represents the number of layering layers of the i-th region of the C-shaped beam, and n represents the total number of regions of the C-shaped beam;
step 3: establishing a finite element simulation model of the bending process of the C-shaped beam;
Step 4: completing parameterized modeling of the bending process of the Double-Double layer composite material C-shaped beam based on a finite element simulation model;
Step 5: taking the total weight Gmin of the C-shaped beam as a target, and establishing a Double-Double layer composite material C-shaped Liang Bianhou-degree optimized mathematical model:
Wherein x i represents the number of layering layers in the ith region of the C-beam; i represents the i-th region of the C-beam, and n represents the total number of regions of the C-beam; g i represents the weight of the ith region of the C-beam as a function of the number of plies x i; g represents the total weight of the C-shaped beam; the constraint condition is that the number of layers of each region is required to be within a feasible region, and the bending rigidity K of the C-shaped beam is not lower than a specified value K 0;
Step 6: and according to the mathematical model and the finite element simulation model, performing variable-thickness optimization on the Double-Double layer composite material C-shaped beam.
The technical scheme of the invention is further described in detail below with reference to the accompanying drawings 1-5.
Referring to fig. 1, the method for optimizing the Double-Double lay-up composite material at the angle of Liang Bianhou degrees in the embodiment comprises the following steps:
s1: equivalently converting the traditional layering composite material into a Double-Double layering composite material based on a tensile stiffness matrix;
In this example, the conventional lay-up composite was [ 90/45/0/45 ] 8t, and to obtain a lay-up angle of the Double-Double lay-up composite equivalent to the tensile stiffness of the conventional lay-up composite, parameters were calculated using the following formula ：
In the formula, a, b and c are the proportions of 0 degree, 45 degrees and 90 degrees in the traditional layering respectively. Will beSubstituted into the following formula:
Obtaining the layering angle of the Double-Double layering composite material Thus the Double-Double lay-up composite is [ 22/67/22/67 ] 8t.
S2: dividing the C-shaped beam into regions along the length direction, namelyI represents the i-th region of the C-beam, x i represents the number of layering layers of the i-th region of the C-beam;
In this embodiment, since the loading position is the middle part of the C-shaped beam, the C-shaped beam should be designed to satisfy symmetry, that is, the thickness of the C-shaped beam is symmetrical along the middle plane of the length direction, so that the C-shaped beam to be optimized is equally divided into 20 regions as shown in fig. 2 along the length direction, each region has a length of 35mm, and the thicknesses of the 20 regions are respectively optimized. Because of symmetry requirement, only 10 areas on the left side of the middle plane are needed to be used as design variables, and the number of layers of the 10 areas is respectively from edge to middle 。
S3: establishing an Abaqus finite element simulation model of the bending process of the C-shaped beam;
In this embodiment, an Abaqus finite element simulation was performed on the bending process of the conventional laid composite C-beam, and the bending load born when the C-deflection reached 1.5mm was 4655.49N.
S4: the parameterization modeling of the bending process of the Double-Double layer composite material C-shaped beam is realized through Python language;
In the present embodiment, the flow of obtaining the script. Py script by parametric modeling is shown in fig. 3. Firstly, establishing a standardized model through Abaqus, and generating a script/py initial script by utilizing a script generated by preprocessing; secondly, simplifying the initial script, reserving effective operation, and deleting invalid operations such as zooming and rotating; setting an import path and a position for parameters to be modified (the thickness and the layering number of each region model) so as to facilitate the subsequent direct call; and debugging the preprocessing script, and finally writing the preprocessing script, and extracting the required data into the designated file.
S5: taking the total weight Gmin of the C-shaped beam as a target, and establishing a Double-Double layer composite material C-shaped Liang Bianhou-degree optimized mathematical model:
Wherein x i represents the number of layering layers in the ith region of the C-beam; i represents the i-th region of the C-beam, and n represents the total number of regions of the C-beam; g i represents the weight of the ith region of the C-beam as a function of the number of plies x i; g represents the total weight of the C-shaped beam; the constraint condition is that the number of layers of each region is required to be within a feasible range, and the bending rigidity K of the C-shaped beam is not lower than a specified value K 0.
In this embodiment, the number of layers in each region of the C-beam should be a multiple of its sub-layering, i.e., a multiple of 4, due to the requirements of the Double-Double layering method. Meanwhile, in order to ensure the uniform change of the thickness of the C-shaped beam in the length direction, the difference of the layers of adjacent areas of the C-shaped beam is at most 4 layers in the optimal design, and the layer number of the area close to the middle is not lower than the thickness of the edge area and the ratio of the layer number of the thickest area to the thinnest area is not more than 2 because the middle area is loaded. In order to further improve the homogenization characteristics of the laminated plate, reduce the solidification deformation, and limit the number of layers of each region of the optimized C-shaped beam to be not less than 20.
According to the above several rules, the feasible region R of the design variable can be expressed as:
wherein m represents the number of sub-plies of each region of the C-beam, each sub-ply comprising 4 plies; j represents the j-th region of the C-beam near the middle region.
S6: performing variable-thickness optimization on the Double-Double layer composite material C-shaped beam according to the established mathematical model and the finite element simulation model;
In this embodiment, as shown in fig. 4, the specific steps in step S6 are as follows:
s601: setting the initial thickness of the C-shaped beam to 40 layers according to the optimization constraint condition of S5 so as to leave an optimization space;
S602: outputting script/py script from the initial thickness, submitting Abaqus to calculate to obtain load deflection curve and node stress-strain data of the C-shaped beam, and judging whether the rigidity of the C-shaped beam meets the requirement; after the satisfaction, the process goes to S603;
s603: and calculating strain energy data of each region, and reducing the thickness of the minimum strain energy region according to requirements. The greater the stored strain energy area, the stronger the load-bearing contribution to the C-beam under flexural loading. Therefore, according to the characteristic that the local thickness of the structure is changed under the condition that the material performance of the Double-Double layer can be kept unchanged, the thickness of the region with large strain energy of the C-shaped beam can be increased, and the thickness of the region with smaller strain energy is reduced, so that the use efficiency of the material performance in the structure is improved, and finally the aim of reducing the weight is fulfilled. The strain energy calculation formula of the ith area of the composite material C-shaped beam is as follows:
Wherein: n i and x i represent the number of grid cells and the number of layers in the ith region of the C-beam, Is positively strained,/>Bending stiffness of p-th grid cell of q-th layer of C-shaped beam,/>Thickness of p-th grid cell of q-th layer of C-shaped beam,/>The area of the p-th grid cell of the q-th layer is represented, and T represents the transpose.
S604: and obtaining all feasible solutions in the feasible region R by adopting an exhaustion method, updating the thickness information of the C-shaped beam according to the feasible solutions, changing the script, re-submitting Abaqus calculation until the rigidity of the C-shaped beam is lower than an optimization target, and outputting the optimization result of the previous step as a final result.
In this embodiment, the bending load experienced by a conventional lay-up composite C-beam when the deflection of the C-beam under bending load reaches 1.5mm is 4655.49N. Therefore, in order to ensure the rigidity of the optimized C-shaped beam, 110% and 100% of the load are used as constraint conditions for optimization, namely, bending acting forces born by the optimized C-shaped beam made of the variable-thickness Double-layer composite material under the same displacement load are not lower than 5121.04N and 4655.49N, and the optimized results obtained by optimizing the two constraint conditions are respectively an optimized result 1 and an optimized result 2, as shown in table 1.
Table 1 optimization results
The total mass of the C-shaped beams before and after optimization is compared, and compared with the non-optimized C-shaped beam with the equal thickness and 32 layers, the weight of the optimized result 1 is reduced by 12.5 percent, and the weight of the optimized result 2 is reduced by 17.5 percent. Simulation analysis is carried out on the optimized variable-thickness Double-layer composite material C-shaped beam, and the bending acting force born by the optimized result 1 under the displacement load of 1.5mm is 5144.32N 10.50% higher than 4655.49N of the traditional layer C-shaped beam; and the bending acting force born by the optimization result 2 is 4771.88N, which is improved by 2.51 percent. The load-deflection curve of the optimized C-shaped beam is shown in fig. 5, and it can be seen that the rigidity of the two optimized results is better than that of the traditional paved C-shaped beam.
The results of calculating the overall mean and overall variance of strain energy for each region of the Double-Double clad composite C-beam before and after optimization are shown in Table 2. As can be seen from the calculation results in table 2, compared with the conventional C-shaped beam made of the composite material of the laying layer before optimization, the average value of the strain energy absorbed by each region of the C-shaped beam after optimization is improved, the optimization result 1 and the optimization result 2 are respectively improved by 11.39% and 4.15%, which means that the C-shaped beam after optimization bears larger bending load during the same deformation, and therefore absorbs more strain energy. In addition, the overall variance of strain energy of each region of the optimized Double-Double composite material C-shaped beam is smaller than that of the conventional composite material C-shaped beam, and according to the data in the table, the overall variances of the optimized result 1 and the optimized result 2 are respectively reduced by 15.37% and 9.61%, so that the strain energy of each region of the optimized C-shaped beam is more uniform. Therefore, the thickness distribution of the optimized C-shaped beam is more reasonable, and the distribution of strain energy in the structure is more uniform and reasonable.
TABLE 2 calculation of Strain energy population means and variances
Although embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the invention, and that variations, modifications, alternatives, and variations may be made in the above embodiments by those skilled in the art without departing from the spirit and principles of the invention.
Claims (10)
1. A Double-Double layer composite material C-type Liang Bianhou DEG optimization method is characterized by comprising the following specific steps:
Equivalently converting the laminated composite material of the C-shaped beam into a Double-Double laminated composite material;
dividing the equivalent converted C-shaped beam into regions along the length direction, namely I represents the i-th region of the C-shaped beam, x i represents the number of layering layers of the i-th region of the C-shaped beam, and n represents the total number of regions of the C-shaped beam;
Establishing a finite element simulation model of the bending process of the C-shaped beam;
completing parameterized modeling of the bending process of the Double-Double layer composite material C-shaped beam based on a finite element simulation model;
Taking the total weight Gmin of the C-shaped beam as a target, and establishing a Double-Double layer composite material C-shaped Liang Bianhou-degree optimized mathematical model:
Wherein x i represents the number of layering layers in the ith region of the C-beam; i represents the i-th region of the C-beam, and n represents the total number of regions of the C-beam; g i represents the weight of the ith region of the C-beam as a function of the number of plies x i; g represents the total weight of the C-shaped beam; the constraint condition is that the number of layers of each region is required to be within a feasible region, and the bending rigidity K of the C-shaped beam is not lower than a specified value K 0;
And according to the mathematical model and the finite element simulation model, performing variable-thickness optimization on the Double-Double layer composite material C-shaped beam.
2. The optimization method of Double-Double layered composite material C-type Liang Bianhou degrees according to claim 1, which is characterized in that: the Double-Double layer composite material of the C-shaped beam is obtained by equivalent conversion of a traditional layer composite material through a tensile rigidity matrix, and the layer angle of the Double-Double layer composite material equivalent to the tensile rigidity of the traditional layer composite material is obtained; the traditional ply composite material comprises a ply composite material with four angles of 0 degree (+/-45 degrees) and 90 degrees.
3. The optimization method of Double-Double layered composite material C-type Liang Bianhou degrees according to claim 2, which is characterized in that: the calculation formula of the equivalent conversion is as follows:
wherein, 、/>Respectively represent the layering angles of Double-Double layering composite materials,/>For calculating the parameters, the formula is as follows:
Wherein a, b and c are the proportions of four angles of 0 degree, plus or minus 45 degrees and 90 degrees in the traditional layered composite material respectively;
The Double-Double layer composite material obtained by equivalent conversion is Where r represents the number of repetitions of the sub-ply and t represents the total laminate.
4. The method for optimizing the type-C Liang Bianhou degree of the Double-Double lay-up composite material according to claim 3, wherein the method comprises the following steps: the parameterized modeling method of the Double-Double laminated composite material C-shaped beam bending process comprises the following steps: firstly, generating script. Py initial script by utilizing script generated by preprocessing on the basis of a finite element simulation model established by Abaqus; secondly, simplifying the initial script, reserving effective operation, and deleting scaling and rotation ineffective operation; setting an import path and a position for parameters to be modified for subsequent direct calling; debugging the initial script after simplifying processing, finally writing the post-processing script, and extracting the needed data into the designated file.
5. The optimization method of Double-Double lay-up composite material C-type Liang Bianhou degrees according to claim 4, wherein the optimization method is characterized by comprising the following steps: constraint conditions of the Double-Double layered composite material C-type Liang Bianhou-degree optimized mathematical model are as follows:
The number of layers in each region of the C-shaped beam should be a multiple of the number of sub-layers, namely a multiple of 4;
the difference between the layers of adjacent areas of the C-shaped beam is at most 4 layers;
The layer number of the C-shaped beam near the middle area is not lower than that of the edge area;
The ratio of the number of layers in the thickest region to the thinnest region of the C-beam is not greater than 2.
6. The optimization method of Double-Double layered composite material C-type Liang Bianhou degrees according to claim 5, wherein the optimization method is characterized by comprising the following steps: the number of layers of each area of the C-shaped beam is not less than 20.
7. The optimization method of Double-Double layered composite material C-type Liang Bianhou degrees according to claim 5, wherein the optimization method is characterized by comprising the following steps: the feasible region is R, expressed as:
wherein m represents the number of sub-plies of each region of the C-beam, each sub-ply comprising 4 plies; j represents the j-th region of the C-beam near the middle region.
8. The optimization method of Double-Double lay-up composite material C-type Liang Bianhou degrees according to claim 7, wherein the optimization method is characterized by comprising the following steps: the variable thickness optimization method for the Double-Double laminated composite material C-shaped beam comprises the following steps:
Setting the initial thickness of a C-shaped beam according to the constraint condition of the Double-Double layered composite material C-shaped Liang Bianhou DEG optimization mathematical model so as to leave an optimization space;
Outputting script/py script from the initial thickness, submitting Abaqus to calculate to obtain load deflection curve and node stress-strain data of the C-shaped beam, and judging whether the rigidity of the C-shaped beam meets the requirement; entering the next step after meeting;
calculating strain energy data of each region of the C-shaped beam, and reducing the thickness of the minimum strain energy region according to the constraint condition;
And obtaining all feasible solutions in the feasible region R by adopting an exhaustion method, updating thickness information of the C-shaped beam according to the feasible solutions, changing a script, re-submitting Abaqus calculation until the bending rigidity of the C-shaped beam is lower than an optimization target, and outputting the optimization result of the previous step as a final result.
9. The optimization method of Double-Double lay-up composite material C-type Liang Bianhou degrees according to claim 8, wherein the optimization method is characterized by comprising the following steps: the strain energy data calculation formula of the ith area of the Double-Double laminated composite material C-shaped beam is as follows:
Wherein: n i and x i represent the number of grid cells and the number of layers of the ith region of the C-beam, Is positively strained,/>Bending stiffness of p-th grid cell of q-th layer of C-shaped beam,/>Thickness of p-th grid cell of q-th layer of C-shaped beam,/>The area of the p-th grid cell of the q-th layer is represented, and T represents the transpose.
10. The optimization method of Double-Double lay-up composite material C-type Liang Bianhou degrees according to claim 8, wherein the optimization method is characterized by comprising the following steps: the rule of reducing the thickness of the strain energy minimum area according to the constraint condition is to take four layers as a unit and reduce four layers at a time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410294681.0A CN117910323B (en) | 2024-03-15 | Double-Double layer composite material C-type Liang Bianhou DEG optimization method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410294681.0A CN117910323B (en) | 2024-03-15 | Double-Double layer composite material C-type Liang Bianhou DEG optimization method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117910323A true CN117910323A (en) | 2024-04-19 |
| CN117910323B CN117910323B (en) | 2024-05-31 |
Family
ID=
Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103573287A (en) * | 2013-11-14 | 2014-02-12 | 山东科技大学 | Method for determining thickness of flexible material of roadside flexible-strength double-layer composite support in gob-side entry retaining |
| CN104699900A (en) * | 2015-03-10 | 2015-06-10 | 西北工业大学 | Design method of composite material laminates with equal thickness |
| US20150299408A1 (en) * | 2012-12-27 | 2015-10-22 | Toray Industries, Inc. | Conductive fiber reinforced polymer composite and multifunctional composite |
| CN105005672A (en) * | 2015-08-04 | 2015-10-28 | 西北工业大学 | Sequential optimized design method of equal thickness laminated composite plate |
| CN105631073A (en) * | 2014-10-31 | 2016-06-01 | 中国航空工业集团公司西安飞机设计研究所 | Method for optimally manufacturing waved beam of composite material |
| CN106874573A (en) * | 2017-01-20 | 2017-06-20 | 西北工业大学 | A kind of method for designing of subregion Varying-thickness composite laminated plate |
| CN107103114A (en) * | 2017-03-23 | 2017-08-29 | 西安飞机工业（集团）有限责任公司 | A kind of composite hatch door optimization method |
| CN107256320A (en) * | 2017-07-28 | 2017-10-17 | 中科机器人科技有限公司 | A kind of small-sized hand based on composite plys throws the finite element method of unmanned plane |
| US20170371980A1 (en) * | 2015-08-11 | 2017-12-28 | Sabic Global Technologies B.V. | Multiple ply layered composite having low areal weight |
| CN108563905A (en) * | 2018-05-02 | 2018-09-21 | 吉林大学 | Automobile B-pillar reinforcement plate carbon fibre reinforced composite optimum design method |
| CN108897965A (en) * | 2018-07-10 | 2018-11-27 | 浙江工业大学 | A kind of design method of multistable composite shell |
| CN109583090A (en) * | 2018-11-30 | 2019-04-05 | 中国运载火箭技术研究院 | A kind of laying optimization method of composite structure |
| CN109583040A (en) * | 2018-11-07 | 2019-04-05 | 中国航空工业集团公司西安飞机设计研究所 | A kind of optimization method considering composite structure parametric continuity |
| CN110362943A (en) * | 2019-09-16 | 2019-10-22 | 成都飞机工业（集团）有限责任公司 | Multi-coordinate splices composite material skin Lay up design method |
| US20220009280A1 (en) * | 2020-04-20 | 2022-01-13 | Jilin University | Variable thickness carbon fiber composite rim structure and ply design method thereof |
| CN114282408A (en) * | 2021-12-15 | 2022-04-05 | 中国航空工业集团公司成都飞机设计研究所 | Layering optimization method for composite material laminated plate with lost layer |
| CN114547791A (en) * | 2022-01-25 | 2022-05-27 | 武汉理工大学 | Multi-stage optimization design method for layering of composite material laminated plate member |
| CN115146375A (en) * | 2022-06-20 | 2022-10-04 | 中国第一汽车股份有限公司 | Simulation analysis method for beam structure of bumper made of carbon fiber composite material |
| CN116136943A (en) * | 2023-04-20 | 2023-05-19 | 太原理工大学 | Composite material pressure vessel layering sequence optimization method |
| CN116882058A (en) * | 2023-07-12 | 2023-10-13 | 河南科技大学 | Composite material lightweight design and optimization method for automobile floor |
| CN117352108A (en) * | 2023-11-03 | 2024-01-05 | 中国飞机强度研究所 | Composite material laminated plate design method based on D-D layering |
Patent Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150299408A1 (en) * | 2012-12-27 | 2015-10-22 | Toray Industries, Inc. | Conductive fiber reinforced polymer composite and multifunctional composite |
| CN103573287A (en) * | 2013-11-14 | 2014-02-12 | 山东科技大学 | Method for determining thickness of flexible material of roadside flexible-strength double-layer composite support in gob-side entry retaining |
| CN105631073A (en) * | 2014-10-31 | 2016-06-01 | 中国航空工业集团公司西安飞机设计研究所 | Method for optimally manufacturing waved beam of composite material |
| CN104699900A (en) * | 2015-03-10 | 2015-06-10 | 西北工业大学 | Design method of composite material laminates with equal thickness |
| CN105005672A (en) * | 2015-08-04 | 2015-10-28 | 西北工业大学 | Sequential optimized design method of equal thickness laminated composite plate |
| US20170371980A1 (en) * | 2015-08-11 | 2017-12-28 | Sabic Global Technologies B.V. | Multiple ply layered composite having low areal weight |
| CN106874573A (en) * | 2017-01-20 | 2017-06-20 | 西北工业大学 | A kind of method for designing of subregion Varying-thickness composite laminated plate |
| CN107103114A (en) * | 2017-03-23 | 2017-08-29 | 西安飞机工业（集团）有限责任公司 | A kind of composite hatch door optimization method |
| CN107256320A (en) * | 2017-07-28 | 2017-10-17 | 中科机器人科技有限公司 | A kind of small-sized hand based on composite plys throws the finite element method of unmanned plane |
| CN108563905A (en) * | 2018-05-02 | 2018-09-21 | 吉林大学 | Automobile B-pillar reinforcement plate carbon fibre reinforced composite optimum design method |
| CN108897965A (en) * | 2018-07-10 | 2018-11-27 | 浙江工业大学 | A kind of design method of multistable composite shell |
| CN109583040A (en) * | 2018-11-07 | 2019-04-05 | 中国航空工业集团公司西安飞机设计研究所 | A kind of optimization method considering composite structure parametric continuity |
| CN109583090A (en) * | 2018-11-30 | 2019-04-05 | 中国运载火箭技术研究院 | A kind of laying optimization method of composite structure |
| CN110362943A (en) * | 2019-09-16 | 2019-10-22 | 成都飞机工业（集团）有限责任公司 | Multi-coordinate splices composite material skin Lay up design method |
| US20220009280A1 (en) * | 2020-04-20 | 2022-01-13 | Jilin University | Variable thickness carbon fiber composite rim structure and ply design method thereof |
| CN114282408A (en) * | 2021-12-15 | 2022-04-05 | 中国航空工业集团公司成都飞机设计研究所 | Layering optimization method for composite material laminated plate with lost layer |
| CN114547791A (en) * | 2022-01-25 | 2022-05-27 | 武汉理工大学 | Multi-stage optimization design method for layering of composite material laminated plate member |
| CN115146375A (en) * | 2022-06-20 | 2022-10-04 | 中国第一汽车股份有限公司 | Simulation analysis method for beam structure of bumper made of carbon fiber composite material |
| CN116136943A (en) * | 2023-04-20 | 2023-05-19 | 太原理工大学 | Composite material pressure vessel layering sequence optimization method |
| CN116882058A (en) * | 2023-07-12 | 2023-10-13 | 河南科技大学 | Composite material lightweight design and optimization method for automobile floor |
| CN117352108A (en) * | 2023-11-03 | 2024-01-05 | 中国飞机强度研究所 | Composite material laminated plate design method based on D-D layering |
Non-Patent Citations (5)
| Title |
|---|
| ZHANG H; ZHANG Z K; BAO J W; WANG Y W: "Study on Stealth Bulletproof Integrated Composite Material Based on Multi-layer Resistance Film", 《2020 6TH GLOBAL ELECTROMAGNETIC COMPATIBILITY CONFERENCE (GEMCCON)》, 22 June 2021 (2021-06-22) * |
| 宋广舒;郑艳萍;赵江铭;: "复合材料沉头搭接强度与渐进损伤研究", 玻璃钢/复合材料, no. 04, 28 April 2017 (2017-04-28) * |
| 彭波;王帅培;: "复合材料无人机弓形起落架有限元仿真与优化", 机械研究与应用, no. 03, 28 June 2020 (2020-06-28) * |
| 杜晨，彭雄奇: "变厚度连续纤维增强复合材料铺层设计优化方法", 《应用数学和力学》, vol. 43, no. 12, 31 December 2022 (2022-12-31) * |
| 汪玮;王志瑾;: "基于等效模型复合材料有效载荷承力结构优化设计", 飞机设计, no. 03, 15 June 2011 (2011-06-15) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kim et al. | Optimum design of composite structures with ply drop using genetic algorithm and expert system shell | |
| Paluch et al. | Combining a finite element programme and a genetic algorithm to optimize composite structures with variable thickness | |
| CN108009323A (en) | A kind of composite structure uncertain optimization method based on multi-state failure load | |
| CN105760619A (en) | Composite material layer optimization design method taking fatigue reliability into account | |
| CN108829914A (en) | A kind of structure of FRP structural member and process integration design method | |
| CN105500867A (en) | Novel metal-free connector composite material sandwiched plate connecting structure and design method thereof | |
| CN111898295A (en) | Finite element modeling method of variable-rigidity composite laminated plate | |
| CN110866353B (en) | Aircraft composite material structure optimization method based on strain neighborhood | |
| Qi-Hua et al. | Crashworthiness design of gradient bionic Al/CFRP hybrid tubes under multiple loading conditions | |
| CN112614014A (en) | Power and economy coordinated development method based on ecological environment view angle | |
| CN117910323B (en) | Double-Double layer composite material C-type Liang Bianhou DEG optimization method | |
| CN117910323A (en) | Double-Double layer composite material C-type Liang Bianhou DEG optimization method | |
| CN109583090B (en) | Layering optimization method of composite material structure | |
| Xu et al. | Analysis of corrugated steel web beam bridges using spatial grid modelling | |
| CN111274624B (en) | Multi-working-condition special-shaped node topology optimization design method based on RBF proxy model | |
| CN111597641B (en) | Material-structure-manufacturing integrated optimization design method suitable for thermoplastic braided composite material space bearing member | |
| CN112926244A (en) | Method for determining ultimate load of perforated part of composite laminated plate | |
| CN110728081B (en) | Composite material layering sequence optimizing system | |
| CN113011014B (en) | Composite material layering optimization method and system | |
| CN111310328A (en) | Point adding updating optimization method of gradient enhanced cooperative agent model | |
| CN115938514A (en) | Laminate press-shear failure analysis method considering buckling and first-layer failure competition | |
| CN112231805B (en) | Combined beam cable-stayed bridge cable force estimation method based on big data | |
| CN110704912B (en) | Method for topological optimization of bridge bracket arm structure under stress constraint | |
| CN113609591A (en) | Defect tolerance-oriented reinforced cylinder shell rapid optimization design method | |
| Yang et al. | Research on Comparative of Multi-Surrogate Models to Optimize Complex Truss Structures |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |