Stress Separation of Interferometrically Measured Isopachics in a Perforated Plate

A method for the stress separation of interferometrically measured isopachics using an Airy stress function is proposed in this study. A Poisson equation that represents the relationship between the sum of principal stresses and an Airy stress function is solved using a finite element method. The Dirichlet boundary condition for solving the Poisson equation is determined by the approximation of an assumed Airy stress function along the boundary of the model. Therefore, the distribution of the Airy stress function is obtained from the measured isopachic contours. Then, the stresses are obtained from the computed Airy stress function. The effectiveness of the proposed method is validated by applying the proposed method to the isopachic contours in a perforated plate obtained by Mach-Zehnder interferometry. Results indicate that stress components around a hole in a plate can be obtained from isopachics by the proposed method.     

刚构连续梁桥墩梁固节点结构的光弹性试验及有限元分析

某铁路特大桥的主桥部分采用预应力混凝土刚构连续梁桥,本文采用光弹性模型冻结应力法对该刚构连续梁桥墩梁固节点结构进行了应力分析,给出了结构边界应力分布和主应力迹线。应用ANSYS软件对墩梁固节点实际结构进行了有限元计算。将实验与有限元分析结果进行了比较,结果表明,实验应力分析法与有限元数值法吻合较好。在有限元分析中,对墩梁固节点结构梁端加载边界的影响区进行了研究,通过合理选择梁的长度可以有效减小局部应力。模型实验和有限元计算的结果为结构的优化设计提供了重要的依据。     

振荡应力奇异性及其强度系数的数值分析方法

本文以具有振荡应力奇异性的平面问题为例,提出了一种利用普通的数值分析结果（由有限地或边界元程序计算得到的应力分量或位移分量）,来确定奇异性次数以及相应的复应力强度系数的数值分析方法。为了验证该方法的有效性,应用平面应变情况下的边界元计算结果,对界面端模型进行了分析。计算结果表明,本方法可以精确地求得振动应力奇异性次数,并且与奇异性对应的复应力强度系数也可以很方便地应用外插法得到。     

边界元法与光弹性实验相结合的应力分析法

本文介绍了边界无法与光弹性实验相结合进行应力分析的方法。叙述了用边界无法求解主应力的基本原理。结合光弹性实验等差线条纹图,可分离出主应力。     

CRACK FACE CONTACT PROBLEM ANALYSIS USING THE SCALED BOUNDARY FINITE ELEMENT METHOD

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

封隔器卡瓦接触应力研究

封隔器上的卡瓦锚定后起到支撑封隔器、锁定胶皮筒的作用,其性能好坏直接影响到油井的产量和生产安全。卡瓦与套管间的咬合力分布决定了封隔器在使用中的成败。咬合力的分布是否合理,将直接关系到套管的损伤程度及卡瓦的寿命。文中同时应用有限元法和三维光弹性技术对封隔器卡瓦进行了接触应力分析,介绍了实验模型的建立和实验步骤以及计算模型的建立和边界条件的考虑,比较了数值计算与三维光弹实验的结果,分析了造成误差的原因。实验应力分析结果为有限元计算模型的建立及边界条件的确定提供了重要的实验保证。研究工作获得了一些重要的结果,为封隔器卡瓦的优化设计提供了相关的参数。     

Load Extraction from Photoelastic Images Using Neural Networks

Photoelastic materials develop colored fringes under white light when subjected to mechanical stresses, which can be viewed through a polariscope. This technique has traditionally been used for stress analysis of loaded components, however, this can also be potentially used in sensing applications where the requirement may be measurement of the stimulating forces causing the generation of fringes. This leads to inverse photoelastic problem where the developed image can be analyzed for the input forces. However, there could be infinite number of possible solutions which cannot be determined by conventional techniques. This paper presents neural networks based approach to solve this problem. Experiments conducted to prove the principle have been verified with theoretical results and finite element analysis of loaded specimens. The developed technique, if generalized, can be implemented for whole-field analysis of the stress patterns involving complex fringes under different loading conditions. This can also provide direct visualization of the stress field, which may find application in a variety of specialized areas including biomedical engineering and robotics.

Load and Boundary Condition Calibration Using Full-field Strain Measurement

For critical load bearing structures, it is often necessary to experimentally determine the load distribution on the structure so that accurate finite element models can be developed for stress and fatigue life predictions. An inverse problem approach is presented here for computing or calibrating the loads and boundary conditions acting on a structure. This enables the creation of more accurate finite element models, especially for structures that have complicated load distribution and compliant boundary conditions. The method presented here involves minimizing the least square error between the strains computed using the finite element model and the strains and displacements obtained experimentally. The nodal loads and the compliance at fixed boundaries are treated as the variables in the optimization problem. The compliance is modeled as springs attached at the nodes that are on the boundary where the structure is restrained. The method is verified by computing the loads and boundary conditions when displacements, maximum shear strain or both are available at large number of points on the surface of the structure. The experimental data set was generated using the luminescent photoelastic coating (LPC) technique.     

