A Numerical Approach to Static Deflection Analysis of an Infinite Beam on a Nonlinear Elastic Foundation : One-Way Spring Model

A numerical procedure proposed by Jang et al. (2011) is applied for the numerical analyzing of static deflection of an infinite beam on a nonlinear elastic foundation. And one-way spring model is used for the modeling of fully nonlinear elastic foundation. The nonlinear procedure involves Green’s function technique and an iterative method using the pseudo spring coefficient. The workability of the numerical procedure is demonstrated through showing the validity of the solution and the convergence test with some external loads.


Introduction
Accurate modeling of nonlinear deflection of an infinite beam on a nonlinear elastic foundation is crucial for material and structural engineering.The research can be applied to strength analysis and practical engineering design application, say, to curved plate manufacturing.Therefore, many theoretical and experimental studies have been carried out on the nonlinear modeling of an infinite beam on a nonlinear elastic foundation.
The closed-form solutions for the static and dynamic response of a uniform beam resting on a linear elastic foundation can be found in several references [1][2][3].Timoshenko [4], Kenney [5], Saito and Murakami [6], and Frýba [7] formulated a closed-form solution using Green's function approach based on a linear assumption.Beaufait and Hoadley [8], Massalas [9], Lakshmanan [10], and Hui [11] proposed the static, dynamic, and elastic stability analysis of a beam resting on a nonlinear elastic foundation.And there are many researches concerning the linear elastic foundation [12][13][14][15].Among the references, Beaufait and Hoadley [8] approximated the relationship of the stress-strain curve to be hyperbolic, but they approximated the bilinear curve to handle the nonlinear problem.The applied nonlinear foundation is active only when the beam is pressing against the foundation, and it is assumed to be inactive in the regions where the beam has been displaced away from the foundation.Soldatos and Selvadurai [16] also applied the hyperbolic-type nonlinear elastic foundation to analyze finite or infinite beams.Lee et al. [17][18][19] developed the exact and semiexact analysis of a nonuniform beam with general elastic end-restraints.Kuo and Lee [20] derived the static deflection of a general elastically end-restrained, nonuniform beam on a nonlinear elastic foundation under axial and transverse forces.
Recently, Jang et al. [21] proposed a new method for assessing the nonlinear deflection of an infinite beam on a nonlinear elastic foundation.They approach the high nonlinear problems using Green's function technique with an iterative method.Jang and Sung [22] proposed a new functional iterative method for static beam deflection, which has a variable cross-section.Choi and Jang [23] proved the existence and uniqueness of the nonlinear deflections of an infinite beam resting on a nonlinear elastic foundation using the Banach fixed point theorem.Jang [24] also proposed a new iterative method for the large deflection of an infinite beam resting on an elastic foundation based on the v. Karman approximation of geometrical nonlinearity.From the existing 0 Deflection  Nonlinear spring force literature, a number of studies have analyzed a beam on an elastic foundation; however, they just use linear plus a nonlinear term of spring force, that is, linear-cubic model.And they are related to the static analysis of nonuniform beams which is resting on a nonlinear elastic foundation, and the recovered solution is not accurate or has many limits.Few studies have fully adopted the nonlinear elastic foundation model, whose spring force is based on one-way spring model, as shown in Figure 1.In the real world, at the steady state, the soil or foundation would not be raised, or they are separable (in Figures 1 and 2).A nonlinear spring force exists when the infinite beam deflects downward but does not exist in case of the other cases.

One way Two way
Although there are many researches, fully nonlinear elastic foundation was not considered.Beaufait and Hoadley [8] and Soldatos and Selvadurai [16] approximated the stressstrain relationship as a bilinear curve.In this paper, one-way spring model is successfully used to examine the real nonlinear elastic foundation, and the nonlinear iterative method proposed by Jang et al. [21] is applied.Some numerical experiments are carried out to report the accuracy of the method, and the convergence of the solution is investigated according to several physical properties of the system.

