Stress tensor determination by modified hydraulic tests on pre-existing fractures:Method and stress constraints

The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determination based on the HTPF method requires at least six tests or a minimum of 14-15 tests(under different conditions)for reliable results.In this study,we modified the HTPF method by considering the shear stress on each pre-existing fracture,which increased the number of equations for the stress tensor determination and decreased the number of tests required.Different shear stresses were attributed to different fractures by random sampling;therefore,the stress tensors were obtained by searching for the optimal solution using the least squares criterion based on the Monte Carlo method.Thereafter,we constrained the stress tensor based on the tensile strength criterion,compressive strength criterion,and vertical stress constraints.The inverted stress tensors were presented and analyzed based on the tensorial nature of the stress using the Euclidean mean stress tensor.Two stress-measurement campaigns in Weifang(Shandong Province,China)and Mercantour road tunnel(France)were implemented to highlight the validity and efficiency of the modified HTPF(M-HTPF)method.The results showed that the M-HTPF method can be applied for stress tensor inversion using only three to four tests on pre-existing fractures,neglecting the stress gradient.The inversion results were confined to relatively small distribution dispersions and were significantly reliable and stable due to the shear stresses on the fractures and the stress constraints employed.The M-HTPF method is highly feasible and efficient for complete stress tensor determination in a single borehole.

Probabilistic-Ellipsoid Hybrid Reliability Multi-Material Topology Optimization Method Based on Stress Constraint

This paper proposes a multi-material topology optimization method based on the hybrid reliability of the probability-ellipsoid model with stress constraint for the stochastic uncertainty and epistemic uncertainty of mechanical loads in optimization design.The probabilistic model is combined with the ellipsoidal model to describe the uncertainty of mechanical loads.The topology optimization formula is combined with the ordered solid isotropic material with penalization(ordered-SIMP)multi-material interpolation model.The stresses of all elements are integrated into a global stress measurement that approximates the maximum stress using the normalized p-norm function.Furthermore,the sequential optimization and reliability assessment(SORA)is applied to transform the original uncertainty optimization problem into an equivalent deterministic topology optimization(DTO)problem.Stochastic response surface and sparse grid technique are combined with SORA to get accurate information on the most probable failure point(MPP).In each cycle,the equivalent topology optimization formula is updated according to the MPP information obtained in the previous cycle.The adjoint variable method is used for deriving the sensitivity of the stress constraint and the moving asymptote method(MMA)is used to update design variables.Finally,the validity and feasibility of the method are verified by the numerical example of L-shape beam design,T-shape structure design,steering knuckle,and 3D T-shaped beam.

APPROACH FOR LAYOUT OPTIMIZATION OF TRUSS STRUCTURES WITH DISCRETE VARIABLES UNDER DYNAMIC STRESS, DISPLACEMENT AND STABILITY CONSTRAINTS

A mathematical model was developed for layout optimization of truss structures with discrete variables subjected to dynamic stress, dynamic displacement and dynamic stability constraints. By using the quasi-static method, the mathematical model of structure optimization under dynamic stress, dynamic displacement and dynamic stability constraints were transformed into one subjected to static stress, displacement and stability constraints. The optimization procedures include two levels, i.e., the topology optimization and the shape optimization. In each level, the comprehensive algorithm was used and the relative difference quotients of two kinds of variables were used to search the optimum solution. A comparison between the optimum results of model with stability constraints and the optimum results of model without stability constraint was given. And that shows the stability constraints have a great effect on the optimum solutions.

A METHOD FOR TOPOLOGICAL OPTIMIZATION OF STRUCTURES WITH DISCRETE VARIABLES UNDER DYNAMIC STRESS AND DISPLACEMENT CONSTRAINTS

A method for topological optimization of structures with discrete variables subjected to dynamic stress and displacement constraints is presented. By using the quasistatic method, the structure optimization problem under dynamic stress and displacement constraints is converted into one subjected to static stress and displacement constraints. The comprehensive algorithm for topological optimization of structures with discrete variables is used to find the optimum solution.

THE TOPOLOGICAL OPTIMIZATION FOR TRUSS STRUCTURES WITH STRESS CONSTRAINTS BASED ON THE EXIST-NULL COMBINED MODEL

