Constitutive Relation of Engineering Material Based on SIR Model and HAM

As an epidemic mathematical model, the SIR model represents the transition of the Susceptible, Infected, and Recovered. The profound implication of the SIR model is viewed as the propagation and dynamic evolutionary process of the different internal components and the characteristics in a complex system subject to external effect. The uniaxial stress-strain curve of engineering material represents the basic constitutive relation, which also represents the damage propagation in the units of the damaged member. Hence, a novel dynamic stress-strain model is established based on the SIR model. The analytical solution and the approximate solution for the proposed model are represented according to the homotopy analysis method (HAM), and the relationship of the solution and the size effect and the strain rate is discussed. In addition, an experiment on the size effect of confined concrete is carried out and the solution of SIR model is suitable for simulation. The results show that the mechanical mechanism of the parameters of the uniaxial stress-strain model proposed in this paper reflects the actual characteristics of the materials.The solution of the SIRmodel can fully and accurately show the change of the mechanical performance and the influence of the size effect and the strain rate.


Introduction
Engineering material means the material used for engineering or the materials used to produce other materials which may be used in engineering.The traditional material in civil engineering includes concrete, steel, soil, wood, and alloy.These materials have different qualities about strength, workability, durability, and resistance against corrosion [1].The materials differ also in their structure, texture, and performance.The strength, toughness, and ductility of construction material are the fundamental guarantee for engineering reliability.With the development of civil engineering, the type and form of civil engineering materials are unceasingly rich and the performance is enhanced also, and the corresponding research and application on experimental techniques, theoretical analysis, numerical simulation, and actual engineering are improved constantly.However, there are still more phenomena to be interpreted and more research issues should be given attention because the engineering material has microcrack, heterogeneity, and evolution in micromechanics and shows elastoplasticity, discreteness, and randomness in the micromechanics.In addition, the mechanical behaviors of engineering materials are also related to the size, the loading rate, the external environment, and other factors.Therefore, the mechanical properties of different engineering materials have obvious differences; the demand about the experiment and numerical simulation for various materials are diverse.Researchers normally investigate the mechanic property of the specified material, and the study on the generalized characteristics for various construction materials is still not intense.
The uniaxial stress-strain curve indicates the simple mechanic property of the material subjected to pure compress or pure tension; representing the basic material constitutive relation and the uniaxial stress-strain relation is also the foundation to establish multiple compound constitutive models.In the elastoplasticity or nonlinear analysis, the uniaxial stress-strain curve plays a decisive role in the accuracy of the numerical results.Furthermore, the computational analysis of the engineering structure subjected to dynamic or cyclic loadings requires the stress-strain models to simulate the response of the structure [2,3].

Journal of Applied Mathematics
There are many types of stress-strain curves models for different material or even one kind of material.As known to all, the stress-strain curve of concrete has a variety of forms or functions because concrete is a type of composite material and the discreteness is obvious.Numerous concrete models have been proposed in the last years.For example, the compressive stress-strain function for concrete consists of five kinds of equations as polynomial, exponential, and trigonometric function, rational fraction, and sectional form and the total number of the corresponding formulas is over 20 [2].
In the macroscopic level, three broad categories can be distinguished: models derived from the theory of elasticity, models based on the theory of plasticity, and models based on the continuum damage theory [3,4].Also, some coupled models based on the association of plasticity and continuum damage theory have been recently developed.Although it has been proved that the models derived from theory of plasticity and continuum damage theory can accurately simulate the observed behavior of concrete, the engineering application of these models is less.This is motivated by the great amount of parameters that are usually needed and the difficulty to obtain them through conventional laboratory tests.From the perspective of another point of view, the structural member with definite shape and materials is a complex system which has huge amounts of units, and the units have the whole process of sustained force, damage, interact, transmitting energy, and propagation, until failure.The characteristics of all the units generally constitute the property of the macromaterial and member.
Hence, taking the damage propagation and transmit in the microunits into consideration, establishing a new and generalized uniaxial constitutive model which can take into account that both the versatility and the dirigibility is necessary.

