Linear and Nonlinear Algorithms for Stress Separation in Photoelasticity

An experimental-numerical hybrid method for the stress separation in photoelasticity is proposed in this study. In the proposed method, boundary conditions for a local finite element model, that is, tractions along boundaries are inversely determined from photoelastic fringes. Two algorithms are proposed for determining the boundary condition. One is a linear algorithm in which the tractions are obtained by the method of linear least-squares from both principal stress difference and principal direction. Another is the nonlinear algorithm in which the tractions are determined only from the principal stress difference. After determining the boundary conditions for the local finite element model, the stresses can be obtained by finite element direct analysis. The effectiveness is demonstrated by applying the proposed method to a perforated plate under tension and contact problems. Results show that the boundary conditions of the local finite element model can be determined from the photoelastic fringes and then the individual stresses can be obtained by the proposed method. Furthermore, the stresses can be evaluated even if the boundary condition is complicated such as at the contact surface. It is expected that the proposed method can be powerful tool for stress analysis.     

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.     

The general stress strain relation of soils involving the rotation of principal stress axes

1.IntroductionThestudyofsoils'constitutiverelationismostlylimitedinsuchStressstatethattherotationofprincipalstressaxesisabletobeignored.Soonlythevalueofprincipalstressneedtobeconsideredandthedirectionofprincipal'stressisregardedtobeunchangeable.Onthiscondition,theprincipalaxesof'stressincrement,strainincrementjstressandstrainarethesame.Therelevantexperimentsandengineeringsshowthattherotationofprincipalstressaxeswillgeneratesignificantplasticdeformationandthenoncoaxalityofstressandstrain.First…     

Phase Shifting Full-Field Interferometric Methods for Determination of In-Plane Tensorial Stress

A new method that combines phase shifting photoelasticity and transmission Coherent Gradient Sensing (CGS) is developed to determine the tensorial stress field in thin plates of photoelastic materials. A six step phase shifting photoelasticity method determines principal stress directions and the difference of principal stresses. The transmission CGS method utilizes a standard four step phase shifting method to measure the x and y first derivatives of the sum of principal stresses. These stress derivatives are numerically integrated using a weighted preconditioned conjugate gradient (PCG) algorithm, which is also used for the phase unwrapping of the photoelastic and CGS phases. With full-field measurement of the sum and difference of principal stresses, the principal stresses may be separated, followed by the Cartesian and polar coordinate stresses using the principal stress directions. The method is demonstrated for a compressed polycarbonate plate with a side V-shaped notch. The experimental stress fields compare well with theoretical stress fields derived from Williams solution for a thin plate with an angular corner.     

基于统一强度理论的非静水压圆形隧道塑性区半径的脆塑性摄动解

基于统一强度理论和弹?脆?塑性模型, 综合考虑围岩强度的中间主应力效应和脆性软化, 采用摄动法建立了非静水压圆形隧道塑性区半径的解析解, 继而探讨所得解析解的适用范围, 并与文献复变函数法、摄动法、数值模拟和总荷载不变法进行对比, 最后分析各因素对隧道塑性区形状和大小的影响规律. 研究结果表明: 所建立的圆形隧道塑性区半径摄动解析解为反映中间主应力效应不同程度的系列解答, 可退化为理想弹塑性模型解答且得到文献中4种方法的正确性和合理性验证, 适用于塑性区完全包围隧道的情况, 对应的隧道塑性区为双轴对称的类椭圆, 具有广泛的理论意义和工程应用价值; 摄动参数对隧道类椭圆形塑性区的大小和长/短轴变化都有明显影响; 隧道塑性区范围随中间主应力效应、围岩峰后强度参数的增加均显著减小, 说明不考虑中间主应力效应的Mohr?Coulomb强度准则解答偏保守, 弹?脆?塑性模型相比理想弹塑性模型更适合隧道塑性分析.      