A new exist-null combined model is proposed for the structural topology optimization. The model is applied to the topology optimization of the truss with stress constraints. Satisfactory computational result can be obtained with more rapid and more stable convergence as compared with the cross-sectional optimization. This work also shows that the presence of independent and continuous topological variable motivates the research of structural topology optimization.

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅰ)--THE ANALYSES OF CONSTRAINT PARAMETERS AND FRACTURE PARAMETERS

In the present paper,three dimensional analyses of some general constraint parameters and fracture parameters near the crack tip of Mode I CT specimens in two different thicknesses are carried out by employing ADINA program.The results reveal that the constraints along the thickness direction are obviously separated into two parts:the keeping similar high constraint field(Z_(1))and rapid reducing constraints one(Z_(2)).The two fields are experimentally confiremed to correspond to the smooth region and the shear lip on the fracture face respectively.So the three dimensional stress structure of Mode I specimens can be derived through discussing the two fields respectively.The distribution of the Crack Tip Opening Displacement(CTOD)along the thickness direction and the three dimensional distribution of the void growth ratio(V_(g))near the crack tip are also obtained.The two fracture parameters are in similar trends along the thickness direction,and both of them can reflect the effect of thickness and that of the loading level to a certain degree.

Dys-psychological Stress Effect on Expressions of P53 and NFκBp65 in Human Ovarian Carcinoma In Vivo

Objective： To investigate the dys-psychological stress effect on the growth of subcutaneous xenotransplanted tumor in nude mice bearing human epithelium ovarian carcinoma, and the influence on P53 and NFκBp65 expressions. Methods： The subcutaneous tumor xenografts were established by implanting human epithelium ovarian carcinoma tissues into nude mice and the dys-psychological stress model was established with restraint. The mice were randomized into the following four treatment groups with each group six mice respectively： tumor group （group A）, normal saline intraperitoneal injection; tumor with stress group （group B）, normal saline intraperitoneal injection; tumor therapy group （group C）, cisplatin intraperitoneal injection; and tumor therapy with stress group （group D）, cisplatin intraperitoneal injection. The expressions of P53 and NFκBp65 in tumor tissues were determined by Western blotting. Results： The expressions of P53 and NFκBp65 in each restraint group were enhanced compared with the control groups （P0.05）. Conclusion： The dys-psychological stress may induce the high expressions of P53 and NFκBp65 proteins and further promote tumor growth.

Three dimensional K-Tz stress fields around the embedded center elliptical crack front in elastic plates

Through detailed three-dimensional （3D） finite element （FE） calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are obtained near the crack front with aspect ratios （a/c） of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. Tz decreases from an approximate value of Poisson ratio v at the crack tip to zero with increasing normalized radial distances （r/a） in the normal plane of the crack front line, and increases gradually when the elliptical parameter angle φ changes from 0° to 90°at the same r/a. With a/c rising to 1.0, Tz is getting nearly independent of φ and is only related to r/a. Based on the present FE calculations for Tz, empirical formulas for Tz are obtained to describe the 3D distribution of Tz for embedded center-elliptical cracks using the least squares method in the range of 0.2 ≤ a/c ≤ 1.0. These Tz results together with the corresponding stress intensity factor K are well suitable for the analysis of the 3D embedded centerelliptical crack from field, and a two-parameter K-Tz principle is proposed.

Modeling of Fatigue Crack Growth Closure Considering the Integrative Effect of Cyclic Stress Ratio,Specimen Thickness and Poisson's Ratio

Key components of large structures in aeronautics industry are required to be made light and have long enough fatigue lives.It is of vital importance to estimate the fatigue life of these structures accurately.Since the FCG process is affected by various factors,no universal model exists due to the complexity of the mechanisms.Most of the existing models are obtained by fitting the experimental data and could hardly describe the integrative effect of most existing factors simultaneously.In order to account for the integrative effect of specimen parameters,material property and loading conditions on FCG process,a new model named integrative influence factor model（IIF） is proposed based on the plasticity-induced crack closure theory.Accordingly to the predictions of crack opening ratio（γ） and effective stress intensity factor range ratio（U） with different material under various loading conditions,predictions of γ and U by the IIF model are completely identical to the theoretical results from the plane stress state to the plane strain state when Poisson＇s ratio equals 1/3.When Poisson＇s ratio equals 0.3,predictions of γ and U by the IIF model are larger than the predictions by the existing model,and more close to the theoretical results.In addition,it describes the influence of R ratios on γ and U effectively in the whole region from-1.0 to 1.0.Moreover,several sets of test data of FCG rates in 5 kinds of aluminum alloys with various specimen thicknesses under different loading conditions are used to validate the IIF model,most of the test data are situated on the predicted curves or between the two curves that represent the specimen with different thicknesses under the same stress ratio.Some of the test data slightly departure from the predictions by the IIF model due to the surface roughness and errors in measurement.Besides,based on the analysis of the physical rule of crack opening ratios,a relative thickness of specimen is defined to describe the influence of material property,specimen thickness and so forth on FCG characteristics conveniently.In conclusion,the relative thickness of specimen simplifies the expression of FCG characteristic and provides a general parameter to analyze the fatigue characteristics of different materials with various thicknesses under different loading conditions.The IIF model describes the integrative effect of existing influence factors explicitly and quantitatively,and provides a helpful tool for fatigue property estimation of practical component and experiment design.