SIR Model and Its Connection between Constitutive Relations of Material
An epidemic represents the sudden outbreak and propagation of a disease, often occurring on a short temporal scare and affecting a significant portion of a population.Epidemics also may exhibit some periodic behavior, as opposed to endemics, which are diseases that are always present to some extent in a population.Epidemiology is the branch of science which essentially deals with the mathematical modeling of propagation of diseases.The first mathematical model of epidemiology was formulated and solved by Daniel Bernoulli in 1760.Kermack and McKendrick [5] illustrated that diseases showed a threshold type of behavior.In other words, if a single person infected by a particular disease passed on the infection to more than one person in turn, an epidemic would occur, while if less than one secondary infection occurred in each primary one, the disease would die out.The study of mathematical epidemiology has grown rapidly, with a large variety of models having been formulated and applied to infectious diseases.The epidemic model can describe the dynamic process of the epidemic, and spread characteristics will be grasped based on the information of the change of population; then the epidemic can be effectively controlled.
The general epidemic model belongs to first order ordinary differential equations, and the most representative is the SIR model.
The SIR epidemic mathematical model is presented based on a system of first order ordinary differential equations and it has been used in the modeling of several infectious diseases, where the parameters need to be estimated by epidemiological data [6].In this model, the variables represent subpopulations of the Susceptible (S) who can catch the disease, the Infected (I) who are infected and can transmit the disease to the Susceptible, and the Removal (R) who had the disease and recovered or died or have developed immunity or have been removed from contact with the other classes.Thus, the model describes the propagation and transformation between different classes.Assume that the total number of population is one unit for the sake of simplicity and the SIR model is written as where () and () denote the Infected and the Susceptible, respectively, and () is Removal and  is time.Hence, () + () + () = 1 is satisfied. is the infectivity coefficient of the typical Lotka-Volterra interaction term, which means the daily contact rate (i.e., the population that each patient effectively contact with the healthy persons every day). is the daily recovery rate or number of patients cured or removed out accounting for the proportion of the total patients. is the increasing rate of Susceptible persons when patients are increasing. is the speed parameter after taking preventive and control measures.The SIR model is subject to the initial conditions (0) =  0 and (0) =  0 , where  0 > 0 and  0 > 0 are given constants.From the introduction above, it is significant that the SIR model is not simply an equation representing the spread of an epidemic, while it actually describes the propagation and dynamic evolution process of the different internal components and the characteristics in a complex system subject to external effect.In general, the SIR model embodies the general characteristics and rules in the similar system with propagation and transformation.The continuous damage and destruction will occur for the material specimen in civil engineering under increasing load, and the unstressed units, the stressed units, and the fractured or invalid units in anyone state should be included.The stress-strain relationship is the typical representation of the propagation and the dynamic evolution.Therefore, an innovative model for the material can be presented referencing to the SIR model, and the two models have similar parameters.To validate this interpretation, a compressed concrete specimen with 60 mm diameter and 120 mm height is selected as an example.The uniaxial failure process of the specimen is simulated based on the three-dimensional mesolevel finite element method [7,8], and the failure charts in different stage are shown in Figure 1.It is clear that there are unstressed units, stressed units, and the fractured or invalid units in all stages, whereas the damage propagation and the proportion of each part are time-variant.The unstressed units are majority and fractured units are few at the initial stage.With the compression being enhanced, the units participate in the resistance gradually and the stressed units boom.After the counteragent in most units decreases, the invalid units increase and the cracks appear, though the stressed units decline and the specimen collapses.In conclusion, the SIR model is suitable for illustrating the failure process of materials, that is, typical damage propagation process in a complex system.
Strain  reflects the deformation capability and damage process of material, corresponding with the time variable  in the SIR model.(), (), and () represent the equivalent stress in unstressed units, the equivalent stress in stressed units, and the equivalent stress in fractured or invalid units, respectively.This model subjects to the initial conditions (0) =  0 and (0) =  0 , where  0 > 0 and  0 > 0 are given constants.In the same manner, () + () + () = 1.In this case,  can be viewed as the transmissibility rate in units, and  represents the failure rate in units. can be expressed as the increase rate of stressed units when the invalid units grow, and  is the reduced speed parameter for invalid units due to stress transmission and distribution.Finally, the stress-strain model of material can be expressed as Equation ( 2) is nonlinear differential equations which represent the coupling and transformation of the three types of units and the full complex process for the material from intact to failure.The system of equations has the unique solution, but the analytical solution cannot be obtained by normal mathematical methods and the results are calculated usually by numerical methods.
Assume the basic condition of the equation above is that  is 1.00,  is 0.15,  is 0.1,  is 0.8, and (0) is 0.02.The consequences of the equivalent stress () under different parameters are solved by numerical method, as shown in Figure 2. It is obvious that the full curves accurately indicate the linear ascent stage, nonlinear yield stage, and the decline stage or strengthening stage after peak stress.Further, the variation of , , or  can obviously affect peak and ductility of the stress-strain curve, and different forms of the curve can be obtained according to the overall adjustment of the corresponding parameters.
Different solutions of SIR model in various conditions can represent the stress-strain curves of diverse types of materials.Figure 3 reveals the stress-strain curves including concrete or rock-soil in compression or tension as well as different types of stress-strain curves of soils and steel obtained by numerical methods.The comparison curves of the real experimental data and numerical solution about the concrete and steel are shown in Figures 4 and 5, respectively.In order to compare the regularity but not the specific data, the peak strain is normalized.By changing the parameter appropriately, the full stress-strain curve obtained from (2) can simulate the peak stress, ductility, and softening stage of concrete with different strength, as shown in Figure 6 relative significance; it can be adjusted to a proportional coefficient to real physical value in practical applications.

