Tikhonov regularization with Incremental Hole-Drilling and the Integral Method in Cross-Ply Composite Laminates

Background: Incremental hole-drilling with the integral method has been extensively used in composite laminates but is sensitive to small measurement errors. Error sensitivity can be reduced by limiting the number of depth increments used in the calculation procedure. This approach is limited if a rapidly varying residual stress distribution exists since the calculated stress in each incremental depth is considered constant. Distortion of stress results can consequently occur due to averaging effects if the depth increments become too large. Tikhonov regularization is usually applied in isotropic materials to smooth the resulting residual stress distribution and reduce stress uncertainties, but has only been applied to composite laminates using the slitting technique. Objective: The intention of this work is to extend the use of Tikhonov regularization to incremental hole-drilling of composite laminates using the integral method. Methods: Finite element modelling is used to calculate the necessary calibration coefficients for unit pulses of uniform stress. Monte Carlo simulation is used to the determine uncertainties in the calculated residual stress distributions. Tikhonov regularization is optimised to reduce the stress uncertainty, while ensuring that the stress solution is not distorted. Results: The method is demonstrated on a GFRP (Glass fibre reinforced plastic) laminate of [0₂/90₂]_s construction and the calculated residual stress field is compared with those obtained by the standard integral method and series expansion. Conclusions: It is found that Tikhonov regularization significantly improves the accuracy of the standard integral method in composite laminates and shows good agreement with the series expansion method.

     

Measurement of the Distribution of Residual Stresses in Layered Thick-Walled GFRP Pipes

The objective of this study is to measure the axial, circumferential, shear and radial residual stress distributions in three thick-walled glass fibre reinforced plastic (GFRP) filament-wound pipes, two of which are layered. The measurement of residual stresses was carried out using a recently published layer removal method which overcomes the limitations of previous techniques and can be applied to layered anisotropic pipes of any wall thickness. Layers of approximately 0.3 mm thickness were incrementally ground from the outer surface of the pipes. The resulting strains were measured on the inner surfaces. A least-squares polynomial was fitted to each measured data set, and used to calculate the corresponding stress distributions. All of the resulting axial, hoop and shear stress distributions adhere to the requirement of self-equilibrium and the radial stress distributions all vanish to zero at the inner and outer surfaces. The radial stresses of the layered pipes showed a tendency to have two peaks, one for each layer, a consequence of the two-stage manufacturing process of these pipes. The measured axial and hoop stresses of all three pipes were similar at the inner surfaces despite significant differences in the stiffnesses in the principal directions arising from different wind angles.     

A general method for obtaining shear stress and normal stress functions from parallel disk rheometry data

The problems of converting the torque and normal force versus rim shear rate data generated by parallel disk rheometers into shear stress and normal stress difference as functions of shear rate are formulated as two independent integral equations of the first kind. Tikhonov regularization is used to obtain approximate solutions of these equations. This way of handling parallel disk rheometer data has the advantage that it is independent of the rheological constitutive equation and noise amplification is kept under control by the user-specified parameter in Tikhonov regularization. If the fluid under test exhibits a yield stress, Tikhonov regularization computation will simultaneously give an estimate of the yield stress. The performance of this method is demonstrated by applying it to a number of data sets taken from the published literature and to laboratory measurements conducted specifically for this investigation.     

Mapping Multiple Components of the Residual Stress Tensor in a Large P91 Steel Pipe Girth Weld Using a Single Contour Cut

The contour method is applied in an innovative manner to measure the distribution of hoop residual stress in a large martensitic-ferritic steel pipe containing a multi-pass girth weld. First, a novel one-step wire electro-discharge machining cut is conducted to divide the pipe lengthways into two halves. The deformation of the cut halves is then measured and analysed in a way that simultaneously gives maps of hoop stress across the wall thickness on both sides of the pipe and automatically accounts for through-thickness hoop bending effects and how they may vary along the pipe. Finally the contour method results are combined with X-ray diffraction residual stress measurements using the principle of superposition to determine the distribution of the axial and radial residual stresses in the pipe. It is thereby demonstrated how the distribution of three direct components of the residual stress tensor in a welded pipe can be readily determined using a “hybrid” contour measurement approach.     

A smoothed inverse eigenstrain method for reconstruction of the regularized residual fields