Topology Optimization with Aperiodic Load Fatigue Constraints Based on Bidirectional Evolutionary Structural Optimization

Because of descriptive nonlinearity and computational inefficiency,topology optimization with fatigue life under aperiodic loads has developed slowly.A fatigue constraint topology optimization method based on bidirectional evolutionary structural optimization(BESO)under an aperiodic load is proposed in this paper.In viewof the severe nonlinearity of fatigue damagewith respect to design variables,effective stress cycles are extracted through transient dynamic analysis.Based on the Miner cumulative damage theory and life requirements,a fatigue constraint is first quantified and then transformed into a stress problem.Then,a normalized termination criterion is proposed by approximatemaximum stress measured by global stress using a P-normaggregation function.Finally,optimization examples show that the proposed algorithm can not only meet the requirements of fatigue life but also obtain a reasonable configuration.

CONSTRAINTS ON THE CRACK TIP PLASTIC FIELDS UNDER LARGE SCALE YIELD

Plane strain elastic-plastic finite element analyses are used tostudy the stress, strain fields around a straight crack in powerhardening plastic material. Center crack panel (CCP), single edgecrack panel (SECP) and double edge crack panel (DECP) tensionspecimens are analyzed with various crack lengths. Two localconstraint parameters, i.e. in-plane stress ratio T_x andout-of-plane constraint T_z are an- alyzed, which are defined astangential stress dividing normal (open) stress and out-of-planestress dividing the sum of tangential stress and normal stressrespectively.

Three-Dimensional Stress Fields in Finite Thickness Plate with Hole Under Shear Load

The theoretical solutions are obtained for the three-dimensional(3-D)stress field in an infinite isotropic elastic plate with a through-the-thickness circular hole subjected to shear load at far field by using Kane and Mindlin′s assumption based on the stress function method.Based on the present solutions,the characteristics of 3-D stress field are analyzed and the emphasis is placed on the effects of the plate thickness and Poisson′s ratio on the deviation of the present 3-D in-plane stress from the related plane stress solutions,the stress concentration and the out-of-plane constraint.The present solutions show that the stress concentration factor reaches its peak value of about 8.9% which is higher than that of the plane stress solutions.As expected,the out-of-plane stress constraint factor can reach 1on the surface of the hole when the plate is a very thick one.

Advanced test methods of material property characterization:high strain-rate testing and experimental simulation of multiaxial stress states

Optimum utilization of the loading capability of engineering materials is an important and active contribution to protect nature's limited resources,and it is the key for economic design methods.In order to make use of the materials' resources,those must be known very well;but conventional test methods will offer only limited informational value.The range of questions raised is as wide as the application of engineering materials,and partially they are very specific.The development of huge computer powers enables numeric modelling to simulate structural behaviour in rather complex loading environments-so the real material behaviour is known under the given loading conditions.Here the art of material testing design starts.To study the material behaviour under very distinct and specific loading conditions makes it necessary to simulate different temperature ranges,loading speeds, environments etc.and mostly there doesn't exist any commonly agreed test standard.In this contribution two popular,non-standard test procedures and test systems will be discussed on the base of their application background,special design features as well as test results and typically gained information:The demand for highspeed tests up to 1000 s^(-1) of strain rate is very specific and originates primarily in the automotive industry and the answers enable CAE analysis of crashworthiness of vehicle structures under crash conditions.The information on the material behaviour under multiaxial loading conditions is a more general one.Multiaxial stress states can be reduced to an equivalent stress,which allows the evaluation of the material's constraint and criticality of stress state.Both discussed examples shall show that the open dialogue between the user and the producer of testing machines allows custom-tailored test solutions.

