热、机械载荷作用下夹杂对应力强度因子的影响

推导了远场应力、热应力耦合作用下含夹杂裂纹体的应力强度因子求解公式,改进了体积力法中的裂纹面合力平衡条件,将应力强度因子的求解归结为解一组积分方程,再将积分方程转化为线性方程组进行数值求解。算例分析结果表明方法正确、有效。将该算法应用于含Al2O3夹杂的FGH95材料应力强度因子分析中,计算结果表明热应力对应力强度因子影响很小。     

On the stress state of a piezoceramic body with a flat crack under symmetric loads

A static-equilibrium problem is solved for an electroelastic transversely isotropic medium with a flat crack of arbitrary shape located in the plane of isotropy. The medium is subjected to symmetric mechanical and electric loads. A relationship is established between the stress intensity factor (SIF) and electric-displacement intensity factor (EDIF) for an infinite piezoceramic body and the SIF for a purely elastic material with a crack of the same shape. This allows us to find the SIF and EDIF for an electroelastic material directly from the corresponding elastic problem, not solving electroelastic problems. As an example, the SIF and EDIF are determined for an elliptical crack in a piezoceramic body assuming linear behavior of the stresses and the normal electric displacement on the crack surface __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 11, pp. 67–77, November 2005.     

基于修正PSO-UKF的SINS/GPS组合导航滤波算法

针对噪声时变特性引起滤波精度下降的问题,提出了一种基于修正粒子群技术( PSO)的自适应UKF算法.为了克服传统粒子群算法过早收敛,容易陷入局部最优的问题,基于粒子的适应值方差提出了一种惯性权值实时修正算法,有效改善了传统PSO算法.在使用新息序列对观测噪声进行实时跟踪的同时,通过构造合理的适应度函数将修正PSO算法和...     

An application of particle swarm optimization to identify the temperature-dependent heat capacity

An inverse analysis of estimating temperature-dependent heat capacity for a heat conduction problem is presented. The heat capacities are assumed to be temperature-dependent of first and second-order polynomial functions, respectively. A fitness function is minimized by implementing the particle swarm optimization algorithm for the inverse analysis. Computation results are provided to demonstrate the performance of the particle swarm optimization algorithm. The comparison to the modified genetic algorithm is also included.     

Equilibrium of an internal transverse crack in a semiinfinite elastic body with thin coating

Static elasticity problems for a half-plane and a strip weakened by a rectilinear transverse crack are studied. In each case, the upper boundary of the body is reinforced by a flexible patch. Various versions of conditions on the lower boundary are considered in the case of the strip. The crack is maintained in the open state by distributed normal forces. The method of generalized integral transforms reduces solving the problem for the equations of equilibriumto solving a singular integral equation of the first kind with the Cauchy kernel with respect to the derivative of the crack opening function. The solutions of the integral equation are constructed by the small parameter and collocation methods for various combinations of the geometric and physical parameters of the problem, and the structure of the solutions is analyzed. The values of the stress intensity factor (SIF) near the crack vertex are obtained.     

Thermostressed state of a piezoelectric bodywith a plane crack under symmetric thermal load

The paper addresses a thermoelectroelastic problem for a piezoelectric body with an arbitrarily shaped plane crack in a plane perpendicular to the polarization axis under a symmetric thermal load. A relationship between the intensity factors for stress (SIF) and electric displacement (EDIF) in an infinite piezoceramic body with a crack under a thermal load and the SIF for a purely elastic body with a crack of the same shape under a mechanical load is established. This makes it possible to find the SIF and EDIF for an electroelastic material from the elastic solution without the need to solve specific problems of thermoelasticity. The SIF and EDIF for a piezoceramic body with an elliptic crack and linear distribution of temperature over the crack surface are found as an example __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 3, pp. 96–108, March 2008.     

基于Kriging代理模型的迭代更新高效反演方法

为提高混凝土坝等大体积结构参数反演效率和精度,减少由于应用有限元进行大量正分析而产生的计算机时,建立了一种结合Kriging代理模型和粒子群优化(PSO)算法的迭代更新反演方法。通过拉丁超立方抽样(LHS)方法确定初始样本点的空间分布,并使用有限元正分析获取对应的响应值,构建粗糙的初始代理模型,结合具有全局寻优能力的PSO算法,反演大体积结构的分区弹性模量,随之再代入有限元模型中,计算获取新的位移响应,并将其作为新样本加入到样本集中,通过迭代更新获得局部更高精度的代理模型。工程实际算例表明,该方法对混凝土坝等大体积结构参数反演精度较高和适用性好,且能大幅减少传统有限元模型反演方法所需消耗的正分析机时,提高反演效率。     

粒子群算法优化的捷联罗经初始对准方法

初始对准是惯导系统的关键技术,罗经法对准是实现捷联惯导系统初始对准的重要手段。罗经对准回路的参数选择直接影响对准结果的好坏。对于不同的捷联惯导系统,罗经回路的最优参数也是不同的。传统的方法是根据经验以及大量的反复试验确定罗经对准参数,不能保证对准参数为最优。针对此问题,提出以水平罗经对准回路阻尼振荡周期T_(d1)和航向罗经对准回路阻尼振荡周期T_(d2)为寻优目标,用粒子群算法对参数(T_(d1),T_(d2))进行寻优的方法,以确定出满足条件的最优对准参数,从而提高捷联罗经初始对准的性能。实验结果表明:粒子群算法能够快速、准确地搜索出罗经对准回路的最优参数,提高捷联罗经对准的性能。将粒子群算法应用到捷联罗经初始对准中是有效的。     

基于荷载形函数的大跨桥梁结构移动荷载识别

大跨桥梁结构中的移动荷载识别是结构健康监测的重要组成部分之一,其作为动力学中的反问题,存在唯一性及稳定性等病态问题。本文首先推导出移动荷载作用下结构响应的卷积形式的离散公式,然后利用有限元理论中的形函数拟合移动荷载,推导出基于荷载形函数的移动荷载识别表达式。将移动荷载的识别转化为荷载形函数拟合系数的识别,降低了需要识别的未知量,减小了逆运算的计算量,消除或减弱逆问题病态性,并提高了计算效率。利用某大跨刚构-连续预应力混凝土桥梁修正后的有限元模型进行数值仿真,考虑路面粗糙度,由模态叠加法计算移动荷载作用下的响应;然后采用荷载形函数方法识别移动荷载,证明该方法在5%以下的噪声时能快速并精确识别复杂桥梁结构的移动荷载。     

A piezoelectric-elastic body with inhomogeneities and a crack under plane electrical and anti-plane mechanical loads

Summary  In this paper, we study a two-dimensional electroelastic problem of an infinite piezoelectric body with two circular piezoelectric inhomogeneities, one of which contains a crack. We formulate the stress intensity factor (SIF) analytically and investigate it numerically. The problem is solved based on Bueckner's principle, and is reduced to a problem of a singular integral equation of the first kind with respect to the distribution function of screw dislocation. The effect of interaction between the two inhomogeneities and the crack on the electroelastic field as well as the control of the SIF by electrical loads is investigated. Received 18 April 2000; accepted for publication 24 October 2000     

On the Path of a Quasi-static Crack in Mode III

A method for finding the path of a quasi-static crack growing in a brittle body is presented. The propagation process is modelled by a sequence of discrete steps optimizing the elastic energy released. An explicit relationship between the optimal growing direction and the parameters defining the local elastic field around the tip is obtained for an anti-plane field. This allows to describe a simple algorithm to compute the crack path. Mathematics Subject Classifications (2000) 74R05, 74B05, 74G70.Gerardo E. Oleaga: Supported by EU-Project Front Singularities University of Leipzig and the Max Planck Institute MIS. Partial support was also provided by the Spanish DGES project BFM2000-0605.     

裂纹体分析的权函数理论与应用: 回顾和展望

断裂力学是工程材料和结构的疲劳与断裂分析、损伤容限设计和结构完整性评定的理论基础. 应力强度因子作为线弹性裂纹尖端奇异场的单一表征参量和裂纹扩展驱动力, 在裂纹体的断裂力学分析中发挥着关键作用. 权函数法为复杂受载裂纹体的应力强度因子求解计算提供了强有力的解析工具, 不但具有远高于各类数值解法的计算效率, 而且精度可靠, 使用方便. 本文结合笔者团队在权函数法方面的长期研究工作, 对该方法自20世纪70年代初提出至今半个世纪以来, 国际断裂界在二维和三维权函数理论与应用方面的主要研究进展作了回顾和评述, 并对其未来发展提出了展望. 主要内容涵盖: 当前国际断裂界广泛应用的3种二维裂纹解析权函数法简介和以格林函数为基准的验证评价; 三维裂纹问题的片条合成权函数法和点载荷权函数法; 权函数法在复杂受载裂纹体的应力强度因子和裂纹张开位移等关键力学参量计算、内聚力/桥连等裂纹模型分析、共线多裂纹权函数理论及其在剩余强度预测等方面的应用, 以及复杂裂纹几何的工程化权函数分析和权函数法的反向应用问题.      

Assessing the effect of residual stresses on the fatigue behavior of integrally stiffened structures

Residual stresses emerge quite often in real structures due to the various manufacturing processes such as, welding, forming, cutting, milling, etc. In such cases, development of cracks at regions influenced by manufacturing operations demand additional attention. In the present work a numerical methodology has been developed, based on three-dimensional Finite Element Analysis, for the calculation of Stress Intensity Factors at cracks in welded components. The residual stress fields, which are used in SIF calculations, have been computed by the numerical simulation of the thermo-mechanical process. A numerical algorithm based on interpolation principles is developed, in order to introduce the three-dimensional field in the computational model of the cracked structure. The SIF calculation methodology is initially validated for the case of a welded plate by comparison of numerical results with existing analytical solutions. A cracked stiffened panel is analysed afterwards and the calculated fatigue crack propagation results are compared to experimentally measured data. Finally, the numerical procedure is applied to study the effect of more complicated residual stress fields on SIF values developing at cracks located in stiffened panels.     

表面裂纹问题瞬态响应的边界元分析

本文将拉氏变换-边界元法用于表面裂纹问题的瞬态响应分析。文中讨论了拉氏反演参数的选择和动态应力强度因子的计算方法。作为程序的考核和离散方案的选择,分别地计算了水平柱体一端固定、另一端受p(t)=poH(t)载荷时的位移响应和具有贯穿裂纹的厚板在两种离散方案时的动态应力强度因子响应。最后,还计算了若干载荷工况的半圆表面裂纹板应力强度因子的瞬态响应,获得了有效的数值结果。     

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.     

粒子群优化算法在传递对准中的应用

给出了一种基于粒子群优化算法的捷联惯导传递对准算法。简单分析了传递对准任务要求和主子惯导惯性器件输出之间的关系,将传递对准问题作为参数优化问题进行求解,给出了基于粒子群优化算法进行传递对准的数学模型。定义了传递对准的优化目标函数,介绍了粒子群优化算法及其应用于传递对准的具体算法设置。用粒子群优化算法求解目标函数的最小值,可获得主子惯导之间的失准角,进行一次校正即可完成传递对准过程。通过计算机仿真对算法进行了验证分析,在仿真条件下（陀螺精度为0.1°/h）,能达到方位0.1°的精度。与其他对准算法一样,算法受载体机动条件的影响较大,一般需要姿态机动来提高陀螺的信噪比。     

STRESS INTENSITY FACTOR OF AN ANTI-PLANE CRACK PARALLEL TO THE WEAK/MICRO-DISCONTINUOUS INTERFACE IN A BI-FGM COMPOSITE

The problem considered is a mode Ⅲ crack lying parallel to the interface of an exponential-type functional graded material （FGM） strip bonded to a linear-type FGM substrate with infinite thickness. By applying the Fourier integral transform, the problem was reduced as a Cauchy singular integral equation with an unknown dislocation density function. The collocation method based on Chebyshev polynomials proposed by Erdogan and Gupta was used to solve the singular integral equation numerically. With the numerical solution, the effects of the geometrical and physical parameters on the stress intensity factor （SIF） were analyzed and the following conclusions were drawn： （a） The region affected by the interface or free surface varies with the material rigidity, and higher material rigidity will lead to bigger affected region. （b） The SIF of the crack in the affected region and parallel to the micro-discontinuous interface is lower than those of the weak discontinuous cases. Reducing the weak-discontinuity of the interface will be beneficial to decrease the SIF of the interface-parallel crack in the region affected by the interface. （c） The effect of the free surface on SIF is more remarkable than that of the interface, and the latter is still more notable than that of the material rigidity. When the effects of the interface and free surface are fixed, increase of the material rigidity will enhance the value of SIF.     

Magnetostrictive/electrostrictive fracture of the piezomagnetic and piezoelectric layers in a multiferroic composite: Anti-plane case

The main purpose of the present work is to study the influences of magnetostriction, electrostriction and piezomagnetic/piezoelectric stiffening on the fracture behavior of a layered multiferroic composite. For comparison, it is assumed that there is a crack, parallel to the interface, in each layer. Methods of cosine transform and Cauchy singular integral equations are used to solve the crack problem. Numerical results of the stress intensity factor (SIF) are provided and the computational accuracy is demonstrated. Discussion on the numerical results indicates that the multiferroic composite consisting of cobalt ferrite and barium titanate layers are more prone to fracture under electric loading than under magnetic loading. In the case of magnetostriction, to increase the shear modulus of the piezomagnetic layer would raise the SIF; but to increase that of the piezoelectric layer would reduce the SIF; in the case of electrostriction, inverse results are obtained. Piezomagnetic stiffening can affect the SIF when the composite is under electrostriction; piezoelectric stiffening can influence the SIF if the composite is under magnetostriction. In addition, it is also revealed that two parallel equal cracks may shield each other even if an interface exists between them.     

The Stress State of a Ferromagnetic with an Elliptic Crack in a Homogeneous Magnetic Field

The problem on the stress–strain state of an infinite isotropic body made of a magnetically soft material and containing an elliptic crack is considered. It is assumed that the body is under an external magnetic field perpendicular to the crack plane. The basic characteristics of the stress–strain state and the magnetic field induced are determined and their singularities near the elliptic crack are studied. Formulas are given for the stress intensity factors for the force and magnetic fields near the crack tip     

A mode III crack crossing the magnetoelectroelastic bimaterial interface under concentrated magnetoelectromechanical loads

A mode III crack cutting perpendicularly across the interface between two dissimilar semi-infinite magnetoelectroelastic solid is studied under the combined loads of a line force, a line electric charge and a line magnetic charge at an arbitrary location. The impermeable conditions are implied on the crack faces. The technique developed in literature for the elastic bimaterial with a crack cutting interface is exploited to treat the magnetoelectroelastic bimaterial. The Riemann-Hilbert problem can be formulated and solved based on complex variable method. Analytical solutions can be obtained for the entire plane. The intensity factors around crack tips can be defined for the elastic, electric and magnetic fields. It shows that, no matter where the load position is, the electric displacement intensity factors (EDIFs), as well as the magnetic induction intensity factors (MIIFs), are identical in magnitude but opposite in sign for both crack tips, on condition that a line force is solely applied. Alternatively, if only a line electric charge is considered, then the stress intensity factors (SIFs) and the MIIFs exhibit the behavior. Likewise, if only a line magnetic charge is applied, it turns to the SIFs and the EDIFs instead. In addition, the dependence of the intensity factors is graphically shown with respect to the location of a line force. It is found that the SIF for a crack tip tends to be infinite if the applied force is approaching the tip itself, but the EDIF, with the complete opposite trend, tends to be vanishing. Finally, focusing on the more practical case of piezoelectric/piezomagnetic bimaterial, variation of the SIF along with the moduli as well as the piezo constitutive coefficients is explored. These analyses may provide some guidance for material selection by minimizing the SIF. It is also believed that the results obtained in this paper can serve as the Green’s function for the dissimilar magnetoelectroelastic semi-infinite bimaterial with a crack cutting the interface under general magnetoelectromechanical loads.