CN103455679A - Finite element analysis method of composite material - Google Patents
Finite element analysis method of composite material Download PDFInfo
- Publication number
- CN103455679A CN103455679A CN2013104030393A CN201310403039A CN103455679A CN 103455679 A CN103455679 A CN 103455679A CN 2013104030393 A CN2013104030393 A CN 2013104030393A CN 201310403039 A CN201310403039 A CN 201310403039A CN 103455679 A CN103455679 A CN 103455679A
- Authority
- CN
- China
- Prior art keywords
- finite element
- model
- layer
- compound substance
- grid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Abstract
The invention relates to a finite element analysis method of a composite material. The method comprises the following steps: 1, building a mathematic model, i.e. building a three-dimensional model by utilizing three-dimensional software; 2, defining material attributes for the built model by utilizing finite element software, performing mesh generation, and converting the model into a finite element model, wherein the material attributes are defined according to the anisotropy of the material, i.e. the winding direction and winding thickness of each layer and the failure criterion of the material are defined; 3, loading boundary conditions on the finite element model, loading conditions on the model according to an actual stress situation, and simulating the actual situation; 4, performing finite element analysis, so as to obtain a stress distribution cloud picture and a displacement situation graph. When the finite element analysis method is used for analysis, the material attributes are defined according to the anisotropy of the material, and calculation is performed respectively according to the winding direction and winding thickness of each layer, so that the winding direction and residual pretightening force of each layer of the composite material can be well simulated.
Description
Technical field
The present invention relates to a kind of finite element method of compound substance.
Background technology
At present, while being calculated for compound substance machinery performance, often utilize traditional mathematical formulae to be calculated, for adopting Wrapping formed compound substance, the winding direction disunity of compound substance and certain pretightning force is arranged when it is wound around, after being wound around, it is layer structure, need to carry out hot setting after being wound around end, in solidification process, the pretightning force of every layer of fiber has release in various degree, the method of traditional mathematics formula is to every layer of winding direction difference, and the compound substance that remaining pretightning force is different is difficult to carry out accurate analysis calculating.
Summary of the invention
The object of the invention is to overcome above defect, a kind of analysis finite element method of compound substance accurately is provided.
Technical scheme of the present invention is: a kind of finite element method of compound substance, and it comprises the following steps:
Step 1: set up mathematical model, according to the size of analytical structure, utilize three-dimensional software to set up three-dimensional model;
Step 2: utilize finite element software to set up model definition material properties and carry out the grid division, convert finite element model to, described definition material properties is defined for the anisotropy according to material, the winding direction of every layer, the failure criteria that is wound around thickness and material are defined, grid is divided can adopt according to the shape of structure different grid dividing mode, can adopt adaptive mesh division, free grid to divide, scan grid division etc.;
Step 3: finite element model is carried out to the boundary condition loading, according to stressed actual conditions, model is carried out to condition loading, simulating actual conditions;
Step 4: above-mentioned model is carried out to finite element analysis, draw stress distribution cloud atlas and displacement condition diagram.
Technique effect of the present invention is, during due to analysis, the definition material properties is for being defined according to the anisotropy of material, by the winding direction of every one deck be wound around thickness and calculate respectively, thereby can well simulate the winding direction of every layer, compound substance, remaining pretightning force.
The accompanying drawing explanation
Fig. 1 is rotor protective sleeve structural representation in the embodiment of the present invention;
Fig. 2 is embodiment of the present invention mathematical model;
Fig. 3 is embodiment of the present invention finite element model.
Embodiment
Below in conjunction with drawings and Examples, the invention will be further described.
With reference to accompanying drawing 1, be a kind of motor rotor construction, it comprises protective sleeve 1, magnet steel 2, rotor core 3, rotating shaft 4.Wherein protective sleeve 1 is compound substance; gross thickness is 2mm; first to four layer is fiber A; winding direction is 4 one-tenth 30 degree angles of rotating shaft, and the thickness of one to four layer is 0.1mm, and the pretightning force of every layer is 500N; layer 5 is fiber B to Floor 12; winding direction is and 4 one-tenth miter angles of rotating shaft, and thickness is 0.2mm, and the pretightning force of every layer is 600N.During the motor operation, protective sleeve 1 can be subject to the centrifugal action of magnet steel 2.To protective sleeve 1, the stressing conditions when the operation of motor carries out sunykatuib analysis in whole finite element analysis, and it comprises the following steps:
Step 1: set up mathematical model; requirement according to analytical structure; utilize three-dimensional software set up three-dimensional model and three-dimensional model is simplified; removal is on the very little structure of result impact; retain the Main Analysis structure; build up mathematical model, this mathematical model comprises rotating shaft 4, rotor core 3, magnet steel 2, protective sleeve 1.Consider the symmetry of structure, choose minimum symmetrical cell and carry out analytical calculation, set up the symmetrical cell model of the minimum shown in Fig. 2;
Step 2: utilize finite element software to set up model definition material properties; the material properties definition of protective sleeve 1 is defined according to anisotropy; protective sleeve 1 winding direction, the failure criteria that is wound around thickness and material of every layer defined, magnet steel 2, rotor core 3, the rotating shaft 4 of isotropic material defined respectively to elastic modulus, Poisson ratio, the density of material itself.Then structure is carried out to the grid division, convert finite element model to, grid is divided can adopt according to the shape of structure different grid dividing mode, as adaptive mesh is divided, free grid is divided, scan grid division etc.Physical dimension is little, easily produces the concentrated position of stress, and grid will guarantee certain precision, to guarantee the accuracy of result of calculation.Set up the finite element model shown in Fig. 3;
Step 3: finite element model is carried out to the boundary condition loading, and what model adopted is minimum symmetrical cell, and therefore symmetrical border will apply symmetrical boundary condition; Rotor core 3 is with magnet steel 2 for normally to contact, and it is right to be arranged to contact; Rotating shaft 4 hypothesis one ends do not have axial displacement, the axial displacement of constraint rotating shaft 4 one end end faces; Protective sleeve 1 is compound substance; certain pretightning force F is arranged during colligation; and with magnet steel 2 close contacts; therefore protective sleeve 1 and magnet steel 2 be arranged to contact right; and four layers of protective sleeve 1 ground floors to the are applied to pretightning force F=500N; layer 5 applies pretightning force F=600N to Floor 12, and by order, ISTRESS applies.Simulated machine operation conditions during calculating; the centrifugal force that magnet steel 2 is applied on protective sleeve 1 need to be realized by block mold being applied to rotating speed; rotating speed is that angular velocity n is rotating speed by formula α=n π/30(α) rotating speed is converted to angular velocity, block mold is applied to angular velocity.It is mainly to simulate actual condition fully that whole model loads, and the centrifugal force of rotating shaft 4 during operation, rotor core 3,2 pairs of protective sleeves 1 of magnet steel calculates and is applied on protective sleeve by software;
Step 4: above-mentioned model is carried out to finite element analysis, draw stress distribution cloud atlas and displacement condition diagram.
Claims (4)
1. the finite element method of a compound substance, is characterized in that, it comprises the following steps:
Step 1: set up mathematical model, according to the size of analytical structure, utilize three-dimensional software to set up three-dimensional model;
Step 2: utilize finite element software to set up model definition material properties and carry out the grid division, convert finite element model to, described definition material properties is for being defined according to the anisotropy of material, the winding direction of every layer, the failure criteria that is wound around thickness and material defined;
Step 3: finite element model is carried out to the boundary condition loading, according to stressed actual conditions, model is carried out to condition loading, simulating actual conditions;
Step 4: above-mentioned model is carried out to finite element analysis, draw stress distribution cloud atlas and displacement condition diagram.
2. the finite element method of a kind of compound substance according to claim 1, is characterized in that, described grid is divided into adaptive mesh to be divided.
3. the finite element method of a kind of compound substance according to claim 1, is characterized in that, described grid is divided into free grid and divides.
4. the finite element method of a kind of compound substance according to claim 1, is characterized in that, described grid is divided into and scans the grid division.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2013104030393A CN103455679A (en) | 2013-09-06 | 2013-09-06 | Finite element analysis method of composite material |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2013104030393A CN103455679A (en) | 2013-09-06 | 2013-09-06 | Finite element analysis method of composite material |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103455679A true CN103455679A (en) | 2013-12-18 |
Family
ID=49738035
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2013104030393A Pending CN103455679A (en) | 2013-09-06 | 2013-09-06 | Finite element analysis method of composite material |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103455679A (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104679943A (en) * | 2015-01-29 | 2015-06-03 | 大连理工大学 | Method for simulating fiber reinforced composite material cutting chip forming |
| CN105205223A (en) * | 2015-08-27 | 2015-12-30 | 湘潭大学 | Finite element modeling method for establishing defective material model |
| CN106096094A (en) * | 2016-05-31 | 2016-11-09 | 北方工业大学 | Method for eliminating anisotropic material cylindrical part drawing lug based on finite element simulation |
| CN106202598A (en) * | 2015-05-07 | 2016-12-07 | 哈尔滨飞机工业集团有限责任公司 | The analysis method of residual compressive strength after the damage of a kind of composite impact |
| CN106354965A (en) * | 2016-09-08 | 2017-01-25 | 中国航空工业集团公司西安飞机设计研究所 | Finite element analysis method of residual strength of laminated board containing layered composite material |
| CN106599364A (en) * | 2016-11-11 | 2017-04-26 | 长春理工大学 | Finite element fluid solid coupling modeling and mode analyzing method of droplet penetrating micro cavity |
| CN106980711A (en) * | 2017-03-06 | 2017-07-25 | 中国电子产品可靠性与环境试验研究所 | The air-tightness analysis method of level Hermetic Package glass insulator based on finite element simulation |
| CN107491627A (en) * | 2017-10-16 | 2017-12-19 | 南京航空航天大学 | Composite drilling process finite element Global sensitivity analysis method |
| CN108090242A (en) * | 2017-08-31 | 2018-05-29 | 镇江春环密封件集团有限公司 | The three-dimensional thermal-structural coupling analysis method of carbon fiber winding composite cylinder |
| CN108491612A (en) * | 2018-03-15 | 2018-09-04 | 东南大学 | The Finite Element Method of scheme of material selection is provided for multiple tube hydraulic bulging process |
| CN111563337A (en) * | 2020-04-03 | 2020-08-21 | 中国航发哈尔滨东安发动机有限公司 | Finite element analysis method for strength of shaft parts |
| CN112417724A (en) * | 2020-11-19 | 2021-02-26 | 贺州学院 | Calcium carbonate filled composite material mechanical property prediction method based on reverse modeling |
| CN112632659A (en) * | 2020-12-08 | 2021-04-09 | 株洲时代新材料科技股份有限公司 | Analysis method for GINA waterstop for immersed tunnel |
| CN113076674A (en) * | 2021-04-09 | 2021-07-06 | 华南理工大学 | Structure modeling and manufacturing method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030154451A1 (en) * | 2002-02-12 | 2003-08-14 | The Boeing Company | Method, system and computer program product for multidisciplinary design analysis of structural components |
| CN101906895A (en) * | 2010-07-23 | 2010-12-08 | 北京玻钢院复合材料有限公司 | Composite electrical cross arm and design method of structure thereof |
| CN102339348A (en) * | 2011-08-11 | 2012-02-01 | 西北工业大学 | Method for detecting micro stress of unidirectional C/SiC (continuous carbon fiber reinforced silicon carbide composites) in oxidation environment |
| CN102517760A (en) * | 2011-12-31 | 2012-06-27 | 机械科学研究总院先进制造技术研究中心 | Laminated weaving formation method for fabricated part made of composite material |
| CN102622472A (en) * | 2012-02-27 | 2012-08-01 | 西北工业大学 | Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins |
-
2013
- 2013-09-06 CN CN2013104030393A patent/CN103455679A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030154451A1 (en) * | 2002-02-12 | 2003-08-14 | The Boeing Company | Method, system and computer program product for multidisciplinary design analysis of structural components |
| CN101906895A (en) * | 2010-07-23 | 2010-12-08 | 北京玻钢院复合材料有限公司 | Composite electrical cross arm and design method of structure thereof |
| CN102339348A (en) * | 2011-08-11 | 2012-02-01 | 西北工业大学 | Method for detecting micro stress of unidirectional C/SiC (continuous carbon fiber reinforced silicon carbide composites) in oxidation environment |
| CN102517760A (en) * | 2011-12-31 | 2012-06-27 | 机械科学研究总院先进制造技术研究中心 | Laminated weaving formation method for fabricated part made of composite material |
| CN102622472A (en) * | 2012-02-27 | 2012-08-01 | 西北工业大学 | Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins |
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104679943B (en) * | 2015-01-29 | 2018-01-26 | 大连理工大学 | A kind of emulation mode cut fibre reinforced composites chip and formed |
| CN104679943A (en) * | 2015-01-29 | 2015-06-03 | 大连理工大学 | Method for simulating fiber reinforced composite material cutting chip forming |
| CN106202598A (en) * | 2015-05-07 | 2016-12-07 | 哈尔滨飞机工业集团有限责任公司 | The analysis method of residual compressive strength after the damage of a kind of composite impact |
| CN106202598B (en) * | 2015-05-07 | 2019-11-15 | 哈尔滨飞机工业集团有限责任公司 | The analysis method of residual compressive strength after a kind of damage of composite impact |
| CN105205223A (en) * | 2015-08-27 | 2015-12-30 | 湘潭大学 | Finite element modeling method for establishing defective material model |
| CN105205223B (en) * | 2015-08-27 | 2018-08-24 | 湘潭大学 | A kind of finite element modeling method for establishing the material model containing defect |
| CN106096094A (en) * | 2016-05-31 | 2016-11-09 | 北方工业大学 | Method for eliminating anisotropic material cylindrical part drawing lug based on finite element simulation |
| CN106354965A (en) * | 2016-09-08 | 2017-01-25 | 中国航空工业集团公司西安飞机设计研究所 | Finite element analysis method of residual strength of laminated board containing layered composite material |
| CN106599364A (en) * | 2016-11-11 | 2017-04-26 | 长春理工大学 | Finite element fluid solid coupling modeling and mode analyzing method of droplet penetrating micro cavity |
| CN106599364B (en) * | 2016-11-11 | 2020-02-21 | 长春理工大学 | Finite element fluid-solid coupling modeling and modal analysis method for micro-cavity liquid drop |
| CN106980711A (en) * | 2017-03-06 | 2017-07-25 | 中国电子产品可靠性与环境试验研究所 | The air-tightness analysis method of level Hermetic Package glass insulator based on finite element simulation |
| CN106980711B (en) * | 2017-03-06 | 2021-02-23 | 中国电子产品可靠性与环境试验研究所 | Finite element simulation-based airtight packaging glass insulator airtightness analysis method |
| CN108090242A (en) * | 2017-08-31 | 2018-05-29 | 镇江春环密封件集团有限公司 | The three-dimensional thermal-structural coupling analysis method of carbon fiber winding composite cylinder |
| CN107491627A (en) * | 2017-10-16 | 2017-12-19 | 南京航空航天大学 | Composite drilling process finite element Global sensitivity analysis method |
| CN107491627B (en) * | 2017-10-16 | 2020-10-09 | 南京航空航天大学 | Finite element global sensitivity analysis method in composite material drilling process |
| CN108491612A (en) * | 2018-03-15 | 2018-09-04 | 东南大学 | The Finite Element Method of scheme of material selection is provided for multiple tube hydraulic bulging process |
| CN111563337A (en) * | 2020-04-03 | 2020-08-21 | 中国航发哈尔滨东安发动机有限公司 | Finite element analysis method for strength of shaft parts |
| CN112417724A (en) * | 2020-11-19 | 2021-02-26 | 贺州学院 | Calcium carbonate filled composite material mechanical property prediction method based on reverse modeling |
| CN112632659A (en) * | 2020-12-08 | 2021-04-09 | 株洲时代新材料科技股份有限公司 | Analysis method for GINA waterstop for immersed tunnel |
| CN113076674A (en) * | 2021-04-09 | 2021-07-06 | 华南理工大学 | Structure modeling and manufacturing method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103455679A (en) | Finite element analysis method of composite material | |
| CN103501088B (en) | A kind of motor lamination core vibration characteristic test system and vibration characteristics preparation method | |
| US9524356B2 (en) | System and methods of generating a computer model of composite component | |
| Rezaei et al. | Modal-based damage identification for the nonlinear model of modern wind turbine blade | |
| WO2012175406A3 (en) | Method for producing a tower segment of a concrete tower of an wind energy installation | |
| CN102401727A (en) | Method for obtaining mechanical joint stiffness | |
| CN104149340A (en) | Method for obtaining properties of metal part based on 3D printing technique | |
| CN105069249A (en) | Method for optimizing cutting path of spacer frame type integral structural member | |
| CN103991556B (en) | A kind of load measurement method | |
| CN104978466A (en) | Three-dimensional definite-number bar arrangement method used for multiple guide wires | |
| CN104573178A (en) | Finite element method for calculating integrated impeller strength | |
| US20150370923A1 (en) | System and Methods of Generating a Computer Model of a Composite Component | |
| CN104089563A (en) | Measuring method for spherical radial clearance of large radial spherical plain bearing | |
| JP2019500665A (en) | How to design a plybook for complex components | |
| CN106442407B (en) | Reliability verification method for additive manufacturing product | |
| CN104268320B (en) | Suitable for the New Type Vibration Isolator Using stiffness matrix rapid Estimation method of satellite sensitive load | |
| CN112307575B (en) | Method for checking qualification of composite material wallboard part for judging assembly result | |
| CN105022890A (en) | Method for backward inferring Young modulus of loudspeaker through resonant frequency of vibrating part of loudspeaker | |
| JP5994269B2 (en) | Motor torque setting method and setting device | |
| US10099433B2 (en) | Datum feature for a composite component | |
| Montano et al. | Bend-twist coupling on rotor blades for wind turbines | |
| CN103008565B (en) | Method for manufacturing sand mould of propeller shaft bracket | |
| CN110569611B (en) | Structural frequency response function sensitivity analysis method based on multi-complex-variable method | |
| CN102860633B (en) | Digital shoe tree girth measurement method based on extended geodesic curve energy optimization | |
| CN104389717B (en) | A kind of method determining Hydropower Unit worm frame compression clearance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20131218 |