考虑湿度应力的深部软岩隧道大变形控制研究

为探明深部高地应力软岩隧道大变形产生机理,制定适应大变形控制措施,以月直山隧道为工程依托,首先明确求解围岩形变压力与松动压力的Kastner与Caquot公式,基于岩体弹塑性力学与连续介质理论建立围岩湿度应力解析解,采用收敛-约束法绘制出3种围岩应力作用下围岩与初期支护特征曲线,分析断面变形过程中围岩与支护结构相互作用规律及3种应力占比演化规律。分别以混凝土与型钢钢架作为二衬支护结构对月直山隧道围岩稳定性进行计算分析,明确考虑湿度应力与松动压力条件下隧道二衬最优支护时机与支护参数,以支护结构安全系数FS评判出最优支护方案并对隧道软岩大变形进行治理。结果表明:当围岩径向位移达到550 mm时,仅考虑形变压力Pi时围岩对支护结构的压力P仅为0.813 MPa,考虑湿度应力Pw与松动压力Pa时P增大为1.372 MPa,湿度应力与松动应力总占比达40.7%,仅考虑形变压力而设计的支护结构不满足围岩稳定性要求;根据“卸压支护”理念,确定以“位移释放峰值upeak=0.325 m”与“间距d=0.7 m”作为增设第二层钢架的最佳支护时机与支护参数,现场二次钢架设置24天后使围岩变形收敛于50.8 mm,围岩大变形得到控制,研究成果可为今后相关隧道工程设计与施工提供参考。

基于物理约束的分布式神经网络三维地应力预测模型

地应力是油气藏全生命周期中井壁稳定分析、水力压裂设计、地层出砂预防、套损预测防治和油气开采措施制定的重要基础参数.针对缺乏三维地震资料情况下的三维地应力场预测问题,以工区已钻井地应力测井解释结果为基础,建立了一种基于物理约束的分布式神经网络三维地应力预测模型.首先,采用测井解释方法获得工区已钻井地应力单井剖面,并通过Kriging插值构建三维地质模型;其次,将测井数据和三维地质模型的三维空间坐标喂入3个并行的全连接神经网络,引入地应力物理约束条件,提出了数据和物理双约束的条件,进而通过三维空间坐标实现三维地应力场的预测;再次,优选了数据和物理双约束的权重参数,与人工神经网络、支持向量回归、随机森林等三种机器学习模型以及克里金插值模型进行了对比;最后,评价了不同机器学习模型预测地应力的效果.研究结果表明:(1)研究区域的地应力满足潜在正断层应力状态,即垂向地应力>最大水平地应力>最小水平地应力;(2)数据和物理双约束权重参数取值对预测结果有显著影响,当权重参数λ=0.2时,三维地应力预测效果最优;(3)与人工神经网络、支持向量回归、随机森林、克里金插值等模型相比,本文模型预测得到的地层三维地应力和孔隙压力更准确,在工区内测试集预测得到的垂向地应力、最大水平地应力、最小水平地应力和地层孔隙压力的最大相对误差分别为0.63%、7.59%、7.16%、3.21%,而且,本文模型能够更加准确捕捉地应力梯度的三维变化特征,水平井段地应力和地层孔隙压力参数的正态分布特征明显.结论认为,物理约束的分布式神经网络模型能有效融合地形构造与地应力之间的联系,提高三维地应力预测的准确性和可解释性,为油田三维地应力场的预测提供了一种新的思路.

约束状态下的有机玻璃零件的加工方法探讨

在生产实践中发现有机玻璃零件处于约束状态时,其所在环境存在较大温差变化,容易产生塑性变形,从而使产品不合格以致报废。该文对处于约束状态下的有机玻璃零件加工工艺进行了研究,重点对托盘组件加工过程中受约束状态下的有机玻璃盘形件加工工艺过程进行了改进、优化和实验验证,得出了约束状态下的有机玻璃零件加工工艺方案,并在实际加工生产中验证了此加工工艺方案的可行性。此方案对有机玻璃零件在自由状态和约束状态下的加工均具有指导意义,具有一定的推广价值。

周向约束对木材压缩吸能特性的影响

