A method for establishing a creep model of low-fat chocolate using fractional order calculation

Abstract

The present invention provides a method for establishing a low-fat chocolate creep model using fractional-order calculation, and relates to the field of food science and technology. The method includes establishing a fractional-order derivative constitutive model based on the Grünwald-Letnikov form according to the creep power-law behavior exhibited by the low-fat chocolate in rheological deformation; establishing a time-discretized fractional-order derivative constitutive equation; establishing an improved Grünwald-Letnikov fractional-order derivative creep model to obtain an approximate discrete numerical solution for creep compliance; performing a uniaxial creep experiment on the low-fat chocolate, fitting the experimental data using the established creep model, and obtaining constitutive parameters; by changing the constitutive parameters, numerically simulating the creep response of the low-fat chocolate predicted by the fractional-order derivative creep model when the constitutive parameters are different, and studying the influence of different constitutive parameters on the fractional-order derivative creep model. The present invention can provide a theoretical basis for the application and optimization of the fractional-order derivative model in the field of food engineering.

Description

Technical Field

[0001] The invention relates to the field of food science and technology, and in particular to a method for establishing a low-fat chocolate creep model by using fractional order calculation. Background Art

[0002] Chocolate is mainly composed of a continuous fat phase with sugar crystals and cocoa solids dispersed in it. However, due to the large amount of saturated fatty acids in fat, excessive intake will be significantly positively correlated with the risk of metabolic syndrome, cardiovascular disease, obesity and diabetes. As a new type of oil structure technology, oleogel shows great potential in partially or completely replacing solid fats containing unhealthy saturated fatty acids in chocolate. However, how to maintain similar rheological properties while ensuring the low-fat and healthy nature of food remains a big challenge. The lack of a concise and effective mathematical framework makes it complicated to simulate the rheological creep response of viscoelastic foods such as low-fat chocolate. Summary of the invention

[0003] To this end, an embodiment of the present invention provides a method for establishing a low-fat chocolate creep model using fractional-order calculation, which is used to solve the problem that existing models in the prior art are difficult to simply and accurately simulate the rheological creep response of viscoelastic foods such as low-fat chocolate.

[0004] In order to solve the above problems, an embodiment of the present invention provides a method for establishing a low-fat chocolate creep model using fractional order calculation, the method comprising:

[0005] Step S1: According to the creep power law behavior of low-fat chocolate in rheological deformation, a fractional derivative constitutive model based on the Grünwald-Letnikov form is established;

[0006] Step S2: Approximately rewrite the definition of the fractional derivative of the Grünwald-Letnikov form in the constitutive model, and establish a time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form;

[0007] Step S3: According to the fractional-order derivative constitutive equation based on the Grünwald-Letnikov form with time discretization, an improved Grünwald-Letnikov fractional-order derivative creep model is established to obtain an approximate discrete numerical solution for the creep compliance;

[0008] Step S4: performing a uniaxial creep experiment on the low-fat chocolate, fitting the experimental data using the established creep model, and obtaining constitutive parameters of the fractional derivative creep model.

[0009] Preferably, the method further comprises:

[0010] Step S5: by changing the constitutive parameters of the fractional derivative creep model, numerically simulating the creep response of low-fat chocolate predicted by the fractional derivative creep model when the constitutive parameters are different, the influence of different constitutive parameters on the fractional derivative creep model is studied.

[0011] Preferably, the fractional-order derivative constitutive model based on the Grünwald-Letnikov form is expressed as:

[0012]

[0013] In the formula, σ(t) and ε(t) refer to the stress and strain of the low-fat chocolate material, respectively; η and G refer to quasi-properties, that is, as intermediate parameters between elasticity and viscosity, they have both the elasticity of the spring and the viscosity of the damper; β and v refer to fractional orders, that is, fractional derivative operators; It means to find the υ-β derivative of σ(t), It means to find the υ-order derivative of ε(t).

[0014] Preferably, the definition of the fractional derivative based on the Grünwald-Letnikov form is as follows:

[0015]

[0016] In the formula, refers to the fractional derivative based on the Grünwald-Letnikov form, dt refers to the time step; m is a variable; F(t) represents an arbitrary function; Refers to the recursive formula, written as and Γ(·) refers to the Eulerian Gamma function; for The value of is calculated as follows: when

