Equations Of Linear Elasticity Research Articles

An Optimal Control Problem for a Rotating Elastic Crane-Trolley-Load System

In this study we present an extension of a model of an elastic crane transporting a load by means of controlling the crane trolley motion and the crane rotation. In addition to the model considered in Kimmerle et al. (2017), we allow for rotations of the crane and include damping of the trolley and moments of inertia as well. We derive a fully coupled system of ordinary differential equations (ODE), representing the trolley and load (modelled as a pendulum), and partial differential equations (PDE), i.e. the linear elasticity equations for the deformed crane beam. The objective to be minimized is a linear combination of the terminal time, the control effort, the kinetic energy of the load, and penalty terms for the terminal conditions. We show the Fréchet-differentiability of the mechanical displacement field with respect to the location of the boundary condition that is moving. This is a crucial point for a further mathematical analysis on the existence of optimal controls and the derivation of necessary optimality conditions. Finally we present first results for the full time-optimal control of the extended model.

Homogenization of periodic micro‐structures with FFT‐based higher order schemes

AbstractA matrix‐free homogenization method making use of the fast Fourier transform (FFT) to work on structured grids was introduced by Moulinec and Suquet in 1994. Due to its simplicity and computational advantages, this so‐called Basic Scheme has provided the basis for many works in the past two decades. An extension presented by Tran, Monchiet and Bonnet allows for the quasi‐static equations of linear elasticity to include effects from higher‐order derivatives by expanding the strain and the displacement as asymptotic series. To compute the solution hierarchically, we derive generalized Lippmann Schwinger‐type equations. Furthermore, we show how existing algorithms have to be adapted for these problems. (© 2017 Wiley‐VCH Verlag GmbH & Co. KGaA, Weinheim)

Weight‐adjusted discontinuous Galerkin methods: Matrix‐valued weights and elastic wave propagation in heterogeneous media

SummaryWeight‐adjusted inner products are easily invertible approximations to weighted L2 inner products. These approximations can be paired with a discontinuous Galerkin (DG) discretization to produce a time‐domain method for wave propagation which is low storage, energy stable, and high‐order accurate for arbitrary heterogeneous media and curvilinear meshes. In this work, we extend weight‐adjusted DG methods to the case of matrix‐valued weights, with the linear elastic wave equation as an application. We present a DG formulation of the symmetric form of the linear elastic wave equation, with upwind‐like dissipation incorporated through simple penalty fluxes. A semidiscrete convergence analysis is given, and numerical results confirm the stability and high‐order accuracy of weight‐adjusted DG for several problems in elastic wave propagation.

Optimal control of an elastic crane-trolley-load system - a case study for optimal control of coupled ODE-PDE systems

ABSTRACTWe present a mathematical model of a crane-trolley-load model, where the crane beam is subject to the partial differential equation (PDE) of static linear elasticity and the motion of the load is described by the dynamics of a pendulum that is fixed to a trolley moving along the crane beam. The resulting problem serves as a case study for optimal control of fully coupled partial and ordinary differential equations (ODEs). This particular type of coupled systems arises from many applications involving mechanical multi-body systems. We motivate the coupled ODE-PDE model, show its analytical well-posedness locally in time and examine the corresponding optimal control problem numerically by means of a projected gradient method with Broyden-Fletcher-Goldfarb-Shanno (BFGS) update.

A Fourier Collocation Approach for Transit-Time Ultrasonic Flowmeter Under Multi-Phase Flow Conditions

A numerical model for a clamp-on transit-time ultrasonic flowmeter (TTUF) under multi-phase flow conditions is presented. The method solves equations of linear elasticity for isotropic heterogeneous materials with background flow where acoustic media are modeled by setting shear modulus to zero. Spatial derivatives are calculated by a Fourier collocation method allowing the use of the fast Fourier transform (FFT) and time derivatives are approximated by a finite difference (FD) scheme. This approach is sometimes referred to as a pseudospectral time-domain method. Perfectly matched layers (PML) are used to avoid wave-wrapping and staggered grids are implemented to improve stability and efficiency. The method is verified against exact analytical solutions and the effect of the time-staggering and associated lowest number of points per minimum wavelengths value is discussed. The method is then employed to model a complete TTUF measurement setup to simulate the effect of a flow profile on the flowmeter accuracy and a study of an impact of inclusions in flowing media on received signals is carried out.

Plasticity Effects in the Hole-Drilling Residual Stress Measurement in Peened Surfaces