通过准静态压缩试验,考察轻木、泡桐、杉木、毛白杨4种木材在周向约束条件时,顺纹和横纹的压缩力学行为、破坏形态、吸能特性,并利用t检验分析其差异显著性。结果表明:①周向约束木材的顺纹压缩应力-应变曲线仍表现为“三阶段”特征,周向约束可抑制木材剪切扩展、大幅屈曲变形等非稳态破坏,曲线平台阶段应力基本与破坏强度持平。②泡桐、杉木、毛白杨顺纹压缩的吸能能力、平台阶段平稳性显著提高,其单位体积吸能值分别提高了15.51%、53.02%、48.15%,理想吸能率增大了35.93%、24.24%、86.00%。③周向约束木材的径向与弦向压缩力学行为、吸能特性相似,木材压头侧边纤维受试样中部拉伸和剪切作用发生竖向剥离,并沿压缩方向产生裂纹。④与无约束压缩相比,周向约束轻木、泡桐横纹压缩应力-应变曲线形状无明显变化,径向压缩单位体积吸能值分别提高了325.51%、182.60%,弦向压缩单位体积吸能值分别提高了430.57%、255.66%;杉木、毛白杨周向约束时横纹压缩应力-应变曲线变为无平台区的“两阶段”形式,理想吸能率均显著降低。

残余应力与拘束对GH4169合金疲劳及蠕变-疲劳裂纹扩展速率影响研究

为深入理解残余应力与拘束在疲劳载荷下的交互作用,以镍基高温合金GH4169为研究对象,选用紧凑拉伸(CT)试样,对不同拘束CT试样的上方施加不同大小的预加载荷从而在裂尖产生不同应力。将该应力作为残余应力,系统研究残余应力和拘束交互作用下GH4169合金的疲劳和蠕变–疲劳裂纹扩展速率。结果表明:随着裂尖残余压应力的增加,不同拘束下GH4169合金的疲劳和蠕变–疲劳裂纹扩展速率均降低。与低拘束试样相比,高拘束试样的疲劳和蠕变–疲劳裂纹扩展速率对残余应力的变化更加敏感,这主要与裂尖Mises应力和垂直于裂纹扩展方向的正应力有关。与疲劳裂纹扩展速率相比,蠕变-疲劳裂纹扩展速率对残余应力的变化更加敏感。

ALTERNATING OPTIMIZATION METHOD FOR ISOGEOMETRIC TOPOLOGY OPTIMIZATION WITH STRESS CONSTRAINTS

Topology optimization (TO) has developed rapidly recently. However, topology optimization with stress constraints still faces many challenges due to its highly non-linear properties which will cause inefficient computation, iterative oscillation, and convergence guarantee problems. At the same time, isogeometric analysis (IGA) is accepted by more and more researchers, and it has become one important tool in the field of topology optimization because of its high fidelity. In this paper, we focus on topology optimization with stress constraints based on isogeometric analysis to improve computation efficiency and stability. A new hybrid solver combining the alternating direction method of multipliers and the method of moving asymptotes (ADMM-MMA) is proposed to solve this problem. We first generate an initial feasible point by alternating direction method of multipliers (ADMM) in virtue of the rapid initial descent property. After that, we adopt the method of moving asymptotes (MMA) to get the final results. Several benchmark examples are used to verify the proposed method, and the results show its feasibility and effectiveness.

圆环非线性恢复力的梁约束模型建模

圆环隔振器是一种以圆环结构为基础的非线性隔振器,圆环结构在压缩变形过程中会受到压力和环面拉伸的耦合作用,从而产生非线性恢复力,对圆环结构非线性恢复力的精确建模是研究隔振器性能的关键。将圆环结构等分为多段曲梁,利用梁约束模型对每段曲梁建立计及几何非线性的力-位移关系模型,结合曲梁间的力传递关系和几何约束关系,建立了圆环整体在压缩过程中的非线性恢复力模型,并计算了圆环变形过程中所有分段点处的正应力。通过电子伺服疲劳试验机对圆环结构在压缩过程中的恢复力进行了测量,验证了梁约束模型的建模精度。研究结果表明,利用梁约束模型可以表征圆环结构非线性恢复力特性,其建模精度与椭圆积分法相当,而模型表达式和求解过程都比椭圆积分法简洁。圆环在压缩量最大时正应力最大,此时最大正应力在圆环的上下端点。梁约束模型的建模精度随着分段数的增加而提高,当分段数大于12时,梁约束模型的恢复力计算误差小于2%。