DIC Challenge 2.0: Developing Images and Guidelines for Evaluating Accuracy and Resolution of 2D Analyses

Experimental Mechanics - The DIC Challenge 2.0 follows on from the work accomplished in the first Digital Image Correlation (DIC) Challenge Reu et al. (Experimental Mechanics 58(7):1067, 1). The...     

Effect of 2D Image Resolution on 3D Stochastic Reconstruction and Developing Petrophysical Trend

Multi-resolution digital rock physics (DRP) makes it possible to up-scale petrophysical properties from micron size to core sample size using two-dimensional (2D) thin section images. Resolution of 3D images and sample size are challenging problems in DRP where high-resolution images are acquired from small samples using inefficient and expensive micro-CT facilities. Three-dimensional stochastic reconstruction is an alternative approach to overcome these challenges. In this paper, we use multi-resolution images and investigate effect of 2D image resolution on 3D stochastic reconstruction and development of petrophysical trends for our two sandstone and carbonate original representative volume elements (RVEs). The proposed method includes three steps. In the first step, the spatial resolution of our original RVEs is decreased synthetically. In the second step, stochastic RVEs are realized for each resolution using two perpendicular images, correlation functions, and phase recovery algorithm. In the reconstruction method, a full set of two-point correlation functions (TPCFs) is extracted from two perpendicular 2D images. Then TPCF vectors are decomposed and averaged to realize 3D stochastic RVEs. In the third step, petrophysical properties like relative and absolute permeability as well as porosity and formation factor are computed. The output is used to develop trends for petrophysical properties in different resolutions. Experimental results illustrate that the proposed method can be used to predict petrophysical properties and reconstruct 3D RVEs for resolutions unavailable in the acquired 2D or 3D data.     

Guidelines for Modeling a 2D Rough Wall Channel Flow

The effects of 2-D roughness elements on the Reynolds stress anisotropy tensor, and of the energy dissipation rate anisotropy tensor of a turbulent channel flow are investigated using data obtained from direct numerical simulations (DNSs). The roughness elements consist of transverse square rods of size , placed on one wall of the channel only. While is kept constant (, is the half-width of the channel), the spacing between the rods is varied from to . The results show that the variation in can dramatically change the structure of the wall region flow. The modeling of the near-wall region needs to reflect the structural changes caused by the variation in . On the basis of the Reynolds stress budgets, attention needs to be given to the turbulent energy and pressure diffusion terms while local isotropy may be a reasonable approximation for the energy dissipation rate, especially over a range of for which the drag is near its maximum.     

二维Euler方程无网格算法的精度分析

在几种典型的点云布点方案中,分析了文献[3]中的无网格算法的计算精度。分析指出影响无网格算法计算精度的主要因素是点云中点的分布是否平衡;增加点云中节点的数目或者采用高阶曲面拟合都不能提高精度;这里的理论分析为无网格算法的布点提供了重要参考。     

Accuracy and performance of numerical wall boundary conditions for steady, 2D,incompressible streamfunction vorticity

Three recent papers have studied fourth-order compact discretizations of the streamfunction vorticity equations. They differed primarily in how the no-slip wall boundary conditions were handled. In this paper, these different formulas are compared to one another, as well as to three newly proposed formulas. Special consideration is paid to the truncation errors; in particular, it is shown that many well-known formulations are actually more accurate by O(h) than previously reported, where h is the mesh size. These new theoretical error rates are confirmed with an analytical model problem. The different formulas are then compared with published driven cavity results, both in terms of accuracy and performance, and the newly proposed high-order Jensen formula is judged to have the marginally best combination of these properties. © 1998 John Wiley & Sons, Ltd.     

An Out-of-Plane Motion Compensation Strategy for Improving Material Parameter Estimation Accuracy with 2D Field Measurements

In-plane surface displacements, when measured with 2D Digital Image Correlation (2D-DIC), are very sensitive to out-of-plane displacement components. Any out-of-plane motion of the surface can pollute the measured field by introducing artificial displacements. These displacements are difficult to separate from the underlying response of the surface and thereby limit the application of 2D-DIC in inverse problems where the test specimen has significant motion in the out-of-plane direction. In the context of inverse problems, we propose to partially relax this condition of no out-of-plane motion in 2D-DIC. With this approach, only the out-of-plane rigid-body motion of the specimen surface, which is initially in-plane, needs to be avoided while the requirement of surface deformations to be primarily in-plane is essentially waived. Compensation, based on the pinhole camera model, for out-of-plane displacements of the surface in response to applied load is included within the error function of the minimization problem. The improvements in material parameter estimation, obtained by using the proposed compensation strategy, are demonstrated by an example. The proposed technique makes it possible to utilize 2D-DIC with a simple conventional lens for an increased number of inverse problems; and in the process avoiding the computational and experimental difficulties associated with 3D measurement methods as well as the high cost and magnification limitations of a telecentric lens.     

Assessment of High Speed Imaging Systems for 2D and 3D Deformation Measurements: Methodology Development and Validation

