Numerical solutions for general inverse heat conduction problem

A family of numerical methods for determining the space-and time-variable heat transfer coefficient, based on experimentally acquired interior temperature-time data, is presented. Newton-type methods are utilized to compute simultaneously the unknown heat transfer coefficient components. To reduce the influence of random errors in the measurement data on the estimated heat transfer coefficients, the noisy data are smoothed using least squares approximation by cubic splines. Three test examples using experimental and random simulated data are used to illustrate the computation efficiency and generality of the present methods.     

Towards High Performance Digital Volume Correlation

We develop speed and efficiency improvements to a three-dimensional (3D) digital volume correlation (DVC) algorithm, which measures displacement and strain fields throughout the interior of a material. Our goal is to perform DVC with resolution comparable to that achieved in 2D digital image correlation, in time that is commensurate with the image acquisition time. This would represent a significant improvement over the current state-of-the-art available in the literature. Using an X-ray micro-CT scanner, we can resolve features at the 5 micron scale, generating 3D images with up to 36 billion voxels. We compute twelve degrees-of-freedom at each correlation point and utilize tricubic spline interpolation to achieve high accuracy. We improve the algorithm’s speed and robustness through an improved coarse search, efficient implementation of spline interpolation, and using smoothing splines to address noisy image data. For DVC, the volume of data, number of correlation points, and work to solve each correlation point grow cubically. We therefore employ parallel computing to handle this tremendous increase in computational and memory requirements. We demonstrate the application of DVC using simulated deformations of 3D micro-CT scans of polymer samples with embedded particles forming an internal pattern.     

Events Maps in a Stick-Slip System

This paper describes a one-dimensional map generated by a two degree-of-freedom mechanical system that undergoes self-sustained oscillations induced by dry friction. The iterated map allows a much simpler representation and a better understanding of some dynamic features of the system. Some applications of the map are illustrated and its behaviour is simulated by means of an analytically defined one-dimensional map. A method of reconstructing one-dimensional maps from experimental data from the system is introduced. The method uses cubic splines to approximate the iterated mappings. From a sequence of such time series the parameter dependent bifurcation behaviour is analysed by interpolating between the defined mappings. Similarities and differences between the bifurcation behaviour of the exact iterated mapping and the reconstructed mapping are discussed.     

Laser surface-contouring and spline data-smoothing for residual stress measurement

We describe non-contact scanning with a confocal laser probe to measure surface contours for application to residual stress measurement. (In the recently introduced contour method, a part is cut in two with a flat cut, and the part deforms by relaxation of the residual stresses. A cross-sectional map of residual stresses is then determined from measurement of the contours of the cut surfaces.) The contour method using laser scanning is validated by comparing measurements on a ferritic steel (BS 4360 grade 50D) weldment with neutron diffraction measurements on an identical specimen. Compared to lower resolution touch probe techniques, laser surface-contouring allows more accurate measurement of residual stresses and/or measurement of smaller parts or parts with lower stress levels. Furthermore, to take full advantage of improved spatial resolution of the laser measurements, a method to smooth the surface contour data using bivariate splines is developed. In contrast to previous methods, the spline method objectively selects the amount of smoothing and estimates the uncertainties in the calculated residual stress map.     

数字图像相关法测量局域变形场中形函数和模板尺寸的影响

数字图像相关法测量均匀变形场已被普遍接受, 其测量结果可与应变片测量结果比较. 然而, 在工程测量中, 针对局域变形场(应变高度集中, 如波特文-勒夏特利埃带、试件缺口附近和裂纹尖端等), 应变片受限于其尺寸, 其测量结果是接触面内的平均应变值. 此时, 采用数字图像相关法能够测量这些局域变形场. 但形函数和模板尺寸等计算参数对计算结果影响很大, 这也导致使用者很难判断计算结果的可靠性. 论文通过对合金拉伸实验获得的不同应变梯度的波特文-勒夏特利埃带和模拟生成的带的计算分析, 发掘了形函数和模板尺寸作用于计算结果的深层机制, 证明了二阶形函数比一阶形函数更适用于高度非均匀的局域变形场. 提出了在局域应变场测量中, 当一阶和二阶形函数计算结果的相对误差小于10% 时, 二阶形函数的结果是可靠的判据.      

Techniques and Applications of the Simulated Pattern Adaptation of Wilkinson’s Method for Advanced Microstructure Analysis and Characterization of Plastic Deformation

