CN109657298A - A kind of carbon fibre composite interface phase mechanical property optimization method - Google Patents
A kind of carbon fibre composite interface phase mechanical property optimization method Download PDFInfo
- Publication number
- CN109657298A CN109657298A CN201811449293.6A CN201811449293A CN109657298A CN 109657298 A CN109657298 A CN 109657298A CN 201811449293 A CN201811449293 A CN 201811449293A CN 109657298 A CN109657298 A CN 109657298A
- Authority
- CN
- China
- Prior art keywords
- carbon fibre
- fibre composite
- mechanical property
- finite element
- element model
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 42
- 239000002131 composite material Substances 0.000 title claims abstract description 35
- 239000000835 fiber Substances 0.000 title claims abstract description 30
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 24
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000013461 design Methods 0.000 claims abstract description 25
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000010586 diagram Methods 0.000 claims description 7
- 230000002068 genetic effect Effects 0.000 claims description 5
- 239000007787 solid Substances 0.000 claims description 5
- 230000004075 alteration Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000012805 post-processing Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 3
- 229920000049 Carbon (fiber) Polymers 0.000 description 3
- 239000004917 carbon fiber Substances 0.000 description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000011157 advanced composite material Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000011347 resin Substances 0.000 description 2
- 229920005989 resin Polymers 0.000 description 2
- NQBHBXSQWRBTIE-VKHMYHEASA-N (2s)-2-(carbamoylamino)pentanediamide Chemical compound NC(=O)CC[C@@H](C(N)=O)NC(N)=O NQBHBXSQWRBTIE-VKHMYHEASA-N 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- Geometry (AREA)
- Game Theory and Decision Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Moulding By Coating Moulds (AREA)
- Reinforced Plastic Materials (AREA)
- Laminated Bodies (AREA)
Abstract
The application belongs to material structure mechanics field, in particular to a kind of carbon fibre composite interface phase mechanical property optimization method, comprising: Step 1: establishing finite element model in ANSYS, the design variable of the finite element model is parameterized;Step 2: reading the finite element model using ISIGHT, and identify the design variable in the finite element model;Step 3: choosing optimization algorithm in the ISIGHT, and define the parameter of the optimization algorithm;Step 4: optimizing iteration to the Optimized model, optimum results are obtained;Step 5: being post-processed to the optimum results.This method confirms that the performance at interface can be controlled by way of bonding, and it is further discovered that, interface performance is different to the various performance effects of composite material, and the interface of composite material is programmable, and can obtain " balance " or optimization of a kind of comprehensive performance.
Description
Technical field
The application belongs to material structure mechanics field, in particular to a kind of carbon fibre composite interface phase mechanical property is excellent
Change method.
Background technique
People have carried out for a long time being related on this critical problem between the microscopical structure and macro property of composite material
Research and exploration, try hard to design composite material with the principle of mesomechanics, push the development and development of advanced composite material (ACM),
And then the inherent law of Material reinforcement, toughening is grasped, especially in the past 30 years, achieving very big research achievement.As a result,
The design work of material is deep on deeper level, i.e. the microstructure design of material.It is right with going deep into for composite material research
In composite micro-structure optimization and design attention by more and more people.But the purpose largely to work is from boundary
The angle of the bonding performance in face is modified interface, and the interface phase micro-structure for lacking carbon fibre composite optimizes
Method.
Thus, it is desirable to have a kind of technical solution overcomes or at least mitigates at least one drawbacks described above of the prior art.
Summary of the invention
There is provided a kind of carbon fibre composite interface phase mechanical property optimization methods for the purpose of the application, with solution pair
The interface phase micro-structure of carbon fibre composite optimizes.
The technical solution of the application is:
A kind of carbon fibre composite interface phase mechanical property optimization method, comprising:
Step 1: establishing finite element model in ANSYS, the design variable of the finite element model is parameterized;
Step 2: reading the finite element model using ISIGHT, and identify the design variable in the finite element model;
Step 3: choosing optimization algorithm in the ISIGHT, and define the parameter of the optimization algorithm;
Step 4: optimizing iteration to the Optimized model, optimum results are obtained;
Step 5: being post-processed to the optimum results.
According at least one embodiment of the application, the design variable packet of the finite element model in the step 1
Include the thickness of the Young's modulus of the mutually each layer in interface, mutually each layer of radius and interface of the mutually each layer in interface.
According at least one embodiment of the application, the optimization algorithm includes neighborhood culture genetic algorithm.
According at least one embodiment of the application, the parameter of the optimization algorithm includes individual UVR exposure length, population
Scale, genetic algebra, crossing-over rate and aberration rate.
According at least one embodiment of the application, the step 3 further include:
Autoexec is created, position and the output optimization of the finite element model are defined in the autoexec
The position of destination file.
According at least one embodiment of the application, the post-processing includes drawing target variable to convert with design variable
Situation map, target variable dominate function and optimal solution with optimization number of run variation diagram, calculating.
It further include drawing the carbon fiber composite after the step 5 according at least one embodiment of the application
Expect the mutually each layer tensile modulus of elasticity figure in interface.
According at least one embodiment of the application, according to the mutually each layer tensile elasticity in the carbon fibre composite interface
Modulus figure draws out the curve graph that tensile modulus of elasticity changes with the interface phase thickness of certain root fiber.
According at least one embodiment of the application, the tensile modulus of elasticity is drawn out with certain root fiber using interpolation
Interface phase thickness variation curve graph.
According at least one embodiment of the application, the finite element model is that the plane of solid matter hexagon volume elements is answered
Power model.
At least there are following advantageous effects in the application:
Phase mechanical property optimization method in carbon fibre composite interface provided by the present application, it was confirmed that the performance at interface can be with
It is controlled by way of bonding, further investigation revealed that, the journey that interface performance influences the various performances of composite material
Degree is different, and the interface of composite material is programmable, and can obtain " balance " or optimization of a kind of comprehensive performance.
Detailed description of the invention
Fig. 1 is the flow chart of the carbon fibre composite interface phase mechanical property optimization method of the application;
Fig. 2 is FEM model schematic diagram in the embodiment of the present application;
Fig. 3 is the stratification state schematic diagram of finite element model in the embodiment of the present application;
Fig. 4 is the figure of changing of the target variable in the embodiment of the present application after finite element model Optimized Iterative;
Fig. 5 be in the embodiment of the present application target variable with optimization number of run variation diagram;
Fig. 6 is the schematic diagram in the embodiment of the present application after finite element model optimization;
Fig. 7 is the mutually each layer tensile modulus of elasticity datagram in the embodiment of the present application median surface;
Fig. 8 is the curve graph of each bed boundary phase thickness variation of tensile modulus of elasticity in the embodiment of the present application.
Specific embodiment
To keep the purposes, technical schemes and advantages of the application implementation clearer, below in conjunction in the embodiment of the present application
Attached drawing, technical solutions in the embodiments of the present application is further described in more detail.In the accompanying drawings, identical from beginning to end or class
As label indicate same or similar element or element with the same or similar functions.Described embodiment is the application
A part of the embodiment, instead of all the embodiments.The embodiments described below with reference to the accompanying drawings are exemplary, it is intended to use
In explanation the application, and it should not be understood as the limitation to the application.Based on the embodiment in the application, ordinary skill people
Member's every other embodiment obtained without creative efforts, shall fall in the protection scope of this application.Under
Face is described in detail embodiments herein in conjunction with attached drawing.
1 to Fig. 8 the application is described in further details with reference to the accompanying drawing.
As shown in Figure 1, a kind of carbon fibre composite interface phase mechanical property optimization method is provided in the present embodiment,
Its step are as follows:
Step 1: as shown in Fig. 2, establish finite element model in ANSYS, by the design variable of the finite element model into
Row parametrization;In this embodiment, the finite element model makees solid matter hexagon volume elements using solid matter hexagon volume elements
It is analyzed for representative volume elements, determines and select Plane stress model when optimizing to fiber distribution, target variable is matter
Amount, equivalent stress and YZ shear stress (axial shear stress), the design variable include the Young's modulus E of the mutually each layer in face, boundary
The thickness d of mutually each layer of radius r and interface of the mutually each layer in face, loads the finite element model, keeps its axial one end solid
Branch, the other end respectively add 1.5% alternating load in the direction Z+, Y-.As shown in figure 3, in the embodiment, the finite element model of building
Interface is mutually divided into five layers.
Step 2: the finite element model is exported as APDL file, the APDL file is read using ISIGHT, and
Identify the design variable in the wherein finite element model;
Step 3: choosing optimization algorithm in the ISIGHT, and define the parameter of the optimization algorithm;In the present embodiment
In, optimization algorithm selection field culture algorithm (NCGA), parameter is chosen as follows:
Individual UVR exposure length (bit) | Population scale | Genetic algebra | Crossing-over rate | Aberration rate |
20 | 20 | 125 | 1.0 | 0.01 |
Using the mutually each layer Young's modulus in interface as design variable, the value interval of design variable is as shown in the table:
Create a batch processing script file, defined in the batch processing script file position of the APDL file with
And the position of output optimum results file;
Step 4: running the batch processing script file in the ISIGHT, and the iteration of the optimization algorithm is carried out,
Obtain the data in the optimum results file;
Step 5: the data in the optimum results file are post-processed, in this embodiment, the post-processing packet
Include draw target variable with design variable change situation figure (as shown in Figure 4), target variable with optimize number of run variation diagram (such as
Shown in Fig. 5), calculate and dominate function and optimal solution.
In the above-described embodiments, finite element model enhances composite material using resin-based carbon fiber, optimizes repeatedly
Dai Hou, quality, equivalent stress and axial shear stress can tend to very small-scale Pareto solution system, the solution system simultaneously to be existed
Range is very small, can be approximately considered three target variables while reach minimum.
Pareto solution system is as shown in the table:
It is found that it is No. 14 design points, the geometry of No. 14 design points that it is maximum, which to dominate functional value, in above-mentioned Pareto solution system
Shape about the mutually each floor modulus in No. 14 design point interfaces as shown in fig. 6, be as follows:
E1 | 185.482822036743 |
E2 | 140.688501167297 |
E3 | 96.6746053695679 |
E4 | 62.1631479263306 |
E5 | 46.7315139770508 |
Wherein, each fiber geometrical characteristic of No. 14 design points such as following table (unit um):
In the present embodiment, the finite element model is optimized, to quality (mass), equivalent stress (seqv) and axis
It is optimized to stress (syz), effect of optimization is as shown in the table:
By above it is found that the finite element model in the present embodiment has significant improvement after optimizing.
After completing the procedure, it is drawn out as shown in Figure 7 for the mutually each floor modulus in No. 14 design point interfaces after optimization
The mutually each layer tensile modulus of elasticity figure in interface draw elasticity using interpolation and by every layer of elasticity modulus value at its center
The curve graph (as shown in Figure 8) of each bed boundary phase thickness variation of modulus, it is more convenient intuitively to see resin-based carbon fiber enhancing
The interface phase mechanical property optimum results of composite material.
The above, the only specific embodiment of the application, but the protection scope of the application is not limited thereto, it is any
Within the technical scope of the present application, any changes or substitutions that can be easily thought of by those familiar with the art, all answers
Cover within the scope of protection of this application.Therefore, the protection scope of the application should be with the scope of protection of the claims
It is quasi-.
Claims (10)
1. a kind of carbon fibre composite interface phase mechanical property optimization method characterized by comprising
Step 1: establishing finite element model in ANSYS, the design variable of the finite element model is parameterized;
Step 2: reading the finite element model using ISIGHT, and identify the design variable in the finite element model;
Step 3: choosing optimization algorithm in the ISIGHT, and define the parameter of the optimization algorithm;
Step 4: optimizing iteration to the Optimized model, optimum results are obtained;
Step 5: being post-processed to the optimum results.
2. phase mechanical property optimization method in carbon fibre composite interface according to claim 1, which is characterized in that described
The design variable of the finite element model in step 1 include the Young's modulus of the mutually each layer in interface, the mutually each layer in interface radius with
And the thickness of the mutually each layer in interface.
3. phase mechanical property optimization method in carbon fibre composite interface according to claim 1, which is characterized in that described
Optimization algorithm includes neighborhood culture genetic algorithm.
4. phase mechanical property optimization method in carbon fibre composite interface according to claim 1, which is characterized in that described
The parameter of optimization algorithm includes individual UVR exposure length, population scale, genetic algebra, crossing-over rate and aberration rate.
5. phase mechanical property optimization method in carbon fibre composite interface according to claim 1, which is characterized in that described
Step 3 further include:
Autoexec is created, position and the output optimum results of the finite element model are defined in the autoexec
The position of file.
6. phase mechanical property optimization method in carbon fibre composite interface according to claim 1, which is characterized in that described
Post-processing includes drawing target variable with design variable change situation figure, target variable with optimization number of run variation diagram, calculating
Dominate function and optimal solution.
7. phase mechanical property optimization method in carbon fibre composite interface according to claim 1, which is characterized in that described
It further include drawing the mutually each layer tensile modulus of elasticity figure in the carbon fibre composite interface after step 5.
8. phase mechanical property optimization method in carbon fibre composite interface according to claim 7, which is characterized in that according to
The mutually each layer tensile modulus of elasticity figure in carbon fibre composite interface draws out tensile modulus of elasticity with the interface of certain root fiber
The curve graph of phase thickness variation.
9. phase mechanical property optimization method in carbon fibre composite interface according to claim 8, which is characterized in that utilize
Interpolation draws out the curve graph that the tensile modulus of elasticity changes with the interface phase thickness of certain root fiber.
10. phase mechanical property optimization method in carbon fibre composite interface according to claim 1, which is characterized in that institute
State the Plane stress model that finite element model is solid matter hexagon volume elements.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811449293.6A CN109657298A (en) | 2018-11-30 | 2018-11-30 | A kind of carbon fibre composite interface phase mechanical property optimization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811449293.6A CN109657298A (en) | 2018-11-30 | 2018-11-30 | A kind of carbon fibre composite interface phase mechanical property optimization method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109657298A true CN109657298A (en) | 2019-04-19 |
Family
ID=66112068
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811449293.6A Pending CN109657298A (en) | 2018-11-30 | 2018-11-30 | A kind of carbon fibre composite interface phase mechanical property optimization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109657298A (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102634183A (en) * | 2012-05-03 | 2012-08-15 | 南京航空航天大学 | Rolling type interface of fiber reinforced resin matrix composite material and preparation method of interface |
CN106777693A (en) * | 2016-11-30 | 2017-05-31 | 中国航空工业集团公司沈阳飞机设计研究所 | A kind of real-time optimization method for designing based on CATIA digital-to-analogues |
CN106815408A (en) * | 2016-12-23 | 2017-06-09 | 上海交通大学 | Long fiber composites random structure is generated and its elastic performance Forecasting Methodology |
-
2018
- 2018-11-30 CN CN201811449293.6A patent/CN109657298A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102634183A (en) * | 2012-05-03 | 2012-08-15 | 南京航空航天大学 | Rolling type interface of fiber reinforced resin matrix composite material and preparation method of interface |
CN106777693A (en) * | 2016-11-30 | 2017-05-31 | 中国航空工业集团公司沈阳飞机设计研究所 | A kind of real-time optimization method for designing based on CATIA digital-to-analogues |
CN106815408A (en) * | 2016-12-23 | 2017-06-09 | 上海交通大学 | Long fiber composites random structure is generated and its elastic performance Forecasting Methodology |
Non-Patent Citations (1)
Title |
---|
单泉: "纤维增强复合材料界面相微结构优化设计", 《中国优秀硕士学位论文全文数据库工程科技I辑》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105512385B (en) | The method that the more material topological optimizations of structure are realized using sequence power function interpolation method | |
Nguyen et al. | IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies | |
CN111046505B (en) | Axial flow pump spoke parameter optimization design method based on response surface model | |
Vyverman et al. | essaMEM: finding maximal exact matches using enhanced sparse suffix arrays | |
Hu et al. | Evolution of RNA interference proteins dicer and argonaute in Basidiomycota | |
CN107065803B (en) | Flexible job shop dynamic dispatching method based on Weight variable scheduling interval | |
Vesting et al. | Development and application of optimisation algorithms for propeller design | |
CN107590325A (en) | A kind of fiber-reinforced composite materials structures optimization method based on Shepard interpolation | |
Chen et al. | Modified differential evolution algorithm using a new diversity maintenance strategy for multi-objective optimization problems | |
CN106599451B (en) | Multi-objective optimization method for RV reducer main bearing | |
CN109657298A (en) | A kind of carbon fibre composite interface phase mechanical property optimization method | |
Gao et al. | Evolution of RAD-and DIV-like genes in plants | |
CN110188468B (en) | Aeroelastic cutting optimization method and system for curved fiber composite material airfoil structure | |
Moezi et al. | 2&3-dimensional optimization of connecting rod with genetic and modified cuckoo optimization algorithms | |
JP2009045627A (en) | Method and program for optimizing press forming condition | |
Maside et al. | Patterns of molecular variation and evolution in Drosophila americana and its relatives | |
CN101582131B (en) | Two dimensional coding and decoding methods of discrete structure genetic optimization | |
WO2014041917A1 (en) | Flow pattern word expression method, word expression device, and program | |
Reutter et al. | Advanced endwall contouring for loss reduction and outflow homogenization for an optimized compressor cascade | |
CN114818369A (en) | Method, system, device and medium for designing continuous transonic wind tunnel section | |
Good et al. | Adaptive protein evolution and regulatory divergence in Drosophila | |
Banichuk et al. | Shape optimization of quasi-brittle axisymmetric shells by genetic algorithm | |
CN108595871B (en) | Hydraulic support optimal design method based on mean difference evolution | |
Woon et al. | Knowledge‐based algorithms in fixed‐grid GA shape optimization | |
CN111859807A (en) | Initial pressure optimizing method, device, equipment and storage medium for steam turbine |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190419 |
|
RJ01 | Rejection of invention patent application after publication |