CN116227155B - Method for researching microscopic mechanical property of propellant - Google Patents
Method for researching microscopic mechanical property of propellant Download PDFInfo
- Publication number
- CN116227155B CN116227155B CN202310032033.3A CN202310032033A CN116227155B CN 116227155 B CN116227155 B CN 116227155B CN 202310032033 A CN202310032033 A CN 202310032033A CN 116227155 B CN116227155 B CN 116227155B
- Authority
- CN
- China
- Prior art keywords
- propellant
- microscopic
- model
- rigidity
- polygonal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Micromachines (AREA)
Abstract
The application discloses a method for researching the micro-mechanical property of a propellant, which belongs to the technical field of micro-mechanics and comprises the following steps: constructing a propellant microscopic particle filling model by using a molecular dynamics method; polygonal meshing is carried out on the propellant microscopic particle filling model by utilizing a meshing tool, and a propellant microscopic polygonal mesh model is obtained; designing a calculation scheme of the rigidity of the elastic particles of the propellant and the rigidity of the viscoelastic matrix; based on the propellant microscopic polygonal grid model, the propellant is subjected to stress relaxation simulation test by combining the propellant elastic particle rigidity and the viscoelastic matrix rigidity calculation scheme, so that the propellant microscopic mechanical property research is realized. The method can effectively cope with the complex microscopic geometric model of the propellant, and realizes the research on the microscopic mechanical properties of the propellant.
Description
Technical Field
The application relates to a method for researching the microscopic mechanical properties of a propellant, and belongs to the technical field of microscopic mechanics.
Background
The development direction of high loading, long service life and large range of solid missiles puts higher demands on the mechanical properties of solid propellants.
Currently, the development of solid propellant is still a semi-empirical semi-theoretical research method, and the development of the mechanical properties of propellant is difficult to develop. On the other hand, the prophetic mechanical property prophetic method of the propellant provides a research thought for the component design of the propellant. However, the conventional finite element method takes a lot of time to perform the meshing in terms of processing the prophetic mechanical property predictions of the propellant, and it is difficult to obtain a good meshing effect.
Disclosure of Invention
The application aims to provide a method for researching the micro-mechanical properties of a propellant, which can efficiently cope with a complex micro-geometric model of the propellant and realize the research on the micro-mechanical properties of the propellant.
In order to achieve the above purpose, the present application provides the following technical solutions:
a method for studying the microscopic mechanical properties of a propellant, comprising:
constructing a propellant microscopic particle filling model by using a molecular dynamics method;
polygonal meshing is carried out on the propellant microscopic particle filling model by utilizing a meshing tool, and a propellant microscopic polygonal mesh model is obtained;
calculating the rigidity of the elastic particles of the propellant and the rigidity of the viscoelastic matrix;
based on the propellant microscopic polygonal grid model, the propellant elastic particle rigidity and the viscoelasticity matrix rigidity are combined to perform stress relaxation simulation test on the propellant, so that the propellant microscopic mechanical property research is realized.
Further, using molecular dynamics methods, constructing a propellant fine particle filling model includes:
determining the number N of propellant grains according to the volume percentage of the propellant;
scattering N particles with initial radius of 0 in the coordinate range of the microscopic size of the propellant;
each of the particles is grown at a predetermined rate until its radius increases to a design value to construct a propellant fine particle filling model.
Further, constructing a space rectangular coordinate system by taking the center of the propellant microscopic particle filling model as an origin, wherein the dimensionless coordinate range of each particle is [ -1,1]Let the boundary vertex coordinates of the ith grain be eta i If eta i >1, let eta i =η i -2, if eta i <-1, let eta i =η i +2。
Further, performing polygonal meshing on the propellant microscopic particle filling model by using a meshing tool, wherein obtaining the propellant microscopic polygonal meshing model comprises:
stretching the propellant microscopic particle filling model from two dimensions to three dimensions to obtain a three-dimensional geometric model;
performing polygon meshing on the three-dimensional geometric model by adopting a ANASYS fluent meshing meshing tool;
and extracting surface information of the polyhedron from the three-dimensional geometric model divided by the polygonal mesh to obtain the propellant microscopic polygonal mesh model.
Further, calculating the propellant elastane particle stiffness and the viscoelastic matrix stiffness comprises:
defining a mapping operator from the virtual cell space to the cell polynomial space;
acquiring an expression of a basis function and a unit stiffness matrix by using the mapping operator;
and calculating the rigidity of the elastic particles of the propellant and the rigidity of the visco-elastic matrix according to the basic function and the expression of the unit rigidity matrix.
Further, the mapping operator has an orthogonal relationship as shown in formula (1) for each function in the virtual cell space:
in the formula (1), P k (E) Is a unit polynomial space, p k Is a unit polynomial space P k (E) Is selected from the group consisting of,to map operator, v h Is a virtual unit space V k (E) E is a cell region.
Further, the expression of the basis function is shown in formula (2):
in the formula (2),as a basis function, i is the ordinal number of the basis function;
the expression of the rigidity matrix is shown in formula (3):
in the formula (3), K E Is a rigidity matrix, I is an identity matrix,g is the second auxiliary matrix.
Further, performing a stress relaxation simulation test on the propellant includes:
restraining the longitudinal displacement and left lateral displacement of the bottom of the propellant microscopic polygonal grid model, applying 5% of longitudinal step strain on the top of the propellant microscopic polygonal grid model, and keeping for 1000s;
calculating the average stress and average strain of the propellant microscopic polygonal grid model according to the rigidity of the propellant elastic particles and the rigidity of the visco-elastic matrix;
calculating the loose modulus of the propellant under different representative volume unit sizes according to the average stress and the average strain of the propellant microscopic polygonal grid model;
reducing the size of the representative volume element until the relative error of the relaxed modulus exceeds 5%;
and according to each relaxation modulus, realizing the research on the microscopic mechanical properties of the propellant.
Further, the average stress of the propellant microscopic polygonal mesh model is the average stress of each node of the propellant microscopic polygonal mesh model, and the average strain of the propellant microscopic polygonal mesh model is the average strain of each node of the propellant microscopic polygonal mesh model.
Further, the representative volume units are micro structures extracted from the composite material as a whole for studying the performance of the composite material.
Compared with the prior art, the application has the beneficial effects that:
compared with the traditional finite element method, the method can effectively cope with the complex microscopic geometric model of the propellant. Under the condition that the total number of degrees of freedom is the same, the application can construct more entity units for propellant microscopic geometric modeling.
Drawings
FIG. 1 is a flow chart of a method for researching the microscopic mechanical properties of a propellant, which is provided by the embodiment of the application;
FIG. 2 is a schematic illustration of a particle according to an embodiment of the present application when a "wire-pressing" event occurs;
FIG. 3 is a schematic illustration of a propellant fine particle filling model provided by an embodiment of the present application;
FIG. 4 is a schematic view of a three-dimensional geometric model provided by an embodiment of the present application;
FIG. 5 is a schematic diagram of a three-dimensional geometric model polygon meshing provided by an embodiment of the present application;
FIG. 6 is a schematic diagram of a propellant microscopic polygon mesh model provided by an embodiment of the present application;
fig. 7 is an axial stress cloud of a propellant mesostructure according to an embodiment of the present application relaxed to a time of 1000 s.
Detailed Description
The technical scheme of the patent is further described in detail below with reference to the specific embodiments.
Embodiments of the present patent are described in detail below, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to like or similar elements or elements having like or similar functions throughout. The embodiments described below by referring to the drawings are exemplary only for explaining the present patent and are not to be construed as limiting the present patent. The embodiments of the present application and technical features in the embodiments may be combined with each other without collision.
Embodiment one:
fig. 1 is a flowchart of a method for researching the micro-mechanical properties of a propellant, which is provided in the first embodiment of the present application, and the flowchart merely shows the logic sequence of the method in this embodiment, and on the premise of not conflicting with each other, in other possible embodiments of the present application, the steps shown or described can be completed in a different sequence from that shown in fig. 1.
Referring to fig. 1, the method of this embodiment specifically includes the following steps:
step one: constructing a propellant microscopic particle filling model by using a molecular dynamics method;
the construction of the propellant microscopic particle filling model by using a molecular dynamics method comprises the following steps:
step A: determining the number N of propellant grains according to the volume percentage of the propellant;
and (B) step (B): scattering N particles with initial radius of 0 in the coordinate range of the microscopic size of the propellant;
the coordinates of the microscopic dimensions of the propellant range from: in the two-dimensional plane, point a (0, 0), point B (L, 0), point C (L, L), and point D (0, L) are four-point lines, where L is the propellant mesoscopic size.
Step C: growing each particle at a predetermined rate until its radius increases to a design value to form a propellant fine particle filling model;
during the radius becoming larger, when particles collide, the particles spring in the opposite direction, and so on, until the radius increases to the design value.
In order to describe the mechanical response of a true mesostructure, the propellant mesogranule filling model must meet the requirements of the periodic geometry and the periodic boundary conditions. For the periodic geometry, if a particle "line-up" event occurs during particle loading, as shown in fig. 2, there is a corresponding particle on the corresponding geometric side (or face in the three-dimensional model) to satisfy the constraint of the periodic geometry. Furthermore, the number of particles and the corresponding geometry are fixed values calculated according to the formulation of the material and cannot be changed by the filling process before the particle filling is carried out.
In order to solve the problems of periodic geometry and periodic boundary conditions, a space rectangular coordinate system is constructed by taking the center of a propellant microscopic particle filling model as an origin, and the dimensionless coordinate range of each particle is [ -1,1]Let the boundary vertex coordinates of the ith grain be eta i If eta i >1, let eta i =η i -2, if eta i <-1, let eta i =η i +2。
By this method, it is possible to have all the particles "pressed" symmetrically on the corresponding geometric border of the model, while the number of particles is increased compared to the number initially calculated, without affecting the total packing fraction. As shown in fig. 3, a propellant fine particle filling model is generated by filling propellant particles by a molecular dynamics method, in fig. 3, circles are particles, and the rest is a matrix.
Step two: polygonal mesh division is carried out on the propellant microscopic particle filling model by utilizing a mesh division tool, and a propellant microscopic polygonal mesh model is obtained;
polygonal meshing is carried out on the propellant microscopic particle filling model by utilizing a meshing tool, and the propellant microscopic polygonal meshing model is obtained, wherein the method comprises the following steps of:
step a: stretching the propellant microscopic particle filling model from two dimensions to three dimensions to obtain a three-dimensional geometric model;
step b: adopting a ANASYS fluent meshing mesh dividing tool to carry out polygonal mesh division on the three-dimensional geometric model;
step c: extracting surface information of a polyhedron from the three-dimensional geometric model divided by the polygonal mesh to obtain a propellant microscopic polygonal mesh model;
at present, polygon meshing is relatively more software, but is a ANASYS fluent meshing meshing tool relatively mature. Since the meshing tool now only supports meshing of three-dimensional geometry, the propellant microscopic particle filling model in fig. 3 needs to be stretched to obtain the three-dimensional geometry model shown in fig. 4, then the obtained three-dimensional geometry model is imported into the ANASYS fluent meshing meshing tool, the three-dimensional geometry model is subjected to polygonal meshing by the ANASYS fluent meshing meshing tool, finally the propellant microscopic polygonal meshing model is obtained by extracting polyhedral surface information, and the obtained propellant microscopic polygonal meshing model is shown in fig. 6.
Step three: calculating the rigidity of the elastic particles of the propellant and the rigidity of the viscoelastic matrix;
the calculation of the propellant elastane particle stiffness and viscoelastic matrix stiffness comprises the steps of:
step I: defining a slave virtual cell space V k (E) To the unit polynomial space P k (E) Mapping operator of (a)
Step II: using mapping operatorsAcquiring a basis function->And cell stiffness matrix K E Is an expression of (2);
step III: root of Chinese characterAccording to the basis functionAnd cell stiffness matrix K E Calculating the propellant elastane particle stiffness and viscoelastic matrix stiffness;
wherein the mapping operatorFor the virtual cell space V k (E) Is a function v h There is an orthogonal relationship as shown in equation (1):
in the formula (1), P k (E) Is a unit polynomial space, p k Is a unit polynomial space P k (E) Is selected from the group consisting of,to map operator, v h Is a virtual unit space V k (E) E is a cell region.
The expression of the basis function is shown in formula (2):
in the formula (2),as a basis function, i is the ordinal number of the basis function;
the expression of the stiffness matrix is shown in formula (3):
in the formula (3), K E Is a rigidity matrix, I is an identity matrix,g is the second auxiliary matrix.
Substituting the elasticity and the viscoelasticity of the propellant into the formula (3), and obtaining the rigidity of the elastic particles and the rigidity of the viscoelasticity matrix of the propellant through calculation.
Step four: based on a propellant microscopic polygonal grid model, combining the rigidity of elastic particles of the propellant and the rigidity of a viscoelastic matrix, performing stress relaxation simulation test on the propellant, and realizing the research on the microscopic mechanical properties of the propellant;
based on a propellant microscopic polygonal grid model, combining the rigidity of propellant elastic particles and the rigidity of a visco-elastic matrix, carrying out stress relaxation simulation test on the propellant, and realizing the research on the microscopic mechanical properties of the propellant, wherein the method comprises the following steps:
step i: restraining longitudinal displacement of the bottom and left lateral displacement of the propellant microscopic polygonal grid model, applying 5% of longitudinal step strain on the top of the propellant microscopic polygonal grid model, and keeping the top of the propellant microscopic polygonal grid model for 1000s;
step ii: calculating average stress and average strain of the propellant microscopic polygonal grid model according to the rigidity of the propellant elastic particles and the rigidity of the visco-elastic matrix;
step iii: calculating the loose modulus of the propellant under different representative volume unit sizes according to the average stress and the average strain of the propellant microscopic polygonal grid model;
step iv: decreasing the size of the representative volume unit until the relative error in the relaxed modulus exceeds 5%;
step v: according to each relaxation modulus, the research on the microscopic mechanical properties of the propellant is realized;
the average stress of the propellant microscopic polygonal mesh model is the average stress of each node of the propellant microscopic polygonal mesh model, and the average strain of the propellant microscopic polygonal mesh model is the average strain of each node of the propellant microscopic polygonal mesh model.
The representative volume units are micro structures extracted from the composite material as a whole for studying the properties of the composite material. The loose moduli of the propellant under different sizes of the representative volume units are calculated respectively, the sizes of the representative volume units are reduced until the relative error of the loose moduli exceeds 5%, and the loose moduli are the loose moduli of the propellant under the components, and as shown in fig. 7, are axial stress cloud diagrams of the mesoscopic structure of the propellant when the loose to 1000 s.
Compared with the traditional finite element method, the method for researching the micro-mechanical properties of the propellant can efficiently cope with complex micro-geometric models of the propellant and realize the research on the micro-mechanical properties of the propellant. Under the condition that the total number of degrees of freedom is the same, the application can construct more entity units for propellant microscopic geometric modeling.
The foregoing is merely a preferred embodiment of the present application, and it should be noted that modifications and variations could be made by those skilled in the art without departing from the technical principles of the present application, and such modifications and variations should also be regarded as being within the scope of the application.
Claims (6)
1. A method for studying the microscopic mechanical properties of a propellant, comprising:
constructing a propellant microscopic particle filling model by using a molecular dynamics method;
polygonal meshing is carried out on the propellant microscopic particle filling model by utilizing a meshing tool, and a propellant microscopic polygonal mesh model is obtained;
calculating the rigidity of the elastic particles of the propellant and the rigidity of the viscoelastic matrix;
based on the propellant microscopic polygonal grid model, combining the propellant elastic particle rigidity and the viscoelasticity matrix rigidity, performing stress relaxation simulation test on the propellant to realize propellant microscopic mechanical property research;
calculating the propellant elastane particle stiffness and viscoelastic matrix stiffness includes:
defining a mapping operator from the virtual cell space to the cell polynomial space;
acquiring an expression of a basis function and a unit stiffness matrix by using the mapping operator;
calculating the rigidity of elastic particles of the propellant and the rigidity of the viscoelastic matrix according to the basic function and the expression of the unit rigidity matrix;
the mapping operator has an orthogonal relationship for each function in the virtual cell space as shown in equation (1):
in the formula (1), P k (E) Is a unit polynomial space, p k Is a unit polynomial space P k (E) Is selected from the group consisting of,to map operator, v h Is a virtual unit space V k (E) E is a unit area;
the expression of the basis function is shown in formula (2):
in the formula (2),as a basis function, i is the ordinal number of the basis function;
the expression of the rigidity matrix is shown in formula (3):
in the formula (3), K E Is a rigidity matrix, I is an identity matrix,g is the second auxiliary matrix.
2. The method of claim 1, wherein constructing a model of propellant microscopic particle packing using molecular dynamics methods comprises:
determining the number N of propellant grains according to the volume percentage of the propellant;
scattering N particles with initial radius of 0 in the coordinate range of the microscopic size of the propellant;
each of the particles is grown at a predetermined rate until its radius increases to a design value to construct a propellant fine particle filling model.
3. The method for studying the microscopic mechanical properties of a propellant according to claim 2, wherein a space rectangular coordinate system is constructed with the center of the propellant microscopic particle filling model as an origin, and the dimensionless coordinate range of each particle is [ -1,1]Let the boundary vertex coordinates of the ith grain be eta i If eta i >1, let eta i =η i -2, if eta i <-1, let eta i =η i +2。
4. The method of claim 1, wherein performing polygonal meshing of the propellant microscopic particle filling model with a meshing tool to obtain a propellant microscopic polygonal mesh model comprises:
stretching the propellant microscopic particle filling model from two dimensions to three dimensions to obtain a three-dimensional geometric model;
performing polygon meshing on the three-dimensional geometric model by adopting a ANASYS fluent meshing meshing tool;
and extracting surface information of the polyhedron from the three-dimensional geometric model divided by the polygonal mesh to obtain the propellant microscopic polygonal mesh model.
5. The method of claim 1, wherein performing a stress relaxation simulation test on the propellant based on the propellant microscopic polygonal mesh model in combination with the propellant elastic particle stiffness and the viscoelastic matrix stiffness, the performing the propellant microscopic mechanical property study comprises:
restraining the longitudinal displacement and left lateral displacement of the bottom of the propellant microscopic polygonal grid model, applying 5% of longitudinal step strain on the top of the propellant microscopic polygonal grid model, and keeping for 1000s;
calculating the average stress and average strain of the propellant microscopic polygonal grid model according to the rigidity of the propellant elastic particles and the rigidity of the visco-elastic matrix;
calculating the loose modulus of the propellant under different representative volume unit sizes according to the average stress and the average strain of the propellant microscopic polygonal grid model;
reducing the size of the representative volume element until the relative error of the relaxed modulus exceeds 5%;
according to each relaxation modulus, realizing the research on the microscopic mechanical properties of the propellant;
the representative volume units are micro structures extracted from the composite material as a whole for studying the performance of the composite material.
6. The method for studying the microscopic mechanical properties of a propellant according to claim 5, wherein the average stress of the microscopic polygonal mesh model of the propellant is an average value of the stress of each node of the microscopic polygonal mesh model of the propellant, and the average strain of the microscopic polygonal mesh model of the propellant is an average value of the strain of each node of the microscopic polygonal mesh model of the propellant.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310032033.3A CN116227155B (en) | 2023-01-10 | 2023-01-10 | Method for researching microscopic mechanical property of propellant |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310032033.3A CN116227155B (en) | 2023-01-10 | 2023-01-10 | Method for researching microscopic mechanical property of propellant |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116227155A CN116227155A (en) | 2023-06-06 |
CN116227155B true CN116227155B (en) | 2023-09-22 |
Family
ID=86579805
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310032033.3A Active CN116227155B (en) | 2023-01-10 | 2023-01-10 | Method for researching microscopic mechanical property of propellant |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116227155B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116913432A (en) * | 2023-07-25 | 2023-10-20 | 中国人民解放军陆军工程大学 | Reconstruction method of propellant simulated filling model based on polygonal grid |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103267679A (en) * | 2013-05-09 | 2013-08-28 | 上海交通大学 | Approximate model technology based composite foamed plastic interface phase mechanical test method |
JP2016184245A (en) * | 2015-03-25 | 2016-10-20 | トヨタ自動車株式会社 | Particle model and data transfer method between mesh models |
CN107391801A (en) * | 2017-06-23 | 2017-11-24 | 中国人民解放军国防科学技术大学 | The structure and methods for using them of propellant Dynamic Thermal viscoelastic Constitutive Model containing dehumidification |
CN108491639A (en) * | 2018-03-26 | 2018-09-04 | 中国石油大学(华东) | Closure shale pore throat simulation method based on nanoparticle impact contact model |
US10460007B1 (en) * | 2015-09-03 | 2019-10-29 | Ansys, Inc. | Systems and methods for solving integral equations |
CN114417659A (en) * | 2021-12-30 | 2022-04-29 | 清华大学 | Rock-fill concrete microscopic three-phase finite element model, and establishment method and device thereof |
-
2023
- 2023-01-10 CN CN202310032033.3A patent/CN116227155B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103267679A (en) * | 2013-05-09 | 2013-08-28 | 上海交通大学 | Approximate model technology based composite foamed plastic interface phase mechanical test method |
JP2016184245A (en) * | 2015-03-25 | 2016-10-20 | トヨタ自動車株式会社 | Particle model and data transfer method between mesh models |
US10460007B1 (en) * | 2015-09-03 | 2019-10-29 | Ansys, Inc. | Systems and methods for solving integral equations |
CN107391801A (en) * | 2017-06-23 | 2017-11-24 | 中国人民解放军国防科学技术大学 | The structure and methods for using them of propellant Dynamic Thermal viscoelastic Constitutive Model containing dehumidification |
CN108491639A (en) * | 2018-03-26 | 2018-09-04 | 中国石油大学(华东) | Closure shale pore throat simulation method based on nanoparticle impact contact model |
CN114417659A (en) * | 2021-12-30 | 2022-04-29 | 清华大学 | Rock-fill concrete microscopic three-phase finite element model, and establishment method and device thereof |
Non-Patent Citations (7)
Title |
---|
"Evaluation of cross-ply laminate stiffness with a non-uniform distribution of transverse matrix cracks by means of a computational meso-mechanic model";Barulich, ND 等;《Composite Structures》;第185卷;第561-572页 * |
"Investigation on Aging Model of Solid Propellant Using the Degree of Crosslinking";Wang, Junqi 等;《International Journal of Aerospace Engineering》;第2022卷;第1--8页 * |
"Numerical Simulation of Three-dimensional Combustion Flows of Pasty Propellant Rocket Motor";Renjie Hu 等;《2022 13th International Conference on Mechanical and Aerospace Engineering (ICMAE)》;第290-295页 * |
"围压下固体推进剂的破坏机理分析";张亮 等;《国防科技大学学报》;第44卷(第6期);第6-13页 * |
"固体推进剂粘弹性泊松比试验与细观预示方法研究";吕轩;《中国优秀硕士学位论文全文数据库(工程科技Ⅰ辑)》;第B017-37页 * |
"复合固体推进剂细观损伤机理及本构模型研究";韩龙;《中国博士学位论文全文数据库 (工程科技Ⅱ辑)》;第C031-17页 * |
"细观结构对HTPB推进剂力学性能影响的数值研究";封涛 等;《功能材料》;第49卷(第1期);第01078-01082页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116227155A (en) | 2023-06-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116227155B (en) | Method for researching microscopic mechanical property of propellant | |
Poole et al. | High-fidelity aerodynamic shape optimization using efficient orthogonal modal design variables with a constrained global optimizer | |
JP5255714B2 (en) | Three-dimensional fluid simulation method | |
CN102609990B (en) | Massive-scene gradually-updating algorithm facing complex three dimensional CAD (Computer-Aided Design) model | |
Ding et al. | An efficient dynamic mesh generation method for complex multi-block structured grid | |
EP2202659A2 (en) | Method of simulating impact events in a multi-processor computer system | |
CN108986220A (en) | A method of accelerating finite element solving object lattice model flexible deformation | |
CN112613243A (en) | Method and device for fluid mechanics simulation and computer readable storage medium | |
CN103886120A (en) | Lightweighting visualization method for digital mockups of large products | |
CN109284529A (en) | A kind of finite element structural analysis method based on VTK | |
CN109960865B9 (en) | GPU acceleration method for dynamic response analysis of thin plate grid-free Galerkin structure | |
Muhamad et al. | Design improvement of steering knuckle component using shape optimization | |
CN113178011B (en) | Cut grid THINC method for solving VOF convection equation | |
CN1828672A (en) | Geometrical data subdividing method based on hexahedron interpolation subdivision | |
CN107145630B (en) | Plate shell structure design and analysis integrated method based on CAD (computer-aided design) cutting curved surface | |
EP2500868B1 (en) | Method for creating finite element model of rubber composite | |
CN109815518A (en) | Vehicle rudder design method based on rotary inertia constraint | |
Barral et al. | Two mesh deformation methods coupled with a changing-connectivity moving mesh method for CFD applications | |
Williams et al. | An accelerated iterative linear solver with GPUs for CFD calculations of unstructured grids | |
US20140129183A1 (en) | Method for producing finite element model | |
CN109948253B (en) | GPU acceleration method for thin-plate meshless Galerkin structure modal analysis | |
CN1680943A (en) | Distributed CAD apparatus | |
CN1822051A (en) | Three dimension net figure simplifying method based on quadrangle folding | |
CN1702656A (en) | Plane layout planning method based on weights | |
Haghi et al. | Mechanical models and finite-element approaches for the structural analysis of photovoltaic composite structures: a comparative study |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |