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.     

Crack identification in elastically restrained vibrating rods

In this paper we consider the problem of identifying an open crack in a longitudinally vibrating rod with smooth variable profile by minimal eigenfrequency data. The crack is assumed to be open during vibration and it is modelled by an elastic spring acting along the rod axis. Most, if not all, the results available in the literature for this inverse problem refer to ideal end conditions, that is the rod is either under free or supported end conditions. As an example of almost optimal result, it is known that the knowledge of the fundamental (positive) natural frequency of the rod under free-free and cantilever end conditions allows for the unique determination of the crack, without any restriction on the damage severity. In this paper we show that the analysis of the analogous crack identification problem for rods under elastically restrained end conditions leads to different results and that, in general, the knowledge of the fundamental frequency belonging to two spectra associated to different end conditions is not sufficient for the uniqueness of the solution. The method we used to solve the inverse problem is of constructive type and it is based on general properties of the eigenfrequencies as functions of the position and severity of the crack. The identification procedure has been tested numerically on rods under various damage scenarios. Numerical results agree well with the theory, even in presence of noisy input data.     

An inversion integral for crack-scattering data

The inverse problem of determining the size, shape and orientation of a flat crack from high-frequency far-field elastic waves scattered by the crack is investigated. The results show that desired information on a crack can be obtained from the first arriving scattered longitudinal waves only. It is shown that an approximate high-frequency solution to the direct problem, based on physical elastodynamics, yields an expression for the scattered far-field of longitudinal motion which suggests a solution to the inverse problem by application of Fourier-type inversion integrals to scattering data. Two kinds of inversion integrals are examined. The inversion problem becomes relatively simple if some a-priori information is available, either on the orientation of the plane of the crack or on a plane of symmetry. The method of inversion is verified for a flat crack of elliptical shape. Some computational technicalities are discussed, and the method is also applied to experimental scattering data.     

Inverse parameter identification of cohesive zone model for simulating mixed-mode crack propagation

Inverse analysis is widely applied to the identification of material properties or model parameters. In order to improve the computational efficiency of the inverse method based on the genetic algorithm, an interpolation scheme upon the response surface constructed by the finite element simulation has been adopted in this paper. Meanwhile, a gradual homogenization treatment scheme has also been presented to improve the convergence of the inverse method based on the Kalman filter algorithm. Both methods are proven effective in dealing with the single-objective inverse problem. However, literature studies show that the adoption of multiple types of experimental information is useful to improve the accuracy of inverse analysis. In this case, it turns into a multiple-objective inverse problem. Our practice proved that the above-mentioned two methods might not yield a proper result if the sensitivity issue of different types of information is not considered. Therefore, another multi-objective inverse method, in combination of the above two optimization algorithms and a weight-estimating scheme that can consider such sensitivity, has been further presented. Finally, by using a mixed-mode crack propagation simulation and two types of experimental information (loading-displacement response curve and crack path profile), the parameters of the cohesive zone model were inversely identified and its simulation results are in good agreement with the experiment.     

The inversion of acoustical impedance profile by methods of characteristics

The inverse problem for the one-dimensional propagation of waves through a medium with an unknown impedance profile is considered. From the reflected pulses caused by an excitation of short duration, this profile is reconstructed by applying the method of characteristics with finite difference approximations. This leads to a simple and fast algorithm, which is demonstrated by the numerical examples in the paper.     

各向异性介质中SH波引起的裂纹扩展

本文利用Green函数法,求解各向异性介质中半无限长裂纹在SH波作用下,以任意速度扩展的问题。首先,利用Laplace变换和Cagniard-de Hoop反演法求解各向异性介质中反平面问题的Green函数,并利用它建立了求解裂纹扩展问题的积分方程。因为方程为Abel型的,所以可得到在SH波作用下,半无限长裂纹扩展问题的解析解。还可求得裂纹端点附近的应力和裂纹表面上位移的表达式。并对裂纹端点附近的奇异性进行讨论。最后讨论了裂纹尖端附近任一点的能量关系。并应用Griffith的能量准则,对裂纹扩展规律进行了讨论。     

Theory of inverse problems for crack bridging stresses determination

A new mathematical procedure has been proposed for the determination of crack bridging stresses from perturbed crack opening displacements (COD) in fiber composites. The problem is an ill-posed inverse problem and the proposed procedure exploited concepts from functional analysis and the theory of inverse problems. The transformation from crack bridging stress into COD has been linearized, where the continuous crack bridging stresses functions of infinite dimensional vector spaces were approximated into finite dimensional subspaces. The coordinate representation of functions and the matrix of the transformation facilitated a numerical solution to the inverse problem. Self consistent direct and inverse problems were solved numerically in context to an example. Results establish that the proposed mathematical procedure produces fairly accurate results for engineering application.     

Permeable cracks between two dissimilar piezoelectric materials

