Analytical solution for the effective elastic properties of rocks with the tilted penny-shaped cracks in the transversely isotropic background

Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with the isotropic background,while the explicit model for the cracked rock with the anisotropic background is rarely investigated in spite of such case being often encountered in the earth.Hence,we first studied dependences of the crack opening displacement tensors on the crack dip angle in the coordinate systems formed by symmetry planes of the crack and the background anisotropy,respectively,by forty groups of numerical experiments.Based on the conclusion from the experiments,the analytical solution was derived for the effective elastic properties of the rock with the inclined penny-shaped cracks in the transversely isotropic background.Further,we comprehensively analyzed,according to the developed model,effects of the crack dip angle,background anisotropy,filling fluid and crack density on the effective elastic properties of the cracked rock.The analysis results indicate that the dip angle and background anisotropy can significantly either enhance or weaken the anisotropy degrees of the P-and SH-wave velocities,whereas they have relatively small effects on the SV-wave velocity anisotropy.Moreover,the filling fluid can increase the stiffness coefficients related to the compressional modulus by reducing crack compliance parameters,while its effects on shear coefficients depend on the crack dip angle.The increasing crack density reduces velocities of the dry rock,and decreasing rates of the velocities are affected by the crack dip angle.By comparing with exact numerical results and experimental data,it was demonstrated that the proposed model can achieve high-precision estimations of stiffness coefficients.Moreover,the assumption of the weakly anisotropic background results in the consistency between the proposed model and Hudson's published theory for the orthorhombic rock.

Analysis of stresses at the center of transversely isotropic Brazilian disk

This article presents the stresses at the center of a Brazilian disk(BD)for transversely isotropic rocks.It is shown that the solution of stresses at the center of an anisotropic disk is a function of the disk radius and the magnitude of applied load,as well as the material orientation with respect to the load axis and two dimensionless ratios with specific physical meanings and limitations.These two dimensionless parameters are the ratios of Young’s modulus and apparent shear modulus,although the ratio of apparent shear modulus will be eliminated if the Saint-Venant assumption is considered.Considerable finite element simulations are carried out to find the stresses at the disk center concerning the material orientation and the two dimensionless parameters.Also,an approximate formula obtained from analytical results,previously proposed in the literature for solving the tensile and compressive stresses at the disk center,is re-written and simplified based on these new definitions.The results of the approximate formula fitted to the analytical results are compared to those obtained from numerical solutions,suggesting a good agreement between the numerical and analytical methods.An approximate equation for the shear stress at the disk center is also formulated based on the numerical results.Finally,the influence of the assumptions for simplification of the proposed formula for the tensile,compressive,and shear stresses at the disk center is discussed,and simple and practical equations are proposed as estimations for the stresses at the center of the BD specimen for low to moderate anisotropic rocks.For highly anisotropic rocks,the reference plots can be used for more accuracy.

Exact simulation for direction-dependent large elastic strain responses of soft fibre-reinforced composites

An explicit form of the elastic strain-energy function for direction-dependent large elastic strain behaviors of soft fiber-reinforced composites is first presented based upon a decoupled approach for simulating complex nonlinear coupling effects.From this form,the exact closed-form solutions are then obtained for the uniaxial tension responses in the fiber and cross-fiber directions.With such exact solutions,the issue of simultaneously simulating strongly coupling nonlinear responses in the fiber and cross-fiber directions may be reduced to the issue of separately treating each decoupled uniaxial stress-strain response,thus bypassing usual complexities and uncertainties involved in identifying a large number of strongly coupled adjustable parameters.The numerical examples given are in good agreement with the experimental data for large strain responses.

Effects of 3D deformation and nonlinear stress–strain relationship on the Brazilian test for a transversely isotropic rock

To improve the accuracy of indirect tensile strength for a transversely isotropic rock in the Brazilian test, this study considered the three-dimensional (3D) deformation and the nonlinear stress–strain relationship. A parametric study of a numerical Brazilian test was performed for a general range of elastic constants, revealing that the 3D modeling evaluated the indirect tensile strength up to 40% higher than the plane stress modeling. For the actual Asan gneiss, the 3D model evaluated the indirect tensile strength up to 10% higher and slightly enhanced the accuracy of deformation estimation compared with the plane stress model. The nonlinearity in stress–strain curve of Asan gneiss under uniaxial compression was then considered, such that the evaluated indirect tensile strength was affected by up to 10% and its anisotropy agreed well with the physical intuition. The estimation of deformation was significantly enhanced. The further validation on the nonlinear model is expected as future research.

Complete solutions for elastic fields induced by point load vector in functionally graded material model with transverse isotropy