A smoothed inverse eigenstrain method is developed for reconstruction of residual field from limited strain measurements. A framework for appropriate choice of shape functions based on the prior knowledge of expected residual distribution is presented which results in stabilized numerical behavior. The analytical method is successfully applied to three case studies where residual stresses are introduced by inelastic beam bending, laser-forming and shot peening. The well-rehearsed advantage of the proposed eigenstrain-based formulation is that it not only minimizes the deviation of measurements from its approximations but also will result in an inverse solution satisfying a full range of continuum mechanics requirements. The smoothed inverse eigenstrain approach allows suppressing fluctuations that are contrary to the physics of the problem. Furthermore, a comprehensive discussion is performed on regularity of the asymptotic solution in the Tikhonov scheme and the regularization parameter is then exactly determined utilizing Morozov discrepancy principle. Gradient iterative regularization method is also examined and shown to have an excellent convergence to the Tikhonov–Morozov regularization results.     

High resolution COD image analysis for health monitoring of reinforced concrete structures through inverse analysis

Direct and inverse problems of a fracture mechanics based RC beam model are solved. Solution of the direct problem that maps crack bridging stresses into crack opening displacements (COD) is straightforward, but the inverse problem is ill-posed, and better solved by the theory of inverse problems. This paper exploits the Tikhonov regularization method to solve the inverse problem, and estimates the force and location of rebar in buried concrete from CODs. Bending tests are carried out on model RC beams in the laboratory to demonstrate the applicability of the method. During the tests, a microscopic camera snaps high resolution digital pictures of cracked concrete surface. The images are analyzed by a software to measure surface CODs that are input into the inverse problem. The practical CODs inevitably include noise due to experimental error, which makes the inverse problem ill-posed, and necessitates regularization. In the current inverse analysis by the Tikhonov regularization method, bridging stress profiles, i.e. variation of the crack bridging stress along the crack length, has been figured out. Results are compared with those from other theoretical methods of analysis as well as with the readings from strain gauges. The method is a suitable non-destructive means for existing structures in cases where the section information is inadequate, or damages/repairs have altered the designed cross-section.     

Using regularization methods for the determination of relaxation and retardation spectra of polymeric liquids

The application of Tikhonov regularization to the determination of relaxation and retardation spectra of viscoelastic fluids is discussed using simulated and real experimental data. It is thereby shown that with this method consistent results for relaxation and retardation spectra can be obtained from experimental data for different material functions. Furthermore, the differences between the most frequently used variants of Tikhonov regularization and the maximum entropy method are discussed. For most calculations the program FTIKREG has been used, which is an especially fast and reliable implementation of Tikhonov's regularization method.     

Mechanics and full-field deformation study of the Ring Hoop Tension Test

We investigate the mechanics of the Ring Hoop Tension Test (RHTT), as a means to assess the properties of anisotropic tubes in the hoop direction. This test involves placing a ring extracted from the tube over two close-fitting D-shaped mandrels that are then parted using a universal testing machine. Since the curvature of the ring does not change during loading, the ring undergoes only stretching. We determine the effects of contact pressure, radial stress, and friction between the tube and mandrels using FEA simulations. The effects of the pre-existing thickness eccentricity and of the specimen-making procedure on the recorded RHTT response are also assessed with a combination of experiments and analysis. We tested tubes from Al-6061-T4 with a diameter/thickness ratio of 20. We found that as the friction increases beyond 0.1–0.15, the state of uniaxial tension is deteriorated, indicating that care must be taken to minimize the tube-mandrel friction. We determined that although these tubes have a relatively large thickness eccentricity (±4% of the nominal thickness), this had no effect on the recorded results. We showed that the tubes should not be turned to remove that eccentricity, as the machining process induces damage that is noticeable in the results. We found that the contact pressure and the contact-induced radial stress cause limited deviations from uniaxial tension, comparable to the case of a tube under axial load and internal pressure which is often used for assessing the material properties in the hoop direction. Central in our analyses is the knowledge of the hoop strain field, which was assessed using 3D Digital Image Correlation. We propose a data reduction procedure for RHTT that accounts for all the above effects. Finally, with all effects accounted for, we establish the anisotropy of the extruded Al-6061-T4 tubes studied.     

