# Numerical Solution of Uncertain Beam Equations Using Double Parametric Form of Fuzzy Numbers

Authors: Smita Tapaswini; S. Chakraverty
Journal: Applied Computational Intelligence and Soft Computing
 (2013)

---

## Abstract

Present paper proposes a new technique to solve uncertain beam equation using double parametric form of fuzzy numbers. Uncertainties appearing in the initial conditions are taken in terms of triangular fuzzy number. Using the single parametric form, the fuzzy beam equation is converted first to an interval-based fuzzy differential equation. Next, this differential equation is transformed to crisp form by applying double parametric form of fuzzy number. Finally, the same is solved by homotopy perturbation method (HPM) to obtain the uncertain responses subject to unit step and impulse loads. Obtained results are depicted in terms of plots to show the efficiency and powerfulness of the methodology.


---

## Body

## 1. Introduction

Recently, theory of differential equations plays a vital role to model physical and engineering problems such as in solid and fluid mechanics, viscoelasticity, biology, physics, and other areas of science. But in actual case, the parameters, variables, and initial conditions involved in the differential equations may be uncertain, or a vague estimation of those are found in general by some observation, experiment, experience, or maintenance  induced error. So, to overcome the uncertainty and vagueness, one may use fuzzy environment in parameters, variables, and initial condition in place of crisp (fixed) ones. So, with these uncertainties the general differential equations turn into fuzzy differential equations (FDEs). In real-life application, it is too complicated to obtain the exact solution to fuzzy differential equations, so one may need a reliable and efficient numerical technique for the solution for fuzzy differential equations.There exist a large number of papers dealing with fuzzy differential equations and its applications in the open literature. Some are reviewed and cited here for better understanding of the present analysis. Chang and Zadeh [1] first introduced the concept of a fuzzy derivative, followed by Dubois and Prade [2] who defined and used the extension principle in their approach. The fuzzy differential equations and fuzzy initial value problems are studied by Kaleva [3, 4] and Seikkala [5]. Various numerical methods for solving fuzzy differential equations are introduced in [6–16]. Ma et al. [6] developed a scheme based on the classical Euler method to solve fuzzy ordinary differential equations. A two-dimensional differential transform method to solve fuzzy partial differential equations (FPDEs) has been studied in [7]. Abbasbandy and Allahviranloo [8] applied Taylor method for the solutions of fuzzy differential equations. Abbasbandy et al. [9] developed a numerical method for solving fuzzy differential inclusions, and in the proposed method, fuzzy reachable set is used to approximate the solution. Variation of constant formula is handled by Khastan et al. [10] to solve first-order fuzzy differential equations. Variational iteration method is discussed by Allahviranloo et al. [11] to obtain the exact solutions of fuzzy wave-like equations with variable coefficients. The concept of generalized H-differentiability is studied by Chalco-Cano and Román-Flores [12] to solve fuzzy differential equations. Mondal and Roy [13] described the solution procedure of a first-order linear non homogeneous ordinary differential equation in fuzzy environment. Palligkinis et al. [14] applied the Runge-Kutta method for more general problems and proved the convergence for s-stage Runge-Kutta methods. Akin et al. [15] developed an algorithm based on α-cut of a fuzzy set for the solution of second-order fuzzy initial value problems. Very recently, Khastan et al. [10] proposed improved Euler method for the solution of fuzzy differential equations.Similarly, many authors studied various other methods to solventh order fuzzy differential equations in [17–21]. Based on the idea of collocation method, Allahviranloo et al. [17] investigated the numerical solution of nth order fuzzy differential equations. Jafari et al. [18] used variational iteration method for solving nth-order fuzzy differential equations recently. Parandin [19] discussed Runge-Kutta method for the numerical solution of fuzzy differential equations of nth-order. Hashemi et al. [20] studied homotopy analysis method for the solution of system of fuzzy differential equations. Zhang and Wang [21] utilized time domain methods for the solutions of nth-order fuzzy differential equations.Recently, homotopy perturbation method is found to be a powerful tool for linear and nonlinear differential equations. The HPM was first developed by He in [22, 23], and many authors applied this method to solve various linear and nonlinear functional equations of scientific and engineering problems [24, 25]. The solution is considered as the sum of infinite series, which converges rapidly to accurate solutions. In the homotopy technique (in topology), a homotopy is constructed with an embedding parameter p∈[0,1], which is considered as a “small parameter”. Very recently homotopy perturbation method has been implemented to a wide class of problems. Few researchers have also investigated the solution of fuzzy differential equations using HPM. Allahviranloo and Hashemzehi [26] successfully applied homotopy perturbation method for the solution of fuzzy-fredholm integral equations. Babolian et al. [27] solved integro-differential equations using homotopy perturbation method. Matinfar and Saeidy [28] handled an application of homotopy perturbation method for fuzzy integral equations. Numerical solution of fuzzy initial value problems under generalized differentiability by HPM is studied by Ghanbari [29]. Recently, Tapaswini and Chakraverty [30] implemented HPM for the solution of nth-order fuzzy linear differential equations.Moreover, Bede [31] described the exact solutions of fuzzy differential equations in his note in an excellent way. Ahmada et al. [32] study analytical and numerical solutions of fuzzy differential equations based on the extension principle. Buckley and Feuring [33] applied two analytical methods for solving nth-order linear differential equations with fuzzy initial conditions. In the first method, they simply fuzzify the crisp solution to obtain a fuzzy function and then check whether it satisfies the differential equation or not. And the second method was just the reverse of the first method.The above literature review reveals that the fuzzy differential equations are always converted to two crisp differential equations. Then corresponding crisp systems are solved to obtain the bounds of fuzzy solution. But here, a new method is proposed based on double parametric form of fuzzy as define in Definition2 where fuzzy differential equation has been converted to a single crisp differential equation. Finally the corresponding crisp (single) differential equation is solved by homotopy perturbation method to obtain the fuzzy solution in double parametric form.This paper is organized as follows. In Section2, some basic preliminaries related to the present investigation are given. Section 3 represents the basic idea of HPM.  Proposed technique is applied with HPM in Section 4 to solve uncertain beam equation. Next, uncertain response analysis is presented followed by numerical results and discussions. Finally, in Section 7 conclusions are drawn.