The paper develops and examines the complete solutions for the elastic field induced by the point load vector in a general functionally graded material(FGM)model with transverse isotropy.The FGMs are approximated with n-layered materials.Each of the n-layered materials is homogeneous and transversely isotropic.The complete solutions of the displacement and stress fields are explicitly expressed in the forms of fifteen classical Hankel transform integrals with ten kernel functions.The ten kernel functions are explicitly expressed in the forms of backward transfer matrices and have clear mathematical properties.The singular terms of the complete solutions are analytically isolated and expressed in exact closed forms in terms of elementary harmonic functions.Numerical results show that the computation of the complete solutions can be achieved with high accuracy and efficiency.

Influence of Anisotropic Permeability and Soret Effect on the Convective Heat and Mass Transfer through a Porous Cavity Saturated by a Non-Newtonian Fluid

In this work, an analytical study is carried out on double-diffusive natural convection through a horizontal anisotropic porous layer saturated with a non-Newtonian fluid by using the Darcy model with the Boussinesq approximations. The horizontal walls of the system are subject to vertical uniform fluxes of heat and mass, whereas the vertical walls are assumed to be adiabatic and impermeable. The Soret effect is taken into consideration. Based on parallel flow approximation theory, the problem is solved in the limit of a thin layer and documented the effects of the physical parameters describing this investigation.

BENDING OF FUNCTIONALLY GRADED PIEZOELECTRIC RECTANGULAR PLATES

By introducing two displacement functions as well as two stressfunctions, two independent state equations with variable coefficientsare derived from the three-dimensional theory equations of piezo-elasticity for transverse isotropy. A laminated approximation is usedto transform the state equations to those with constant coefficientsin each sub-layer. The bending problem of a functionally gradedrectangular plate is then analyzed based on the state equations.Numerical results are presented and the effect of material gradi- entindex is discussed.

OPTIMAL DESIGN OF A 6-AXIS FORCE TRANSDUCER BASED ON STEWART PLATFORM RELATED TO SENSITIVITY ISOTROPY

The design method of a 6-axis force robot's transducer based on the Stewartplatform is detailed. For this purpose, the sensitivity isotropy evaluation criteria of thetransducer are defined, and by the aid of computer, the relationships between the criteria and theparameters of all the transducers based on the Stewart platform are investigated within thegeometric model of the solution space, which can provide the theoretical background for the optimalconstruction design of the 6-axis force transducer related to the sensitivity isotropy.

Some recent advances in 3D crack and contact analysis of elastic solids with transverse isotropy and multifield coupling

Significant progress has been made in mixed boundary-value problems associated with three-dimensional（3D） crack and contact analyses of advanced materials featuring more complexities compared to the conventional isotropic elastic materials.These include material anisotropy and multifield coupling,two typical characteristics of most current multifunctional materials.In this paper we try to present a state-of-the-art description of 3D exact/analytical solutions derived for crack and contact problems of elastic solids with both transverse isotropy and multifield coupling in the latest decade by the potential theory method in the spirit of V.I.Fabrikant.whose ingenious breakthrough brings new vigor and vitality to the old research subject of classical potential theory.We are particularly interested in crack and contact problems with certain nonlinear features.Emphasis is also placed on the coupling between the temperature field（or the like） and other physical fields（e.g.,elastic,electric,and magnetic fields）.We further highlight the practical significance of 3D contact solutions,in particular in applications related to modern scanning probe microscopes.

Interplay between Micro-Anisotropy and Macro-Isotropy on Evolution of Non-Equilibrium Morphology

Non-equilibrium morphology has received much attention from both scientific and engineering points of view for its intricate pattern selection mechanisms and useful industrial application. Most study of non-equilibrium is about the metal, alloy and other simple system. The complex silicate system is rarely involved. However, silicate is very important in geosciences and ceramic industry. In this paper, two kinds of non-equilibrium crystal morphologies of silicate: dendrite of diopside and spherulite of plagioclase, were introduced. Combining with the other kinds of non-equilibrium morphologies, the characteristics of micro-macro and anisotropy-isotropy of the non-equilibrium morphologies were discussed. Dendrite of diopside is micro- and macro-anisotropic, spherulite of plagioclase is micro-anisotropic, but macro-isotropic, fractal of NH4Cl is also micro-anisotropic, but macro-isotropic, dense-branching morphology (DBM) formed in non-crystalline system is micro-and macro-isotropic. Based on the micro-macro interplay on the pattern formation, it is proposed that the interplay between micro-anisotropy of crystal structure vs macro-isotropy of undercooling in crystal growth system will control the morphological evolution. The nucleation rate related to the anisotropy for the morphological evolution was also discussed. The fact that diopside develops dendrite and plagioclase develop spherulite in our experiment is due to their structural anisotropy difference.

A local homogenization approach for simulating the reinforcement effect of the fully grouted bolt in deep underground openings

Fully grouted bolts are a key component of the support system for underground openings.Although considerable effort has been made in the simulation of the reinforcement effect of the fully grouted bolts on the rock masses surrounding underground openings,most of the work has limited significance since the structural element approach is used.This study proposes a local homogenization approach(L-H approach)that integrates elastoplastic mechanics,composite mechanics,and analytical approaches with numerical simulation to effectively simulate the reinforcement effect of the fully grouted bolt on deep surrounding rock masses.In the L-H approach,the representative volume of bolted rock mass(RVBRM)with a fully grouted bolt is established based on the original mesh model utilized in the rockbolt element approach.The RVBRM is a regular quadrangular prism with a cross-sectional size equal to the bolt spacing and a length equal to the bolt length.The RVBRM is homogenized by the L-H approach from a unidirectional bolt-reinforced composite into a homogeneous transversely isotropic medium whose mechanical properties are described by a new transversely isotropic elastoplastic model.The L-H parameters for the RVBRM are obtained using analytical approaches,composite mechanics,and known parameters of the rock mass and bolt.Using the L-H approach,the reinforcement effect of the fully grouted bolt on the bolted rock specimen and the surrounding rock mass in Jinping II Diversion Tunnel#2 with a depth greater than 2000 m is simulated.The results show that the predictions of the L-H approach are more in agreement with the physical model results of bolted rock specimen and provide a more realistic response of the bolted surrounding rock mass.The L-H approach demonstrates that fully grouted bolts with common bolt spacings and diameters substantially enhance the elastic modulus,shear strength,and tensile strength of the rock mass in the direction of the bolt axis.

Force/Moment Isotropy of 8/4-4 Parallel Six-Axis Force Sensor Based on Performance Atlases

A six-axis force sensor with parallel 8/4-4 structure is introduced and its measurement principle is analyzed.Based on condition numbers of Jacobian matrix spectral norm of the sensor,the relationship between the force and moment isotropy and some structural parameters is deduced.Orthogonal test methods are used to determine the degree of primary and secondary factors that have significant effect on sensor characteristics.Furthermore,the relationship between each performance index and the structural parameters of the sensor is analyzed by the method of the atlas,which lays a foundation for structural optimization design of the force sensor.

The reflection and transmission of elastic waves at a interface between a transversely isotropi celastic soils and saturated soils

According to Biot′s wave equation of transversely isotropic saturated soil, this paper deduces the general equation of the reflection coefficients and transmission coefficients when qP 1 wave goes through from saturated soil to elastic media. The effects of anisotropies and boundary drainage condition on reflection coefficients and transmission coefficients are analyzed by numerical method. The idea of this paper can be applied to the case when qSV wave or qP 2 wave goes through from saturated soils to elastic soils.

Modeling shale with consideration of bedding planes by cohesive finite element method

Shale contains distributed directional bedding planes,which make the shale transverse isotropic.To model shale with consideration of the bedding planes,a cohesive finite element method(CFEM)is developed based on the randomized triangular mesh.The interface orientation generated from such mesh tends to be uniformly distributed with the element number increasing.To represent the bedding plane,the interfaces aligned with the bedding plane are assigned the cohesive law that characterizes the bedding plane while the other interfaces are assigned the cohesive law that characterizes the matrix.By this means,the anisotropy characteristics of the stiffness and the strength of shale are well represented.The simulation examples demonstrate that the bedding plane has a significant influence on the fracture trajectory,which is consistent with the observation in the experiment.It is suggested that this modeling method of shale is feasible.It provides an alternative approach to fracture simulation in shale.

NOVEL 6-SPS PARALLEL 3-DIMENSIONAL PLATFORM MANIPULATOR AND ITS FORCE/MOTION TRANSMISSION ANALYSIS

The unique design for a novel 6-SPS parallel 3-dimensional platformmanipulator with an orthogonal configuration is investigated. The layout feature of the parallelmanipulator is described. Its force/motion transmission capability, evaluation criteria arepresented. At the orthogonal configuration, the criteria and the relationships between the criteriaand the link lengths are analyzed, which is important since it can provide designer a piece ofvaluable information about how to choose the linear actuators. From the analysis of the results itis shown that the force/motion transmission capabilities of the parallel manipulator arecharacterized by isotropy at the orthogonal configuration. The manipulator is particularly suitablefor certain applications in 6-DOF micromanipulators and 6-axis force/moment transducers.

ON THE BENDING, VIBRATION AND STABILITY OF LAMINATED RECTANGULAR PLATES WITH TRANSVERSELY ISOTROPIC LAYERS