A STRESS VECTOR-BASED CONSTITUTIVE MODEL FOR COHESIONLESS SOIL (Ⅱ)-APPLICATION

ＩｎｔｒｏｄｕｃｔｉｏｎＡｃｏｍｐａｎｉｏｎｐａｐｅｒ [1 ]ｄｅｓｃｒｉｂｅｄａｓｔｒｅｓｓｖｅｃｔｏｒ＿ｂａｓｅｄｃｏｎｓｔｉｔｕｔｉｖｅｍｏｄｅｌｆｏｒｃｏｈｅｓｉｏｎｌｅｓｓｓｏｉｌ.Ｔｈｅｍｏｄｅｌｉｓｂａｓｅｄｏｎｔｈｅｅｑｕｉｖａｌｅｎｔｄｅｃｏｍｐｏｓｉｔｉｏｎｏｆｔｈｅａｃｔｉｏｎｅｆｆｅｃｔｏｆｓｔｒｅｓｓｖｅｃｔｏｒｆｏｒｔｈｅｍａｉｎ .Ｔｈａｔｉｓ ,ｕｎｄｅｒｐｌａｎｅｓｔｒａｉｎｃｏｎｄｉｔｉｏｎ ,ｔｈｅａｃｔｉｏｎｅｆｆｅｃｔｏｆｓｔｒｅｓｓｖｅｃｔｏｒｃａｎｂｅｅｑｕｉ…     

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.     

Stress-strain state of an inclined elliptical defect in a plate under biaxial loading

A mathematical model for the stress-strain state of a plate with an inclined elliptical defect under biaxial loading is considered. Exact formulas for stresses in polar coordinates, displacements, principal stresses, maximum shear stress, and stress intensity in the case of a plane stress state of the plate were obtained by the Kolosov-Muskhelishvili method. Simulation results are compared with experimental data obtained by holographic interferometry.     

含主应力轴旋转的土体本构关系研究进展

主应力轴旋转对岩土工程的影响日益受到人们的重视．本文从纯主应力轴旋转、纯应力洛德角变化、多种因素变化、排水与不排水等方面较全面地描述了含主应力轴旋转情况下土体的基本变形特性．对当前较有影响的含主应力轴旋转的土体本构模型（基于一般应力空间的土体本构模型、运动硬化模型、边界面模型、土体弹塑性应力应变关系的完全应力增量表述等等）进行了较为系统的评述．提出了合主应力轴旋转的实验研究中存在的核心问题与建立合主应力轴旋转的土体本构关系的根本途径．      

Data analysis in dynamic photoelasticity

A method to separate the principal stresses and to obtain displacements from dynamic photoelastic fringe patterns is presented. The method utilizes geometric characteristics of either the model or the stress wave to provide the additional equation necessary for separation. The method is illustrated in ten specific cases.     

A STRESS VECTOR-BASED CONSTITUTIVE MODEL FOR COHESIONLESS SOIL (Ⅱ)-APPLICATION

The stress vector-based constitutive model for cohesionless soil, proposed by SHI Hong-yan et al., was applied to analyze the deformation behaviors of materials subjected to various stress paths. The result of analysis shows that the constitutive model can capture well the main deformation behavior of cohesionless soil, such as stress-strain nonlinearity, hardening property, dilatancy, stress path dependency, non-coaxiality between the principal stress and the principal strain increment directions, and the coupling of mean effective and deviatoric stress with deformation. In addition, the model can also take into account the rotation of principal stress axes and the influence of intermediate principal stress on deformation and strength of soil simultaneously. The excellent agreement between the predicted and measured behavior indicates the comprehensive applicability of the model. Biography: SHI Hong-yan (1962-), Associate Professor, Doctor     

高地应力岩体多孔爆破破岩机制