Analytical Solution and Approximate Solution of SIR Model
The numerical solution of the SIR model is accurate but it must be calculated for any specific case, not representing the panorama and rule of the solution.Therefore, the analytical solution of the SIR model is necessary and it has important mathematical and mechanical meaning.Many corresponding studies of the analytical solution of the SIR model have been carried out, but the problem is not solved perfectly because the model has strong nonlinearity.There are some classic nonlinearity techniques that can be considered, such as Adomian decomposition method [9,10], delta expansion method, and perturbation method [11,12].All these methods have some limitations, such that these techniques do not provide a convenient way to adjust and control the convergence region and rate of approximation series.
To overcome the mentioned limitations, Liao has proposed the homotopy analysis method (HAM) for nonlinear problems and then modified it step by step [13], and many nonlinear problems in different fields have been successfully solved by HAM.
Many problems such as boundary layer similarity solution for forced, natural, and mix convection in porous medium, heat transfer, and fluid mechanic problems are nonlinear inherently.The analytical solutions of nonlinear ordinary differential equations can be solved by HAM [14][15][16].
HAM is based on homotopy, a concept in topology.The notion of equivalent maps or processes, where one can be deformed into the other, is the fundamental structure of homotopy.Two maps ,  :  →  of topological spaces are homotopic if there exists a map  :  ×  →  such that (, 0) = () and (, 1) = () for  ∈ .Here,  ×  denotes the product of  with unit interval [0, 1] of real numbers.The map  is called the homotopy between  and .The basic idea of HAM is shown in the nonlinear differential equation as follows: where  is nonlinear operator, ⃗  is a vector of spatial variables,  denotes time, and ( ⃗ , ) is an unknown function.Boundary or initial conditions can be treated in a similar manner, which we avoid here just for simplicity [17].
Generalizing the concept of traditional homotopy, the socalled zero-order deformation equation is established as where  ∈ [0, 1] is the embedding parameter, ( ⃗ , ; ) is an unknown function, ( ⃗ , ) is an auxiliary function,  is a linear operator,  0 () is the initial guess, and ℎ is a convergence-control parameter.As  increases from 0 to 1, ( ⃗ , ; ) vary from initial trial  0 () to the exact solution ( ⃗ , ).If this variation is smooth enough, we construct the Maclaurin series of ( ⃗ , ; ) at  = 0 and the coefficients of all the higher terms can be obtained from the higherorder deformation equations.Differentiating the zerothorder deformation equation above  times with respect to the embedding parameter , then setting  = 0, and finally dividing by !, we have the so-called th-order deformation equations: with   = (  [( ⃗ , ; )]/(! ))| =0 = !⋅   ().In this way, a nonlinear equation is transformed into a series of linear equations.The exact solution of ( ⃗ , ) is then approximated by the summation of the Maclaurin series at  = 1.One has great freedom for the choice of linear operator .For example, one can even select linear operator of different order as compared to the original nonlinear problem.To simplify the applications of HAM, Liao has suggested some rules, that is, the rule of solution expression, the rule of solution existence, and the rule of ergodicity for coefficients of homotopy series solution [13].
In HAM, there is sufficient space for the convergence of the approximation by the introduction of ℎ which greatly improves the early homotopy analysis method.It provides us with a simple way to ensure the convergence of the series solutions of nonlinear problems.This is an obvious advantage of HAM over homotopy perturbation method (HPM).In fact, HPM is just a special case of HAM when ℎ = −1, as pointed out by Liao and Abbasbandy and, in general, proved by Sajid and Hayat [18,19].Liao also has pointed out that HAM logically contains other nonperturbation methods such as Adomian's decomposition method [9], the -expansion method, and Lyapunov's artificial small parameter method [12].
By choosing the auxiliary linear operators and from deformation equations and deformation derivative condition, we have where  , and  , are coefficients.Substituting the above expressions into (10) to (11), the recurrence formulas are where  and  are parameters to meet deformation equation.Finally, we get analytic solution of expressions in (2): At this point,  =  − (∞) ≈ .According to the above conclusions and characteristics of solution, the approximate analytical solution of ( 17) can be expressed as where  belongs to even number and   =  +1 when  is an odd number.The above approximate analytical solution can usually achieve adequate satisfactory results when  is 2 or 4. As a simplified format, the equation above can also be transformed as follows: The comparison of the real concrete strain-stress curve, the numerical solution curve, and the approximate analytical solution of the SIR model is shown as Figure 4, and the comparison of steel is shown in Figure 5.It can be seen that both the numerical solution and the approximate analytical solution fit the original value precisely.It is worth noting that the approximate analytical solution reflects the connotation and characteristics of the analytical solutions, but it nevertheless has a certain similarity, and the specific value of   in the solution should be adjusted appropriately according to actual condition and the numerical solution in order to obtain better accuracy.In addition, Umemura and Aoyama [20] have presented an exponential constitutive relation for concrete, which is similar with the model as (18) in  equals 2. However, the original exponential relation is determined by experimental data fitting, and the proposed solution in this paper has the theoretical basis and generality; the results can be verified with each other.The effect of various parameters in (18) on solutions is shown in Figure 7.
The stress-strain curve of compression concrete with different strength is shown as Figure 8, and the simulation curves are similar with the experimental data in shapes and the variation rules, indicating that the approximate analytical solution put forward in this paper can also embody the intrinsic characteristics and variation of the SIR model.