## 2. Preliminaries

In this section, we present some notations, definitions, and preliminaries which are used further in this paper.Definition 1 (single parametric form of fuzzy numbers).
The triangular fuzzy numberU~=(a,b,c) can be represented with an ordered pair of functions through α-cut approach; namely, [u_(α),u-(α)]=[(b-a)α+a,-(c-b)α+c] where, α∈[0,1].
Theα-cut form is known as parametric form or single parametric form of fuzzy numbers.
It may be noted that the lower and upper bounds of the fuzzy numbers satisfy the following requirements.(i)
u_(α) is a bounded left continuous nondecreasing function over  [0,1].(ii)
u-(α) is a  bounded right continuous nonincreasing function over [0,1].(iii)
u_(α)≤u-(α),0≤α≤1.Definition 2 (double parametric form of fuzzy number).
Using the parametric form as discussed in Definition1 one has  U~=[u_(α),u-(α)]. Now one may write this as crisp number with double parametric form as U~(α,β)=β(u-(α)-u_(α))+u_(α), where α  and β∈[0,1].


## 3. Homotopy Perturbation Method [22, 23]

To illustrate the basic idea of this method, we consider the following nonlinear differential equation of the form:(1)A(u)-f(r)=0,r∈Ω,
with the boundary condition
(2)B(u,∂u∂n)=0,r∈Γ,
where A is a general differential operator, B a boundary operator, f(r) a known analytical function, and Γ is the boundary of the domain Ω. A can be divided  into two parts which are L and N, where L is linear and N is nonlinear. Equation (1) can therefore be written as follows:
(3)L(u)+N(u)-f(r)=0,r∈Ω.
By the homotopy technique, we construct a homotopy U(r,p):Ω×[0,1]→R, which satisfies
(4)H(U,p)=(1-p)[L(U)-L(v0)]+p[A(U)-f(r)]=0,p∈[0,1],r∈Ω
or
(5)H(U,p)=L(U)-L(v0)+pL(v0)+p[N(U)-f(r)]=0,
where r∈Ω, p∈[0,1] is an embedding parameter, and v0 is an initial approximation of (1). Hence, it is obvious that
(6)H(U,0)=L(U)-L(v0)=0,H(U,1)=A(U)-f(r)=0,
and the changing process of p from 0 to 1 is just that of H(U,p) from L(U)-L(v0) to A(U)-f(r). In topology, this is called deformation, L(U)-L(v0) and A(U)-f(r) is called homotopic. Applying the perturbation technique [22, 23], due to the fact that 0≤p≤1 can be considered as a small parameter, we can assume that the solution of (4) or (5) can be expressed as a series in p as follows:
(7)U=u0+pu1+p2u2+p3u3+⋯,
when p→1, (4) or (5) corresponds to (3); Then (7) becomes the approximate solution of (3), that is,
(8)u=limp→1U=u0+u1+u2+u3+⋯.
The convergence of the series (8) has been proved in [22, 23].