Electron Backscatter Diffraction (EBSD) based Orientation Imaging Microscopy (OIM) is used routinely at ~500 materials laboratories worldwide for the characterization and development of diverse crystalline materials. Statistically significant data sets (~10⁷ individual EBSD measurements) can be collected and analyzed within time periods of acceptable beam stability (~10⁵s). However, limitations in angular and spatial resolution have motivated a continued search for more robust EBSD-based methods. Herein is a gathered presentation of advanced techniques in use, intended as a guide to researchers in selecting the most appropriate method for their work. Wilkinson’s method has been shown to increase angular resolution nearly two orders of magnitude to ±0.006°, facilitating measurement of elastic strain, lattice curvature, and dislocation density. A simulated pattern adaptation of Wilkinson’s method extends these measurement capabilities to polycrystalline materials, by avoiding the need for an experimental strain free reference pattern. The angular resolution limit obtained is ~0.04°. Accurate pattern center calibration, essential to the high resolution methods, is accomplished by parallelization of band edges projected onto a sphere centered at the interaction volume. FFT powered cross-correlation functions improve the spatial resolution near grain boundaries and correct for measurement inaccuracies induced by overlapping patterns. To corroborate these claims, exemplary results taken from a wedge-indented nickel single crystal, cold-worked copper polycrystal, and rolled nickel polycrystal are shown.     

On the Accuracy of Elastic Strain Field Measurements by Laue Microdiffraction and High-Resolution EBSD: a Cross-Validation Experiment

Determining the accuracy of elastic strain measurements in plastically deformed alloys is an experimental challenge. To develop a novel cross-validation procedure, a controlled elasto-plastic strain gradient was created in a stainless steel single crystal by four point bending deformation. The corresponding elastic strain field was probed, with an intragranular spatial resolution, in-situ by Laue microdiffraction and ex-situ by High Resolution EBSD. Good agreement is found for the two independent measurements and the predictions of a mechanical model, at plastic strains below 0.5 %. The accuracy of the measurements is estimated at 3.2 × 10^{? 4}.     

Spline collocation approach to boundary value problems

A well-known Stokes problem is discussed by a cubic spline collocation method. Two consecutive cubic splines are obtained for the problem. The results by this method are compared with those of an orthogonal collocation method. The selection of the length of the subintervals of the range of the boundary value problem is also justified. The results obtained by these two methods are compared with the analytic solution. The methods involve simple algebra, and hence the calculations do not require the help of a computer. Necessary error analysis has been carried out.     

An optimisation method for separating and rebuilding one-dimensional dispersive waves from multi-point measurements. Application to elastic or viscoelastic bars

When using a classical SHPB (split Hopkinson pressure bar) set-up, the useful measuring time is limited by the length of the bars, so that the maximum strain which can be measured in material testing applications is also limited. In this paper, a new method with no time limits is presented for measuring the force and displacement at any station on a bar from strain or velocity measurements performed at various places on the bar. The method takes the wave dispersion into account, as must inevitably be done when making long time measurements. It can be applied to one-dimensional and single-mode waves of all kinds propagating through a medium (flexural waves in beams, acoustic waves in wave guides, etc.). With bars of usual sizes, the measuring time can be up to 50 times longer than the time available with classical methods. An analysis of the sensitivity of the results to the accuracy of the experimental data and to the quality of the wave propagation modelling was also carried out. Experimental results are given which show the efficiency of the method.     

金属表面改性中界面力学性能的实验研究

提出了适用于微区测量的两步数字图像相关法,采用二阶泰勒多项式子区位移模式有效地降低了截断误差;采用B-样条插值函数提高了亚像素重建的精度;采用拟牛顿迭代法提高了计算速度等,使微区内应变测量精度控制在100~150με。在残余应变测量中通过数字标记点相关技术实现了试件的精确复位,提高了其测量精度。文中利用上述方法测量了镀层改性试件跨界面的位移场和应变场,分析了裂纹对镀层的影响以及测出了基材中引起的残余应变;同时对镀层改性界面结合性能作了定量比较。     

基于余弦样条函数的有限点阵法

应用节点影响域的概念,提出了基于余弦样条函数的有限点阵方法.利用余弦样条函数的性质,通过张量积的形式构造余弦样条函数正则解空间,用于逼近场函数,余弦样条基的线性组合使得边界条件处理如同有限元法一样方便.余弦样条与边界型方法结合,可用于求解不规则域问题.数值实例的计算结果表明,文中方法避免了高阶多项式构造形函数所带来的数值振荡,解的连续性不受限制,为改进计算精度而加密点阵所导致的计算量增加量较少,计算收敛快.     