Size Effect and Strain Rate Based on SIR Solution
For real engineering materials, the mechanical properties also related to size, load mode, and external environment besides their own composition and characteristics and mainly include size effects, strain rate, and multiaxis loading.In this paper, the uniaxial stress-strain curve on the influence of size effect and strain rate of are mainly discussed.The size effect of material strength refers to phenomenon that large size member in brittle material usually fractures under a lower nominal stress than geometrically similar small-size member (the nominal stress being defined as the load divided by the characteristic cross section area).Early researchers suppose that any observed size effect should be described by extreme value statistics prevailed in structural engineering.In the mid 1970s, the fact that there exists a purely deterministic size effect, caused by energy release associated with stress redistribution prior to failure and that this energetic size effect usually dominates in the so-called quasibrittle structures (i.e., structures in which fracture propagation is preceded by a relatively large fracture process zone which, in contrast to brittle ductile fracture of metals, exhibits almost no plastic deformations but undergoes progressive softening due to microcracks) gradually emerged.Bažant and Chen [21] and Bažant and Planas [22] summarized six main causations for size effect which include boundary layer, diffusion phenomena, hydration heat, randomness of material strength, energy release, and the fractal character of the crack surface.Recent research focuses on three main types of size effects, namely, the statistical size effect due to randomness of strength, the energy release size effect, and the possible size effect due to fractality of fracture or microcracks.
Recent research has shown that the elastic modulus and the peak stress of brittle material will gradually decrease with the increase of the member size and depth-width ratio; at the same time, the strain value at the peak stress changes a little, and the descent rate of stress and fragility at the softening section will also reduce [23], as shown in Figures 2 and 7.In some cases, the curves of different sizes in softening section can even intersect as shown in Figures 6 and 8.
In fact, for the large size member with the same axial force, the total number and size of microcracks are bigger and the domain that the sustained damage occurs is larger.Hence, the equivalent strain energy for damage needs fewer paths and the damage propagation needs shorter paths, which lead to the probability of regional brittle failure; that is, the failure rate in unit and its growth rate both increase.
Considering the above characteristic of size effect, the SIR model can simulate the size effect of brittle material and the whole process of failure.The advantage of the application of the SIR model is the physical significance being explicit and the variation can be realized by enhancing the failure rate in unit  and the reduced speed parameter for invalid elements due to stress transmission and distribution , referring to Figure 2.
The materials such as concrete and rocks are typical rate sensitive for their strength, ductility, and failure mode will significantly change in different strain rates.Existing research shows that the elastic modulus and the ultimate strength of rate sensitive material will enhance with increasing strain rate, and the strain value at the peak stress changes a little [23].In different strain rates, the stress-strain full curve of concrete is basically consistent with the whole curve under static load in shape; however, the ductility increases slightly, as shown in Figures 2 and 7.
For the materials subjected to the same axial force, if the strain rate increased, the internal microcracks in the mortar substrate are late to fully extent, but the quantity and degree of the damaged aggregate relatively increase which lead to the enhancement of the failure strength.At the same time, the unit failure rate increased, and the increase rate of the stressed unit sustainably grows.
Therefore, the SIR model can show the influence of the strain rate from the physical sense by properly adjusting parameters, and the simulation is realized by enhancing the failure rate  and the increase rate of stressed units when the invalid units grow , referring to Figure 2.
The above discussion illustrates the direct relationship between size effect, strain rate, and the SIR model with its numerical solution from mechanical principle and propagation characteristic.In practical applications, the approximate analytic solution of the model is more convenient, so the study below will discuss the adjustment of coefficient   in the approximate analytic solution for the purpose of representing size effect and strain rate.
The most widely used theory is the size effect law proposed by Bažant and Chen [21] based on a large number of experiments, and this method put forward the size effect unified formula in certain extent according to plasticity theory or elastic theory: where   is the nominal stress when material is damaged,  is dimensionless parameter,   is the strength of the quasi brittle material,  is the characteristic length of the structure, and  0 is the constant related to structural shape.It can be seen from Figure 6 that the consideration about size effect can be realized by modulating   and   properly in the approximate analytical solution of SIR.
In the usual study, the strain rate effect coefficient of compression concrete is described as the following exponential form according to European standard [24]: where ε is the current material strain rate, ε 0 is the quasi static strain rate, taken as 3 × 10 −5 /s, and  cs and  cd are the static and dynamic compression strength of concrete, respectively.One has  = 1/(5 + 0.9  cs ) and log   = 6.156  − 2.
For tensioned concrete where ε is the current material strain rate and ε 0 is the quasi static strain rate, taken as 3 × 10 −5 /s, and  cs and  cd are the static and dynamic tension strength of concrete, respectively.One has  = 1/(10 + 0.6)  ts and log  = 7.11 − 2.33.
The strain rate effect coefficient of steel I is as follows: where ε is the current material strain rate and ε 0 is the quasi static strain rate, taken as 3 × 10 −4 /s, and  ys and  yd are the static and dynamic strength, respectively. 1 and  2 are the test parameters obtained by regression analysis method [25,26].
Referring to Figure 7, the strain rate effects on concrete can be simulated by modulating   and   properly for the approximate analytical solution of SIR model.In addition, if   is taken as a small negative number when  is an odd, the characteristics of strengthening after yielding for metal materials can also be simulated, and the stress-strain curve of metal with high strain rate or high strength can be obtained by enlarging the absolute value of   or   , as shown in Figure 9.
In conclusion, it is significant that the coefficient   of the approximate analytical solution implies the multiple parameters interaction in SIR model, and the size effect and strain rate in the whole load process can be achieved by regulating   .The factor   on size effect can be achieved from (19) and the factor   on strain rate can be achieved from (21) to (23).Thus, the correction coefficient about   is where   and   are adjustment coefficient on size effect and strain rate, respectively, and the more precise value should be determined by experimental data.