局部有限元模型边界条件的直接传递方法

提出一种将整体分析得到的节点力或节点位移直接传递到精细化局部有限元模型的方法,即部分混合单元法。沿精细化局部有限元模型周边建立一组过渡单元,该组过渡单元采用与整体模型一致的单元类型和模拟方式,其外侧边界上的节点与整体模型节点的相对坐标对应,内侧边界与精细化局部有限元模型采用基于面约束的方式连接。在外侧边界上根据节点坐标对应施加整体分析获得的节点力或节点位移,过渡单元就可直接将边界条件传递到精细化局部有限元模型。通过贵州红水河特大桥钢-混结合段的精细化有限元分析,验证了本文方法的实用性和有效性。     

V形切口应力强度因子的一种边界元分析方法

将V形切口结构分成围绕切口尖端的小扇形和剩余结构两部分. 尖端处扇形域应力场表示成关于尖端距离$\rho$的渐近级数展开式,从线弹性理论方程推导出了一组分析平面V形切口奇异性的常微分方程特征值问题,通过求解特征方程,得到前若干个奇性指数和相应的特征向量. 再将切口尖端的位移和应力表示为有限个奇性阶和特征向量的组合. 然后用边界元法分析挖去小扇形后的剩余结构. 将位移和应力的线性组合与边界积分方程联立,求解获得切口根部区域的应力场、应力幅值系数和整体结构的位移和应力. 从而准确计算出平面V形切口的奇异应力场和应力强度因子.      

Separation of principal stresses in dynamic photoelasticity

A new and effective method used to separate the transient principal stresses for dynamic photoelasticity is proposed. This is a hybrid method combining the optical method of dynamic caustics and the boundary element numerical method. Firstly, a modified Cranz-Schardin spark camera is used to record simultaneously the isochromatic fringe patterns of photoelasticity and the shadow spot patterns in the dynamic process. By means of the isochromatic fringe patterns, the difference between transient principal stresses in the whole domain and the principal stresses along the free boundary can be solved. In addition, the method of caustics is a very powerful technique for measuring the concentrative load. Then, the sum of the principal stresses is calculated by the boundary integral equation obtained from the Laplace integral transform of the wave equation. So, the transient principal stresses can be determined from the experimental and numerical results. As an example, the transient principal stresses in a polycarbonate disk under an impact load are resolved. Concurrently published in the Chinese Edition of Acta Mechanica Sinica, Vol. 26, No. 1, 1994     

Extended Kantorovich method for local stresses in composite laminates upon polynomial stress functions

The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress functions. The stress fields are initially assumed by means of the Lekhnitskii stress functions under the plane strain state. Applying the principle of complementary virtual work, the coupled ordinary differential equations are obtained in which the solutions can be obtained by solving a generalized eigenvalue problem. Then an iterative procedure is estab-lished to achieve convergent stress distributions. It should be noted that the stress function based extended Kantorovich method can satisfy both the traction-free and free edge stress boundary conditions during the iterative processes. The stress components near the free edges and in the interior regions are calculated and compared with those obtained results by finite element method (FEM). The convergent stresses have good agreements with those results obtained by three dimensional (3D) FEM. For generality, various layup configurations are considered for the numerical analysis. The results show that the proposed polynomial stress function based extended Kan-torovich method is accurate and efficient in predicting the local stresses in composite laminates and computationally much more efficient than the 3D FEM.     

Contact stress analysis of an elastic half-plane containing multiple inclusions

This paper presents a two-dimensional contact stress analysis to investigate the effects of multiple inclusions on the contact pressure and subsurface stresses in an elastic half-plane. The boundary element method is used to analyze the contact problem where a set of integral equations is derived on the contact region and the matrix–inclusion interfaces. As the contact region is unknown a priori, an iterative procedure is implemented to determine the actual contact region and the contact pressure, and the tractions and displacements on the matrix–inclusion interfaces are obtained by solving the integral equations numerically. Numerical results show that the inclusions near contact surface could cause significant alterations in the contact pressure distribution. The stiff inclusions could toughen the surrounding material and reduce the internal stresses while the soft inclusions could increase the subsurface stresses.     

Determining Individual Stresses Throughout a Pinned Aluminum Joint by Reflective Photoelasticity

Pinned (bolted) joints are an extremely important, but difficult to analyze, structural or mechanical element. They are a class of inverse problems in which the stresses at the pin/hole interface are typically unknown. Moreover, stresses vary non-linearly with applied load. Failures of mechanical or structural systems frequently initiate at connections. Although almost always present, many stress analyses of such mechanical connections ignore friction for simplicity. The stresses are evaluated here in an aluminum connector using a series solution of an Airy stress function, the coefficients being evaluated from known boundary tractions (near, but not including the contact region on the hole) and photoelastically measured data obtained from a bonded birefringent coating. Both friction and pin/hole clearance are accounted for, and individual stresses are evaluated full-field, including on the contact boundary of the hole.     