Parametric study on single shot peening by dimensional analysis method incorporated with finite element method

Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and resistance to corrosion fatigue,cracking,etc.Compressive residual stress and dent profile are important factors to evaluate the effectiveness of shot peening process.In this paper,the influence of dimensionless parameters on maximum compressive residual stress and maximum depth of the dent were investigated.Firstly,dimensionless relations of processing parameters that affect the maximum compressive residual stress and the maximum depth of the dent were deduced by dimensional analysis method.Secondly,the influence of each dimensionless parameter on dimensionless variables was investigated by the finite element method.Furthermore,related empirical formulas were given for each dimensionless parameter based on the simulation results.Finally,comparison was made and good agreement was found between the simulation results and the empirical formula,which shows that a useful approach is provided in this paper for analyzing the influence of each individual parameter.     

Intralaboratory Repeatability of Residual Stress Determined by the Slitting Method

This paper presents repeated slitting method measurements of the residual stress versus depth profile through the thickness of identically prepared samples, which were made to assess repeatability of the method. Measurements were made in five 17.8 mm thick blocks cut from a single plate of 316L stainless steel which had been uniformly laser peened to induce a deep residual stress field. Typical slitting method techniques were employed with a single metallic foil strain gage on the back face of the coupon and incremental cutting by wire EDM. Measured residual stress profiles were analyzed to assess variability of residual stress as a function of depth from the surface. The average depth profile had a maximum magnitude of −668 MPa at the peened surface. The maximum variability also occurred at the surface and had a standard deviation of 15 MPa and an absolute maximum deviation of 26 MPa. Since measured residual stress exceeded yield strength of the untreated plate, microhardness versus depth profiling and elastic–plastic finite element analysis were combined to bound measurement error from inelastic deformation.     

Determination of the relaxation time spectrum from dynamic moduli using an edge preserving regularization method

For the numerical determination of a continuous relaxation time spectrum from rheological data so-called regularization methods are necessary. The regularization method based on Tikhonov regularization and commonly used in rheology assumes the relaxation time spectrum to be smooth. Accordingly, this method is not able to properly resolve spectra containing edges or at least large curvatures as proposed in the literature. In order to reconstruct such spectra a novel edge preserving regularization (EPR) method was developed which extends the properties of the commonly used method. In order to demonstrate the efficiency of EPR it is compared with the commonly used regularization method by means of simulated and experimental data. This reveals that EPR is able to prove whether or not the theoretically predicted edge containing spectra for high-molecular-weight, monodisperse polymers are compatible with experimental data. Received: 11 June 1999/Accepted: 1 September 1999     

Application of X-ray diffraction, micromagnetic and hole drilling methods for residual stress determination in a ball bearing steel ring

A basic understanding of distortion problems requires the analysis of a complete manufacturing process including an almost complete overview of residual stress states in the component during each production step. To reduce the measurement time in the future, three measurements methods (X-ray diffraction, micromagnetic and blind hole drilling methods) have been used to analyze residual stress states in machined AISI 52100 ball bearing rings. X-ray diffraction was used as a state-of-the-art method for machining induced residual stresses with pronounced gradients. The ring exhibited a complex residual stress state with high tensile residual stresses at the surface, a strong gradient in depth, and also showed some variation along the outer circumference due to a superimposition of machining induced residual stresses and effects from the clamping device process. Due to this surface state, micromagnetic signals depend on the analyzing frequency. A calibration of the signals was only possible with the X-ray diffraction data. The results of the three different measurement methods correlate reasonably well.     

Investigation of the Dynamic Behavior of a Griffith Crack in a Piezoelectric Material Strip Subjected to the Harmonic Elastic Anti-plane Shear Waves by Use of the Non-local Theory

In this paper, the dynamic behavior of a Griffith crack in a piezoelectric material strip subjected to the harmonic anti-plane shear waves is investigated by use of the non-local theory for impermeable crack surface conditions. To overcome the mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress and the electric displacement near at the crack tip. By means of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations. These equations are solved using the Schmidt method. Contrary to the classical solution, it is found that no stress and electric displacement singularity is present near the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress near the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the thickness of the strip, the circular frequency of incident wave and the lattice parameter.     

Identification of a moving boundary for a heat conduction problem in a multilayer medium

In this paper, we give a uniqueness theorem for the moving boundary of a heat problem in a composite medium. Through solving the Cauchy problem of heat equation in each subdomain, we finally find an approximation to the moving boundary for one-dimensional heat conduction problem in a multilayer medium. The numerical scheme is based on the use of the method of fundamental solutions and a discrete Tikhonov regularization technique with the generalized cross-validation choice rule for a regularization parameter. Numerical experiments for five examples show that our proposed method is effective and stable.     