[0017] Preferably, the definition of the fractional-order derivative in the Grünwald-Letnikov form in the constitutive model is approximately rewritten to establish a time discretized fractional-order derivative constitutive equation based on the Grünwald-Letnikov form, specifically including:

[0018] First, consider discretizing time t, that is, writing it as t = {t1, t2, …, t N-1 ,t N}, and t N =t N-1 +dt, formula (2) can be further rewritten as:

[0019]

[0020] According to formula (3), the time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form is given:

[0021]

[0022] Further expand formula (4):

[0023]

[0024] Reintegrating formula (5), we get:

[0025]

[0026] Where N represents the number of time t.

[0027] Preferably, according to the fractional-order derivative constitutive equation based on the Grünwald-Letnikov form of time discretization, an improved Grünwald-Letnikov fractional-order derivative creep model is established to obtain an approximate discrete numerical solution for the creep compliance, specifically including:

[0028] In order to analyze the creep behavior of low-fat chocolate, the creep compliance J(t) is first given:

[0029]

[0030] Then the fractional derivative formula of the constant σ0 is introduced as follows:

[0031]

[0032] Then substitute formula (8) into formula (6) to obtain:

[0033]

[0034] Then, formula (9) is rearranged to obtain the improved Grünwald-Letnikov fractional-order derivative discretization strain formula:

[0035]

[0036] Then the improved Grünwald-Letniko fractional derivative creep model is derived:

[0037]

[0038] In the formula, σ0 refers to constant stress and is a constant.

[0039] Based on the same inventive concept, an embodiment of the present invention further provides a system for establishing a low-fat chocolate creep model using fractional order calculation, and the system is used to implement the method for establishing a low-fat chocolate creep model using fractional order calculation as described above, specifically comprising:

[0040] The constitutive model building module is used to establish a fractional derivative constitutive model based on the Grünwald-Letnikov form according to the creep power law behavior of low-fat chocolate in rheological deformation;

[0041] The module for rewriting the definition of fractional derivatives is used to approximately rewrite the definition of fractional derivatives in the Grünwald-Letnikov form in the constitutive model and establish a time-discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form.

[0042] The creep model building module is used to establish an improved Grünwald-Letnikov fractional derivative creep model based on the time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form, and obtain an approximate discrete numerical solution for the creep compliance;

[0043] The experimental verification and parameter acquisition module is used to conduct uniaxial creep experiments on low-fat chocolate, fit the experimental data using the established creep model, and obtain the constitutive parameters of the fractional derivative creep model.

[0044] Preferably, the system further comprises:

[0045] The numerical simulation and parameter analysis module is used to change the constitutive parameters of the fractional derivative creep model, numerically simulate the creep response of low-fat chocolate predicted by the fractional derivative creep model when the constitutive parameters are different, and study the influence of different constitutive parameters on the fractional derivative creep model.

[0046] An embodiment of the present invention also provides an electronic device, which includes a processor, a memory and a bus system, wherein the processor and the memory are connected via the bus system, the memory is used to store instructions, and the processor is used to execute the instructions stored in the memory to implement the above-mentioned method of establishing a low-fat chocolate creep model using fractional-order calculation.

[0047] An embodiment of the present invention further provides a computer storage medium storing a computer software product. The computer software product includes several instructions for enabling a computer device to execute the above-mentioned method of establishing a low-fat chocolate creep model using fractional-order calculation.

[0048] It can be seen from the above technical solutions that the present invention has the following beneficial effects:

[0049] 1) The present invention uses an improved Grünwald-Letnikov fractional derivative creep model to simulate the creep power-law behavior of low-fat chocolate materials. The model has fewer parameters and can reflect the change process of material properties during the deformation of the chocolate material under constant stress.

[0050] 2) Based on the difference method, the present invention provides a novel Grünwald-Letnikov fractional-order derivative creep model by approximately rewriting the fractional-order derivative in the Grünwald-Letnikov form, which is beneficial to improving the accuracy of modeling.

[0051] 3) The creep model established by the present invention can be verified by uniaxial creep experiments, and the constitutive parameters can be obtained by fitting the experimental data. This experimental verification method not only enhances the reliability of the model, but also makes the acquisition of parameters more direct and accurate. At the same time, by changing the constitutive parameters, the influence of different parameters on the prediction results of the creep model can be conveniently studied, which provides convenience for parameter adjustment and optimization in practical applications.