Test on the Size Effect of Confined Concrete
As mentioned above, the uniaxial construction relationship of single material based on SIR model and HAM method is presented.However, the normal member of civil engineering consists of both concrete and steel, and the constructive relationship of reinforced concrete and confined concrete and corresponding properties such as size effect need more research.
In order to verify the feasibility of the construction relationship for confined concrete based on the SIR model, six reinforced concrete prism specimens confined by square stirrups were made.Each type of specimen contains two same members, and the volumetric percentage of stirrups is 1.26% for all.The configurations of stirrups are shown in Figure 10.The design parameters of the specimens are listed in Table 1.For the concrete used, the normal prismatic compressive strength is 42.67 N/mm 2 , the ultimate compression strain is 0.0022, and Young's modulus is 3.08 × 10 4 N/mm 2 .For the steel bars, the average yield strength is 480 N/mm 2 , the ultimate strength is 665 N/mm 2 , and Young's modulus is 2.05 × 10 5 N/mm 2 .
The electrohydraulic servo testing machine with 4 × 10 4 kN maximum range was used as the load device and the continuous axial monotonic load was applied, as shown in Figure 11.When the force is less than the ultimate bearing    compression is measured by the displacement meter with 200 mm measuring range and the gauge length is 2/3 of the specimen height.The average stress-strain relation of cross section curves of the specimens with different sizes is shown in Figure 13.The abscissa  is obtained by the measured displacement divided by the respective gauge length.Ordinate  is the relative stress, equal to the axial load value divided by the cross-sectional area of the specimen.According to the stress-strain relation of the specimens with different sizes, the approximate solution based on nonlinear data fitting and ( 19) is presented as shown in Figure 14.The failure rate  is assumed as 0.1, and the parameters in (19) are listed in Table 2.It is significant that the curves of the approximate solution fit well with the original data and the variation trend and the variance rule among the curves are clearly revealed.The SIR model and the corresponding approximate solution are suitable for confined concrete and complex composite material with size effect.

Conclusion
Though there are various constitutive relation models of engineering material at present, the united constitutive relation model is rare and it is necessary to establish a model which can present many properties in different materials by concise form.
Under the external load, the member composed of engineering materials is damaged inevitably and the damage (a) Initial stage (b) Stressed stage (c) Failure stage Unstressed units (susceptible) Stressed units (infected)Fractured or invalid units (removal) Cracks

Figure 1 :
Figure 1: Damage propagation process and different units based on microlevel simulation.

Figure 2 :
Figure 2: Numerical solutions of SIR for different parameter.

Figure 3 :
Figure 3: Constitutive curves for different engineering material based on SIR model.

Figure 6 :
Figure 6: Constitutive curves for concrete based on SIR model.

Figure 7 :
Figure 7: Approximate solutions based on SIR model with different parameter.

Figure 8 :
Figure 8: Constitutive relation of concrete with different strength based on approximate solutions.

2 )Figure 9 :Figure 10 :
Figure 9: Constitutive relation of steels with different strain rate based on approximate solutions.

Figure 12 :
Figure 12: Final failure mode of specimens.

Table 1 :Figure 13 :
Figure 13: Stress-strain curve of gross area for all prism specimens.

Table 2 :
Parameters of the approximate solution.