Use of 3-D Digital Image Correlation to Characterize the Mechanical Behavior of a Fiber Reinforced Refractory Castable

Refractory castables exhibit very low fracture strain levels when subjected to tension or bending. The main objective of this work is to show that 3-D digital image correlation (3-D DIC) allows such low strain levels to be measured. Compared to mechanical extensometer measurements, 3-D DIC makes it possible to reach similar strain resolution levels and to avoid the problem of position dependance related to the heterogeneous nature of the strain and to strain localization phenomena. First, the 3-D DIC method and the experimental set-up are presented. Secondly, an analysis of the 3-D DIC method is performed in order to evaluate the resolution, the standard uncertainty and the spatial resolution for both displacement and strain measurements. An optimized compromise between strain spatial resolution and standard uncertainty is reached for the configuration of the experimental bending test. Finally, the macroscopic mechanical behavior of a fiber reinforced refractory castable (FRRC) is studied using mechanical extensometry and 3-D DIC in the case of tensile and four-point bending tests. It is shown that similar results are obtained with both methods. Furthermore, in the case of bending tests on damaged castable, 3-D DIC results demonstrate the ability to determine Young’s modulus from heterogeneous strain fields better than by using classical beam deflection measurements.     

The accuracy of residual stress measurement by the hole-drilling method

This paper describes the verification of the accuracy of residual stress measurement by the hole-drilling method. The strain measurement is simulated by the use of the indirect fictitious-boundary integral method. As an example, a finite rectangular plate subjected to initial stress is treated, and a simulated measurement of the residual stress is made using the strain relieved during hole drilling. The accuracy of residual stress measurement is estimated by comparing the simulated measured residual stress with the actual residual stress, i.e., the given initial stress. The results are shown for various distances and angles of strain gages. Also, the influences of the eccentricity of the hole from the center of the strain gages and the effect of a boundary near the hole are examined.     

Strain in transition joints measured by high-resolution moiré photography

The strain distribution produced by heating in a plate cut from a welded austenitic-ferritic steel transition joint has been measured using high resolution moiré photography Two methods of doing this are discussed one using a random pattern on the specimen surface and the other using a regular, orthogonal grating stencilled on to the surface. The strains parallel and perpendicular to the interface are presented and show strain peaks near the interface. The effect on strain of introducing a simulated crack along the interface was also examinec. The measurements are compared with the strains computed from a finite-element analysis.     

Sensitivity of in-Plane Strain Measurement to Calibration Parameter for out-of-Plane Specimen Rotations

In practice, out-of-plane motions usually are not avoidable during experiments. Since 2D–DIC measurements are vulnerable to parasitic deformations due to out-of-plane specimen motions, three-dimensional digital image correlation (StereoDIC or 3D–DIC) oftentimes is employed. The StereoDIC method is known to be capable of accurate deformation measurements for specimens subjected to general three-dimensional motions, including out-of-plane rotations and displacements. As a result, there has been limited study of the deformation measurements obtained when using StereoDIC to measure the displacement and strain fields for a specimen subjected only to out-of-plane rotation. To assess the accuracy of strain measurements obtained using stereovision systems and StereoDIC when a specimen undergoes appreciable out of plane rotation, rigid body out-of-plane rotation experiments are performed in the range ?40⁰?≤?θ?≤?40⁰ using a two-camera stereovision system. Results indicate that (a) for what would normally be considered “small angle” calibration processes, the measured normal strain in the foreshortened specimen direction due to specimen rotation increases in a non-linear manner with rotation angle, with measurement errors exceeding ±1400με and (b) for what would normally be considered “large angle” calibration processes, the magnitude of the errors in the strain are reduced to ±300με. To theoretically assess the effect of calibration parameters on the measurements, two separate analyses are performed. First, theoretical strains due to out-of-plane rigid body rotation are determined using a pinhole camera model to project a series of three-dimensional object points into the image plane using large angle calibration parameters and then re-project the corresponding sensor plane coordinates back into the plane using small angle calibration parameters. Secondly, the entire imaging process is also simulated in order to remove experimental error sources and to further validate the theory. Results from both approaches confirmed the same strain error trends as the experimental strain measurements, providing confidence that the source of the errors is the calibration process. Finally, variance based sensitivity analyses show that inaccuracy in the calibrated stereo angle parameter is the most significant factor affecting the accuracy of the measured strain.     