Euler-Bernoulli's Beam on a Nonlinear Elastic Foundation:
A Nonlinear Spring Model.In this paper, the nonlinear spring force is fully analyzed by the one-way spring model instead of the conventional mathematical form of the two-way spring model [25].
The well-known classical Euler-Bernoulli's beam theory is considered for the solution procedure which is a simplification of elasticity which provides a means of calculating the load-carrying and deflection characteristics of beams.The governing equation for the linear deflection of an infinite beam on an elastic foundation that satisfies the fourth-order differential equation is as follows (the weight of the beam is neglected): And the reaction force (), and , , , () and () are Young's modulus, the mass moment of inertia, a linear spring coefficient, a nonlinear part of spring force, and external load, respectively.The boundary and Therefore, ( 1) and ( 3) together form a well-defined boundary value problem.Timoshenko [4], Kenney [5], Saito and Murakami [6], and Frýba [7] derived the general linear solutions neglecting the nonlinear part, (), in (2) as follows: where Green's function  can be defined as  1.
The loading condition is assumed to be localized, so () in (4) satisfies the boundary conditions in (3).Deriving the linear solution in (4), a uniformly upward, nonlinear spring force depends on the beam deflection .
In this study, the main idea of the present study is proposed by Jang et al. [21].The nonlinear spring model is used for the formulation of a realistic nonlinear elastic foundation.A pseudo linear spring coefficient   and a (real) spring force () are used.Therefore, the fourth-order differential equation in ( 1) is equivalent to the following equation: where the nonlinear spring force  depends on the deflection : or Table 1: Three cases of the exact solution.
Case Exact solution () In (8),    is the pseudolinear spring force term and is finally compensated, so it does not affect the nonlinear solution.Therefore, (4) shows that (8) must be equivalent to the following equation: where Green's function with pseudo linear spring coefficient can be expressed as Regarding (8), the nonlinear relation for  can be derived as follows: where the function  can be written as Equation ( 11) is a nonlinear Fredholm integral equation of the second kind for .
To examine the nonlinear iterative procedure in ( 13) by simulation, it should be discretized as follows: where   denotes the weights for the integration rule.The number  in the summation of ( 14) denotes the total segments of the interval (−, ), and  is a sufficiently large value satisfying (3).

Numerical Experiments
In this section, numerical experiments are performed to determine the validity of the iterative method.This study assumes a nonlinear spring force () in (7) and examines how accurate the solution converges to an exact solution.The convergence of the method is investigated with some external loading conditions.

Nonlinear Spring Model.
For simplicity, the present study considers an infinite beam on a nonlinear elastic foundation, whose spring force is derived as follows: where () in ( 7) is chosen immediately as a cubic form

Comparison with Exact Solution.
To determine if the iterative method converges to an exact solution, the 3 cases of exact solutions listed in Table 1 are first assumed and are illustrated in Figure 3.These cases are chosen to be infinitely differentiable and satisfy the conditions in (3).The following are assumed:  =  =  = 1,   = 3, and  = 0.2, and the initial guess of the deflection is  0 = 0. Simpson's integration rule is applied to the numerical integration.Three cases of external loads, shown in Figure 4 are obtained by substituting the exact solutions in Table 1 to (6).Under the loading condition, the nonlinear iterative method in Section 2 is applied to obtain the solution.Figure 5 compares the exact and numerical solutions, and Figure 6 shows the convergence behavior of the solutions in Figure 5.The errors of the solutions at the th iteration are defined as follows: Three cases of Error () are plotted as a function of the iteration number in Figure 7.The solutions converge at 200-300 iterations.

3.3.
Convergence of the Procedure.The accuracy of the applied iterative method for the nonlinear spring model is proven in Section 3.2.In this section, 2 cases of loading conditions are investigated to show the convergence of the solutions.
The locally distributed rectangular-type loadings in Figure 8 are taken The number of iterations ()  Figure 9 presents the convergence behavior of the solutions according to the 3 cases of   :  =  =  = 1 and  = 0.2.As specified in Section 2.1,   does not affect the converged iterative solutions, but the convergence characteristics changed.The solution converged to a steady state faster when   is small.In Figure 10, two cases of converged solutions are compared: the numerical results using the one-way and two-way spring models.The loading condition is  normal (single) () in (18).They are different only when the beam deflected upward; that is, the deflections near  = ±5 are larger when the beam is separable (or free) from the foundation, whereas the oneway spring model is not.Of course, at  = 0, the solutions are equivalent because the physical system is the same in Figure 1 for  > 0. Therefore, Figure 10 shows the validity of the applied iterative method for the nonlinear spring model.
Another loading condition in Figure 11 is taken as follows: Other numerical experiments are also conducted according to the 3 cases of a linear spring coefficient , whereas the properties ( =  =  = 1 and  = 0.2) are fixed, and the solution converged (Figure 12).Finally, the validity and accuracy of the applied iterative procedure are investigated.The nonlinear spring force is considered, and the iterative procedure is applied successfully for the solution.Convergence of the deflections according to the external loading conditions is observed for the effects of the pseudo   and real .

Conclusion
In this work, we succeeded in applying the numerical method proposed by Jang et al. [21] to find the static deflection of an infinite beam on a full nonlinear elastic foundation.For that, one-way spring model is considered for the formulation of the nonlinear elastic one.Since the problem concerns one-way spring force, the governing equation for the static beam deflection is transformed as in (11).For the solution, an iterative procedure is applied for the calculation of the high nonlinearity.Some numerical experiments are carried out for showing the validity and the fast convergence of the applied numerical method.And we can also find that the results converge to the solutions fast for certain external loads.

Figure 1 :Figure 2 :
Figure 1: Nonlinear elastic foundation model: a nonlinear spring model (one way) and a conventional one (two way).

Figure 5 :
Figure 5: Numerical solutions compared to the exact solutions.

Figure 6 :
Figure 6: Convergence behaviors of the iterative solutions.

Figure 7 :
Figure 7: Errors between the exact solutions and iterative solutions.

Figure 13 :
Figure 13: Validity of the solution in Figure 12(a).