Benchmark analytical solutions from beams with shared eigenpair

Structures with governing equations having identical inertial terms but somewhat differing stiffness terms can be termed flexurally analogous. An example of such a structure includes an axially loaded non-uniform beam and an unloaded uniform beam, for which an exact solution exists. We find that there exist shared eigenpairs (frequency and mode shapes) for a particular mode between such structures. Non-uniform beams with uniform axial loads, gravity loaded beams and rotating beams are considered and shared eigenpairs with uniform beams are found. In general, the derived flexural stiffness functions (FSF's) for the non-uniform beams required for the existence of shared eigenpair have internal singularities, but some of the singularities can be removed by an appropriate selection of integration constants using the theory of limits. The derived functions yield an insight into the relationship between the axial load and flexural stiffness of axially loaded beam structures. The derived functions can serve as benchmark solutions for numerical methods. (C) 2016 Elsevier Ltd. All rights reserved.

Advancing Edge Speeds of Epithelial Monolayers Depend on Their Initial Confining Geometry

Collective cell migrations are essential in several physiological processes and are driven by both chemical and mechanical cues. The roles of substrate stiffness and confinement on collective migrations have been investigated in recent years, however few studies have addressed how geometric shapes influence collective cell migrations. Here, we address the hypothesis that the relative position of a cell within the confinement influences its motility. Monolayers of two types of epithelial cells-MCF7, a breast epithelial cancer cell line, and MDCK, a control epithelial cell line-were confined within circular, square, and cross-shaped stencils and their migration velocities were quantified upon release of the constraint using particle image velocimetry. The choice of stencil geometry allowed us to investigate individual cell motility within convex, straight and concave boundaries. Cells located in sharp, convex boundaries migrated at slower rates than those in concave or straight edges in both cell types. The overall cluster migration occurred in three phases: an initial linear increase with time, followed by a plateau region and a subsequent decrease in cluster speeds. An acto-myosin contractile ring, present in the MDCK but absent in MCF7 monolayer, was a prominent feature in the emergence of leader cells from the MDCK clusters which occurred every similar to 125 mu m from the vertex of the cross. Further, coordinated cell movements displayed vorticity patterns in MDCK which were absent in MCF7 clusters. We also used cytoskeletal inhibitors to show the importance of acto-myosin bounding cables in collective migrations through translation of local movements to create long range coordinated movements and the creation of leader cells within ensembles. To our knowledge, this is the first demonstration of how bounding shapes influence long-term migratory behaviours of epithelial cell monolayers. These results are important for tissue engineering and may also enhance our understanding of cell movements during developmental patterning and cancer metastasis.

Boundary knot method for 2D and 3D Helmholtz and convection-diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration-free, boundary-type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non-singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection-diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed.

Numerical Analysis of Theoretical Model of the Rf MEMs Switches

An improved electromechanical model of the RF MEMS (radio frequency microelectromechanical systems) switches is introduced, in which the effects of intrinsic residual stress from fabrication processes, axial stress due to stretching of beam, and fringing field are taken into account. Four dimensionless numbers are derived from the governing equation of the developed model. A semi-analytical method is developed to calculate the behavior of the RF MEMS switches. Subsequently the influence of the material and geometry parameters on the behavior of the structure is analyzed and compared, and the corresponding analysis with the dimensionless numbers is conducted too. The quantitative relationship between the presented parameters and the critical pull-in voltage is obtained, and the relative importance of those parameters is given.

An Analytical Solution for Three-Dimensional Diffraction of Plane P-Waves by a Hemispherical Alluvial Valley with Saturated Soil Deposits

An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion technique. Unlike previous studies, in which the saturated soil deposits were simulated with the single-phase elastic theory, in this paper, they are simulated with Biot's dynamic theory for saturated porous media, and the half space is assumed as a single-phase elastic medium. The effects of the dimensionless frequency, the incidence angle of P-wave and the porosity of soil deposits on the surface displacement magnifications of the hemispherical alluvial valley are investigated. Numerical results show that the existence of a saturated hemispherical alluvial valley has much influence on the surface displacement magnifications. It is more reasonable to simulate soil deposits with Biot's dynamic theory when evaluating the displacement responses of a hemispherical alluvial valley with an incidence of P-waves.

A Numerical Study of Indentation Using Indenters of Different Geometry

Finite element simulation of the Berkovich, Vickers, Knoop, and cone indenters was carried out for the indentation of elastic-plastic material. To fix the semiapex angle of the cone, several rules of equivalence were used and examined. Despite the asymmetry and differences in the stress and strain fields, it was established that for the Berkovich and Vickers indenters, the load-displacement relation can closely be simulated by a single cone indenter having a semiapex angle equal to 70.3degrees in accordance with the rule of the volume equivalence. On the other hand, none of the rules is applicable to the Knoop indenter owing to its great asymmetry. The finite element method developed here is also applicable to layered or gradient materials with slight modifications.

Analytical Interaction of the Acoustic Wave and Turbulent Flame

A modified resonance model of a weakly turbulent flame in a high-frequency acoustic wave is derived analytically. Under the mechanism of Darrieus-Landau instability, the amplitude of flame wrinkles, which is as functions of the expansion coefficient and the perturbation wave number, increases greatly independent of the 'stationary' turbulence. The high perturbation wave number makes the resonance easier to be triggered but weakened with respect to the extra acoustic wave. In a closed burning chamber with the acoustic wave induced by the flame itself, the high perturbation wave number is to restrain the resonance for a realistic flame.

Numerical and Analytical Study on the Pull-in Instability of Micro-Structure under Electrostatic Loading

The one-mode analysis method on the pull-in instability of micro-structure under electrostatic loading is presented. Taylor series are used to expand the electrostatic loading term in the one-mode analysis method, which makes analytical solution available. The one-mode analysis is the combination of Galerkin method and Cardan solution of cubic equation. The one-mode analysis offers a direct computation method on the pull-in voltage and displacement. In low axial loading range, it shows little difference with the established multi-mode analysis on predicting the pull-in voltages for three different structures (cantilever, clamped-clamped beams and the plate with four edges simply-supported) studied here. For numerical multi-mode analysis, we also show that using the structural symmetry to select the symmetric mode can greatly reduce both the computation effort and the numerical fluctuation.

Analytical Study of Idealized Two-Dimensional Cellular Detonations

In this study, the idealized two-dimensional detonation cells were decomposed into the primary units referred to as sub-cells. Based on the theory of oblique shock waves, an analytical formula was derived to describe the relation between the Mach number ratio through triple-shock collision and the geometric properties of the cell. By applying a modified blast wave theory, an analytical model was developed to predict the propagation of detonation waves along the cell. The calculated results show that detonation wave is, first, strengthened at the beginning of the cell after triple-shock collision, and then decays till reaching the cell end. The analytical results were compared with experimental data and previous numerical results; the agreement between them appears to be good, in general.

Modeling Nonlinear Peeling of Ductile Thin Films-Critical Assessment of Analytical Bending Models Using FE Simulations

Three analytical double-parameter criteria based on a bending model and a two-dimensional finite element analysis model are presented for the modeling of ductile thin film undergoing a nonlinear peeling process. The bending model is based on different governing parameters: (1) the interfacial fracture toughness and the separation strength, (2) the interfacial fracture toughness and the crack tip slope angle, and (3) the interfacial fracture toughness and the critical Mises effective strain of the delaminated thin film at the crack tip. Thin film nonlinear peeling under steady-state condition is solved with the different governing parameters. In addition, the peeling test problem is simulated by using the elastic-plastic finite element analysis model. A critical assessment of the three analytical bending models is made by comparison of the bending model solutions with the finite element analysis model solutions. Furthermore, through analyses and comparisons for solutions based on both the bending model and the finite element analysis model, some connections between the bending model and the finite element analysis model are developed. Moreover, in the present research, the effect of different selections for cohesive zone shape on the ductile film peeling solutions is discussed.

Seismic Responses of a Hemispherical Alluvial Valley to SV Waves: A Three-Dimensional Analytical Approximation

An analytical solution to the three-dimensional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier-Bessel series expansion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot's dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn: (1) there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; (2) the normalized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; (3) with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.

Using internal fibre geometry to improve the performance of pin-loaded holes in composite materials

The anisotropic nature of fibre reinforced composites leads to large stress concentrations around pin-loaded holes through standard weave cloths. Proper understanding of how this anisotropic nature affects the load distribution around holes can be utilised to reduce these con-centrations if sufficient thought is given to the internal fibre geometry near to the hole. Such local reinforcements need not be highly complex and can be readily produced without excessive effort, producing significant improvements in performance. © 1996 Kluwer Academic Publishers.

The effect of inclined soil layers on surface vibration from underground railways using a semi-analytical approach

Ground vibration due to underground railways is a significant source of disturbance for people living or working near the subways. The numerical models used to predict vibration levels have inherent uncertainty which must be understood to give confidence in the predictions. A semi-analytical approach is developed herein to investigate the effect of soil layering on the surface vibration of a halfspace where both soil properties and layer inclination angles are varied. The study suggests that both material properties and inclination angle of the layers have significant effect ( ± 10dB) on the surface vibration response. © 2009 IOP Publishing Ltd.