Showing papers on "Displacement field published in 2012"

Automated identification and indexing of dislocations in crystal interfaces

Abstract: We present a computational method for identifying partial and interfacial dislocations in atomistic models of crystals with defects. Our automated algorithm is based on a discrete Burgers circuit integral over the elastic displacement field and is not limited to specific lattices or dislocation types. Dislocations in grain boundaries and other interfaces are identified by mapping atomic bonds from the dislocated interface to an ideal template configuration of the coherent interface to reveal incompatible displacements induced by dislocations and to determine their Burgers vectors. In addition, the algorithm generates a continuous line representation of each dislocation segment in the crystal and also identifies dislocation junctions.

A Continuum Phase Field Model for Fracture

Abstract: A variational free-discontinuity formulation of brittle fracture was given by Francfortand Marigo [1], where the total energy is minimized with respect to the crackgeometry and the displacement field simultaneously. The entire evolution of cracksincluding their initiation and branching is determined by this minimization principlerequiring no further criterion. However, a direct numerical discretization of themodel faces considerable difficulties as the displacement field is discontinuous inthe presence of cracks.

A new higher order shear deformation theory for sandwich and composite laminated plates

Abstract: A new shear deformation theory for sandwich and composite plates is developed. The proposed displacement field, which is “m” parameter dependent, is assessed by performing several computations of the plate governing equations. Therefore, the present theory, which gives accurate results, is relatively close to 3D elasticity bending solutions. The theory accounts for adequate distribution of the transverse shear strains through the plate thickness and tangential stress-free boundary conditions on the plate boundary surface, thus a shear correction factor is not required. Plate governing equations and boundary conditions are derived by employing the principle of virtual work. The Navier-type exact solutions for static bending analysis are presented for sinusoidally and uniformly distributed loads. The accuracy of the present theory is ascertained by comparing it with various available results in the literature.

Refined beam elements with only displacement variables and plate/shell capabilities

Abstract: This paper proposes a refined beam formulation with displacement variables only. Lagrange-type polynomials, in fact, are used to interpolate the displacement field over the beam cross-section. Three- (L3), four- (L4), and nine-point (L9) polynomials are considered which lead to linear, quasi-linear (bilinear), and quadratic displacement field approximations over the beam cross-section. Finite elements are obtained by employing the principle of virtual displacements in conjunction with the Unified Formulation (UF). With UF application the finite element matrices and vectors are expressed in terms of fundamental nuclei whose forms do not depend on the assumptions made (L3, L4, or L9). Additional refined beam models are implemented by introducing further discretizations over the beam cross-section in terms of the implemented L3, L4, and L9 elements. A number of numerical problems have been solved and compared with results given by classical beam theories (Euler-Bernoulli and Timoshenko), refined beam theories based on the use of Taylor-type expansions in the neighborhood of the beam axis, and solid element models from commercial codes. Poisson locking correction is analyzed. Applications to compact, thin-walled open/closed sections are discussed. The investigation conducted shows that: (1) the proposed formulation is very suitable to increase accuracy when localized effects have to be detected; (2) it leads to shell-like results in case of thin-walled closed cross-section analysis as well as in open cross-section analysis; (3) it allows us to modify the boundary conditions over the cross-section easily by introducing localized constraints; (4) it allows us to introduce geometrical boundary conditions along the beam axis which lead to plate/shell-like cases.

Continuum phase field modeling of dynamic fracture: variational principles and staggered FE implementation

Abstract: The modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations of complex crack topologies including branching. This drawback can be overcome by a diffusive crack modeling based on the introduction of a crack phase field as proposed in Miehe et al. (Comput Methods Appl Mech Eng 19:2765–2778, 2010a; Int J Numer Meth Eng 83:1273–1311, 2010b), Hofacker and Miehe (Int J Numer Meth Eng, 2012). In this work, we summarize basic ingredients of a thermodynamically consistent, variational-based model of diffusive crack propagation under quasi-static and dynamic conditions. It is shown that all coupled field equations, in particular the balance of momentum and the gradient-type evolution equation for the crack phase field, follow as the Euler equations of a mixed rate-type variational principle that includes the fracture driving force as the mixed field variable. This principle makes the proposed formulation extremely compact and provides a perfect basis for the finite element implementation. We then introduce a local history field that contains a maximum energetic crack source obtained in the deformation history. It drives the evolution of the crack phase field. This allows for the construction of an extremely robust operator split scheme that updates in a typical time step the history field, the crack phase field and finally the displacement field. We demonstrate the performance of the phase field formulation of fracture by means of representative numerical examples, which show the evolution of complex crack patterns under dynamic loading.

Correlation of multi-temporal ground-based optical images for landslide monitoring: Application, potential and limitations

Abstract: The objective of this work is to present a low-cost methodology to monitor the displacement of continuously active landslides from ground-based optical images analyzed with a normalized image correlation technique. The performance of the method is evaluated on a series of images acquired on the Super-Sauze landslide (South French Alps) over the period 2008–2009. The image monitoring system consists of a high resolution optical camera installed on a concrete pillar located on a stable crest in front of the landslide and controlled by a datalogger. The data are processed with a cross-correlation algorithm applied to the full resolution images in the acquisition geometry. Then, the calculated 2D displacement field is orthorectified with a back projection technique using a high resolution DEM interpolated from Airborne Laser Scanning (ALS) data. The heterogeneous displacement field of the landslide is thus characterized in time and space. The performance of the technique is assessed using differential GPS surveys as reference. The sources of error affecting the results are then discussed. The strongest limitations for the application of the technique are related to the meteorological, illumination and ground surface conditions inducing partial or complete loss of coherence among the images. Small movements of the camera and the use of a mono-temporal DEM are the most important factors affecting the accuracy of the ortho-rectification of the displacement field. As the proposed methodology can be routinely and automatically applied, it offers promising perspectives for operational applications like, for instance, in early warning systems.

A refined shear deformation theory for free vibration of functionally graded plates on elastic foundation

Abstract: A refined shear deformation theory for free vibration of functionally graded plates on elastic foundation is developed. The displacement field is chosen based on assumptions that the in-plane and transverse displacements consist of bending and shear components, and the shear components of in-plane displacements give rise to the parabolic variation of shear strain through the thickness in such a way that shear stresses vanish on the plate surfaces. Therefore, there is no need to use shear correction factor. Material properties of functionally graded plate are assumed to vary according to power law distribution of the volume fraction of the constituents. The elastic foundation is modeled as Pasternak foundation. Equations of motion are derived using Hamilton’s principle. Closed-form solution of rectangular plates is derived, and the obtained results are compared well with three-dimensional elasticity solutions and third-order shear deformation theory solutions. Finally, the influences of power law index, thickness ratio, foundation parameter, and boundary condition on the natural frequency of plates have been investigated.

Analysis of Sandwich Beams With a Compliant Core and With In-Plane Rigidity—Extended High-Order Sandwich Panel Theory Versus Elasticity

Abstract: A new one-dimensional high-order theory for orthotropic elastic sandwich beams is formulated. This new theory is an extension of the high-order sandwich panel theory (HSAPT) and includes the in-plane rigidity of the core. In this theory, in which the compressibility of the soft core in the transverse direction is also considered, the displacement field of the core has the same functional structure as in the high-order sandwich panel theory. Hence, the transverse displacement in the core is of second order in the transverse coordinate and the in-plane displacements are of third order in the transverse coordinate. The novelty of this theory is that it allows for three generalized coordinates in the core (the axial and transverse displacements at the centroid of the core and the rotation at the centroid of the core) instead of just one (midpoint transverse displacement) commonly adopted in other available theories. It is proven, by comparison to the elasticity solution, that this approach results in superior accuracy, especially for the cases of stiffer cores, for which cases the other available sandwich computational models cannot predict correctly the stress fields involved. Thus, this theory, referred to as the “extended high-order sandwich panel theory” (EHSAPT), can be used with any combinations of core and face sheets and not only the very “soft” cores that the other theories demand. The theory is derived so that all core=face sheet displacement continuity conditions are fulfilled. The governing equations as well as the boundary conditions are derived via a variational principle. The solution procedure is outlined and numerical results for the simply supported case of transverse distributed loading are produced for several typical sandwich configurations. These results are compared with the corresponding ones from the elasticity solution. Furthermore, the results using the classical sandwich model without shear, the first-order shear, and the earlier HSAPT are also presented for completeness. The comparison among these numerical results shows that the solution from the current theory is very close to that of the elasticity in terms of both the displacements and stress or strains, especially the shear stress distributions in the core for a wide range of cores. Finally, it should be noted that the theory is formulated for sandwich panels with a generally asymmetric geometric layout. [DOI: 10.1115/1.4005550]

Multiresolution MR elastography using nonlinear inversion

Abstract: Purpose: Nonlinear inversion (NLI) in MR elastography requires discretization of the displacement field for a finite element(FE) solution of the “forward problem”, and discretization of the unknown mechanical property field for the iterative solution of the “inverse problem”. The resolution requirements for these two discretizations are different: the forward problem requires sufficient resolution of the displacement FE mesh to ensure convergence, whereas lowering the mechanical property resolution in the inverse problem stabilizes the mechanical property estimates in the presence of measurement noise. Previous NLI implementations use the same FE mesh to support the displacement and property fields, requiring a trade-off between the competing resolution requirements. Methods: This work implements and evaluates multiresolution FE meshes for NLI elastography, allowing independent discretizations of the displacements and each mechanical property parameter to be estimated. The displacement resolution can then be selected to ensure mesh convergence, and the resolution of the property meshes can be independently manipulated to control the stability of the inversion. Results: Phantom experiments indicate that eight nodes per wavelength (NPW) are sufficient for accurate mechanical property recovery, whereas mechanical property estimation from 50 Hzin vivobrain data stabilizes once the displacement resolution reaches 1.7 mm (approximately 19 NPW). Viscoelasticmechanical property estimates of in vivobrain tissue show that subsampling the loss modulus while holding the storage modulus resolution constant does not substantially alter the storage modulusimages. Controlling the ratio of the number of measurements to unknown mechanical properties by subsampling the mechanical property distributions (relative to the data resolution) improves the repeatability of the property estimates, at a cost of modestly decreased spatial resolution. Conclusions: Multiresolution NLI elastography provides a more flexible framework for mechanical property estimation compared to previous single mesh implementations.

Morphable displacement field based image matching for face recognition across pose

Abstract: Fully automatic Face Recognition Across Pose (FRAP) is one of the most desirable techniques, however, also one of the most challenging tasks in face recognition field. Matching a pair of face images in different poses can be converted into matching their pixels corresponding to the same semantic facial point. Following this idea, given two images G and P in different poses, we propose a novel method, named Morphable Displacement Field (MDF), to match G with P's virtual view under G's pose. By formulating MDF as a convex combination of a number of template displacement fields generated from a 3D face database, our model satisfies both global conformity and local consistency. We further present an approximate but effective solution of the proposed MDF model, named implicit Morphable Displacement Field (iMDF), which synthesizes virtual view implicitly via an MDF by minimizing matching residual. This formulation not only avoids intractable optimization of the high-dimensional displacement field but also facilitates a constrained quadratic optimization. The proposed method can work well even when only 2 facial landmarks are labeled, which makes it especially suitable for fully automatic FRAP system. Extensive evaluations on FERET, PIE and Multi-PIE databases show considerable improvement over state-of-the-art FRAP algorithms in both semi-automatic and fully automatic evaluation protocols.

Identification of plastic constitutive parameters at large deformations from three dimensional displacement fields