The incremental hole-drilling technique (IHD) is a widely established and accepted technique to determine residual stresses in peened surfaces. However, high residual stresses can lead to local yielding, due to the stress concentration around the drilled hole, affecting the standard residual stress evaluation, which is based on linear elastic equations. This so-called plasticity effect can be quantified by means of a plasticity factor, which measures the residual stress magnitude with respect to the approximate onset of plasticity. The observed resultant overestimation of IHD residual stresses depends on various factors, such as the residual stress state, the stress gradients and the material’s strain hardening. In peened surfaces, equibiaxial stresses are often found. For this case, the combined effect of the local yielding and stress gradients is numerically and experimentally analyzed in detail in this work. In addition, a new plasticity factor is proposed for the evaluation of the onset of yielding around drilled holes in peened surface layers. This new factor is able to explain the agreement and disagreement found between the IHD residual stresses and those determined by X-ray diffraction in shot-peened steel surfaces.

Near conserving energy numerical schemes for two-dimensional coupled seismic wave equations

Two-dimensional coupled seismic waves, satisfying the equations of linear isotropic elasticity, on a rectangular domain with initial conditions and periodic boundary conditions, are considered. A quantity conserved by the solution of the continuous problem is used to check the numerical solution of the problem. Second order spatial derivatives, in the x direction, in the y direction and mixed derivative, are approximated by finite differences on a uniform grid. The ordinary second order in time system obtained is transformed into a first order in time system in the displacement and velocity vectors. For the time integration of this system, second order and fourth order exponential splitting methods, which are geometric integrators, are proposed. These explicit splitting methods are not unconditionally stable and the stability condition for time step and space step ratio is deduced. Numerical experiments displaying the good behavior in the long time integration and the efficiency of the numerical solution are provided.

Dispersion Properties of Lamb Waves in an Elastic Layer–Ideal Liquid Half-Space System

The problem of the propagation of normal waves in an elastic layer interacting with a compressible ideal liquid half-space is stated and solved using the three-dimensional linear equations of classical elasticity for the solid and the three-dimensional linearized Euler equations for the compressible ideal fluid. The dispersion curves are plotted for quasi-Lamb modes over a wide frequency range. The effect of the ideal compressible fluid and the elastic layer thickness on the dispersion of the phase velocities of quasi-Lamb modes in hydroelastic waveguides is analyzed. The localization properties of these modes in hydroelastic waveguides are studied. The numerical results are presented in the form of plots and analyzed.

A spectral-Tchebychev solution for three-dimensional dynamics of curved beams under mixed boundary conditions

This paper presents the application of the spectral-Tchebychev (ST) technique for solution of three-dimensional dynamics of curved beams/structures having variable and arbitrary cross-section under mixed boundary conditions. To accurately capture the vibrational behavior of curved structures, a three-dimensional (3D) solution approach is required since these structures generally exhibit coupled motions. In this study, the integral boundary value problem (IBVP) governing the dynamics of the curved structures is found using extended Hamilton's principle where the strain energy is expressed using 3D linear elasticity equation. To solve the IBVP numerically, the 3D spectral Tchebychev (3D-ST) approach is used. To evaluate the integral and derivative operations defined by the IBVP and to render the complex geometry into an equivalent straight beam with rectangular cross-section, a series of coordinate transformations are applied. To validate and assess the performance of the presented solution approach, two case studies are performed: (i) curved beam with rectangular cross-section, (ii) curved and pretwisted beam with airfoil cross-section. In both cases, the results (natural frequencies and mode shapes) are also found using a finite element (FE) solution approach. It is shown that the difference in predicted natural frequencies are less than 1%, and the mode shapes are in excellent agreement based on the modal assurance criteria (MAC) analyses; however, the presented spectral-Tchebychev solution approach significantly reduces the computational burden. Therefore, it can be concluded that the presented solution approach can capture the 3D vibrational behavior of curved beams as accurately as an FE solution, but for a fraction of the computational cost.

Nonlocal theory of curved rods. 2-D, high order, Timoshenko’s and Euler-Bernoulli models

Abstract New models for plane curved rods based on linear nonlocal theory of elasticity have been developed. The 2-D theory is developed from general 2-D equations of linear nonlocal elasticity using a special curvilinear system of coordinates related to the middle line of the rod along with special hypothesis based on assumptions that take into account the fact that the rod is thin. High order theory is based on the expansion of the equations of the theory of elasticity into Fourier series in terms of Legendre polynomials. First, stress and strain tensors, vectors of displacements and body forces have been expanded into Fourier series in terms of Legendre polynomials with respect to a thickness coordinate. Thereby, all equations of elasticity including nonlocal constitutive relations have been transformed to the corresponding equations for Fourier coefficients. Then, in the same way as in the theory of local elasticity, a system of differential equations in terms of displacements for Fourier coefficients has been obtained. First and second order approximations have been considered in detail. Timoshenko’s and Euler-Bernoulli theories are based on the classical hypothesis and the 2-D equations of linear nonlocal theory of elasticity which are considered in a special curvilinear system of coordinates related to the middle line of the rod. The obtained equations can be used to calculate stress-strain and to model thin walled structures in micro- and nanoscales when taking into account size dependent and nonlocal effects.