## 4. Proposed Method

Here, we first convert the fuzzy differential equation to interval-based fuzzy differential equation using single parametric form. Then by using double parametric form, interval-based fuzzy differential equation is reduced to crisp differential equation. Next, we apply homotopy perturbation method to solve the corresponding differential equation. As per the titled problem, let us now consider the fuzzy beam equation(9)ρA∂2v~∂t2+c∂v~∂t+EI∂4v~∂x4=F(x,t).
This may be written as
(10)∂2v~∂t2+cρA∂v~∂t+EIρA∂4v~∂x4=F(x,t)ρA,
where ρ, A, c, E, and I represent the mass, density, cross sectional area, damping coefficients per unit length, Young’s modulus of elasticity, and moment of inertia of the beam. F(x,t) is the externally applied force and v~(x,t;α) is the transverse fuzzy displacement. Fuzzy initial conditions are considered as v~(0)=v~′(0)=(-0.1,0,0.1). It may be noted that if the initial condition is crisp, then we may have the initial condition as v(0)=v′(0)=0. Here, the initial condition has been taken as fuzzy with an idea that the condition may actually be uncertain; namely, it may be due to error in observation or experiment and so forth, where we may take the error or uncertainty as left and right spread as −0.1 and 0.1, respectively in terms of fuzzy triangular membership function. As such this will force the governing differential equation as a whole as uncertain. So naturally, the outcome or the output (result) must be uncertain. This way we may have the actual essence of the uncertainty in response which may benefit the engineers to understand the safety of the system in a better way. So, we need to have efficient methods to handle these problems.As per the single parametric form we may write the above fuzzy differential equation (10) as
(11)[∂2v_(t;α)∂t2,∂2v-(t;α)∂t2]+cρA[∂v_(t;α)∂t,∂v-(t;α)∂t]+EIρA[∂4v_(t;α)∂x4,∂4v-(t;α)∂x4]=F(x,t)ρA
subject to fuzzy initial condition
(12)[v_(0;α),v-(0;α)]=[v_′(0;α),v-′(0;α)]=[0.1α-0.1,0.1-0.1α],
where α∈[0,1].Next, using the double parametric form (as discussed in Definition2) above (11), it can be expressed as
(13){β(∂2v-(t;α)∂t2-∂2v_(t;α)∂t2)+∂2v_(t;α)∂t2}+cρA{β(∂v-(t;α)∂t-∂v_(t;α)∂t)+∂v_(t;α)∂t}+EIρA{β(∂4v-(t;α)∂x4-∂4v_(t;α)∂x4)+∂4v_(t;α)∂x4}=F(x,t)ρA
subject to the fuzzy initial conditions
(14){β(v-(0;α)-v_(0;α))+v_(0;α)}={β(v-′(0;α)-v_′(0;α))+v_′(0;α)}={β(0.2-0.2α)+(0.1α-0.1)},
where β∈[0,1].Let us now denote(15){β(∂2v-(t;α)∂t2-∂2v_(t;α)∂t2)+∂2v_(t;α)∂t2}=∂2v(t;α,β)∂t2,{β(∂v-(t;α)∂t-∂v_(t;α)∂t)+∂v_(t;α)∂t}=∂v(t;α,β)∂t,{β(∂4v-(t;α)∂x4-∂4v_(t;α)∂x4)+∂4v_(t;α)∂x4}=∂4v(t;α,β)∂x4,{β(v-(0;α)-v_(0;α))+v_(0;α)}=v(0;α,β),{β(v-′(0;α)-v_′(0;α))+v_′(0;α)}=v′(0;α,β).
Substituting these values in (13), we get
(16)∂2v(t;α,β)∂t2+cρA∂v(t;α,β)∂t+EIρA∂4v(t;α,β)∂x4=F(x,t)ρA,
with initial conditions
(17)v(0;α,β)=v′(0;α,β)={β(0.2-0.2α)+(0.1α-0.1)}.Hence, solving the corresponding crisp differential equation, one may get the solution asv(t;α,β). To obtain the lower and upper bound of the solution in single parametric form, we may put β=0 and 1, respectively. This may be represented as
(18)v(t;α,0)=v_(t,α),v(t,α,1)=v-(t,α).Now, to solve (16), we have applied homotopy perturbation method. According to HPM, we may construct a simple homotopy for an embedding parameter p∈[0,1] as follows:
(19)(1-p)∂2v(t;α,β)∂t2+p[∂2v(t;α,β)∂t2+cρA∂v(t;α,β)∂t+EIρA∂4v(t;α,β)∂x4-F(x,t)ρA]=0
or
(20)∂2v(t;α,β)∂t2+p[cρA∂v(t;α,β)∂t+EIρA∂4v(t;α,β)∂x4-F(x,t)ρA]=0.
Here, p is considered as a small homotopy parameter 0≤p≤1. For p=0,  (19) and (20) become a linear equation; that is, ∂2v(t;α,β)/∂t2=0, which is easy to solve. For p=1, (19) and (20) turn out to be same as the original equation (16). This is called deformation in topology. ∂2v(t;α,β)/∂t2 and (c/ρA)(∂v(t;α,β)/∂t)+(EI/ρA)(∂4v(t;α,β)/∂x4)-(F(x,t)/ρA) are called homotopic. Next, we can assume the solution of  (19) or (20) as a power series expansion in p as
(21)v(t;α,β)=v0(t;α,β)+pv1(t;α,β)+p2v2(t;α,β)+p3v3(t;α,β)+⋯,
where vi(t;α,β) for i=0,1,2,3,… are functions yet to be determined. Substituting (21) into  (19) or (20) and equating the terms with the identical powers of p, we have
(22)p0:∂2v0(t;α,β)∂t2=0,p1:∂2v1(t;α,β)∂t2+cρA∂v0(t;α,β)∂t+EIρA∂4v0(t;α,β)∂x4-F(x,t)ρA=0,p2:∂2v2(t;α,β)∂t2+cρA∂v1(t;α,β)∂t+EIρA∂4v1(t;α,β)∂x4=0,p3:∂2v3(t;α,β)∂t2+cρA∂v2(t;α,β)∂t+EIρA∂4v2(t;α,β)∂x4=0,
and so on.Choosing initial approximationv(0;α,β)=v′(0;α,β) and applying the inverse operator Ltt-1 (which is the inverse of the operator Ltt=∂2/∂t2) on both sides of every equation in (22)  one may obtain the following equations:
(23)v0(t;α,β)=v′(0;α,β)t+v(0;α,β),v1(t;α,β)=Ltt-1(-cρA∂v0(t;α,β)∂t-EIρA∂4v0(t;α,β)∂x4+F(x,t)ρA),v2(t;α,β)=Ltt-1(-cρA∂v1(t;α,β)∂t-EIρA∂4v1(t;α,β)∂x4),v3(t;α,β)=Ltt-1(-cρA∂v2(t;α,β)∂t-EIρA∂4v2(t;α,β)∂x4),
and so on.Now, substituting these terms in (21) with p→1, one may get the approximate solution of (16) as
(24)v(t;α,β)=v0(t;α,β)+v1(t;α,β)+v2(t;α,β)+v3(t;α,β)+⋯.
The solution series converge very rapidly. Proof of convergence of the above series may be found in [22, 23]. The rapid convergence means that only few terms are required to get the approximate solutions.