Abstract: The aim of this paper is to provide a general procedure to extract the constitutive parameters of a plasticity model starting from displacement measurements and using the Virtual Fields Method. This is a classical inverse problem which has been already investigated in the literature, however several new features are developed here. First of all the procedure applies to a general three-dimensional displacement field which leads to large plastic deformations, no assumptions are made such as plane stress or plane strain although only pressure-independent plasticity is considered. Moreover the equilibrium equation is written in terms of the deviatoric stress tensor that can be directly computed from the strain field without iterations. Thanks to this, the identification routine is much faster compared to other inverse methods such as finite element updating. The proposed method can be a valid tool to study complex phenomena which involve severe plastic deformation and where the state of stress is completely triaxial, e.g. strain localization or necking occurrence. The procedure has been validated using a three dimensional displacement field obtained from a simulated experiment. The main potentialities as well as a first sensitivity study on the influence of measurement errors are illustrated.

Evaluation of generalized continuum substitution models for heterogeneous materials

Abstract: Several extensions of standard homogenization methods for composite materials have been proposed in the literature that rely on the use of polynomial boundary conditions enhancing the classical affine conditions on the unit cell. Depending on the choice of the polynomial, overall Cosserat, second gradient, or micromorphic homogeneous substitution media are obtained. They can be used to compute the response of the composite when the characteristic length associated with the variation of the applied loading conditions becomes of the order of the size of the material inhomogeneities. A significant difference between the available methods is the nature of the fluctuation field added to the polynomial expansion of the displacement field in the unit cell, which results in different definitions of the overall stress and strain measures and higher order elastic moduli. The overall higher order elastic moduli obtained from some of these methods are compared in the present contribution in the case of a specific periodic two-phase composite material. The performance of the obtained overall substitution media is evaluated for a chosen boundary value problem at the macroscopic scale for which a reference finite element solution is available. Several unsatisfactory features of the available theories are pointed out, even though some model predictions turn out to be highly relevant. Improvement of the prediction can be obtained by a precise estimation of the fluctuation at the boundary of the unit cell.

A time‐domain recursive method to analyse transient wave propagation across rock joints

Abstract: The present investigation is concerned with transient wave propagation in a rock mass with a set of parallel joints by using a recursive method. According to the displacement field of a rock mass with a set of parallel joints, the interaction between four plane waves (two longitudinal-waves and two transverse-waves) and a joint is analysed first. With the displacement discontinuity model and the time shifting function, the wave propagation equation based on the recursive method in time domain for obliquely longitudinal-(P) or transverse-(S) waves across a set of parallel joints is established. The joints are assumed linearly elastic. The analytical solution obtained by the proposed method is compared with the existing results for some special cases, including oblique incidence across a single joint and normal incidence across a set of parallel joints. By verification, it is found that the solutions by the proposed method match very well with the existing methods. The applicability and limitations of the new method are then discussed for incident waves with different properties.

Resummed propagators in multicomponent cosmic fluids with the eikonal approximation

Abstract: We introduce the eikonal approximation to study the effect of the large-scale motion of cosmic fluids on their small-scale evolution. This approach consists in collecting the impact of the long-wavelength displacement field into a single or finite number of random variables, whose statistical properties can be computed from the initial conditions. For a single dark matter fluid, we show that we can recover the nonlinear propagators of renormalized perturbation theory. These are obtained with no need to assume that the displacement field follows the linear theory. Then we extend the eikonal approximation to many fluids. In particular, we study the case of two nonrelativistic components and we derive their resummed propagators in the presence of isodensity modes. Unlike the adiabatic case, where only the phase of small-scale modes is affected by the large-scale advection field, the isodensity modes change also the amplitude on small scales. We explicitly solve the case of cold dark matter-baryon mixing and find that the isodensity modes induce only very small corrections to the resummed propagators.

Lagrangian perturbations and the matter bispectrum I: fourth-order model for non-linear clustering