[0052] 4) The creep model established by the present invention is not only applicable to materials such as low-fat chocolate, but can also be extended to other viscoelastic food materials with similar creep behaviors. The versatility and extensibility of this model make it have broad application prospects in the field of food engineering. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] In order to more clearly illustrate the implementation cases of the present invention or the technical solutions in the prior art, the following is a brief description of the drawings required for use in the embodiments. By referring to the drawings, the features and advantages of the present invention will be more clearly understood. The drawings are schematic and should not be understood as limiting the present invention in any way. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work. Among them:

[0054] Figure 1 A flow chart of a method for establishing a low-fat chocolate creep model using fractional order calculation provided in an embodiment;

[0055] Figure 2 Schematic diagram of a fractional derivative model in an embodiment;

[0056] Figure 3 It is a comparison diagram of creep test and simulation results at a stress level of 200Pa in the embodiment;

[0057] Figure 4 This is a creep prediction data diagram obtained by changing the model constitutive parameter υ in the embodiment;

[0058] Figure 5This is a creep prediction data diagram obtained by changing the model constitutive parameter β in the embodiment;

[0059] Figure 6 The block diagram of a system for establishing a low-fat chocolate creep model using fractional order calculation provided in an embodiment. DETAILED DESCRIPTION

[0060] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0061] Embodiment 1

[0062] like Figure 1 As shown, an embodiment of the present invention proposes a method for establishing a low-fat chocolate creep model using fractional order calculation, the method comprising:

[0063] Step S1: According to the creep power law behavior of low-fat chocolate in rheological deformation, a fractional derivative constitutive model based on the Grünwald-Letnikov form is established;

[0064] Step S2: Approximately rewrite the definition of the fractional derivative of the Grünwald-Letnikov form in the constitutive model, and establish a time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form;

[0065] Step S3: According to the fractional-order derivative constitutive equation based on the Grünwald-Letnikov form with time discretization, an improved Grünwald-Letnikov fractional-order derivative creep model is established to obtain an approximate discrete numerical solution for the creep compliance;

[0066] Step S4: performing a uniaxial creep experiment on the low-fat chocolate, fitting the experimental data using the established creep model, and obtaining constitutive parameters of the fractional derivative creep model;

[0067] Step S5: by changing the constitutive parameters of the fractional derivative creep model, numerically simulating the creep response of low-fat chocolate predicted by the fractional derivative creep model when the constitutive parameters are different, the influence of different constitutive parameters on the fractional derivative creep model is studied.

[0068] It can be seen from the above technical scheme that the present invention proposes a method for establishing a creep model of low-fat chocolate using fractional-order calculations. By using fractional-order derivatives in the form of Grünwald-Letnikov to establish a constitutive model, the creep power-law behavior of low-fat chocolate in rheological deformation can be more accurately described. By approximately rewriting the fractional-order derivatives in the form of Grünwald-Letnikov, a time-discretized constitutive equation is established, which greatly simplifies the numerical calculation process. The creep model established by the present invention can be verified by uniaxial creep experiments, and the constitutive parameters are obtained by fitting the experimental data. The creep model established by the present invention is not only applicable to materials such as low-fat chocolate, but can also be extended to other viscoelastic food materials with similar creep behaviors. The versatility and extensibility of this model make it have broad application prospects in the field of food engineering.

[0069] In this embodiment, in step S1, based on the creep power-law behavior of low-fat chocolate in rheological deformation, the fractional derivative model can be used to accurately describe the power-law response of chocolate with fewer constitutive parameters. Figure 2 The figure shows a fractional derivative model formed by connecting two Scott-Blair models in series. The fractional derivative constitutive model based on the Grünwald-Letnikov form is expressed as:

[0070]

[0071] In the formula, σ(t) and ε(t) refer to the stress and strain of the low-fat chocolate material, respectively; η and G refer to quasi-properties, that is, as intermediate parameters between elasticity and viscosity, they have both the elasticity of the spring and the viscosity of the damper; β and υ refer to fractional orders, that is, fractional derivative operators; It means to find the v-β derivative of σ(t), It means to take the v-th order derivative of ε(t).