An interfacial crack with electrically permeable surfaces between two dissimilar piezoelectric ceramics under electromechanical loading is investigated. An exact expression for singular stress and electric fields near the tip of a permeable crack between two dissimilar anisotropic piezoelectric media are obtained. The interfacial crack-tip fields are shown to consist of both an inverse square root singularity and a pair of oscillatory singularities. It is found that the singular fields near the permeable interfacial crack tip are uniquely characterized by the real valued stress intensity factors proposed in this paper. The energy release rate is obtained in terms of the stress intensity factors. The exact solution of stress and electric fields for a finite interfacial crack problem is also derived.     

Parameter identification for cracks in elastic plate structures based on remote strain fields

A method for the detection of cracks in plate structures is presented. In contrast to most of the common monitoring concepts taking advantage of the reflection of elastic waves at crack faces, the presented approach is based on the strain measured at different locations on the surface of the structure. This allows both the identification of crack position parameters, such as length, location and angles with respect to a reference coordinate system and the calculation of stress intensity factors (SIF). The solution of the direct problem is performed on the basis of the BFM (body force method). The inverse problem is solved applying the particle swarm optimization (PSO) algorithm. The BFM is based on the principle of linear superposition which allows the calculation of the strain field in a cracked body. The strain at an arbitrary point in the structure is replaced by the strain provided by body force doublets in the uncracked structure. The doublets as well as external loads are parameters which have to be determined solving the inverse problem by minimizing a fitness function, which is defined by a square sum of residuals between measured strain distributions and computed ones for an assumed crack. The PSO algorithm applied to the fitness function operates on the basis of a swarm of candidate solutions. Once knowing loading and crack parameters, the SIF can be determined.     

Finite element analysis of a bimaterial interface crack

In this investigation, the enriched element method developed by Benzley was extended to treat the stress analysis problem involving a bimaterial interface crack. Unlike crack problems in isotropic elasticity, where the stress singularity at the crack tip is of the inverse square root type, the interface crack contains an additional oscillatory singularity. Although the effect of this oscillatory characteristic is confined to a region very close to the crak tip, it nevertheless requires proper treatment in order to obtain accurate predictions on the stress intensity factors. Using appropriate crack tip stress and displacement expressions, the enriched element method can model the stress singularity for an interface crack exactly. The finite element implementation of this method has been made on the code APES. Stress intensity factor results predicted by the modified APES program compare favorably with those available in the literature. This indicates tha the enriched element technique provides an accurate and efficient numerical tool for the analysis of bimaterial interface crack problems.     

LONGITUDINAL SHEAR PROBLEMS OF COLLINEAR RIGID LINE INCLUSIONS IN ANISOTROPIC MATERIALS

ＬＯＮＧＩＴＵＤＩＮＡＬＳＨＥＡＲＰＲＯＢＬＥＭＳＯＦＣＯＬＬＩＮＥＡＲＲＩＧＩＤＬＩＮＥＩＮＣＬＵＳＩＯＮＳＩＮＡＮＩＳＯＴＲＯＰＩＣＭＡＴＥＲＩＡＬＳＪｉａｎｇｃｎｉ－ｐｉｎｇ（蒋持平）（ＤｅｐａｒｉｔｍｅｎｔｏｆＦｌｉｇｈｔＶｅｈｉｃｌｅＤｅｓ...     

Solution of multiple crack problem in a finite plate using coupled integral equations

This paper studies a numerical solution of multiple crack problem in a finite plate using coupled integral equations. After using the principle of superposition, the multiple crack problem in a finite plate can be converted into two problems: (a) the multiple crack problem in an infinite plate and (b) a usual boundary value problem for the finite plate. For the former problem, the Fredholm integral equation is used. For the latter problem, a BIE based on complex variable is suggested in which a Cauchy singular kernel exists. For the proposed BIE, after using the inverse matrix technique, the dependence of the traction at a domain point from the boundary tractions is formulated indirectly. This is a particular advantage of the present study. Several numerical examples are provided and the computed results for stress intensity factor and T-stress at crack tips are given.     

Two non-elliptical inhomogeneities with internal uniform stresses interacting with a mode-III crack

We use the method of Green's functions to analyze an inverse problem in which we aim to identify the shapes of two non-elliptical elastic inhomogeneities, embedded in an infinite matrix subjected to uniform remote stress, which enclose uniform stress distributions despite their interaction with a finite mode-III crack. The problem is reduced to an equivalent Cauchy singular integral equation, which is solved numerically using the Gauss–Chebyshev integration formula. The shapes of the two inhomogeneities and the corresponding location of the crack can then be determined by identifying a conformal mapping composed in part of a real density function obtained from the solution of the aforementioned singular integral equation. Several examples are given to demonstrate the solution.     

A numerical analysis of cracks emanating from a surface elliptical hole in infinite body in tension

This paper deals with such a kind of surface crack problem with an approximately same depth, which is called a liked-plane crack problem. Based on the previous investigations on internal rectangular crack and surface rectangular crack in infinite solid in tension and a hybrid displacement discontinuity method (a boundary element method) proposed recently by Yan, a numerical approach for the liked-plane crack problem in hand is presented. Numerical examples are given to illustrate the numerical approach is simple, yet accurate for calculating the SIFs of a liked-plane crack. Specifically, a pair of cracks emanating from a surface elliptical hole in infinite body in tension are investigated in detail.     