Ultra high-speed and moderate speed image acquisition platforms have been characterized, with special emphasis on the variability and accuracy of the measurements obtained when employed in either 2D or 3D computer vision systems for deformation and shape measurements. Specifically, the type of image distortions present in both single channel cameras (HS-CMOS) and multi-channel image intensified cameras (UHS-ICCD) are quantified as part of the overall study, and their effect on the accuracy of experimental measurements obtained using digital image correlation have been determined. Results indicate that established methods for noise suppression and recently developed models for distortion correction can be used effectively in situations where the primary intensity noise components are characterized by minimal cross-talk and stationary spatial distortions. Baseline uniaxial tension experiments demonstrate that image correlation measurements using high speed imaging systems are unbiased and consistent with independent deformation measurements over the same length scale, with point-to-point strain variations that are similar to results obtained from translation experiments. In this study, the point-to-point variability in strain using the image intensified system is on the order of 0.001, whereas the non-intensified system had variability of 0.0001. Results confirm that high speed imaging systems can be utilized for full field two and three-dimensional measurements using digital image correlation methods.

A segregated implicit solution algorithm for 2D and 3D laminar incompressible flows

A segregated algorithm for the solution of laminar incompressible, two- and three-dimensional flow problems is presented. This algorithm employs the successive solution of the momentum and continuity equations by means of a decoupled implicit solution method. The inversion of the coefficient matrix which is common for all momentum equations is carried out through an approximate factorization in upper and lower triangular matrices. The divergence-free velocity constraint is satisfied by formulating and solving a pressure correction equation. For the latter a combined application of a preconditioning technique and a Krylov subspace method is employed and proved more effecient than the approximate factorization method. The method exhibits a monotonic convergence, it is not costly in CPU time per iteration and provides accurate solutions which are independent of the underrelaxation parameter used in the momentum equations. Results are presented in two- and three-dimensional flow problems.     

Integral Identities for a Semi-infinite Interfacial Crack in 2D and 3D Elasticity

The paper is concerned with the problem of a semi-infinite crack at the interface between two dissimilar elastic half-spaces, loaded by a general asymmetrical system of forces distributed along the crack faces. On the basis of the weight function approach and the fundamental reciprocal identity (Betti formula), we formulate the elasticity problem in terms of singular integral equations relating the applied loading and the resulting crack opening. Such formulation is fundamental in the theory of elasticity and extensively used to solve several problems in linear elastic fracture mechanics (for instance various classic crack problems in homogeneous and heterogeneous media).     

On variational formulations with rigid-body constraints for finite elasticity: Applications to 2D and 3D finite element simulations

We investigate a recently proposed variational principle with rigid-body constraints and present an extension of its implementation in three dimensional finite elasticity problems. Through numerical examples, we illustrate that the proposed variational principle with rigid-body constraints is applicable to both single field and mixed finite elements of arbitrary order and geometry, e.g. triangular/tetrahedral and quadrilateral/hexagonal elements, in two and three dimensions. Moreover, we demonstrate that, as compared to the commonly adopted approach of discretizing the rigid domains with standard finite elements, the proposed formulation requires neither discretization nor numerical integration in the interior of each rigid domain. As a comparative result, the variational formulation may reduce the total number of degrees of freedom of the resulting finite element system and provide improved accuracy.     

Uniaxial experimental study of the acoustic emission and deformation behavior of composite rock based on 3D digital image correlation (DIC)

In this paper, uniaxial compression tests were carried out on a series of composite rock specimens with different dip angles, which were made from two types of rock-like material with different strength. The acoustic emission technique was used to monitor the acoustic signal characteristics of composite rock specimens during the entire loading process. At the same time, an optical non-contact 3 D digital image correlation technique was used to study the evolution of axial strain field and the maximal strain field before and after the peak strength at different stress levels during the loading process. The effect of bedding plane inclination on the deformation and strength during uniaxial loading was analyzed. The methods of solving the elastic constants of hard and weak rock were described. The damage evolution process, deformation and failure mechanism, and failure mode during uniaxial loading were fully determined. The experimental results show that the θ = 0?–45?specimens had obvious plastic deformation during loading, and the brittleness of the θ = 60?–90?specimens gradually increased during the loading process. When the anisotropic angle θincreased from 0?to 90?, the peak strength, peak strain,and apparent elastic modulus all decreased initially and then increased. The failure mode of the composite rock specimen during uniaxial loading can be divided into three categories:tensile fracture across the discontinuities(θ = 0?–30?), slid-ing failure along the discontinuities(θ = 45?–75?), and tensile-split along the discontinuities(θ = 90?). The axial strain of the weak and hard rock layers in the composite rock specimen during the loading process was significantly different from that of the θ = 0?–45?specimens and was almost the same as that of the θ = 60?–90?specimens. As for the strain localization highlighted in the maximum principal strain field, the θ = 0?–30?specimens appeared in the rock matrix approximately parallel to the loading direction,while in the θ = 45?–90?specimens it appeared at the hard and weak rock layer interface.     