[0072] The definition of fractional derivatives based on the Grünwald-Letnikov form is as follows:

[0073]

[0074] In the formula, refers to the fractional derivative based on the Grünwald-Letnikov form, dt refers to the time step; m is a variable; F(t) represents an arbitrary function; Refers to the recursive formula, written as and Γ(·) refers to the Eulerian Gamma function; for The value of is calculated as follows: when It should be noted that the coefficient in front of the function F(t) has a smaller and smaller influence as m gradually increases, and F(t) multiplied by the previous coefficient will eventually tend to 0. This calculation method reflects the memory effect of the fractional derivative model, has global correlation, and can accurately describe the deformation behavior of viscoelastic materials.

[0075] In this embodiment, in step S2, in order to obtain a numerical solution for the creep response of low-fat chocolate, the definition of the fractional derivative of the Grünwald-Letnikov form in the constitutive model is approximately rewritten, and a time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form is established, specifically including:

[0076] First, consider discretizing time t, that is, writing it as t = {t1, t2, …, t N-1 ,t N}, and t N =t N-1 +dt, formula (2) can be further rewritten as:

[0077]

[0078] According to formula (3), the time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form is given:

[0079]

[0080] Further expand formula (4):

[0081]

[0082] Reintegrating formula (5), we get:

[0083]

[0084] Where N represents the number of time t.

[0085] Formula (6) is the specific form of the time discretization based on the fractional derivative constitutive equation of Grünwald-Letnikov form. According to formula (6), solving the current value of ε is essentially to take the past historical displacement into account and weight ε N-m Then, the numerical approximate solution of ε in the discretized form in the time domain is obtained.

[0086] In this embodiment, in step S3, an improved Grünwald-Letnikov fractional-order derivative creep model is established according to the time discretized fractional-order derivative constitutive equation based on the Grünwald-Letnikov form, and an approximate discrete numerical solution for creep compliance is obtained, which specifically includes:

[0087] In order to analyze the creep behavior of low-fat chocolate, the creep compliance J(t) is first given:

[0088]

[0089] In the formula, σ0 refers to the constant stress, which is a constant. Therefore, the fractional derivative formula introducing the constant σ0 is:

[0090]

[0091] Then substitute formula (8) into formula (6) to obtain:

[0092]

[0093] Then, formula (9) is rearranged to obtain the improved Grünwald-Letnikov fractional-order derivative discretization strain formula:

[0094]

[0095] Then the improved Grünwald-Letniko fractional derivative creep model is derived:

[0096]

[0097] In this embodiment, in step S4, a uniaxial creep experiment is performed on the low-fat chocolate, and the experimental data are fitted using the established creep model to obtain the constitutive parameters of the fractional derivative creep model.

[0098] like Figure 3 The creep compliance test data obtained through uniaxial creep experiments are shown. The creep response of five low-fat chocolates 1, 2, 3, 4 and 5 with different oil gel contents are predicted by using the improved Grünwald-Letnikov fractional derivative creep model, and the creep simulation results are obtained. It can be seen from the fitting results shown in Table 1 below that the improved Grünwald-Letnikov fractional derivative creep model can accurately capture the power-law creep behavior of low-fat chocolate materials.

[0099] Table 1 Parameter fitting values ​​of chocolate samples 1, 2, 3, 4, 5 for uniaxial creep test

[0100]

[0101] In this embodiment, in step S5, by changing the constitutive parameters of the fractional derivative creep model, the creep response of the low-fat chocolate predicted by the fractional derivative creep model with different constitutive parameters is numerically simulated to study the influence of different constitutive parameters on the fractional derivative creep model.

[0102] Specifically, the material parameter η is given as 9×10 4 Pas υ and G = 3.2 × 10 4 Pas β , change the model constitutive parameters υ and β, take β = 0.2, and take different values ​​of υ, such as Figure 4 As shown. Take υ = 0.7, β takes different values, such as Figure 5 The creep response of the viscoelastic material low-fat chocolate predicted by the fractional-order derivative creep model with different constitutive parameters was numerically simulated to study the influence of different constitutive parameters on the fractional-order derivative creep model.