Direct determination of bulk strain soliton parameters in solid polymeric waveguides

Bulk strain solitons in solids attract considerable attention in applications due to their very small decay, permanent bell shape and the wave parameter dependence on waveguide elasticity and geometry. One of the reasons of an evident gap between theory and numerical simulation of strain solitons propagation in various waveguides and comparatively rare experimental verification of rigorous results lies in extensive variability of physical constants data available for polymers. We show how a dramatic improvement of experimental setup provides new opportunities in solitary wave observation and its parameters measurements.Digital holography based on high-speed registration cameras allowed us to refine the accuracy of measurements, and precise pulse synchronization provided direct measurements of bulk strain soliton velocity with proper accuracy. It confirms the fact that the soliton parameters depend rather on the waveguide geometry and material, not on an initial pulse, which power provides either a single soliton or a soliton train.     

Nonlinear Eulerian thermoelasticity for anisotropic crystals

A complete continuum thermoelastic theory for large deformation of crystals of arbitrary symmetry is developed. The theory incorporates as a fundamental state variable in the thermodynamic potentials what is termed an Eulerian strain tensor (in material coordinates) constructed from the inverse of the deformation gradient. Thermodynamic identities and relationships among Eulerian and the usual Lagrangian material coefficients are derived, significantly extending previous literature that focused on materials with cubic or hexagonal symmetry and hydrostatic loading conditions. Analytical solutions for homogeneous deformations of ideal cubic crystals are studied over a prescribed range of elastic coefficients; stress states and intrinsic stability measures are compared. For realistic coefficients, Eulerian theory is shown to predict more physically realistic behavior than Lagrangian theory under large compression and shear. Analytical solutions for shock compression of anisotropic single crystals are derived for internal energy functions quartic in Lagrangian or Eulerian strain and linear in entropy; results are analyzed for quartz, sapphire, and diamond. When elastic constants of up to order four are included, both Lagrangian and Eulerian theories are capable of matching Hugoniot data. When only the second-order elastic constant is known, an alternative theory incorporating a mixed Eulerian–Lagrangian strain tensor provides a reasonable approximation of experimental data.     

Distortion Correction Protocol for Digital Image Correlation after Scanning Electron Microscopy: Emphasis on Long Duration and Ex-Situ Experiments

Recently, scanning electron microscopy (SEM) has been used for digital image correlation (DIC), as micrographs can be acquired with high magnification, providing improved resolution to quantify strain heterogeneities. However, it has been shown by researchers that SEM images can be problematic due to inherent electromagnetic distortions that are not present in optical images. Drift, spatial distortions, and magnification uncertainties are the main issues that can seriously affect the accuracy of localized strain measurements. The present work focuses on long duration experiments, for which images are taken days or weeks apart. We have proposed a systematic procedure to reduce drift, correct spatial distortion, and account for magnification variations between pairs of acquired images. Additionally, SEM parameters are discussed and chosen to increase the signal-to-noise ratio and improve the accuracy of the DIC measurements. The spatial distortion correction increases the repeatability of the correlated values and the precision of the measurements. The implementation for this type of correction is done by applying the measured distortion gradient of a certified grid onto the DIC strain field. The magnification adjustment increases the reliability of the strain maps, ensuring the measurements are in agreement with the actual strain induced during the experiment. We have presented a systematic protocol for ex-situ DIC experiments within the SEM and some basic cross-check procedures that can be performed to evaluate the reliability of the reference grid and the precision of the final strain map.     

Analysis of strain localization during tensile tests by digital image correlation

This paper presents an imaging technique developed to study the strain localization phenomena that occur during the tension of thin, flat steel samples. The data are processed using digital speckle image correlation to derive the two in-plane components of the displacement vectors. The authors observe that the calculation of the intercorrelation function reveals a systematic error and propose a numerical method to limit its influence. Plastic incompressibility and thin-sheet assumptions are used to derive the third displacement component and, hence, the various strain and strain rate components. Numerous checks are presented at each step in processing the data to determine the final accuracy of the strain measurements. It is estimated that this accuracy is quite sufficient to track the inception and the development of localization. Examples of possible application are presented for mild steels whose strain localization mechanisms appear to be precocious and gradual.