Formation of a 2D vortex pair and its 3D breakup: an experimental study

 Vortex pairs are studied using a dye tracing visualisation technique and a particle tracking velocimetry system. The vortex pairs are produced by gravity induced inlets of water issued through a uniform gap. The inlet Reynolds number is Re=Ud/ν≈875 in all tests (d being the gap width and U the cross sectional mean velocity), i.e. the flow is in the laminar regime. Initially, the dipolar vortex structure is two-dimensional, but after travelling a distance of a few times its own width, the flow structure becomes unstable, breaks up and changes into a three-dimensional flow structure. The breakup appears to be caused by an axial flow in the core centres of each vortex of the dipolar structure. These axial flows are induced by boundary effects related to the von Karman viscous pump. After the breakup, it is believed that a vortex ring is formed through reconnection of rudiments from the dipolar structure mediated by the wall induced vorticity. Received: 20 November 1995/Accepted: 14 November 1996     

Synthesis of highly convergent 2D and 3D finite elements for short acoustic wave simulation

New higher-order finite elements of enhanced convergence properties for acoustic wave simulation are presented in the paper. The element matrices are obtained by combining modal synthesis and optimization techniques in order to achieve minimum errors of higher modes of the computational domain. As a result, simulation models of propagating wave pulses require a smaller number of finite element divisions per wavelength compared to the conventional element model thus significantly reducing computational costs. Though finite element matrices are obtained in optimization, the resulting patterns of the matrices are versatile and further can be used in any wave propagation model. The mass matrices of the elements are diagonal, so explicit time integration schemes are applicable. The usage of new elements is especially efficient in situations where wavelengths of the simulated signal are much shorter than the dimensions of the computational domain. This is referred to as short wave propagation analysis. The results of wave propagation simulation for ultrasonic measurements are presented as application examples. The B-scans and computed dispersion curves are provided for visual interpretation of the results.     

Comparison of finite element techniques for 2D and 3D crack analysis under impact loading

For several technical applications the dynamic aspect in fracture mechanics cannot be neglected. When the reliability of components with macroscopic cracks has to be assessed, the consideration of dynamic effects may lead to much higher stress intensity factors than under static conditions. In this paper three different methods to calculate the dynamic stress intensity factor for the mode-I loading of stationary cracks are compared. Based on two- and three-dimensional finite element simulations, the dynamic stress intensity factor is computed with the dynamic J-integral, the modified crack closure integral and the displacement interpretation method. The theoretical fundamentals of all three methods are summarized in the paper and the numerical implementation is explained briefly. Results for different models are shown and compared to findings in the literature.     

The Domain-Boundary Element Method (DBEM) for hyperelastic and elastoplastic finite deformation: axisymmetric and 2D/3D problems

Summary  This paper presents the solution of geometrically nonlinear problems in solid mechanics by the Domain-Boundary Element Method. Because of the Total-Lagrange approach, the arising domain and boundary integrals are evaluated in the undeformed configuration. Therefore, the system matrices remain unchanged during the solution procedure, and their time-consuming computation needs to be performed only once. While the integral equations for axisymmetric finite deformation problems will be derived in detail, the basic ideas of the formulation in two and three dimensions can be found in [1]. The present formulation includes torsional problems with finite deformations, where additional terms arise due to the curvilinear coordinate system. A Newton–Raphson scheme is used to solve the nonlinear set of equations. This involves the solution of a large system of linear equations, which has been a very time-consuming task in former implementations, [1, 2]. In this work, an iterative solver, i.e. the generalized minimum residual method, is used within the Newton–Raphson algorithm, which leads to a significant reduction of the computation time. Finally, numerical examples will be given for axisymmetric and two/three-dimensional problems. Received 29 August 2000; accepted for publication 10 October 2000     

压电介质边界元法及奇异性处理

本文从压电材料的基本方程出发,利用功的互等原理推导了边界积分方程,并详细地讨论了边界元的计算步骤,利用等参变换,着重研究了在边界元计算中基本解的奇异性问题,对各种情况讨论了系数矩阵H和G的算法,并给出院具体的表达式,作为算例,选取了均匀薄板和开孔薄板PZY-4压电材料,计算结果表明,本文提出的边界元的计算格式和奇异性的处理方法相当有效。     

3D Digital Volume Correlation of Synchrotron Radiation Laminography Images of Ductile Crack Initiation: An Initial Feasibility Study

A feasibility study of measuring 3D displacement fields in the bulk during ductile crack initiation via combined Synchrotron Radiation Computed Laminography (SRCL) and Digital Volume Correlation (DVC) is performed. In contrast to tomography, SRCL is a technique that is particularly adapted to obtain three-dimensional (3D) reconstructed volumes of objects that are laterally extended (i.e., in 2 directions) and thin in the third direction, i.e. sheet-like objects. In-situ laminography data of an initiating crack ahead of a machined notch are used with a voxel size of 0.7 μm. The natural contrast of the observed 2XXX Al-alloy caused by intermetallic particles and initial porosity is used to measure displacement fields via a global DVC technique assuming a continuous displacement field. An initial performance study is carried out on data of the same undeformed material but after a substantial shift of the laminography rotation axis with respect to the imaged specimen. Volume correlations between different loading steps provide displacement fields that are qualitatively consistent with the remote loading conditions. Computed strain fields display a strain concentration close to the notch tip.