[0103] like Figure 4 and Figure 5 The figure shows that the creep compliance J(t) of chocolate material predicted by the fractional derivative creep model is obtained by changing the constitutive parameters υ and β in the fractional derivative creep model. By studying the influence of the constitutive parameters on the modeling of the fractional derivative creep model, the wide applicability of the model is judged. Since different chocolate materials have different constitutive parameters, it is proved that changing the parameter values ​​of the model can predict the creep response of different viscoelastic food materials.

[0104] The improved Grünwald-Letnikov fractional derivative creep model proposed in the present invention can accurately simulate the power-law creep response of low-fat chocolate, and can provide an important theoretical basis for the formulation and design of chocolate food.

[0105] Embodiment 2

[0106] like Figure 6 As shown, the present invention provides a system for establishing a low-fat chocolate creep model by using fractional order calculation, and the system is used to implement the method for establishing a low-fat chocolate creep model by using fractional order calculation in the above embodiment 1, specifically comprising:

[0107] The constitutive model establishment module 100 is used to approximately rewrite the definition of the fractional-order derivative in the Grünwald-Letnikov form in the constitutive model, and establish a time-discretized fractional-order derivative constitutive equation based on the Grünwald-Letnikov form;

[0108] The fractional derivative definition rewriting module 200 is used to approximately rewrite the definition of the fractional derivative in the Grünwald-Letnikov form in the constitutive model, and establish a time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form;

[0109] A creep model building module 300 is used to build an improved Grünwald-Letnikov fractional-order derivative creep model according to a time-discretized fractional-order derivative constitutive equation based on the Grünwald-Letnikov form, and obtain an approximate discrete numerical solution for creep compliance;

[0110] The experimental verification and parameter acquisition module 400 is used to perform a uniaxial creep experiment on low-fat chocolate, fit the experimental data using the established creep model, and obtain the constitutive parameters of the fractional derivative creep model.

[0111] The numerical simulation and parameter analysis module 500 is used to change the constitutive parameters of the fractional derivative creep model, numerically simulate the creep response of low-fat chocolate predicted by the fractional derivative creep model when the constitutive parameters are different, and study the influence of different constitutive parameters on the fractional derivative creep model.

[0112] The system of this embodiment for establishing a low-fat chocolate creep model by using fractional-order calculation is used to implement the aforementioned method for establishing a low-fat chocolate creep model by using fractional-order calculation. Therefore, the specific implementation method of the system for establishing a low-fat chocolate creep model by using fractional-order calculation can be seen in the embodiment of the method for establishing a low-fat chocolate creep model by using fractional-order calculation. For example, the constitutive model establishment module 100, the fractional-order derivative definition rewriting module 200, the creep model establishment module 300, the experimental verification and parameter acquisition module 400, and the numerical simulation and parameter analysis module 500 are respectively used to implement steps S1, S2, S3, S4, and S5 in the aforementioned method for establishing a low-fat chocolate creep model by using fractional-order calculation. Therefore, its specific implementation method can refer to the description of the corresponding embodiments of each part. In order to avoid redundancy, it will not be repeated here.

[0113] Embodiment 3

[0114] An embodiment of the present invention also provides an electronic device, which includes a processor, a memory and a bus system, wherein the processor and the memory are connected via the bus system, the memory is used to store instructions, and the processor is used to execute the instructions stored in the memory to implement the above-mentioned method of establishing a low-fat chocolate creep model using fractional-order calculation.

[0115] Embodiment 4

[0116] An embodiment of the present invention further provides a computer storage medium storing a computer software product. The computer software product includes several instructions for enabling a computer device to execute the above-mentioned method of establishing a low-fat chocolate creep model using fractional-order calculation.

[0117] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0118] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0119] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable device provide for implementing the process in the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0120] Obviously, the above embodiments are merely examples for clear explanation and are not intended to limit the implementation methods. For those skilled in the art, other different forms of changes or modifications can be made based on the above description. It is not necessary and impossible to list all the implementation methods here. The obvious changes or modifications derived from these are still within the protection scope of the invention.

Claims

