CN115270479A - Theoretical calculation method for tensile property of composite propellant - Google Patents
Theoretical calculation method for tensile property of composite propellant Download PDFInfo
- Publication number
- CN115270479A CN115270479A CN202210922297.1A CN202210922297A CN115270479A CN 115270479 A CN115270479 A CN 115270479A CN 202210922297 A CN202210922297 A CN 202210922297A CN 115270479 A CN115270479 A CN 115270479A
- Authority
- CN
- China
- Prior art keywords
- epsilon
- tensile strength
- phase
- calculation method
- maximum
- 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
- 239000002131 composite material Substances 0.000 title claims abstract description 38
- 238000004364 calculation method Methods 0.000 title claims abstract description 34
- 239000003380 propellant Substances 0.000 title claims abstract description 17
- 239000002245 particle Substances 0.000 claims abstract description 44
- 239000011159 matrix material Substances 0.000 claims abstract description 22
- 238000002474 experimental method Methods 0.000 claims abstract description 4
- 239000000853 adhesive Substances 0.000 claims description 16
- 230000001070 adhesive effect Effects 0.000 claims description 16
- 239000000758 substrate Substances 0.000 claims description 12
- 235000015220 hamburgers Nutrition 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 239000000463 material Substances 0.000 abstract description 8
- 238000004904 shortening Methods 0.000 abstract description 2
- 239000004449 solid propellant Substances 0.000 description 6
- 238000000034 method Methods 0.000 description 3
- 238000009864 tensile test Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011230 binding agent Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
Images
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a theoretical calculation method for tensile property of a composite propellant, which comprises the following steps: s1: measuring the specific surface area S and the surface area average particle size D of the filling particles by using a particle size spectrometer; s2: inquiring the empirical binding coefficient B and the particle density of the particles on the surface of the matrix, and performing a tensile experiment on the matrix according to an experimental formula to obtain the maximum tensile strength of the matrix; s3: substituting the formula of the established empirical composite material, and calculating the maximum tensile strength of the composite material containing the particles. The invention adopts a theoretical calculation method, thereby shortening the total test time, reducing the test amount, reducing the cost of manpower and material resources, solving the problems of long test time consumption, large test amount and large consumption of manpower and material resources in the prior art, and guiding the prediction test by comparing the maximum tensile strength and the maximum elongation calculated by theory with the test.
Description
Technical Field
The invention relates to the technical field of composite solid propellants, in particular to a theoretical calculation method for the maximum tensile strength of a composite solid propellant.
Background
The composite solid propellant is a composite material which takes a binder as a matrix and is filled with an energetic solid filler, and the mechanical property of the solid propellant is one of the most important physical properties of the propellant, which means that the propellant is subjected to various exotic response characteristics. Similar to other materials, these response characteristics are typically described in terms of modulus, strength, deformation, poisson's ratio, etc. parameters that occur when a load is applied. In the design of the charge, parameters related to the failure process, such as yield stress and elongation, breaking strength and elongation, are also used for description.
The application of the breaking strength of the composite material is successfully the theoretical prediction of the tensile strength and the maximum tensile rate of the fiber-filled composite material and the short fiber rubber composite material, and indicates that the strength of the composite material depends on the strength of a matrix, the volume of a filler, the particle size, the bonding strength of an adhesive and particles and the like.
As mentioned above, the factors influencing the mechanical properties of the composite energy-containing solid propellant are very complex, the tensile strength of the composite propellant is mainly determined by tests at present, few reports about the theoretical calculation method of the tensile strength and the maximum elongation of the composite solid propellant at home and abroad are reported at present, generally, after the formula is determined, a large number of repeated tests are needed to stretch the composite material, a large amount of manpower and material resources are needed, in addition, the test period is long, a large amount of 6-8 days are needed for one round of test, the test result can have certain fluctuation, and the test accuracy is influenced by a plurality of factors such as the test times, the charge process and the like.
Disclosure of Invention
The invention aims to solve the defects in the prior art and provides a theoretical calculation method for tensile property of a composite propellant.
In order to achieve the purpose, the invention adopts the following technical scheme:
a theoretical calculation method for tensile property of a composite propellant comprises the following steps:
s1: determining the specific surface area S, the surface area average particle diameter D and the volume fraction phi f of the filling particles by using a particle size spectrometer;
s2: inquiring the empirical binding coefficient B of the particles on the surface of the substrate and the particle density phi ρ Carrying out a tensile test on the substrate according to the test formula to obtain the maximum tensile strength of the substrate;
s3: substituting the established empirical formula of the composite material, and calculating the maximum tensile strength of the composite material containing the particles;
s4: introducing the related viscoelastic performance of each phase in three-phase viscoelasticity, establishing a three-phase viscoelasticity model, and performing stress-strain analysis on the model;
s5: the maximum tensile strength is obtained by calculation through the correlation of the time of the maximum tensile strength and the time of the maximum elongation, and the maximum elongation epsilon is obtained by solving the maximum tensile strength yc And finishing the calculation.
Further, for the step S3, the specific calculation method is as follows:
further, the specific manner of stress-strain analysis used in step S4 is as follows: sigma 1 =E 3 ε 2 ;
σ=σ 1 +σ 2 =ε 1 E 1 ;
ε=ε 1 +ε 2 +ε 3 ;
Wherein E is 1 : particle phase elastic modulus;
E 3 : adhesive matrix elastic modulus;
η 2 : middle sticky tapeThe viscosity coefficient of the mixture;
η 3 : adhesive matrix adhesive viscosity coefficient;
ε 1 : the phase change experienced by the particle phase;
ε 2 : the phase change experienced by the mesophase;
ε 3 : the phase change suffered by the matrix phase;
σ 1 : the base spring portion is strained;
σ 2 : stress on the pot-sticking part of the substrate;
σ=σ 1 +σ 2 : the stress to which the whole is subjected;
ε=ε 1 +ε 2 +ε 3 : the strain experienced overall.
Further, step S5 specifically includes the following steps:
s501: on the basis of the step S4, canceling the constitutive equations of epsilon 1, epsilon 2 and epsilon 3 to obtain the following formulas:
and eliminating epsilon 1, epsilon 2 and epsilon 3 to obtain a Burgers model constitutive equation:
wherein:
ε=C*t;
wherein: t =0, e =0, σ =0
Then, solving for:
in the formula, s is a transformation parameter;
Subjecting the formula (2) to inverse Laplace transform, and finishing to obtain a stretching speed v 0 The tensile stress σ (t) at direct drawing over time t is:
s502: maximum tensile strength to be solved, and E 1 ,E 3 ,η 2 ,η 3 In the calculation mode of the formula (1), solving to obtain the maximum elongation rate epsilon yc 。
Compared with the prior art, the invention has the beneficial effects that:
the invention adopts a theoretical calculation method, thereby shortening the total test time, reducing the test amount, reducing the cost of manpower and material resources, and solving the problems of long test time consumption, large test amount and large consumption of manpower and material resources in the prior art;
dividing the composite into a matrix phase, an interface phase, and a particle phase, thereby calculating tensile properties;
the tensile strength of the steel is calculated by combining a compound empirical formula with test data
Solving the maximum stretching rate through the built burger constitutive model;
the maximum tensile strength and the maximum elongation are calculated theoretically and compared with the test, so that the prediction test is guided.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention.
FIG. 1 is a viscoelastic constitutive model diagram of a theoretical calculation method for tensile property of a composite propellant, which is provided by the invention;
FIG. 2 is a three-phase tensile strength schematic diagram of a theoretical calculation method for tensile properties of a composite propellant.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments.
Referring to fig. 1-2, the theoretical calculation method of the tensile property of the composite propellant comprises the following steps:
s1: determining the specific surface area S, the surface area average particle diameter D and the volume fraction phi f of the filling particles by using a particle size spectrometer;
s2: inquiring the empirical binding coefficient B of the particles on the surface of the substrate and the particle density phi ρ Carrying out a tensile test on the substrate according to the test formula to obtain the maximum tensile strength of the substrate;
s3: substituting the formula of the established empirical composite material, and calculating the maximum tensile strength of the composite material containing the particles;
s4: introducing the related viscoelastic performance of each phase in three-phase viscoelasticity, establishing a three-phase viscoelasticity model, and performing stress-strain analysis on the model;
specifically, the stress-strain analysis is performed in the following manner: sigma 1 =E 3 ε 2 ;
σ=σ 1 +σ 2 =ε 1 E 1 ;
ε=ε 1 +ε 2 +ε 3 ;
Wherein E is 1 : particle phase elastic modulus;
E 3 : adhesive matrix elastic modulus;
η 2 : an intermediate adhesive viscosity coefficient;
η 3 : adhesive matrix adhesive viscosity coefficient;
ε 1 : the phase change experienced by the particle phase;
ε 2 : the phase change experienced by the mesophase;
ε 3 : the phase change suffered by the matrix phase;
σ 1 : the base spring portion is strained;
σ 2 : stress on the part of the substrate stuck to the pot;
σ=σ 1 +σ 2 : the overall stress experienced;
ε=ε 1 +ε 2 +ε 3 : the strain experienced by the population;
s5: the maximum tensile strength is obtained by calculation through the correlation of the time of the maximum tensile strength and the time of the maximum elongation, and the maximum elongation epsilon is obtained by solving the maximum tensile strength yc And completing the calculation.
The step S5 specifically includes the following steps:
s501: on the basis of the step S4, eliminating the constitutive equations of epsilon 1, epsilon 2 and epsilon 3 to obtain the following formulas:
and eliminating epsilon 1, epsilon 2 and epsilon 3 to obtain a Burgers model constitutive equation:
wherein:
ε=C*t;
wherein: t =0, e =0, σ =0
Then, solving, according to the boundary conditions:
in the formula, s is a transformation parameter;
Subjecting the formula (2) to inverse Laplace transform, and finishing to obtain a stretching speed v 0 The tensile stress σ (t) at the time of direct drawing with time t is:
s502: maximum tensile strength to be solved, and E 1 ,E 3 ,η 2 ,η 3 In the calculation mode of the formula (1), the maximum elongation epsilon is obtained by solving yc 。
In order to better understand the technical solution of the present invention, the following description is further explained with reference to data and drawings.
Preferred embodiments of the same invention
The formulation (mass fraction) of the composite propellant table is shown in table 1.
TABLE 1
A theoretical calculation method for tensile property of a composite propellant comprises the following steps:
the method comprises the following steps: determining the specific surface area S, the surface area average particle diameter D and the volume fraction phi f of the filling particles by using a particle size spectrometer;
wherein, the physical signs of the granules are shown in Table 2.
Name of particle | Specific surface area S (m) 2 /g) | Surface area average particle diameter D (. Mu.m) |
AP (100-140 mesh) | 2.3 | 64.184 |
AL powder (13 μm) | 345.6 | 11.9 |
TABLE 2
Step two: inquiring the empirical binding coefficient B of the particles on the surface of the substrate and the particle density phi ρ According to the experimental formula, carrying out a tensile experiment on the matrix to obtain the maximum tensile strength of the matrix;
the particle interfacial properties are shown in Table 3.
TABLE 3
Step three: substituting the formula of the established empirical composite material, and calculating the maximum tensile strength of the composite material containing the particles, wherein m is 0.67;
the specific calculation method is as follows:
as shown in fig. 1, step four: establishing three-phase viscoelasticity, and carrying out stress-strain analysis on the model;
specifically, the stress-strain analysis is performed in the following manner: sigma 1 =E 3 ε 2 ;
σ=σ 1 +σ 2 =ε 1 E 1 ;
ε=ε 1 +ε 2 +ε 3 ;
Wherein E is 1 : particle phase elastic modulus;
E 3 : adhesive matrix elastic modulus;
η 2 : an intermediate adhesive viscosity coefficient;
η 3 : adhesive matrix adhesive viscosity coefficient;
ε 1 : the phase change experienced by the particle phase;
ε 2 : the phase change experienced by the mesophase;
ε 3 : the phase change suffered by the matrix phase;
σ 1 : the base spring portion is strained;
σ 2 : stress on the part of the substrate stuck to the pot;
σ=σ 1 +σ 2 : the overall stress experienced;
ε=ε 1 +ε 2 +ε 3 : the strain experienced overall.
Step five: on the basis of the fourth step, eliminating the constitutive equation of epsilon 1, epsilon 2 and epsilon 3 to obtain the following formula:
elimination of epsilon 1 ,ε 2 And ε 3 Get the constitutive equation of Burgers model
Wherein,
ε=C*t;
where t =0, e =0, σ =0,
the original equation can be solved according to the boundary conditions
In the formula, s is a transformation parameter;is the laplace transform of stress. Subjecting the formula (2) to inverse Laplace transform, and finishing to obtain a stretching speed v 0 The tensile stress sigma (t) in direct drawing is in relation to the time t
α=(p1+sqrt(p1^2-4p2))/2*p2;
β=(p1-sqrt(p1^2-4p2)/2*p2。
Step six: the maximum tensile strength solved for in the third step, and E 1 ,E 3 ,η 2 ,η 3 Substituting the formula (1), solving the equation to obtain the maximum elongation epsilon yc And finishing the calculation.
Wherein, the theoretical calculation is compared with the experiment, and the specific result is shown in the table 4.
TABLE 4
Compared with the experimental data, the error is in a small range, the calculation result can predict the test, the number of times of the test is minimized, the test efficiency is improved, unnecessary non-technical cost is avoided, powerful guidance and reference are provided for the current composite material tensile test, and a theoretical calculation method is adopted, so that the total test time is shortened, the test amount is reduced, the labor and material cost is reduced, and the problems of long test time consumption, large test amount and large consumption of manpower and material resources in the prior art are solved;
the maximum tensile strength and the maximum elongation are calculated by theory and compared with the test, thereby guiding the prediction test.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art should be considered as the technical solutions and the inventive concepts of the present invention within the technical scope of the present invention.
Claims (4)
1. A theoretical calculation method for tensile property of a composite propellant is characterized by comprising the following steps:
s1: measuring the specific surface area S and the surface area average particle size D of the filling particles by using a particle size spectrometer;
s2: inquiring the empirical binding coefficient B of the particles on the surface of the substrate and the particle density phi ρ According to the experimental formula, carrying out a tensile experiment on the matrix to obtain the maximum tensile strength of the matrix;
s3: substituting the established empirical formula of the composite material, and calculating the maximum tensile strength of the composite material containing the particles;
s4: introducing the related viscoelastic performance of each phase in three-phase viscoelasticity, establishing a three-phase viscoelasticity model, and performing stress-strain analysis on the model;
s5: the maximum tensile strength is obtained by calculation through the correlation of the time of the maximum tensile strength and the time of the maximum elongation, and the maximum elongation epsilon is obtained by solving the maximum tensile strength yc And finishing the calculation.
3. the theoretical calculation method for tensile property of composite propellant according to claim 2, wherein the stress-strain analysis in step S4 is performed in the following manner:
σ 1 =E 3 ε 2 ;
σ=σ 1 +σ 2 =ε 1 E 1 ;
9=ε 1 +9 2 +9 3 ;
wherein E is 1 : particle phase elastic modulus;
E 3 : adhesive matrix elastic modulus;
η 2 : an intermediate adhesive viscosity coefficient;
η 3 : adhesive matrix adhesive viscosity coefficient;
ε 1 : the phase change experienced by the particle phase;
ε 2 : the phase change experienced by the mesophase;
ε 3 : the phase change suffered by the matrix phase;
σ 1 : the base spring portion is strained;
σ 2 : stress on the part of the substrate stuck to the pot;
σ=σ 1 +σ 2 : the overall stress experienced;
ε=ε 1 +ε 2 +ε 3 : the overall strain experienced.
4. The theoretical calculation method for tensile property of composite propellant according to claim 3, wherein the step S5 further comprises the following steps:
s501: on the basis of the step S4, canceling the constitutive equations of epsilon 1, epsilon 2 and epsilon 3 to obtain the following formulas:
and eliminating epsilon 1, epsilon 2 and epsilon 3 to obtain a Burgers model constitutive equation:
wherein:
ε=C*t;
wherein: t =0, e =0, σ =0
Then, solving for:
in the formula, s is a transformation parameter;
Inverse Laplace transformation of formula (2)Alternatively, it is finished to a drawing rate v 0 The tensile stress σ (t) at the time of direct drawing with time t is:
s502: maximum tensile strength to be solved, and E 1 ,E 3 ,η 2 ,η 3 In the calculation mode of the formula (1), solving to obtain the maximum elongation rate epsilon yc 。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210922297.1A CN115270479A (en) | 2022-08-02 | 2022-08-02 | Theoretical calculation method for tensile property of composite propellant |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210922297.1A CN115270479A (en) | 2022-08-02 | 2022-08-02 | Theoretical calculation method for tensile property of composite propellant |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115270479A true CN115270479A (en) | 2022-11-01 |
Family
ID=83746684
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210922297.1A Pending CN115270479A (en) | 2022-08-02 | 2022-08-02 | Theoretical calculation method for tensile property of composite propellant |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115270479A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115906316A (en) * | 2022-11-21 | 2023-04-04 | 江苏科技大学 | Core-shell carrier design and preparation method without negative mechanical influence on matrix |
-
2022
- 2022-08-02 CN CN202210922297.1A patent/CN115270479A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115906316A (en) * | 2022-11-21 | 2023-04-04 | 江苏科技大学 | Core-shell carrier design and preparation method without negative mechanical influence on matrix |
CN115906316B (en) * | 2022-11-21 | 2023-09-22 | 江苏科技大学 | Design and preparation method of core-shell carrier without negative mechanical influence on matrix |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108897900B (en) | Method for predicting fatigue short crack propagation life under multi-axis variable amplitude loading | |
CN111554362B (en) | Method for establishing HBPRC dynamic damage constitutive model | |
Steeves et al. | Collapse mechanisms of sandwich beams with composite faces and a foam core, loaded in three-point bending. Part II: experimental investigation and numerical modelling | |
CN108627328B (en) | A kind of Notched specimen Multiaxial Fatigue Life Prediction method of coupling stress gradient modification | |
Elber | Failure mechanics in low-velocity impacts on thin composite plates | |
CN113420376B (en) | Multi-scale-based impact mechanical property simulation method for carbon fiber composite material | |
CN115270479A (en) | Theoretical calculation method for tensile property of composite propellant | |
CN111144016B (en) | Multi-axis stress state concrete variable four-parameter damage structure | |
Walters et al. | Cohesive finite element modeling of the delamination of HTPB binder and HMX crystals under tensile loading | |
CN110274826A (en) | A kind of hard metal material multiaxis high cycle fatigue failure prediction method based on single shaft fatigue S-N curve | |
CN103942418A (en) | Method for determining specific fatigue strength rated value under multi-axial loading condition | |
Naik et al. | Fabrication and strength analysis of hybrid jute-glass-silk fiber polymer composites based on hybrid Taguchi-WASPAS method | |
Charkaluk et al. | Revisiting the dang van criterion | |
Praveenkumar et al. | Selecting optimized mix proportion of bagasse ash blended high performance concrete using analytical hierarchy process (AHP) | |
Lina et al. | FEM analysis of spring-backs in age forming of aluminum alloy plates | |
CN112287474A (en) | Algorithm, equipment and storage medium for simulating mechanical behavior of polyvinyl alcohol fiber concrete | |
Jagannatha et al. | Influence of carbon & glass fiber reinforcements on flexural strength of epoxy matrix polymer hybrid composites | |
Harba et al. | Numerical analysis of reinforced concrete circular columns strengthening with CFRP under concentric and eccentric loadings | |
Tvergaard | Debonding of short fibres among particulates in a metal matrix composite | |
Khelifa et al. | Fatigue study of E-glass fiber reinforced polyester composite under fully reversed loading and spectrum loading | |
Goldberg et al. | A study of influencing factors on the tensile response of a titanium matrix composite with weak interfacial bonding | |
Paredes et al. | Multi-response optimization of mechanical properties of hybrid (fiberglass/abaca woven) in polyester matrix using desirability function based on DOE | |
Ueda et al. | Quasi-visco-elasto-plastic constitutive model of concrete for fatigue simulation | |
Kusmiran et al. | Numerical analysis of Composite with Natural Fiber Reinforcement using Finite Element Method: Leaf Spring Composite Application | |
Naseem et al. | Effect of ageing on mechanical properties of composite solid propellants |
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 |