The differential elastic deflection line of a beam axis is widespread in engineering practice. Non-uniform settlement of a tank yields additional stresses in its structure, possibly leading to failure, that is deterioration or damage, triggering a significant threat to petroleum storage tanks. In order to capture these phenomena, the tank bottom plate is modelled by a simply supported beam resting either on stiff bedding [7] or on elastic bedding [18]. Other models are applied to assess internal forces in the storage tank structure [7], [17], [18]. The US Standard [3] makes use of working out design criteria of these models for large storage tanks.

Deflection of a rectangular plate, assumed a preliminary small cylindrical curvature [16], is expressed by a differential equation:

Equations of a similar sort are applied to model the intercomponent interaction of a vehicle/railway track/track bed/subsoil system. In railway engineering domain, the beam models on elastic foundation are applied in the range stated above, see [4], [5], [13], [14], and in railway track geometry issues [11]. The differential equation (1) is also applied in quantum mechanics to deal with the particle vibration in potential well analysis. The relevance of this equation is available.

A wide application range makes it reasonable to propose a new solution method of the equations, keeping in mind critical solution aspects, especially in numerical approach.

The modern design engineering practice requires a fluent use of CAD (Computer-aided design) methods in computational mechanics. Contemporary constitutive theories bring about non-linearity in the system of differential equations defining the initial boundary value problem. The equivalent non-linear equation systems require complex numerical procedures.

A widespread use of CAD methods linked with finite element method (FEM) or finite difference method (FDM) provides a vehicle to develop commercial software. The users, however, are frequently banned from numerical procedures; thus their interaction with the software is limited to operate the so-called black box, that is, they enter the input data and receive the final solution. The problem appears when the results of the computing process are not the ones expected by the designer.

It is a generally non-trivial task to relevantly interpret the computational, discretisation-based results. Whilst attention is paid to the finite element dimensions, marking the FE (finite element) mesh sensibility, the computational results may strongly differ from the observations, for example, in the steel plate cyclic load test, involving shear degradation [2].

Solution ambiguity may appear in the simplest possible cases, for example, beam deflection or other engineering issues incorporating similar mathematical models.

Most boundary problems make it extremely difficult to fulfil the conditions of a well-stated boundary problem, that is, to fulfil Hadamard's conditions. These conditions concern solution existence, uniqueness and stability, the second one is the focus of this article.

Non-classical operational calculus is basically aimed at creating uniform mathematical methods to solve specified forms of differential, difference or integral equations. It naturally generalises classified mathematical issues of a Heaviside idea or Mikusinski operators.

In classification of Mathematical Reviews^{®} Mathematical Subject Classification (MSC2010), they are described under No. 44 Integral Transforms, Operational Calculus, and then in detail under No. 44A40 – Calculus of Mikusinski and other operational calculi.

The non-classical operational calculus allows to uniformly describe and analyse a variety of engineering (mechanical) problems concerning continuum or discrete fields, including methods to solve differential or difference equations and their systems. Three linear operations, _{q} and _{q}, and two linear spaces, ^{1} and ^{0}, are defined, assuming ^{1} ⊂ ^{0}, to read as

The following conditions hold:
_{q} and _{q} are included in [8, 9].

Note that many operations exist to fulfil the conditions stated before; their examples are shown in [8]. In the article, it is sufficient to mention two representations only.

In a continuous field case, the following conditions hold:

In a discrete field, that is, the space _{k}}=(_{0}, _{1}, _{2}, ...), it can be considered as:

Whilst there exist an unambiguously defined solution to the equation,
^{1} and ^{0} commutative with operations _{q}, the solution may be called a generalised exponential function, which is defined as

Whilst the model features an operation

In the case of endomorphism
_{q} = |_{t}_{=0}, where ^{0} and ^{1} are the function spaces of relevant classes [9], the character of a generalised exponential function is

It is easy to check that the function
_{q}_{q}_{q}

Given a discrete field with the operations,
_{0}}, the exponential function is a series defined as

Figure 1 shows the diagram of series (15) for different values of number

The shape and characteristics of a generalised exponential function ^{Rtq}

Equation (13) leads to the following:
_{0}}={1}.

The diagrams in Figures 1 and 2 mark sequential solutions to some equations, possibly to represent stable and unstable problems on both modelling and solution phases, including the use of approximate methods [1].

The equation for the deflection of a beam of length

Equations (20) and (21) correspond to equations (18) and (19), respectively. Whilst variable stiffness is assumed, equation (18) is analysed by operational methods described in [11] by applying a chain-type connection of two generalised inertial parts of the first order.

Owing to a specific supporting variant, the equations are complemented by relevant boundary conditions. In the simply supported case, the following conditions hold:

It is a considerable task to analytically find a deflection function of the beam, thus approximative methods are open, however, excluding a ‘blind user’ routine.