1. A method for establishing a low-fat chocolate creep model using fractional order calculation, characterized in that: include: Step S1: According to the creep power law behavior of low-fat chocolate in rheological deformation, a fractional-order derivative constitutive model based on the Grünwald-Letnikov form is established, wherein the fractional-order derivative constitutive model based on the Grünwald-Letnikov form is expressed as: In the formula, σ(t) and ε(t) refer to the stress and strain of the low-fat chocolate material, respectively; η and G refer to quasi-properties, that is, as intermediate parameters between elasticity and viscosity, they have both the elasticity of the spring and the viscosity of the damper; β and υ refer to fractional orders, that is, fractional derivative operators; It means to find the υ-β derivative of σ(t), Refers to the υ-order derivative of ε(t); The definition of fractional derivatives based on the Grünwald-Letnikov form is as follows: In the formula, refers to the fractional derivative based on the Grünwald-Letnikov form, dt refers to the time step; m is a variable; F(t) represents an arbitrary function; Refers to the recursive formula, written as and Γ(·) refers to the Eulerian Gamma function; for The value of is calculated as follows: When m→∞, Step S2: Approximately rewrite the definition of the fractional derivative in the Grünwald-Letnikov form in the constitutive model, and establish a time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form, specifically including: First, consider discretizing time t, that is, writing it as t = {t1, t2, …, t N-1 ,t N }, and t N =t N-1 +dt, formula (2) can be further rewritten as: According to formula (3), the time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form is given: Further expand formula (4): Reintegrating formula (5), we get: Where N represents the number of time t; Step S3: According to the time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form, an improved Grünwald-Letnikov fractional derivative creep model is established to obtain an approximate discrete numerical solution for the creep compliance, which specifically includes: In order to analyze the creep behavior of low-fat chocolate, the creep compliance J(t) is first given: Then the fractional derivative formula of the constant σ0 is introduced as follows: Then substitute formula (8) into formula (6) to obtain: Then, formula (9) is rearranged to obtain the improved Grünwald-Letnikov fractional-order derivative discretization strain formula: Then the improved Grünwald-Letniko fractional derivative creep model is derived: In the formula, σ0 refers to the constant stress and is a constant; Step S4: performing a uniaxial creep experiment on the low-fat chocolate, fitting the experimental data using the established creep model, and obtaining constitutive parameters of the fractional derivative creep model.

2. The method for establishing a low-fat chocolate creep model using fractional order calculation according to claim 1, characterized in that: The method further comprises: Step S5: by changing the constitutive parameters of the fractional derivative creep model, numerically simulating the creep response of low-fat chocolate predicted by the fractional derivative creep model when the constitutive parameters are different, the influence of different constitutive parameters on the fractional derivative creep model is studied.

3. A system for establishing a creep model of low-fat chocolate using fractional order calculation, characterized in that: The system is used to implement the method for establishing a low-fat chocolate creep model by using fractional order calculation as described in any one of claims 1 to 2, specifically comprising: The constitutive model building module is used to establish a fractional derivative constitutive model based on the Grünwald-Letnikov form according to the creep power law behavior of low-fat chocolate in rheological deformation; The module for rewriting the definition of fractional derivatives is used to approximately rewrite the definition of fractional derivatives in the Grünwald-Letnikov form in the constitutive model and establish a time-discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form. The creep model building module is used to establish an improved Grünwald-Letnikov fractional derivative creep model based on the time discretized fractional derivative constitutive equation based on the Grünwald-Letnikov form, and obtain an approximate discrete numerical solution for the creep compliance; The experimental verification and parameter acquisition module is used to conduct uniaxial creep experiments on low-fat chocolate, fit the experimental data using the established creep model, and obtain the constitutive parameters of the fractional derivative creep model.

4. The system for establishing a low-fat chocolate creep model using fractional order calculation according to claim 3, characterized in that: The system further comprises: The numerical simulation and parameter analysis module is used to change the constitutive parameters of the fractional derivative creep model, numerically simulate the creep response of low-fat chocolate predicted by the fractional derivative creep model when the constitutive parameters are different, and study the influence of different constitutive parameters on the fractional derivative creep model.

5. An electronic device, characterized in that: The electronic device includes a processor, a memory and a bus system, the processor and the memory are connected through the bus system, the memory is used to store instructions, and the processor is used to execute the instructions stored in the memory to implement the method for establishing a low-fat chocolate creep model using fractional order calculation as described in any one of claims 1 to 2.

6. A computer storage medium, characterized in that: The computer storage medium stores a computer software product, which includes several instructions for enabling a computer device to execute the method for establishing a low-fat chocolate creep model using fractional-order calculation as described in any one of claims 1 to 2.

Images