深部岩体爆破破岩是爆炸荷载与高地应力共同作用的结果。基于一些简化假设，建立了一个高地应力岩体双孔爆破计算模型，采用光滑粒子流体力学-有限元方法耦合数值模拟方法，研究了高地应力作用下炮孔间裂纹的传播及贯通过程，分析了炮孔周围应力场动态演化过程与分布特征。研究结果表明:爆破引起的岩体开裂主要是环向动拉应力所致，地应力对岩体的压缩降低了炮孔周围环向动拉应力、缩短了环向动拉应力的作用时间，因而对爆炸致裂起抑制作用；静水地应力条件下多孔爆破时，垂直于炮孔连线方向传播的爆生裂纹更易受到地应力的抑制；对于高地应力岩体爆破，炮孔间的裂纹扩展长度随地应力水平的提高而减小，裂纹主要沿最大主应力方向扩展，因此沿最大主应力方向布置炮孔、缩短炮孔间距有利于炮孔间裂纹的连接贯通，形成良好的爆破开挖面。     

斜拉桥塔锚固区光弹性应力分析与等倾线图像处理

本文采用光弹性冻结应力法，对独塔斜拉桥塔锚固区进行了应力分析，给出了有关截面的边界应力分布及主应力迹线，从而为优化工程设计提供了重要的参考依据。本文还提出了一种新的提取等倾线的方法──图像对数增强相减法，实现了等倾线与等差线条纹的分离，同时使条纹得到了细化，提高了处理等倾线的速度与准确度。编制了自动绘制主应力迹线的程序，利用它绘制了桥塔截面的次主应力迹线。     

页岩储层体积压裂的地应力变化研究