A related case concerns neutral axis definition in the cross-sectional plane of a compressed bar, for example, a concrete post covered with steel or composite coating (reinforced by carbon or glass fibres). The neutral axis passes through the cross section of the beam (post); it is the locus of zero normal stresses and axial strains, that is, the undeflected region. Whilst no axial forces act in a non-curved bar made of homogeneous material, the neutral axis at each cross section intercepts the centroid. The presence of cross-sectional axial forces makes the neutral axis shift. The location of neutral axis is linked with structural state, possibly detected in every time instant based on the deformation measurement.

The application of numerical method to assess deflections of bars brings about a number of disadvantages. Whilst a differential equation is represented by its approximation, for example, a difference ratio, it is probable to get quantitatively and qualitatively incorrect results even after entering a simple input differential equation form.

Equation (18) or its equivalent form (20) may be transformed, assuming

Here the assumptions are still general. Whilst in operational calculus, the operation

The derivative
_{h} operation, that is,

Stating _{k}, the differential equation (24) turns into difference equation (recurrence equation); thus a continuous–discrete transformation is conducted here:

The above equation may be expressed by means of operation _{0} is the assumed natural number, fulfilling condition _{0}

Let _{0}_{0}, where _{0} is assumed. Then the equation (24) considered in continuous field with conditions (22) has one solution, and the solution is equal to zero.

The solution of the recurrence equation (27) with condition _{0}=0 according to (17) can be written as

Introducing the second condition _{k0}=0 leads to

If it happens,

If condition (31) holds, the difference equation approximates the differential equation with conditions (28) and shows an infinite solution range, keeping in mind that the differential equation shows exactly a single zero solution.

The presented numerical methodology brings an essential quantitative and qualitative problem reclassification, that is, the solution passes from a differential unique form to an ambiguous sort. The variation between the exact and the approximate solutions (the latter presents values at nodes) is larger whilst taking a large factor

The work [8] shows the operational calculus model, detecting the form ^{2}

Both Kerr's model analysis and Winkler–Pasternak foundation concept regarded in [15] showed another practical application of (27); it used a shear layer (soil) with geosynthetic membranes. The presented model highly resembles the real geosynthetic-soil systems, possible road or railway flexible pavements [10].

Kerr introduced a model, even a large spectrum of models made of Winkler, Pasternak and Héteny layers. The equation of Kerr's foundation model is similar to the form presented in [10] and Figure 3 whilst not regarding a geosynthetic layer, that is, _{p}=0 and _{t}=_{b}=_{t}, _{s} are upper and lower elastic foundation constants, and

In a one-dimensional case, the model equation (32) modifies the Laplace operator ∇^{2}, which becomes a common second derivative

Sequences (15) are applied to determine responses of discrete two-dimensional systems [9]. Whilst discretising a continuous system of distributed parameters, the sequences can be applied to approximate the system response. The systems cover the performance of civil engineering plate and shell structures, wave propagation in soils and so on [12].

The presented discrete time signals are relevant to approximately determine the response of some class of differential problems adjusted to mechanics, for example, axially compressed columns. The signals may also be used to approximate harmonic and inharmonic signals and modulated signals.

The application of properly generalised exponential and trigonometric function makes it possible to work numerically stable or unstable procedures. A continuous–discrete engineering problem transformation causes a hazard of its quantitative and qualitative re-classification. In the course of continuous–discrete re-modelling of an engineering system, the initial unique solution turns into ambiguous solution or no solution at all. The numerical procedure user should be aware of such situations and take them into account in the design process.

Even a highly selected constitutive model is not enough for the considered boundary problem to relevantly reflect the real analysed case. The static or dynamic load feature possibly complements the differential equation by inertial effects, effects of elastic and shock wave propagation, loading rate and so on. Numerical modelling is also highly influenced by the specific material instability or imperfections, that is, density and temperature variations, microcracks showcasing ductile fracture patterns out of yield zones.

All the enlisted issues consequently increase the algorithm sophistication and computational effort [6]. On the other hand, they trigger more complex mathematical forms, making the errors in numerical procedures more difficult to detect.

The design quality progress is chiefly an effect of progress in numerical domain. However, the engineers should stay critical to numerical results of the ‘black box’ software, excluding the algorithm tracing.

Reviews on Finite Element Modeling Practices of Stone Columns for Soft Soil Stabilization Beneath an Embankment Dam Proposal of concept for structural modelling of hybrid beams Resonance of a structure with soil elastic waves released in non-linear hysteretic soil upon unloading Settlement Analysis of a Sandy Clay Soil Reinforced with Stone Columns The evolution of the shape of composite dowels Parametric study of the earth dam's behaviour subjected to earthquake Impact of longwall mining on slope stability – A case study