CN114462146B - Method for constructing propellant creep constitutive model considering aging damage and finite element application - Google Patents
Method for constructing propellant creep constitutive model considering aging damage and finite element application Download PDFInfo
- Publication number
- CN114462146B CN114462146B CN202210107098.5A CN202210107098A CN114462146B CN 114462146 B CN114462146 B CN 114462146B CN 202210107098 A CN202210107098 A CN 202210107098A CN 114462146 B CN114462146 B CN 114462146B
- Authority
- CN
- China
- Prior art keywords
- propellant
- aging
- creep
- damage
- tensor
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- 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/04—Ageing analysis or optimisation against ageing
-
- 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/08—Thermal analysis or thermal optimisation
-
- 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)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Analytical Chemistry (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Crystallography & Structural Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention discloses a propellant creep constitutive model construction and finite element application method considering aging damage, wherein the construction method comprises the following steps: acquiring creep parameters of the propellant; acquiring an aging parameter of the propellant; presetting relative variation of effective bearing area; determining an aging development equation according to the obtained aging parameters; constructing a creep constitutive model of the propellant according to the creep parameters; determining a damage variable according to the relative variation of the effective bearing area, and determining a damage development equation according to the damage variable; establishing a propellant creep constitutive model considering aging damage according to the determined aging development equation, the damage development equation and the constructed propellant creep constitutive model; the finite element application method comprises the following steps: and decomposing the propellant creep damage constitutive model which is constructed by the method and takes the aging effect into consideration, and deducing an aging increment form and a consistent tangential stiffness array. The method can calculate the creep damage mechanical behavior and the aging mechanical behavior of the propellant grain structure in the storage stage.
Description
Technical Field
The invention relates to a method for constructing a propellant creep constitutive model considering aging damage and a method for applying a finite element, belonging to the technical field of energetic material constitutive model research.
Background
The solid rocket engine is a power device of a remote rocket projectile, and the solid propellant is an important component of the solid rocket engine. During long-distance rocket projectile storage for a long time, the grain structure can generate an aging effect besides a creep and damage effect. In order to accurately simulate the structural response of an engine during long-term storage, it is first necessary to establish an accurate model of propellant creep damage that accounts for aging effects.
At present, the aging constitutive relation of the propellant is mostly established on the basis of the relaxed constitutive structure, and for the creep constitutive structure, the research on the aging effect of the propellant on the basis of the aging constitutive relation is not reported. In order to realize the refined modeling of the aging creep damage effect inside the grain structure during long-term storage, it is necessary to construct a propellant creep damage constitutive model considering the aging effect.
At present, the existing finite element simulation tool does not have the function of creep damage of the structure considering the aging effect, and is difficult to be directly utilized. Therefore, a secondary development interface of commercial software is needed to be adopted to realize development and application of the propellant creep damage organization considering the aging effect.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, provides a method for constructing a propellant creep constitutive model considering aging damage and applying a finite element, and can calculate the creep damage mechanical behavior and the aging mechanical behavior of a propellant grain structure in the storage stage.
In order to achieve the purpose, the invention adopts the following technical scheme:
in one aspect, the invention provides a method for constructing a propellant creep constitutive model considering aging damage, which comprises the following steps:
acquiring creep parameters of the propellant;
acquiring an aging parameter of the propellant;
presetting relative variation of effective bearing area;
determining an aging development equation according to the obtained aging parameters;
constructing a creep constitutive model of the propellant according to the creep parameters;
determining a damage variable according to the relative variation of the effective bearing area, and determining a damage development equation according to the damage variable;
and establishing a propellant creep constitutive model considering aging damage according to the determined aging development equation, the damage development equation and the constructed propellant creep constitutive model.
Further, the determining an aging development equation according to the obtained aging parameters includes the following steps:
the aging development parameters of the propellant are obtained through formulas (1) to (2), and are specifically as follows:
wherein, T 0 For storage temperature, k B Is a function of the Boltzmann constant,aging reaction activation energy without stress strain field effect, t' aging time, h P Is a Planck constant, A is an aging reaction rate constant, upsilon m At maximum crosslink density, α c The undetermined coefficient of the chemical reaction equation,beta is the initial ageing development parameter of the propellant c Aging hair for propellantUnfolding parameters;
the relative change in the degree of crosslinking of the propellant is obtained by the formulae (3) to (4), as follows:
wherein v () is a function of the degree of crosslinking of the propellant, T' is the aging temperature, B is a parameter of the aging development of the propellant, v 0 Is the initial crosslink density;
the aging development equation is determined by the formula (5), which is as follows:
wherein J () is a creep compliance function taking into account aging effects, J 0 () Is the creep modulus without aging, J 1 () Is the change in creep compliance caused by aging.T is the loading time as a function of the relative change in the degree of crosslinking of the propellant.
Further, the injury development equation includes equations (6) - (9), which are as follows:
wherein D is a damage variable, D 0 For initial damage of the propellant, D △V To extend the new lesions formed, σ th The method comprises the following steps of setting a damage stress threshold value, wherein sigma is mechanical stress before propellant damage, gamma is a first parameter of propellant damage, K is a second parameter of propellant damage, beta is a third parameter of propellant damage, h () is a stress state function, and x is an input parameter of the stress state function.
Further, the propellant creep constitutive model comprises formulas (10) - (12), which are as follows:
wherein epsilon ij For mechanical strain, J ijkl () Is a creep tensor function, xi is the conversion time of time t, xi' is the conversion time of the first integrand tau, sigma kl τ is a first integrand function, T is a loading time, T () is a temperature, tt' is a second integrand function, a T Is a temperature shift factor, C 1 For the first parameter of the propellant WLF equation, C 2 For the second parameter of the WLF equation for the propellant, T is the current temperature, T r Is the reference temperature.
Further, the establishing of the propellant creep constitutive model considering the aging damage comprises updating the mechanical strain by the formula (13), which is specifically as follows:
wherein, the first and the second end of the pipe are connected with each other,as a function of the effective stress tensor of the propellant after damage.
In another aspect, the present invention provides a finite element application method of a propellant creep damage constitutive model considering aging effect, comprising the following steps:
decomposing the propellant creep damage constitutive model which is constructed by the method and takes the aging effect into consideration, and obtaining a partial strain tensor and a spherical strain tensor;
dispersing the partial strain tensor and the spherical strain tensor, and deducing an aging increment form of the propellant creep damage constitutive model according to the dispersed partial strain tensor and the spherical strain tensor;
deriving a consistent tangential stiffness array according to the derived aging increment form;
programming the aging increment form and the consistent tangential stiffness array to obtain a subprogram, and calling the subprogram by adopting finite element software.
Further, decomposing the propellant creep damage constitutive model considering the aging effect to obtain the partial strain tensor and the spherical strain tensor including the formula (14) as follows:
wherein the content of the first and second substances,as a function of the effective stress deflection tensor,as a function of the effective stress sphere tensor, e ij () In order to be a function of the partial strain,is the spherical strain tensor, J () is the propellant creep compliance function, and ν is the propellant Poisson ratio;
the creep compliance function of the propellant comprises the formula (15), which is as follows:
wherein N is J The term for the creep compliance level,to balance the creep compliance without considering the aging effect,to account for the initial creep compliance of the aging effect,to account for the nth creep compliance coefficient of aging effects,to account for the change in equilibrium creep compliance of the aging effect,in order to consider the nth creep compliance coefficient of the aging effect, t is loading time;
the effective stress deflection tensor function includes the formula (16), which is specifically as follows:
wherein S is ij () Is a biased stress tensor function;
the effective stress sphere tensor function includes the following formula (17):
wherein σ kk () As a function of the ball stress tensor.
Further, the aging increment form of the propellant creep damage constitutive model comprises a formula (18), and the aging increment form comprises the following specific formula:
wherein < DELTA > ∈ ij (t m+1 ) Is t m+1 Increment of the strain tensor at time, Δ e ij (t m+1 ) Is t m+1 The increment of the moment strain offset is increased,is t m+1 Increment of the global strain tensor of time, delta ij Is the fourth kronecker symbol.
Further, the consistent tangential stiffness array comprises formulas (19) - (20), which are as follows:
wherein, C ijkl () Is a tangential stiffness tensor function, delta sigma ij (t m+1 ) Is t m+1 Stress increment at time,. DELTA.. Epsilon kl (t m+1 ) Is t m+1 Increment of the moment strain tensor, C 2222 (t m+1 ) Is t m+1 First component of the tangential stiffness tensor at time, C 2233 (t m+1 ) Is t m+1 The second component of the tangential stiffness tensor at time, C 2323 (t m+1 ) Is t m+1 Third component of the tangential stiffness tensor at time, C 3333 (t m+1 ) Is t m+1 The fourth component of the tangential stiffness tensor at time, C 1133 (t m+1 ) Is t m+1 Time of dayThe fifth component of the tangential stiffness tensor, C 1313 (t m+1 ) Is t m+1 The sixth component of the tangential stiffness tensor at time, C 1111 (t m+1 ) Is t m+1 The seventh component of the tangential stiffness tensor at time, C 1122 (t m+1 ) Is t m+1 The eighth component of the tangential stiffness tensor at time, C 1212 (t m+1 ) Is t m+1 A ninth component of the tangential stiffness tensor at time;
the tangential stiffness tensor function t m+1 Of all components of the tangential stiffness tensor at the time, other components are zero except the first component to the ninth component.
Compared with the prior art, the invention has the following beneficial effects:
according to the invention, the development relation of aging and damage is respectively determined by considering aging and damage effects, so that the constructed propellant creep deformation constitutive model containing the aging damage not only can reflect the creep deformation rule of the propellant grain structure, but also can accurately calculate the damage mechanical behavior and the aging mechanical behavior of the propellant grain structure in the creep deformation stage.
According to the method, the aging increment form and the consistent tangential stiffness array are deduced by decomposing the propellant creep constitutive model containing the aging damage, the application program is obtained by programming the aging increment form and the consistent tangential stiffness array, and the application program is called by using finite element software, so that a theoretical basis and an implementation means can be provided for the fine structural integrity analysis of the remote rocket projectile charge column structure in the storage stage.
Drawings
FIG. 1 is a flow chart of a method for constructing a creep constitutive model of a propellant considering aging damage according to the present invention;
FIG. 2 is a flow chart of a finite element application method of the propellant creep constitutive model considering aging damage according to the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
Example 1
The embodiment provides a method for constructing a creep constitutive model of a propellant considering aging damage, and with reference to fig. 1, the method comprises the following steps:
acquiring creep parameters of the propellant;
acquiring an aging parameter of the propellant;
presetting relative variation of effective bearing area;
determining an aging development equation according to the obtained aging parameters;
constructing a propellant creep constitutive model according to creep parameters;
determining a damage variable according to the relative variation of the effective bearing area, and determining a damage development equation according to the damage variable;
and establishing a propellant creep constitutive model considering aging damage according to the determined aging development equation, the damage development equation and the constructed propellant creep constitutive model.
According to the invention, the development relation of aging and damage is respectively determined by considering aging and damage effects, so that the constructed propellant creep deformation constitutive model containing the aging damage not only can reflect the creep deformation rule of the propellant grain structure, but also can accurately calculate the damage mechanical behavior and the aging mechanical behavior of the propellant grain structure in the creep deformation stage.
Example 2
On the basis of example 1, this example details the method of determining the aging development equation, the damage development equation, the propellant creep constitutive model, and the propellant creep constitutive model taking into account the aging damage.
Determining an aging development equation
Determining an aging development equation according to the obtained aging parameters, comprising the following steps:
the aging development parameter beta of the propellant is obtained by the formulas (1) to (2) c The method comprises the following steps:
wherein, T 0 For storage temperature, k B Is the Boltzmann constant and is,aging reaction activation energy without stress strain field effect, t' aging time, h P Is the Planck constant, A is the aging reaction rate constant, upsilon m At maximum crosslink density, α c The undetermined coefficient of the chemical reaction equation,is the initial aging development parameter of the propellant.
S12, obtaining the relative change of the crosslinking degree of the propellant through formulas (3) to (4), wherein the relative change is as follows:
wherein v () is a function of the degree of crosslinking of the propellant, T' is the aging temperature, B is a parameter of the aging development of the propellant, v 0 Is the initial crosslink density;
the aging development equation is determined by the formula (5), which is as follows:
wherein J () is a creep compliance function taking into account aging effects, J 0 () Is the creep modulus without aging, J 1 () Is the change in creep compliance caused by aging.T is the loading time as a function of the relative change in the degree of crosslinking of the propellant.
Equation of development of lesions
The injury development equation includes equations (6) - (9), which are as follows:
wherein D is a damage variable, D 0 For initial damage of the propellant, D △V To spread the new lesions formed, sigma th The method comprises the following steps of setting a damage stress threshold value, wherein sigma is mechanical stress before propellant damage, gamma is a first parameter of propellant damage, K is a second parameter of propellant damage, beta is a third parameter of propellant damage, h () is a stress state function, and x is an input parameter of the stress state function.
(III) propellant creep constitutive model
The propellant creep deformation constitutive model comprises formulas (10) to (12), and specifically comprises the following formulas:
wherein epsilon ij For mechanical strain, J ijkl () Is a creep tensor function, xi is the conversion time of time t, xi' is the conversion time of the first integrand tau, sigma kl τ is a first integrand, T () is temperature, tt' is a second integrand, a T Is a temperature shift factor, C 1 For the first parameter of the propellant WLF equation, C 2 For the second parameter of the propellant WLF equation, T is the current temperature, T r Is the reference temperature.
(IV) propellant creep constitutive model considering aging damage
Establishing a propellant creep constitutive model considering aging damage comprises updating mechanical strain through a formula (13), which specifically comprises the following steps:
wherein the content of the first and second substances,as a function of the effective stress tensor of the propellant after damage.
Example 3
The embodiment provides a finite element application method of a propellant creep damage constitutive model considering aging effect, and with reference to fig. 2, the method comprises the following steps:
decomposing the propellant creep damage constitutive model which is constructed in the embodiment 1 or 2 and takes the aging effect into consideration, and obtaining a partial strain tensor and a spherical strain tensor;
dispersing the partial strain tensor and the spherical strain tensor, and deducing an aging increment form of the propellant creep damage constitutive model according to the dispersed partial strain tensor and the spherical strain tensor;
deriving a consistent tangential stiffness array according to the derived aging increment form;
programming the aging increment form and the consistent tangential stiffness array to obtain a subprogram, and calling the subprogram by adopting finite element software.
According to the method, the aging increment form and the consistent tangential stiffness array are deduced by decomposing the propellant creep constitutive model containing the aging damage, the application program is obtained by programming the aging increment form and the consistent tangential stiffness array, and the application program is called by using finite element software, so that a theoretical basis and an implementation means can be provided for the fine structural integrity analysis of the remote rocket projectile charge column structure in the storage stage.
Example 4
On the basis of embodiment 3, this embodiment details a method for decomposing a propellant creep damage constitutive model considering an aging effect, a method for deriving an aging increment form of the propellant creep damage constitutive model, and a method for deriving a consistent tangential stiffness matrix.
Decomposition of propellant creep damage constitutive model considering aging effect
Decomposing a propellant creep damage constitutive model considering an aging effect to obtain a partial strain tensor and a spherical strain tensor, and comprising the following steps of:
s21 decomposes the strain tensor function into a partial strain tensor function and a spherical strain tensor function by the decomposition expression (21), which is specifically as follows:
wherein e is ij (t) is the partial strain tensor function, ε ij () As a function of the strain tensor, S ij () Is a biased stress tensor function, σ ij () As a function of the spherical stress tensor, Y 1 () As a function of the first creep of the propellant, Y 2 () As a function of the second creep of the propellant,is a spherical strain tensor function;
in application, the spherical strain tensor is updated by the formula (211), which is specifically as follows:
wherein alpha is T Is the coefficient of thermal expansion, theta is the change in temperature, epsilon kk Is the mechanical sphere strain tensor;
when the temperature change is applied, the temperature change comprises a formula (212), which is specifically as follows:
Θ=T-T 0 (212)
wherein T is the current temperature of the propellant, T 0 The initial temperature of the propellant.
S22, shearing the first creep function of the propellant and the second creep function of the propellant one by decomposing the formulas (22) to (24) to obtain a creep modulus function and a volume creep modulus function of the propellant creep constitutive model considering the aging damage, wherein the creep modulus function and the volume creep modulus function are as follows:
χ(t)=2(1+ν)J(t,t′) (23)
B(t)=3(1-2ν)J(t,t′) (24)
wherein x () is the creep modulus function of the damage-considered creep constitutive model of the propellant, and B () is the bulk creep modulus function.
S23 decomposes the polarization strain tensor function by expression (25), and decomposes the spherical strain tensor function by expression (26), as follows:
wherein, the first and the second end of the pipe are connected with each other,as a function of the effective stress deflection tensor,as a function of the effective stress sphere tensor.
S24, updating the bias strain tensor function and the spherical strain tensor function by the equation (14) according to the Stieltjes convolution integral definition, specifically as follows:
wherein the content of the first and second substances,as a function of the effective stress deflection tensor,as a function of the effective stress sphere tensor, e ij () In order to be a function of the partial strain,is a spherical strain tensor function;
in application, the creep compliance function of the propellant comprises the formula (15), which is as follows:
wherein N is J The term for the creep compliance level,to balance the creep compliance without considering the aging effect,to account for the initial creep compliance of the aging effect,to account for the nth creep compliance coefficient of aging effects,to account for changes in equilibrium creep compliance due to aging effects,the nth creep compliance coefficient for considering the aging effect, and t is loading time;
in application, the effective stress deflection tensor function includes the formula (16), which is as follows:
wherein S is ij () Is a biased stress tensor function;
in this embodiment, the effective stress sphere tensor function includes the following equation (17):
wherein σ kk () As a function of the spherical stress tensor.
(II) deducing aging increment form of propellant creep damage constitutive model
The aging increment form of the propellant creep damage constitutive model comprises a formula (18), and the aging increment form comprises the following specific formula:
wherein < DELTA > ∈ ij (t m+1 ) Is t m+1 Increment of the strain tensor at time, Δ e ij (t m+1 ) Is t m+1 The increment of the moment strain offset is increased,is t m+1 Increment of the global strain tensor of time, delta ij Is a fourth kronecker symbol.
(III) deriving a consistent tangential stiffness matrix
Deriving a consistent tangential stiffness matrix from the derived incremental form of aging, comprising the steps of:
s31 simplifies the aging increment form by equations (28) to (29), specifically as follows:
wherein the content of the first and second substances,is t m+1 The increase in time relative to the crosslink density,is t m+1 The increment of the equivalent stress offset at that time,is t m+1 Increment of the ball stress tensor at time, gamma 0,J () Is a function of a first intermediate variable, gamma 1,J () In order to construct the second intermediate variable function,to construct the third intermediate variable function,is the fourth intermediate variable function of the constitutive.
In application, the compound is obtained by the formula (30)The method comprises the following specific steps:
wherein the content of the first and second substances,in order to construct the fifth intermediate variable,is the sixth intermediate variable of the structure.
Obtaining gamma by the formulae (31) to (32) 0,J The method comprises the following steps:
wherein, the first and the second end of the pipe are connected with each other,is a seventh intermediate variable function of the constitutive model;
obtaining gamma by the formula (33) 1,J The method comprises the following steps:
s32, deducing a strain increment comprising the formulas (35) - (37) through a simplified aging increment form, wherein the specific expression is as follows:
wherein, the first and the second end of the pipe are connected with each other,is the eighth intermediate variable of the structure.
S33, deriving a consistent tangential stiffness array according to the strain increment, wherein the consistent tangential stiffness array comprises the following formulas (19) - (20):
wherein, C ijkl () Is a tangential stiffness tensor function, delta sigma ij (t m+1 ) Is t m+1 Stress increment at time,. DELTA.. Epsilon kl (t m+1 ) Is t m+1 Increment of strain tensor at time, C 2222 (t m+1 ) Is t m+1 First component of the tangential stiffness tensor at moment, C 2233 (t m+1 ) Is t m+1 The second component of the tangential stiffness tensor at time, C 2323 (t m+1 ) Is t m+1 Third component of the tangential stiffness tensor at time, C 3333 (t m+1 ) Is t m+1 The fourth component of the tangential stiffness tensor at time, C 1133 (t m+1 ) Is t m+1 The fifth component of the tangential stiffness tensor at time, C 1313 (t m+1 ) Is t m+1 The sixth component of the tangential stiffness tensor at time, C 1111 (t m+1 ) Is t m+1 The seventh component of the tangential stiffness tensor at time, C 1122 (t m+1 ) Is t m+1 The eighth component of the tangential stiffness tensor at time, C 1212 (t m+1 ) Is t m+1 A ninth component of the tangential stiffness tensor at time;
furthermore, the tangential stiffness tensor function t m+1 Of all components of the tangential stiffness tensor at the time, other components are zero except the first component to the ninth component.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and so forth) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While the present invention has been described with reference to the particular illustrative embodiments, it is to be understood that the invention is not limited to the disclosed embodiments, but is intended to cover various modifications, equivalent arrangements, and equivalents thereof, which may be made by those skilled in the art without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (7)
1. A method for constructing a propellant creep constitutive model considering aging damage is characterized by comprising the following steps:
acquiring creep parameters of the propellant;
acquiring an aging parameter of the propellant;
presetting relative variation of effective bearing area;
determining an aging development equation according to the obtained aging parameters;
constructing a creep constitutive model of the propellant according to the creep parameters;
determining a damage variable according to the relative variation of the effective bearing area, and determining a damage development equation according to the damage variable;
establishing a propellant creep constitutive model considering aging damage according to the determined aging development equation, the damage development equation and the constructed propellant creep constitutive model;
the injury development equation comprises formulas (6) - (9), which are as follows:
wherein D is a damage variable, D 0 For initial damage of the propellant, D ΔV To extend the new lesions formed, σ th The method comprises the following steps of (1) setting a damage stress threshold value, wherein sigma is mechanical stress before propellant damage, gamma is a first parameter of propellant damage, K is a second parameter of propellant damage, beta is a third parameter of propellant damage, h () is a stress state function, and x is an input parameter of the stress state function;
the propellant creep deformation constitutive model comprises formulas (10) to (12), and specifically comprises the following steps:
wherein epsilon ij For mechanical strain, J ijkl () As a function of the creep tensor, ξ isThe reduced time, ξ' is the reduced time of the first integrand τ, σ kl τ is a first integrand, T () is temperature, tt' is a second integrand, a T Is a temperature shift factor, C 1 For the first parameter of the propellant WLF equation, C 2 For the second parameter of the WLF equation for the propellant, T is the current temperature, T r Is the reference temperature.
2. The method for constructing the propellant creep constitutive model considering the aging damage as claimed in claim 1, wherein the determining the aging development equation according to the obtained aging parameters comprises the following steps:
the aging development parameters of the propellant are obtained through formulas (1) to (2), and are specifically as follows:
wherein, T 0 For storage temperature, k B Is the Boltzmann constant and is,aging reaction activation energy without stress strain field effect, t' aging time, h P Is the Planck constant, A is the aging reaction rate constant, upsilon m At maximum crosslink density, α c The undetermined coefficient of the chemical reaction equation,for the initial ageing development parameter of the propellant, beta c Is a propellant aging development parameter;
the relative change in the degree of crosslinking of the propellant is obtained by the formulae (3) to (4), as follows:
wherein v () is a function of the degree of crosslinking of the propellant, T' is the aging temperature, B is a parameter of the aging development of the propellant, v 0 Is the initial crosslink density;
the aging development equation is determined by the formula (5), which is as follows:
3. The method for constructing the aging damage-considered propellant creep constitutive model according to claim 2, wherein the establishing of the aging damage-considered propellant creep constitutive model comprises updating mechanical strain by a formula (13), and specifically comprises the following steps:
4. A finite element application method of a propellant creep damage constitutive model considering an aging effect is characterized by comprising the following steps:
decomposing the aging effect-considered propellant creep damage constitutive model constructed according to any one of claims 1 to 3 to obtain a partial strain tensor and a spherical strain tensor;
dispersing the partial strain tensor and the spherical strain tensor, and deducing an aging increment form of the propellant creep damage constitutive model according to the dispersed partial strain tensor and the spherical strain tensor;
deriving a consistent tangential stiffness array according to the derived aging increment form;
programming the aging increment form and the consistent tangential stiffness array to obtain a subprogram, and calling the subprogram by adopting finite element software.
5. The finite element application method of the aging effect considered propellant creep damage constitutive model as claimed in claim 4, wherein the decomposing of the aging effect considered propellant creep damage constitutive model to obtain the partial strain tensor and the spherical strain tensor comprises the following formula (14):
wherein the content of the first and second substances,as a function of the effective stress deflection tensor,as a function of the effective stress sphere tensor, e ij () In order to be a function of the partial strain,is the spherical strain tensor, and J () is the propellant creep complianceAs a function, ν is the propellant poisson's ratio;
the creep compliance function of the propellant comprises the formula (15), which is as follows:
wherein N is J The number of terms for the creep compliance level,to balance the creep compliance without considering the aging effect,to account for the initial creep compliance of the aging effect,to account for the nth creep compliance coefficient of aging effects,to account for the change in equilibrium creep compliance of the aging effect,the nth creep compliance coefficient for considering aging effect; t is loading time;
the effective stress deflection tensor function comprises a formula (16), which is as follows:
wherein S is ij () Is a biased stress tensor function;
the effective stress sphere tensor function includes the following formula (17):
wherein σ kk () As a function of the ball stress tensor.
6. The finite element application method of the propellant creep damage constitutive model considering the aging effect as claimed in claim 4, wherein the aging increment form of the propellant creep damage constitutive model comprises an expression (18), and the aging increment form comprises the following specific expression:
wherein, delta epsilon ij (t m+1 ) Is t m+1 Increment of the time-of-day strain tensor, Δ e ij (t m+1 ) Is t m+1 The increment of the offset in strain at the moment,is t m+1 Increment of the ball strain tensor at time, Δ e ij (t m ) Is t m The increment of the moment strain offset is increased,is t m Increment of the strain tensor of the sphere at time, delta ij Is a fourth kronecker symbol.
7. The finite element application method of a propellant creep damage constitutive model considering the aging effect as claimed in claim 4, wherein the consistent tangential stiffness matrix comprises the following equations (19) - (20):
wherein, C ijkl () As a function of the tangential stiffness tensor, Δ σ ij (t m+1 ) Is t m+1 Stress increment of time, delta epsilon kl (t m+1 ) Is t m+1 Increment of the moment strain tensor, C 2222 (t m+1 ) Is t m+1 First component of the tangential stiffness tensor at moment, C 2233 (t m+1 ) Is t m+1 The second component of the tangential stiffness tensor at time, C 2323 (t m+1 ) Is t m+1 Third component of the tangential stiffness tensor at time, C 3333 (t m+1 ) Is t m+1 The fourth component of the tangential stiffness tensor at time, C 1133 (t m+1 ) Is t m+1 The fifth component of the tangential stiffness tensor at time, C 1313 (t m+1 ) Is t m+1 The sixth component of the tangential stiffness tensor at time, C 1111 (t m+1 ) Is t m+1 The seventh component of the tangential stiffness tensor at time, C 1122 (t m+1 ) Is t m+1 The eighth component of the tangential stiffness tensor at time, C 1212 (t m+1 ) Is t m+1 A ninth component of the moment tangent stiffness tensor;
the tangential stiffness tensor function t m+1 Of all components of the tangential stiffness tensor at the time, other components are zero except the first component to the ninth component.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210107098.5A CN114462146B (en) | 2022-01-28 | 2022-01-28 | Method for constructing propellant creep constitutive model considering aging damage and finite element application |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210107098.5A CN114462146B (en) | 2022-01-28 | 2022-01-28 | Method for constructing propellant creep constitutive model considering aging damage and finite element application |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114462146A CN114462146A (en) | 2022-05-10 |
CN114462146B true CN114462146B (en) | 2023-02-03 |
Family
ID=81411998
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210107098.5A Active CN114462146B (en) | 2022-01-28 | 2022-01-28 | Method for constructing propellant creep constitutive model considering aging damage and finite element application |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114462146B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115204003B (en) * | 2022-06-15 | 2023-09-29 | 中国人民解放军火箭军工程大学 | Method for constructing nonlinear creep-type viscoelastic constitutive model of polymer material |
CN115859624B (en) * | 2022-12-05 | 2023-06-13 | 中国人民解放军陆军工程大学 | Propellant fine constitutive relation quantification method |
CN115906583B (en) * | 2022-12-16 | 2023-08-01 | 中国人民解放军陆军工程大学 | Grain structure integrity simulation analysis method and system based on virtual unit method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101509855A (en) * | 2009-03-17 | 2009-08-19 | 华东理工大学 | On-line prediction method for high-temperature pipe damage and longevity |
CN109783985A (en) * | 2019-02-22 | 2019-05-21 | 湖南大学 | The building and application method of a kind of inertia device with epoxy glue constitutive model |
CN109829231A (en) * | 2019-01-24 | 2019-05-31 | 北京理工大学 | A kind of propellant mechanics prediction technique based on CMDB propellant damage process |
WO2020143284A1 (en) * | 2019-01-11 | 2020-07-16 | 华东理工大学 | Abaqus-based multiaxial creep fatigue prediction method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107391801B (en) * | 2017-06-23 | 2020-06-26 | 中国人民解放军国防科学技术大学 | Construction and application method of dehumidifying propellant dynamic thermal bonding elastic constitutive model |
CN111709176B (en) * | 2020-07-27 | 2023-09-12 | 上海圣之尧智能科技有限公司 | Finite element simulation method and system considering plasticity and damage constitutive relation |
-
2022
- 2022-01-28 CN CN202210107098.5A patent/CN114462146B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101509855A (en) * | 2009-03-17 | 2009-08-19 | 华东理工大学 | On-line prediction method for high-temperature pipe damage and longevity |
WO2020143284A1 (en) * | 2019-01-11 | 2020-07-16 | 华东理工大学 | Abaqus-based multiaxial creep fatigue prediction method |
CN109829231A (en) * | 2019-01-24 | 2019-05-31 | 北京理工大学 | A kind of propellant mechanics prediction technique based on CMDB propellant damage process |
CN109783985A (en) * | 2019-02-22 | 2019-05-21 | 湖南大学 | The building and application method of a kind of inertia device with epoxy glue constitutive model |
Non-Patent Citations (7)
Title |
---|
A new constitutive equation for solid propellant with the effects of aging and viscoelastic Poisson’s ratio;Cui,HR 等;《MECCANICA》;20180731;第53卷(第9期);第2393-2410页 * |
HTPB复合固体推进剂含损伤和老化本构研究;阳建红等;《推进技术》;20021228;第23卷(第06期);第509-512页 * |
固体推进剂黏弹性本构模型及其有限元应用研究;崔辉如等;《国防科技大学学报》;20181028;第40卷(第05期);第27-32页 * |
基于改进Harris函数的岩石蠕变损伤模型及二次开发;李云等;《水电能源科学》;20200425;第38卷(第04期);第112-116页 * |
复合推进剂微裂纹损伤本构模型研究;彭威等;《固体火箭技术》;20030625;第26卷(第02期);第33-37页 * |
有限应变下蠕变-弹塑性-损伤耦合本构模型的数值模拟;李锡夔;《计算力学学报》;19970515;第14卷(第02期);第134-142页 * |
考虑粘弹性泊松比的固体推进剂本构模型研究及应用;崔辉如;《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》;20190115(第1期);C031-808 * |
Also Published As
Publication number | Publication date |
---|---|
CN114462146A (en) | 2022-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114462146B (en) | Method for constructing propellant creep constitutive model considering aging damage and finite element application | |
Miller | An inelastic constitutive model for monotonic, cyclic, and creep deformation: Part I—Equations development and analytical procedures | |
Zocher et al. | A three‐dimensional finite element formulation for thermoviscoelastic orthotropic media | |
CN112364535B (en) | Creep fatigue life prediction method based on crystal plasticity | |
CN114462147B (en) | Method for constructing damage-containing propellant creep deformation constitutive model and method for applying finite element | |
CN106777457B (en) | Reliability assessment software system for solid engine grain structure | |
CN111090957A (en) | High-temperature structure dangerous point stress-strain calculation method | |
Ramshaw et al. | Accelerated artificial compressibility method for steady-state incompressible flow calculations | |
Rizzi et al. | Comparison of nonlinear random response using equivalent linearization and numerical simulation | |
Morjaria et al. | Finite element analysis of time‐dependent inelastic deformation in the presence of transient thermal stresses | |
JP2007309704A (en) | Numerical analytical method, analyzer, and program for lead material | |
Liu et al. | A one-dimensional model for the aging and damage effects on the response of composite double-base propellants | |
CN115455593A (en) | Real-time hybrid test method based on multi-task loading, electronic equipment and storage medium | |
CN115146347A (en) | Method and equipment for calculating concrete shrinkage creep in bridge structure and storage medium | |
CN109829222B (en) | Compression mechanics constitutive method of composite propellant with wide strain rate and wide temperature range | |
Mueller | Real-time hybrid simulation with online model updating | |
Tulkki et al. | Viscoelastic modelling of Zircaloy cladding in-pile transient creep | |
Ben Naceur et al. | Micromechanical modeling of the ratcheting behavior of 304 stainless steel | |
CN117828929A (en) | Construction of propellant creep mechanism considering damage softening and finite element application method | |
Ye et al. | Optimal sensor placement for substructural response reconstruction | |
CN115859624B (en) | Propellant fine constitutive relation quantification method | |
Bass et al. | Numerical solution of the evolution equations of damage and rate-dependent plasticity | |
Fukasawa et al. | Development on Rubber Bearings for Sodium-Cooled Fast Reactor: Part 6—Proposal of New Type of Hysteresis Model for Ultimate Behavior | |
Heczko | Numerical modelling of oxidative ageing in deformed state | |
CN111814274A (en) | Polynomial dimension decomposition multi-degree-of-freedom rotor power system statistical moment analysis method |
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 |