Abstract: We investigate the Lagrangian perturbation theory of a homogeneous and isotropic universe in the non-relativistic limit, and derive the solutions up to the fourth order. These solutions are needed for example for the next-to-leading order correction of the (resummed) Lagrangian matter bispectrum, which we study in an accompanying paper. We focus on flat cosmologies with a vanishing cosmological constant, and provide an in-depth description of two complementary approaches used in the current literature. Both approaches are solved with two different sets of initial conditions — both appropriate for modelling the large-scale structure. Afterwards we consider only the fastest growing mode solution, which is not affected by either of these choices of initial conditions. Under the reasonable approximation that the linear density contrast is evaluated at the initial Lagrangian position of the fluid particle, we obtain the nth-order displacement field in the so-called initial position limit: the nth order displacement field consists of 3(n-1) integrals over n linear density contrasts, and obeys self-similarity. Then, we find exact relations between the series in Lagrangian and Eulerian perturbation theory, leading to identical predictions for the density contrast and the peculiar-velocity divergence up to the fourth order.

Displacement fields from point cloud data: Application of particle imaging velocimetry to landslide geodesy

Abstract: [1] Acquiring spatially continuous ground-surface displacement fields from Terrestrial Laser Scanners (TLS) will allow better understanding of the physical processes governing landslide motion at detailed spatial and temporal scales. Problems arise, however, when estimating continuous displacement fields from TLS point-clouds because reflecting points from sequential scans of moving ground are not defined uniquely, thus repeat TLS surveys typically do not track individual reflectors. Here, we implemented the cross-correlation-based Particle Image Velocimetry (PIV) method to derive a surface deformation field using TLS point-cloud data. We estimated associated errors using the shape of the cross-correlation function and tested the method's performance with synthetic displacements applied to a TLS point cloud. We applied the method to the toe of the episodically active Cleveland Corral Landslide in northern California using TLS data acquired in June 2005–January 2007 and January–May 2010. Estimated displacements ranged from decimeters to several meters and they agreed well with independent measurements at better than 9% root mean squared (RMS) error. For each of the time periods, the method provided a smooth, nearly continuous displacement field that coincides with independently mapped boundaries of the slide and permits further kinematic and mechanical inference. For the 2010 data set, for instance, the PIV-derived displacement field identified a diffuse zone of displacement that preceded by over a month the development of a new lateral shear zone. Additionally, the upslope and downslope displacement gradients delineated by the dense PIV field elucidated the non-rigid behavior of the slide.

Simulation-Based Joint Estimation of Body Deformation and Elasticity Parameters for Medical Image Analysis

Abstract: Estimation of tissue stiffness is an important means of noninvasive cancer detection. Existing elasticity reconstruction methods usually depend on a dense displacement field (inferred from ultrasound or MR images) and known external forces. Many imaging modalities, however, cannot provide details within an organ and therefore cannot provide such a displacement field. Furthermore, force exertion and measurement can be difficult for some internal organs, making boundary forces another missing parameter. We propose a general method for estimating elasticity and boundary forces automatically using an iterative optimization framework, given the desired (target) output surface. During the optimization, the input model is deformed by the simulator, and an objective function based on the distance between the deformed surface and the target surface is minimized numerically. The optimization framework does not depend on a particular simulation method and is therefore suitable for different physical models. We show a positive correlation between clinical prostate cancer stage (a clinical measure of severity) and the recovered elasticity of the organ. Since the surface correspondence is established, our method also provides a nonrigid image registration, where the quality of the deformation fields is guaranteed, as they are computed using a physics-based simulation.

Acoustofluidics 3: Continuum mechanics for ultrasonic particle manipulation.

Abstract: Although ultrasonic particle manipulation is performed in fluids, the adjacent solid material and its modelling play an important part in the design of such devices, as the impedance difference between liquids and solids is not very large. The modelling of the solid is done in the framework of linear elastodynamics, which is described by the displacement field, the strain field and the stress field. These quantities are related by kinematical relations, local linear momentum conservation and constitutive laws. Also material damping is important and can be modelled in the framework of linear viscoelasticity. Because of the finite size of the devices, resonant modes are important to analyse. In addition some of the elements behave as mechanical structures obeying specific equations that take into account the structural boundaries. A practical example is given, which shows the full complexity of device modelling, which is preferably done using the Finite Element Method.