页岩储层属于致密超低渗透储层,需改造形成复杂缝网才有经济产能.体积压裂是页岩储层增产改造的主要措施,而地应力场特别是水平主应力差值是体积压裂的关键控制因素. 理论研究表明:(1)当初始两向水平主应力差较小时,容易形成缝网,反之不易产生缝网;(2)人工裂缝的形成能够改变地层初始应力场. 因此应在前人研究的基础上优化设计压裂方式,以克服和翻转初始水平主应力差值,产生体积缝网.基于此,建立了页岩气藏水平井体积压裂数值模型,模型中采用多孔介质流固耦合单元模拟页岩基质的行为,采用带有孔压的"cohesive"单元描述水力裂缝的性质,模型对"Texas Two-Step" 压裂方法进行了数值模拟,模拟结果得到了压裂过程中地层应力场的分布及其变化,模拟结果和解析公式计算结果吻合良好.模拟结果表明:(1)裂缝的产生减弱了地层应力场的各向异性;(2 对于低水平应力差页岩储层,采用"Texas Two-Step"压裂方法可以产生缝网. 对于采用"Texas Two-Step"压裂方法无法产生缝网的高应力差页岩储层,提出了三次应力"共振" 和四次应力"共振" 压裂方法并进行了数值模拟,模拟结果得到了压裂过程中页岩储层应力场的分布及其变化,得到了缝网形成的区域,模拟结果表明:(1)对于高应力差页岩储层,采用"Texas Two-Step" 压裂方法无法产生缝网;(2)对于高应力差页岩储层,三次应力"共振" 和四次应力"共振"压裂方法是有效的体积压裂缝网形成的方法.      

主应力迹线的自动绘制

主应力迹线是代表主应力方向的曲线簇。在土木建筑工程中主应力迹线是布筋方向的依据,但是在工程中通常都是徒手或用曲线板近似绘制主应力迹线。作者提出一种新的主应力迹线自动的绘制方法。该方法是根据主应力迹线满足过任意三条等倾线的主应力都可以用二次曲线来描述这一条件,再利用主应力迹线本身所具有的性质列出方程组解算出描述主应力迹线的一元二次方程,然后根据方程描绘出所求的主应力迹线。实验证明这种方法对于工程上进行主应力迹线的绘制很有效。     

主应力轴旋转下黏土累积变形的分数阶元件模型

在长期交通载荷作用下土体塑性累积变形本构模型对路基沉降计算至为关键.元件组合模型可以计算岩土体循环累积应变,但现有的各类元件模型未能反映饱和软黏土的主应力轴循环旋转现象.在对饱和软黏土进行等向固结条件下的主应力轴循环旋转加载试验及非等向固结下的循环扭剪试验基础上,将Abel黏壶代替Burgers模型中的Newton黏壶,得到分数阶Burgers模型;利用遗传算法优化循环塑性累积应变的Burgers模型和分数阶Burgers模型的参数,通过对比两组模型的计算值与试验值,发现分数阶模型更适合模拟计算循环载荷下饱和软黏土的累积变形.     

Photoplastic stress analysis considering rate effect

The possibilities of cellulose acetate as a photoplastic model material used in the uniaxial tension test were presented by the authors in Ref. 1. Empirical formulas were proposed for photoplastic stress analysis considering the effect of strain rate and temperature. In this paper, the stress-strain-optic laws derived from the uniaxial tension test were confirmed precisely under the biaxial stress field. Consequently, the biaxial stress-strain-optic laws of cellulose acetate can be represented by simply substituting the difference of the principal stress and the difference of the principal strain into the terms of the stress and the strain in the empirical formulas for the uniaxial stress state, respectively. Therefore, the empirical formulas considering the influence of strain rate were obtained for the photo-viscoelastoplastic plane-stress analysis during the usual static loading. This method was applied to an elastoplastic problem of a finite plate with a circular hole under uniaxial tension. The distribution of stress and strain, the stress-concentration factors in the elastoplastic region and the development of the plastic region were obtained experimentally taking the rate effect into account. As a result, the stress on a minimum section in a model increased as head speed increased. On the other hand, the distribution of strain, the development of the plastic region and the stress-concentration factors were almost never influenced by head speed.     

饱和砂土的三剪弹塑性边界面模型研究(一)——模型理论

基于边界面模型方法研究了不同密实度和含泥量的饱和砂土在静、动荷载作用下的力学和变形特性。其中,初始边界面方程是通过改进修正剑桥模型椭圆形屈服方程得到的,后继边界面方程则是在前次边界面方程的基础上采用中心点映射法获得。针对原模型破坏应力比为定值,不能反映全应力状态和黏聚力在土体受力过程中的作用这一问题,基于三剪统一强度准则,通过等量代换法、坐标平移法对破坏应力比、相变应力比进行了修正,使修正后的椭圆形边界面适用于不同密实度和含泥量的饱和砂土。在此基础上应用边界面理论,研究了饱和砂土的三剪弹塑性边界面本构模型及有限元算式。该模型能够描述密实砂土在受力时表现的剪胀特性,同时可以反映中间主应力效应、应力区间效应、拉压差效应和黏聚力对应力-应变特性的影响。     

Stress-strain state and stress intensity factors of an inclined elliptic defect in a plate under biaxial loading

We consider a mathematical model of the stress-strain state of a plate with an inclined elliptic defect. We obtain approximate formulas for the stress tensor, the displacements, and the principal stresses near the defect vertex. The obtained formulas are compared with the results obtained by the holographic photoelasticity method.     

NONLINEAR MICRO－MECHANICAL MODEL FOR PLAIN WOVEN FABRIC

The warp yarns and weft yarns of plain woven fabric which, being the principal axes of material of fabric, are orthogonal in the original configuration, but are obliquely crossed in the deformed configuration in general. The orthotropic constitutive model is unsuitable for fabric. In the oblique principal axes system the relations between loaded stress vectors and stress tensor are investigated, the stress fields of micro-weaving structures of fabric due to pure shear are carefully studied and, finally, a nonlinear micro-mechanical model for plain woven fabric is proposed. This model can accurately describe the nonlinear mechanical behavior of fabric observed in experiments. Under the assumption of small deformation and linearity of mechanical properties of fabric the model will degenerate into the existing linear model.