Strong-form approach to elasticity: Hybrid finite difference-meshless collocation method (FDMCM)

• FDMCM combines the finite difference and the strong form meshless method. • FDMCM scheme uses a uniform Cartesian grid embedded in complex geometry. • The method is suited for problems with geometries known imperfectly from images. • FDMCM method combines the advantages of both FD and meshless methods. We propose a numerical method that combines the finite difference (FD) and strong form (collocation) meshless method (MM) for solving linear elasticity equations. We call this new method FDMCM. The FDMCM scheme uses a uniform Cartesian grid embedded in complex geometries and applies both methods to calculate spatial derivatives. The spatial domain is represented by a set of nodes categorized as (i) boundary and near boundary nodes, and (ii) interior nodes. For boundary and near boundary nodes, where the finite difference stencil cannot be defined, the Discretization Corrected Particle Strength Exchange (DC PSE) scheme is used for derivative evaluation, while for interior nodes standard second order finite differences are used. FDMCM method combines the advantages of both FD and DC PSE methods. It supports a fast and simple generation of grids and provides convergence rates comparable to weak formulations. We demonstrate the appropriateness and robustness of the proposed scheme through various benchmark problems in 2D and 3D. Numerical results show good accuracy and h -convergence properties. The ease of computational grid generation makes the method particularly suited for problems where geometries are very complicated and known only imperfectly from images, frequently occurring in e.g. geomechanics and patient-specific biomechanics, where the proposed FDMCM method, after its extension to non-linear regime, appears to be a promising alternative to the traditional weak form-based numerical schemes used in the field.

HDG methods for elastodynamics

We derive and analyze a hybridizable discontinuous Galerkin (HDG) method for approximating weak solutions to the equations of time-harmonic linear elasticity on a bounded Lipschitz domain in three dimensions. The real symmetry of the stress tensor is strongly enforced and its coefficients as well as those of the displacement vector field are approximated simultaneously at optimal convergence with respect to the choice of approximating spaces, wavenumber, and mesh size. Sufficient conditions are given so that the system is indeed transferable onto a global hybrid variable that, for larger polynomial degrees, may be approximated via a smaller-dimensional space than the original variables. We construct several variants of this method and discuss their advantages and disadvantages, and give a systematic approach to the error analysis for these methods. We touch briefly on the application of this error analysis to the time-dependent problem, and finally, we examine two different implementations of the method over various polynomial degrees and numerically demonstrate the convergence properties proven herein.

On the existence, uniqueness and stability of solutions for semi-linear generalized elasticity equation with general damping term

In this paper, we consider a semi-linear generalized hyperbolic boundary value problem associated to the linear elastic equations with general damping term and nonlinearities of variable exponent type. Under suitable conditions, local and global existence theorems are proved. The uniqueness of the solution have been gotten by eliminating some hypotheses that have been imposed by other authors for different particular problems. We show that any solution with nontrivial initial datum becomes stable.

Shear Stiffness of Homogeneous, Orthotropic, Prismatic Beams

Exact analytical expressions for shear stiffness constants are provided herein for homogeneous, orthotropic, prismatic beams with arbitrary cross-sectional shapes. The analytical expressions are extracted, without invoking any approximations, from exact solutions of the linear equations of three-dimensional elasticity for flexure. Owing to the elasticity solutions being in terms of a stress function for cross sections with arbitrary shapes, the shear stiffness constants are also provided in terms of the same stress function. For elliptical and rectangular cross sections, this stress function is known so that shear stiffness constants are provided in closed form. These closed-form expressions constitute a standard with which the accuracy of two-dimensional beam cross-sectional analyses to calculate shear stiffness in the presence of orthotropic materials may be assessed. The calculated shear stiffness constants, from two such cross-sectional analyses, are successfully validated in this manner.

Coupled variational formulations of linear elasticity and the DPG methodology

This article presents a general approach akin to domain-decomposition methods to solve a single linear PDE, but where each subdomain of a partitioned domain is associated to a distinct variational formulation coming from a mutually well-posed family of broken variational formulations of the original PDE. It can be exploited to solve challenging problems in a variety of physical scenarios where stability or a particular mode of convergence is desired in a part of the domain. The linear elasticity equations are solved in this work, but the approach can be applied to other equations as well. The broken variational formulations, which are essentially extensions of more standard formulations, are characterized by the presence of mesh-dependent broken test spaces and interface trial variables at the boundaries of the elements of the mesh. This allows necessary information to be naturally transmitted between adjacent subdomains, resulting in coupled variational formulations which are then proved to be globally well-posed. They are solved numerically using the DPG methodology, which is especially crafted to produce stable discretizations of broken formulations. Finally, expected convergence rates are verified in two different and illustrative examples.