Effect of diffusion-induced bending on diffusion-induced stress near the end faces of an elastic hollow cylinder

Diffusion-induced stress plays an important role in determining structural integrity of mechanical structures used in lithium-ion batteries and microelectromechanical devices. Incorporating the diffusion-induced bending in the analysis of the diffusion-induced stress in an elastic hollow cylinder, analytical forms of the diffusion-induced resultant axial stress and hoop stress have been formulated for the traction-free condition and the built-in condition at the end faces of the cylinder. Using these results, the evolution of the diffusion-induced stress at the end faces of a hollow, elastic electrode due to the insertion of lithium is discussed under the potentiostatic operation. The end faces of the electrode experience compressive hoop stress through the thickness in contrast to the stress state in the hollow cylinder far away from the end faces. The magnitude of the resultant hoop stress decreases with increasing the diffusion time for the traction-free end faces; it increases with increasing the diffusion time near the inner surface for the built-in end faces.     

刚粘塑性圆环受法向冲击时的动力屈曲

本文研究刚粘塑性圆环受法向冲击时的动力屈曲，文中考虑了冲击过程中圆环厚度的变化，环向力的变化以及基本变形速度的变化。通过用数值方法求解扰动运动微分方程，结果表明，当规定了临界放大函数后，仅当冲击速度位于某一范围之内时，才会发生动力屈曲。     

爆破地震波作用下城市隧道结构动力响应的敏感性

选取纵波作为研究对象,采用波函数展开法,推导了无限弹性介质中复合衬砌结构隧道在爆破地震波作用下动力响应问题的解析解,结合南京红山南路隧道群开挖工程开展了隧道结构动应力集中因数的敏感性分析,并通过曲线拟合分别得到了隧道围岩及内衬环向动应力集中因数随隧道结构动力学及几何参数变化的回归方程。分析结果表明:隧道围岩和二次衬砌层的弹性模量及围岩泊松比对考察点处环向动应力的影响较大,而初期支护层弹性模量的影响几乎可以忽略不计;衬砌层厚度变化对围岩环向动应力的影响要远大于内衬,在隧道施工设计时可以通过增加衬砌层的厚度平衡隧道结构的受力状态,但增加初期支护层的效果不如二次衬砌层。     

The MFS for the Cauchy problem in two-dimensional steady-state linear thermoelasticity

We study the reconstruction of the missing thermal and mechanical fields on an inaccessible part of the boundary for two-dimensional linear isotropic thermoelastic materials from over-prescribed noisy (Cauchy) data on the remaining accessible boundary. This problem is solved with the method of fundamental solutions (MFS) together with the method of particular solutions (MPS) via the MFS-based particular solution for two-dimensional problems in uncoupled thermoelasticity developed in Marin and Karageorghis, 2012a, Marin and Karageorghis, 2013. The stabilisation/regularization of this inverse problem is achieved by using the Tikhonov regularization method (Tikhonov and Arsenin, 1986), whilst the optimal value of the regularization parameter is selected by employing Hansen’s L-curve method (Hansen, 1998).     

网络RTK参考站间模糊度动态解算中病态方程的一种解算方法

在网络RTK参考站间的模糊度估计中,若误差方程严重病态,将导致模糊度解与其准确值偏差较大或整周模糊度无法固定,因此提出了一种适于网络RTK模糊度动态解算的新方案：1）法方程病态性的判断;2）Tikhonov正则化解算病态方程;3）LAMBDA方法搜索固定整周模糊度。同时,深入研究了Tikhonov正则化矩阵的构造方法和正则化参数的选取准则。最后以实例验证了采用此方案解算病态方程是可行的,通过选取合适的正则化参数可以解得准确的整周模糊度;详细讨论了选择不同的正则化参数对模糊度解算结果的影响。     

Obtaining the shear stress shear rate relationship and yield stress of liquid foods from Couette viscometry data

A method based on Tikhonov regularisation is used to convert the Couette viscometry data of a number of liquid foods into shear stress vs shear rate relationships. For liquid foods that have a yield stress and if the viscometry data cover the appropriate range of shear stress, Tikhonov regularisation is used to compute simultaneously the yield stress. The versatility of Tikhonov regularisation is demonstrated by applying it to process data from wide as well as narrow gap Couette viscometers. The results generated by Tikhonov regularisation are compared against that obtained by methods currently employed by rheologists.