## 5. Uncertain Response Analysis

Let us consider the external applied forceF(x,t) as
(25)F(x,t)=f(x)g(t),
where f(x) is a specified space dependent deterministic function and g(t) is time dependent process. In the following paragraph, we will examine the uncertain response of the dynamic system (5) subject to two different loading conditions.


### 5.1. Unit Step Function Response

Now, the response of the beam subject to a unit step load has been considered of the formg(t)=Bu(t), where u(t) is the Heaviside function and B is a constant. By using HPM, we have
(26)v0(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}(1+t),v1(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(-cρAt22)+fBt22ρA,v2(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(c2ρ2A2t3Γ(4))-fBcρ2A2t3Γ(4)-EIBf4ρ2A2t4Γ(5),v3(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(-c3ρ3A3t4Γ(5))+fBc2ρ3A3t4Γ(5)+2EIBcf4ρ3A3t5Γ(6)+E2I2Bf8ρ3A3t6Γ(7),
and so on, where f(i)=∂if/∂xi.In a similar manner, the rest of the components can be obtained. Therefore, the solution can be written in a general form as(27)v(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2(r+1)×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+3).
To obtain the solution bound in single parametric form we may put β=0 and 1, respectively, to get the lower and upper bound of the solution as
(28)v_(t;α,0)=v_(t;α)=(0.1α-0.1){(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2(r+1)×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+3),v-(t;α,1)=v-(t;α)=(0.1-0.1α)×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2(r+1)×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+3).


### 5.2. Impulse Response

Next, let us consider the response of the beam subject to a unit impulse load of the formg(t)=δ(t), where δ(t) is the unit impulse function. Using HPM in this case again, one may have
(29)v0(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}(1+t),v1(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(-cρAt22)+ftρA,v2(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(c2ρ2A2t3Γ(4))-fcρ2A2t2Γ(3)-EIf4ρ2A2t3Γ(4),v3(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(-c3ρ3A3t4Γ(5))+fc2ρ3A3t3Γ(4)+2EIcf4ρ3A3t4Γ(5)+E2I2f8ρ3A3t5Γ(6),
and so on, where f(i)=∂if/∂xi.In a similar manner, the rest of the components can be obtained. Therefore, the solution can be written in a general form as(30)v(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2r+1×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+2).
Again, to obtain the solution bound for this case in single parametric form, one may put β=0 and 1, respectively, to get the lower and upper bound of the solution as
(31)v_(t;α,0)=v_(t;α)=(0.1α-0.1)×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2r+1×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+2),v-(t;α,1)=v-(t;α)=(0.1-0.1α)×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2r+1×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+2).


## 5.1. Unit Step Function Response

Now, the response of the beam subject to a unit step load has been considered of the formg(t)=Bu(t), where u(t) is the Heaviside function and B is a constant. By using HPM, we have
(26)v0(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}(1+t),v1(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(-cρAt22)+fBt22ρA,v2(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(c2ρ2A2t3Γ(4))-fBcρ2A2t3Γ(4)-EIBf4ρ2A2t4Γ(5),v3(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(-c3ρ3A3t4Γ(5))+fBc2ρ3A3t4Γ(5)+2EIBcf4ρ3A3t5Γ(6)+E2I2Bf8ρ3A3t6Γ(7),
and so on, where f(i)=∂if/∂xi.In a similar manner, the rest of the components can be obtained. Therefore, the solution can be written in a general form as(27)v(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2(r+1)×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+3).
To obtain the solution bound in single parametric form we may put β=0 and 1, respectively, to get the lower and upper bound of the solution as
(28)v_(t;α,0)=v_(t;α)=(0.1α-0.1){(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2(r+1)×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+3),v-(t;α,1)=v-(t;α)=(0.1-0.1α)×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2(r+1)×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+3).


## 5.2. Impulse Response

Next, let us consider the response of the beam subject to a unit impulse load of the formg(t)=δ(t), where δ(t) is the unit impulse function. Using HPM in this case again, one may have
(29)v0(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}(1+t),v1(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(-cρAt22)+ftρA,v2(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(c2ρ2A2t3Γ(4))-fcρ2A2t2Γ(3)-EIf4ρ2A2t3Γ(4),v3(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×(-c3ρ3A3t4Γ(5))+fc2ρ3A3t3Γ(4)+2EIcf4ρ3A3t4Γ(5)+E2I2f8ρ3A3t5Γ(6),
and so on, where f(i)=∂if/∂xi.In a similar manner, the rest of the components can be obtained. Therefore, the solution can be written in a general form as(30)v(t;α,β)={β(0.2-0.2α)+(0.1α-0.1)}×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2r+1×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+2).
Again, to obtain the solution bound for this case in single parametric form, one may put β=0 and 1, respectively, to get the lower and upper bound of the solution as
(31)v_(t;α,0)=v_(t;α)=(0.1α-0.1)×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2r+1×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+2),v-(t;α,1)=v-(t;α)=(0.1-0.1α)×{(1+ρAc)-(ρAce(-c/ρA)t)}+1ρA∑r=0∞(-1)rr!(EIρA)rf(4r)t2r+1×∑j=0∞(-cρA)j(j+r)!tjj!Γ(j+2r+2).


## 6. Numerical Results and Discussions

In this section, we present the numerical solution of uncertain beam equation using HPM. It is a gigantic task to include all the results with respect to various parameters involved in the corresponding equation here. As discussed above, uncertain responses subject to unit step and impulse responses have been considered. Computed results are depicted in terms of plots.  For the numerical result, we have considered the simply supported beam and(32)f(x)=sin(πxL).
Here, the numerical computations have been done by truncating the infinite series in (28) and (31) to a finite number of terms. For numerical simulations, let us denote c/m and EI/ρA,  respectively, as 2ηω3/2 and ω2, where ω is the natural frequency and η is the damping ratio. The values of the parameters are taken as B=1, ρA=1, L=π,  x=1/2,  and m=1.Next, by takingx=1/2 with natural frequency ω=5 rad/s and 10 rad/s and damping ratio η=0.05 as constant and varying both the t and α, the obtained fuzzy results for unit and impulse responses are depicted in Figures 1 and 2, respectively.Fuzzy responses atx=1/2 with natural frequency (a) ω=5 rad/s and (b) ω=10 rad/s and damping ratio η=0.05 for unit step response.
(a)
(b)Fuzzy responses atx=1/2 with natural frequency (a) ω=5 rad/s and (b) ω=10 rad/s and damping ratio η=0.05 for unit impulse response.
(a)
(b)As discussed above, for both cases withα=1, the fuzzy initial condition converts into crisp initial condition. It is interesting to note that for both the responses, lower and upper bounds of the fuzzy solutions are same for α=1.


## 7. Conclusions

In this paper, homotopy perturbation method has successfully been applied to the solution of an uncertain simply supported beam using double parametric form of fuzzy numbers. The proposed double parametric form approach is found to be easy and straightforward. Here, performance of the method is shown by using triangular fuzzy number. It is interesting to note forα=1, in both the responses lower response is equal to the upper response. Though the solution by HPM is of the form of an infinite series, it can be written in a closed form in some cases. The main advantage of HPM is the capability to achieve exact solution and rapid convergences with few terms.



---

*Source: 764871-2013-10-20.xml*