基于压电振动能量俘获的弯曲结构损伤监测研究

建立了含裂纹损伤的曲梁压电能量俘获系统在强迫振动下的动力学模型. 基于Prescott型压电曲梁力电耦合振动方程的解析解和裂纹截面处的连续性条件, 求解了含裂纹损伤的压电曲梁的格林函数. 根据线性叠加原理, 对含裂纹的力电耦合模型的系统方程解耦, 得到强迫振动下含裂纹损伤的曲梁压电俘能器的输出电压. 在得到模型的强迫振动解析解后, 提出逆方法检测结构中的裂纹损伤, 这一检测方法适用于处于振动状态下的结构. 在数值计算中, 令裂纹深度为零, 通过对比本文的解析解与现有文献中的解析解, 验证了本文解的有效性. 分别分析了含裂纹损伤的压电曲梁的电压响应与裂纹深度、裂纹位置、材料的几何参数以及阻尼之间的关系. 研究结果表明: 裂纹的存在对曲梁式压电俘能器的影响比直梁式更加复杂; 裂纹出现时, 损伤曲梁在健康曲梁的一阶频率值处一定会出现波动并被激励出二阶频率, 此时的二阶频率是开路中健康压电曲梁的一阶频率值; 通过对电压响应的检测可以确定的损伤裂纹的深度和在结构中出现的位置范围; 利用振动问题的解来检测压电曲梁的健康状况是可行且准确的.      

Discontinuous deformation gradients away from the tip of a crack in anti-plane shear

Summary This paper examines the finite elastostatic problem for a crack under anti-plane shear conditions. A class of incompressible, homogeneous, isotropic, elastic materials is considered for which the governing displacement equation of equilibrium is elliptic at both sufficiently small and sufficiently large shearing strains but is hyperbolic at an intermediate level of strains. The small-scale nonlinear crack problem is considered and an exact solution is found through use of the hodograph method.The principal feature of the solution is the presence of two pairs of curves issuing from points on the crack-faces, different from the crack-tips, across which the displacement gradients and stresses suffer jump discontinuities.The results communicated in this paper were obtained in the course of an investigation supported by the NSF Materials Research Laboratory at Brown University.     

THE EVALUATION OF STRESS INTENSITY FACTORS OF PLANE CRACK FOR ORTHOTROPIC PLATE WITH EQUAL PARAMETER BY F2LFEM

In this paper, the evaluation of stress intensity factor of plane crack problems for orthotropic plate of equal-parameter is investigated using a fractal two-level finite element method (F2LFEM). The general solution of an orthotropic crack problem is obtained by assimilating the problem with isotropic crack problem, and is employed as the global interpolation function in F2LFEM. In the neighborhood of crack tip of the crack plate, the fractal geometry concept is introduced to achieve the similar meshes having similarity ratio less than one and generate an infinitesimal mesh so that the relationship between the stiffness matrices of two adjacent layers is equal. A large number of degrees of freedom around the crack tip are transformed to a small set of generalized coordinates. Numerical examples show that this method is efficient and accurate in evaluating the stress intensity factor (SIF).     

Hydrodynamic methods in the inverse problem of gravitational prospecting

In [1, 2], a dynamical method is proposed for solving stationary inverse problems of potential theory, including the inverse problem of gravitational prospecting. It is based on analogy with the problem of establishing the interface of two immiscible fluids flowing in a porous medium. In the present paper, a system of two functional equations is derived from which one can obtain, as special cases, an equation corresponding to the method of [1, 2], and also a system of equations that enables one to propose a new and different method for solving the inverse problem of gravitational prospecting. Equations are derived in polar coordinates for plane Cauchy problems corresponding to both methods, and the results are also given of the solution of some model problems by these methods. Finally, ways of generating new methods of solution of the inverse problem of gravitational prospecting are considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 63–71, July–August, 1980.     

裂纹技术的基本原理与应用研究进展

本文介绍裂纹技术的基本原理及其应用。裂纹技术通过使裂纹沿给定路线扩展去达到切割材料的目的,它由魏庆同和郎福元在文[2]中首次提出。本文还概要讨论了裂纹技术中四个典型的力学问题,特别是断裂力学的逆问题。笔者希望这篇文章能够引起国内外学者对裂纹技术及其问题的关注。      

The dynamic stress intensity factor analysis of adhesively bonded material interface crack with damage under shear loading

This paper studies the dynamic stress intensity factor (DSIF) at the interface in an adhesive joint under shear loading. Material damage is considered. By introducing the dislocation density function and using the integral transform, the problem is reduced to algebraic equations and can be solved with the collocation dots method in the Laplace domain. Time response of DSIF is calculated with the inverse Laplace integral transform. The results show that the mode Ⅱ DSIF increases with the shear relaxation parameter, shear module and Poisson ratio, while decreases with the swell relaxation parameter. Damage shielding only occurs at the initial stage of crack propagation. The singular index of crack tip is -0.5 and independent on the material parameters, damage conditions of materials, and time. The oscillatory index is controlled by viscoelastic material parameters.