A method based on newly presented state space formulations is developed for analyzing the bending, vibration and stability of laminated transversely isotropic rectangular plates with simply supported edges. By introducing two displacement functions and two stress functions, two independent state equations were constructed based on the three-dimensional elasticity equations for transverse isotropy. The original differential equations are thus decoupled with the order reduced that will facilitate obtaining solutions of various problems. For the simply supported rectangular plate, two relations between the state variables at the tap and bottom surfaces were established. In particular, for the free vibration ( stability) problem, it is found that there exist two independent classes: One corresponds to the pure in-plane vibration (stability) and the other to the general bending vibration (stability). Numerical examples are finally presented and the effects of same parameters are discussed.

垂直发育裂隙介质中PP波扰动法近似反射系数研究

为了实现裂隙介质中的地震AVO(amplitude variation with offset)反演,基于Christoffel方程和边界条件,推导了EDA介质中PP波、PS1波、PS2波精确反射系数表达式,并采用扰动法推导了极端弱各向异性介质中PP波的近似反射系数,通过弹性系数简化得到HTI、裂隙EDA介质中PP波的近似反射系数。模型计算结果表明:在弱各向异性条件下,HTI(horizontal transverse isotropy)介质中地震波的近似反射系数公式计算精度较高,相对误差在4%以内;当裂隙介质在弱各向异性(各向异性系数<0.1)的情况下,当入射角在40°以内时,PP波近似反射系数与精确反射系数的绝对误差在10^(-4)以内,相对误差小于4%,但随着入射角的加大其计算误差有所增加。对各向异性较强烈(各向异性系数达到0.2)的裂隙介质,当入射角相对较小时(小于45°时),裂隙介质中PP波近似反射系数计算式对于强各向异性的介质仍然成立。通过对反射系数的近似研究,可以将裂隙介质中反射系数的非线性问题转为线性问题,进而利用这些特性进行参数反演,有利于提高反演速度。

Water inflow forecasting for tunnel considering nonlinear variation of permeability coefficient

To assess the water inflow which is more suitable to the actual conditions of tunnel,an empirical correlation about the permeability coefficient changing with depth is introduced.Supposing that the surrounding rock is heterogeneous isotropy,the formula for calculating water inflow of tunnel with the nonlinear variation of permeability coefficient is deduced.By the contrast analysis with the existing formulas,the presented method has the similar value to them;moreover,the presented method has more simple form and easy to use.Due to parameter analysis,the water inflow decreases after considering the nonlinear variation of permeability coefficient.When the attenuation coefficient a>0,the water inflow increases first till reaches the maximum at a certain depth,then decreases and is close to 0 finally if deep enough.Thus,it is better to keep away from the certain depth where it is with the maximum water inflow for safe operation and economical construction,and reduce the water damage.Based on the analysis,the radius of tunnel has less impact on the amount of water inflow,and the water inflow just increases by 6.7% when the radius of tunnel increases by 1 m.

The Dosimetric Characteristics and Potential Limitation in Clinical Application of a Low Energy Photon Intra-Operative Radiotherapy System

Purpose: To investigate the dosimetric characteristics of a low energy photon intra-operative radiotherapy (IORT) system and explore its potential limitation in clinical application. Methods: A special water phantom, a parallel-plate ionization chamber and an electrometer were used to measure the depth dose rate, isotropy of dose distribution in X/Y plane, dosimetry reproducibility of bare probe and spherical applicators of different size which were used in comparison with the system data. Results: The difference in depth dose rate between the measurement and system data for bare probe is -2.16% ± 1.36%, the range of the relative deviation for isotropy in the X/Y plane is between -1.9% and 2.1%. The difference in depth dose rate, transfer coefficient, isotropy in X/Y plane between the measurement and system data for the whole set of spherical applicators is -10.0% - 2.3%, -8.9% - 4.2% and -1.6% - 2.6%, respectively. Higher surface dose rate and steeper gradient depth dose are observed in smaller spherical applicators. The depth dose rate and isotropy for bare probe and spherical applicators have been shown good reproducibility. The uncertainty of measurement is associated with the positioning accuracy, energy response, noise current and correction function f’(R). Conclusions: Thorough commissioning of the low energy photon IORT system helps us better understand the dosimetry characteristics, verify the system data, obtain adequate data for clinical application and routine quality assurance. The steep gradient depth dose and limited treatment range may restrain its potential in clinical application.

Mathematical osteon model for examining poroelastic behaviors

An extended and reasonable stress boundary condition at an osteon exte- rior wall is presented to solve the model proposed by Remond and Naili. The obtained pressure and fluid velocity solutions are used to investigate the osteonal poroelastic behaviors. The following results are obtained. （i） Both the fluid pressure and the velocity amplitudes are proportional to the strain amplitude and the loading frequency. （ii） In the physiological loading state, the key role governing the poroelastic behaviors of the osteon is the strain rate. （iii） At the osteon scale, the pressure is strongly affected by the permeability variations, whereas the fluid velocity is not.