An arbitrary Lagrangian‐Eulerian framework with exact mass conservation for the numerical simulation of 2D rising bubble problem

SummaryAn arbitrary Lagrangian‐Eulerian framework, which combines the advantages of both Lagrangian and Eulerian methods, is presented to solve incompressible multiphase flow problems. The incompressible Navier‐Stokes equations are discretized using the side‐centered unstructured finite volume method, where the velocity vector components are defined at the midpoint of each cell face, while the pressure term is defined at element centroids. The pressure field is treated to be discontinuous across the interface with the discontinuous treatment of density and viscosity. The surface tension term at the interface is treated as a force tangent to the interface and computed with several different approaches including the use of Legendre polynomials. In addition, the several different discretizations of interface kinematic boundary conditions are investigated. For the application of the interface kinematic boundary condition, a special attention is given to satisfy both local and global discrete geometric conservation law to conserve the total mass of both species at machine precision. The mesh vertices are deformed by solving the linear elasticity equations due to the normal displacement of interface. The resulting algebraic equations are solved in a fully coupled manner, and a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned subdomain is used for the resulting fully coupled system. The method is validated by simulating the classical benchmark problem of a single rising bubble in a viscous fluid due to buoyancy. The results of numerical simulations are found out to be in an excellent agreement with the earlier results in the literature. The mass of the bubble is conserved, and discontinuous pressure field is obtained to avoid errors due to the incompressibility condition in the vicinity of the interface, where the density and viscosity jumps occur.

FIRST PRINCIPLES DERIVATION OF A STRESS FUNCTION FOR AXIALLY SYMMETRIC ELASTICITY PROBLEMS, AND APPLICATION TO BOUSSINESQ PROBLEM

In this work, a stress function is derived from first principles to describe the behaviour of three dimensional axially symmetric elasticity problems involving linear elastic, isotropic homogeneous materials. In the process, the fifteen governing partial differential equations of linear isotropic elasticity were reduced to the solution of the biharmonic problem involving the stress function. thus simplifying the solution process. The stress function derived was found to be identical with the Love stress function. The stress function was then applied to solve the axially symmetric problem of finding the stress fields, strain fields and displacement fields in the semi-infinite linear elastic, isotropic homogeneous medium subject to a point load P acting at the origin of coordinates also called the Boussinesq problem. The results obtained in this study for the stresses and displacements were exactly identical with those from literature, as obtained by Boussinesq.http://dx.doi.org/10.4314/njt.v36i3.15

FIRST PRINCIPLES DERIVATION OF A STRESS FUNCTION FOR AXIALLY SYMMETRIC ELASTICITY PROBLEMS, AND APPLICATION TO BOUSSINESQ PROBLEM

In this work, a stress function is derived from first principles to describe the behaviour of three dimensional axially symmetric elasticity problems involving linear elastic, isotropic homogeneous materials. In the process, the fifteen governing partial differential equations of linear isotropic elasticity were reduced to the solution of the biharmonic problem involving the stress function. thus simplifying the solution process. The stress function derived was found to be identical with the Love stress function. The stress function was then applied to solve the axially symmetric problem of finding the stress fields, strain fields and displacement fields in the semi-infinite linear elastic, isotropic homogeneous medium subject to a point load P acting at the origin of coordinates also called the Boussinesq problem. The results obtained in this study for the stresses and displacements were exactly identical with those from literature, as obtained by Boussinesq. http://dx.doi.org/10.4314/njt.v36i3.15

A High-Resolution Finite Volume Seismic Model to Generate Seafloor Deformation for Tsunami Modeling

A high-resolution finite volume method approach to incorporating time-dependent slip across rectangular subfaults when modeling general fault geometry is presented. The fault slip is induced by a modification of the Riemann problem to the linear elasticity equations across cell interfaces aligned with the subfaults. This is illustrated in the context of the high-resolution wave-propagation algorithms that are implemented in the open source Clawpack software (www.clawpack.org), but this approach could be easily incorporated into other Riemann solver based numerical methods. Surface deformation results are obtained in both two and three dimensions and compared to those given by the steady-state, homogeneous half-space Okada solution.

Space–time L2 estimates for elastic waves and applications

We establish a kind of weighted space–time L2 estimate, which belongs to Keel–Smith–Sogge type estimates, for perturbed linear elastic wave equations in Rn(n≥3). We also apply this estimate to give an alternative and short proof for John's almost global existence theorem for 3-D nonlinear elastic waves (John 1988, [10]), and show the global existence of small smooth solutions for nonlinear elastic waves in high space dimensions n≥4.