Solution of the connection problems between a finite holed plate and a stiffener by using the partitioning concept of the generalized variational method

In this paper a new finite element method is presented, in which complex functions are chosen to be the finite element model and the partitioning concept of the generalized variational method is utilized. The stress concentration factors for a finite holed plate welded by a stiffener are calculated and the analytical solutions in series form are obtained. From some computer trials it is demonstrated that the problem of displacement compatibility and continuity of tractions between the holed plate and the stiffener is successfully analysed by using this method. Since only three elements need to be formulated, relatively less storage is required than the usual finite element methods. Furthermore, the accuracy of solutions is improved and the computer time requirements are considerably reduced. Numerical results of stress concentration factors and stresses along the welded-line which may be referential to engineers are shown in tables.     

Integrated-plane photoelastic method—Application of photoelastic isodynes

It has been shown by Pindera and Mazurkiewicz that a new type of scattered-light modulation in the plane of a two-dimensional photoelastic object can be obtained when the stationary integrated photoelastic method developed by Pindera and Straka is applied in a scanning mode and when the transfer function of the photoelastic system satisfies certain conditions. The new type of light modulation, called field of isodynes by the authors, carries information on stress components normal to the direction of propagation of primary beam, and on corresponding total-force component. The points where this stress component is equal to zero can be easily determined. The classical scattered-light modulation along a chosen line represents a cross section of a corresponding isodynes field. It is shown that these features of the method of isodynes make it possible to easily determine the distribution and values of normal stress components at any arbitrary rectilinear cross section, and to check immediately the accuracy of measurements. The experimental determination of contact stresses and contact regions using the method of isodynes is especially simple and elegant.     

An Experimental-Numerical Hybrid Method for Elastic Contact Problems*

ABSTRACT

An experimental-numerical hybrid technique for determining the contact stress distribution between two elastic bodies having both frictionless as well as bonded contact is discussed in this paper. The hybrid method makes use of experimental data collected at a section far from the contact surface and the numerically generated influence coefficients, in terms of the applied unit normal and shear stresses. The experimental data, i.e., the differences in normal stresses and the shear stress, are obtained using photoelastic analysis for the examples illustrated in this paper. When substituted into equations corresponding to the unit normal and shear stress applied in the contact region, this results in a set of algebraic equations which, when solved, allow the contact stress distribution to be obtained. This method is illustrated with examples involving simple and complex geometries of the contacting bodies.     

Photoelastic model measurement with rotated principal axes by scattered-light photoelasticity

The method of measurement of a three-dimensional photoelastic model with rotated principal axes has not yet been fully experimentally established. It is known that a three-dimensional photoelastic model can be reduced to an optically equivalent model. In this paper, the optically equivalent model is realized from a stratified model consisting of two layers of the frozen stress model. The secondary principal stress direction and the relative phase retardation of the frozen disk model in the stratified frozen stress model are determined for the entire field from Stokes parameters obtained by scattered-light photoelasticity using unpolarized light. The accuracy of these values is confirmed by a comparison with results only from the frozen stress disk model.     

A HYBRID METHOD FOR ANALYZING THE DYNAMIC RESPONSES OF CAVITIES OR SHELLS BURIED IN AN ELASTIC HALF-PLANE

In this paper, based on a variational formalism which originally proposed by Mei [1] for infinite elastic medium and extended by Yeh, et al. [2,3] for elastic half-plane, a hybrid method which combines the finite element and series expansion method is implemented to solve the diffraction of plane waves by a cavity buried in an elastic half-plane. The finite domain which encloses all inhomogeneities including the cavity can be easily formulated by finite element methods. The unknown boundary data obtained by subtracting the known free fields from the total fields which include the boundary nodal displacements and tractions at the interface between the finite domain and the surrounding elastic half-plane are not independent of each other and can be correlated through aseries repre sentation. Due to the continuity condition at the interface, the same series representation is still valid for the exterior elastic half-plane to represents the scattered wave. The unknown coefficients of this series are treated as generalized coordinates and can be easily formulated by the same variational principle. The expansion function of the series is composed of basis function. Each basis function is constructed from the basis function for an infinite plane by superimposing an additional homogeneous reflective term to satisfy both traction free conditions at ground surface and radiation conditions at infinity. The numerical results are made against those obtained by boundary element methods, and good agreements are found.     

Improving fringe definition by pseudo-solarization

A method is described for obtaining well-defined fringes as well as fractional fringes in photoelastic stress analysis. It is based on an edge effect which occurs during the pseudo-solarization of films during development. Sharp fractional fringes are obtained throughout the field of the model with a high degree of accuracy. The technique can be used with all types of polariscopes.