Generalized response displacement methods for seismic analysis of underground structures with complex cross section

The response displacement method(RDM)is recommended for the seismic analysis of underground structures in the transverse direction for many codes,including bases for design of structures-seismic actions for designing geotechnical works(ISO 23469)and code for seismic design of urban rail transit structures(GB 50909-2014).However,there are some obvious limitations in the application of RDM.Springs and the shear stress of the soil could be approximately evaluated for the structures having a simple cross section,such as rectangular and circular structures.It is necessary to propose simplified seismic analysis methods for structures with complex cross sections.This paper refers to the idea of RDM and proposes three generalized response displacement methods(GRDM).In GRDM1,a part of the soil surrounding a structure is selected to generate a generalized underground structure with a rectangular cross section,and the same analysis model as RDM is applied to analyze the responses of the structure.In GRDM2,a hollow soil model without a generalized structure is used to compute the equivalent load caused by the relative displacement of the soil,and the soil-structure interaction model is applied to calculate the responses of the structure.In GRDM3,a continuous soil model is applied to compute the equivalent load caused by the relative displacement and shear stress of the soil,and the soil-structure interaction model is applied to analyze the responses of the structure,which is the same as the model used in GRDM2.The time-history analysis method(THAM)is used to evaluate the accuracy of the proposed simplified methods.Results show that the error of GRDM1 is about 20%,while the error is only 5%for GRDM2 and GRDM3.Among the three proposed methods,GRDM3 has obvious advantages regarding calculation efficiency and accuracy.Therefore,it is recommended to use GRDM3 for the seismic response analysis of underground structures that have conventional simple or complex cross sections.

Mechanical mechanism analysis of tension type anchor based on shear displacement method

Based on the fact that the shear stress along anchorage segment is neither linearly nor uniformly distributed, the load transfer mechanism of the tension type anchor was studied and the mechanical characteristic of anchorage segment was analyzed. Shear stress?strain relationship of soil surrounding anchorage body was simplified into three-folding-lines model consisting of elastic phase, elasto-plastic phase and residual phase considering its softening characteristic. Meanwhile, shear displacement method that has been extensively used in the analysis of pile foundation was introduced. Based on elasto-plastic theory, the distributions of displacement, shear stress and axial force along the anchorage segment of tension type anchor were obtained, and the formula for calculating the elastic limit load was also developed accordingly. Finally, an example was given to discuss the variation of stress and displacement in the anchorage segment with the loads exerted on the anchor, and a program was worked out to calculate the anchor maximum bearing capacity. The influence of some parameters on the anchor bearing capacity was discussed, and effective anchorage length was obtained simultaneously. The results show that the shear stress first increases and then decreases and finally trends to the residual strength with increase of distance from bottom of the anchorage body, the displacement increases all the time with the increase of distance from bottom of the anchorage body, and the increase of velocity gradually becomes greater.

An extended displacement discontinuity method for analysis of stress wave propagation in viscoelastic rock mass

An extended displacement discontinuity method (EDDM) is proposed to analyze the stress wave propagation in jointed viscoelastic rock mass (VRM).The discontinuities in a rock mass are divided into two groups.The primary group with an average geometrical size larger than or in the same order of magnitude of wavelength of a concerned stress wave is defined as 'macro-joints',while the secondary group with a high density and relatively small geometrical size compared to the wavelength is known as 'micro-defects'.The rock mass with micro-defects is modeled as an equivalent viscoelastic medium while the macro-joints in the rock mass are modeled explicitly as physical discontinuities.Viscoelastic properties of a micro-defected sedimentary rock are obtained by longitudinally impacting a cored long sedimentary rod with a pendulum.Wave propagation coefficient and dynamic viscoelastic modulus are measured.The EDDM is then successfully employed to analyze the wave propagation across macro-joint in VRM.The effect of the rock viscosity on the stress wave propagation is evaluated by comparing the results of VRM from the presented EDDM with those of an elastic rock mass (ERM) from the conventional displacement discontinuity method (CDDM).The CDDM is a special case of the EDDM under the condition that the rock viscosity is ignored.Comparison of the reflected and transmitted waves shows that the essential rock viscosity has a significant effect on stress wave attenuation.When a short propagation distance of a stress wave is considered,the results obtained from the CDDM approximate to the EDDM solutions,however,when the propagation distance is sufficiently long relative to the wavelength,the effect of rock viscosity on the stress wave propagation cannot be ignored.

Propagation simulation and dilatancy analysis of rock joint using displacement discontinuity method

A revised displacement discontinuity method(DDM) program is developed for the simulation of rock joint propagation and dilatancy analysis. The non-linear joint model used in the program adopts Barton-Bandis normal deformation model, Kulhaway shear deformation model and Mohr-Coulomb criterion. The joint propagation criterion is based on the equivalent stress intensity factor which can be obtained by regression analysis. The simulated rock joint propagation accords well with the existing knowledge. The closure and opening of joint is investigated by DDM, and it is shown that if the opening volume of propagated joint is larger than closure volume of the old joint, the joint dilatancy occurs. The dilatancy condition is mainly controlled by the normal stiffness of the rock joint. When the normal stiffness is larger than the critical value, joint dilatancy occurs. The critical normal stiffness of rock joint changes with the joint-load angle, and joint dilatancy is most possible to occur at 30°.

Analysis of Three-dimensional Crack Propagation by Using Displacement Discontinuity Method

A technique for modelling of three-dimensional(3D)quasi-statically propagating cracks in elastic bodies by the displacement discontinuity method(DDM)was described.When the crack is closed,the Mohr-coulomb rule on the two contacted surfaces of the crack must be satisfied.A simple iterative method was adopted in order to consider three different states of cracks.Under the assumption that the advance of the point on the crack front would occur only in the normal plane which is through this edge point,the maximum energy release rate criterion is modified to be used as the criterion for the crack growth.With discretization,the process of crack propagation can be seen as the advance of the vertices of the crack front.The program MCP3D was developed based on these theories to simulate the 3D quasi-static crack propagation.A numerical example of a penny-shaped crack subject to tension and compression in an infinite elastic media was analyzed with MCP3D,and the results in comparison with others' show that the present method for 3D crack propagation is effective.

Simulation of crack coalescence mechanism underneath single and double disc cutters by higher order displacement discontinuity method

The present research is focused on the numerical crack coalescence analysis of the micro-cracks and cracks produced during the cutting action of TBM disc cutters. The linear elastic fracture mechanics(LEFM) concepts and the maximum tangential stress criterion are used to investigate the micro crack propagation and its direction underneath the excavating discs. A higher order displacement discontinuity method with quadratic displacement discontinuity elements is used to estimate the stress intensity factors near the crack tips. Rock cutting mechanisms under single and double type discs are simulated by the proposed numerical method.The main purposes of the present modeling are to simulate the chip formation process of indented rocks by single and double discs.The effects of specific disc parameters(except speed) on the thrust force Ft, the rolling force Fr, and the specific energy ES are investigated. It has been shown that the specific energy(energy required to cut through a unit volume of rock) of the double disc is less than that of the single disc. Crack propagation in rocks under disc cutters is numerically modeled and the optimum ratio of disc spacing S to penetration depth Pd(i.e. S/Pd ratio) of about 10 is obtained, which is in good agreement with the theoretical and experimental results cited in the literature.

On the accuracy of higher order displacement discontinuity method(HODDM) in the solution of linear elastic fracture mechanics problems

The higher order displacement discontinuity method(HODDM) utilizing special crack tip elements has been used in the solution of linear elastic fracture mechanics(LEFM) problems. The paper has selected several example problems from the fracture mechanics literature(with available analytical solutions) including center slant crack in an infinite and finite body, single and double edge cracks, cracks emanating from a circular hole. The numerical values of Mode Ⅰ and Mode Ⅱ SIFs for these problems using HODDM are in excellent agreement with analytical results(reaching up to 0.001% deviation from their analytical results). The HODDM is also compared with the XFEM and a modified XFEM results. The results show that the HODDM needs a considerably lower computational effort(with less than 400 nodes) than the XFEM and the modified XFEM(which needs more than 10000 nodes) to reach a much higher accuracy. The proposed HODDM offers higher accuracy and lower computation effort for a wide range of problems in LEFM.

A method to calculate displacement factors using SVM

In order to improve the precision of mining subsidence prediction, a mathematical model using Support Vector Machine （SVM） was established to calculate the displacement factor. The study is based on a comprehensive analysis of factors affecting the displacement factor, such as mechanical properties of the cover rock, the ratio of mining depth to seam thickness, dip angle of the coal seam and the thickness of loose layer. Data of 63 typical observation stations were used as a training and testing sample set. A SVM regression model of the displacement factor and the factors affecting it was established with a kernel function, an insensitive loss factor and a properly selected penalty factor. Given an accurate calculation algorithm for testing and analysis, the results show that an SVM regression model can calcu- late displacement factor precisely and reliable precision can be obtained which meets engineering requirements. The experimental results show that the method to calculation of the displacement factor, based on the SVM method, is feasible. The many factors affecting the displacement factor can be consid- ered with this method. The research provides an efficient and accurate approach for the calculation of displacement in mining subsidence orediction.

Powell dynamic identification of displacement parameters of indeterminate thin-walled curve box based on FCSE theory

The FCSE controlling equation of pinned thinwalled curve box was derived and the indeterminate problem of continuous thin-walled curve box with diaphragm was solved based on flexibility theory. With Bayesian statistical theory,dynamic Bayesian error function of displacement parameters of indeterminate curve box was founded. The corresponding formulas of dynamic Bayesian expectation and variance were deduced. Combined with one-dimensional Fibonacci automatic search scheme of optimal step size,the Powell optimization theory was utilized to research the stochastic identification of displacement parameters of indeterminate thin-walled curve box. Then the identification steps were presented in detail and the corresponding calculation procedure was compiled. Through some classic examples,it is obtained that stochastic performances of systematic parameters and systematic responses are simultaneously deliberated in dynamic Bayesian error function. The one-dimensional optimization problem of the optimal step size is solved by adopting Fibonacci search method. And the Powell identification of displacement parameters of indeterminate thin-walled curve box has satisfied numerical stability and convergence,which demonstrates that the presented method and the compiled procedure are correct and reliable.During parameters鈥？iterative processes,the Powell theory is irrelevant with the calculation of finite curve strip element（FCSE） partial differentiation,which proves high computation effciency of the studied method.

Velocity plus displacement equivalent force control for real-time substructure testing

This paper employs a velocity plus displacement（V＋D）-based equivalent force control（EFC） method to solve the velocity/displacement difference equation in a real-time substructure test. This method uses type 2 feedback control loops to replace mathematical iteration to solve the nonlinear dynamic equation. A spectral radius analysis of the amplification matrix shows that the type 2 EFC-explicit, Newmark-β method has beneficial numerical characteristics for this method. Its stability limit of Ω = 2 remains unchanged regardless of the system damping because the velocity is achieved with very high accuracy during simulation. In contrast, the stability limits of the central difference method using direct velocity prediction and the EFC-average acceleration method with linear interpolation are shown to decrease with an increase in system damping. In fact, the EFC-average acceleration method is shown to change from unconditionally stable to conditionally stable. We also show that if an over-damped system with a damping ratio of 1.05 is considered, the stability limit is reduced to Ω =1.45. Finally, the results from an experiment with a single-degree-of-freedom structure installed with a magneto-rheological（MR） damper are presented. The results demonstrate that the proposed method is able to follow both displacement and velocity commands with moderate accuracy, resulting in improved test performance and accuracy for structures that are sensitive to both velocity and displacement inputs. Although the findings of the study are promising, additional test data and several further improvements will be required to draw general conclusions.

General displacement arch-cantilever element method for stress analysis of arch dam

Based on the general displacement method and the basic hypothesis of the trial load method, a new advanced trial load method, the general displacement arch-cantilever element method, was proposed to derive the transformation relation of displacements and loads between the surface nodes and middle plane nodes. This method considers the nodes on upstream and downstream surfaces of the arch dam to be exit nodes （master nodes）, and the middle plane nodes to be slave nodes. According to the derived displacement and load transformation matrices, the equilibrium equation treating the displacement of middle plane nodes as a basic unknown variable is transformed into one that treats the displacement of upstream and downstream nodes as a basic unknown variable. Because the surface nodes have only three degrees of freedom （DOF）, this method can be directly coupled with the finite element method （FEM）, which is used for foundation simulation to analyze the stress of the arch dam with consideration of dam-foundation interaction. Moreover, using the FEM, the nodal load of the arch dam can be easily obtained. Case studies of a typical cylindrical arch dam and the Wudongde arch dam demonstrate the robustness and feasibility of the proposed method.

A Study on Propagation Mechanism of Fracture Systems in Rock Masses by Discontinuity Displacement Method

The main task of fracture mechanics of rock masses is the study on the propagating mechanism of fractures in rock masses , which can be efficiently conducted by discontinuty displacement (DD) numerical evaluation . Firstly ,the element stress and displacement are analysed and the principle and steps of the numerical calculation of stress intensity factor and fracture extension force are introduced .The numerical results of parallel and echelon fracture systems ,which are compared with real field fractures .are presented. Finally . a simple engineering application example is presented .

A Breakthrough of Primary Structure in Indeterminate Analysi--Irregular Force Method and Irregular Displacement Method

As the strict limitation of primary structure in traditional force method and displacement method in indeterminate analysis may lead to complicated high-order linear equations, a breakthrough of the limitation, i.e., the application of irregular force method and irregular displacement method, would be introduced in this paper to ease the difficulty of hand computations. By using hyperstatic primary structures and partly chained primary structures, the primary structures of force method and displacement method are reformed, and the order of the system is decreased. The technique is explained through examples. The significance of the new method is summarized.

Numerical Simulation of Non-linear Joint Using Displacement Discontinuity Method

An iterative algorithm for modeling of non-linear joint by the displacement discontinuity method (DDM) was described, and the effect of the non-linear joint on the in-situ stress field was investigated in this paper. The Barton-Bandis (BB) non-linear joint model and failure criterion were adopted in the new DDM program. Using this program, the stress field around the non-linear joint was obtained, the parameters analysis of the joint was carried out, and the deformation and stress distribution of the joint were studied. The simulation results show that: (1)the in-situ stress is significantly affected by the joint; (2)the increase of stiffness, friction angle, and thickness of the joint affect the stress concentration in different ways; (3)the influence distance of the joint changes with the angle of the joint; (4)the deformation and stress of the joint change with the point position.

APPLICATION OF THE METHOD OF THE RECIPROCAL THEOREM TO FINDING DISPLACEMENT SOLUTIONS OF CUBES

In this paper the method of reciprocal theorem is extended to find solutions of three-D problems of elasticity.First we give the basic solution of the cube with six surfaces fixed as the basic system and then using the reciprocal theorem between the basic system acted on by unit concentrated loads and the actual system with prescribed surface displacements, we find displacement solution of the actual system.

Analysis of fracture propagation and shale gas production by intensive volume fracturing

This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation model takes into account the interaction between hydraulic fracture and natural fracture by means of the displacement discontinuity method(DDM)and the Picard iterative method.The shale gas flow considers multiple transport mechanisms,and the flow in the fracture network is handled by the embedded discrete fracture model(EDFM).A series of numerical simulations are conducted to analyze the effects of the cluster number,stage spacing,stress difference coefficient,and natural fracture distribution on the stimulated fracture area,fractal dimension,and cumulative gas production,and their correlation coefficients are obtained.The results show that the most influential factors to the stimulated fracture area are the stress difference ratio,stage spacing,and natural fracture density,while those to the cumulative gas production are the stress difference ratio,natural fracture density,and cluster number.This indicates that the stress condition dominates the gas production,and employing intensive volume fracturing(by properly increasing the cluster number)is beneficial for improving the final cumulative gas production.

Adaptive Soft-sensor Modeling Algorithm Based on FCMISVM and Its Application in PX Adsorption Separation Process

To overcome the problem that soft sensor models cannot be updated with the process changes, a soft sensor modeling algorithm based on hybrid fuzzy c-means （FCM） algorithm and incremental support vector machines （ISVM） is proposed. This hybrid algorithm FCMISVM includes three parts： samples clustering based on FCM algorithm, learning algorithm based on ISVM, and heuristic sample displacement method. In the training process, the training samples are first clustered by the FCM algorithm, and then by training each clustering with the SVM algorithm, a sub-model is built to each clustering. In the predicting process, when an incremental sample that represents new operation information is introduced in the model, the fuzzy membership function of the sample to each clustering is first computed by the FCM algorithm. Then, a corresponding SVM sub-model of the clustering with the largest fuzzy membership function is used to predict and perform incremental learning so the model can be updated on-line. An old sample chosen by heuristic sample displacement method is then discarded from the sub-model to control the size of the working set. The proposed method is applied to predict the p-xylene （PX） purity in the adsorption separation process. Simulation results indicate that the proposed method actually increases the model＇s adaptive abilities to various operation conditions and improves its generalization capability.

DDM regression analysis of the in-situ stress field in a non-linear fault zone

A multivariable regression analysis of the in-situ stress field, which considers the non-linear deformation behavior of faults in practical projects, is presented based on a newly developed three-dimensional displacement discontinuity method （DDM） program. The Bar- ton-Bandis model and the Kulhaway model are adopted as the normal and the tangential deformation model of faults, respectively, where the Mohr-Coulomb failure criterion is satisfied. In practical projects, the values of the mechanical parameters of rock and faults are restricted in a bounded range for in-situ test, and the optimal mechanical parameters are obtained from this range by a loop. Comparing with the traditional finite element method （FEM）, the DDM regression results are more accurate.

Settlement calculation for long-short composite piled raft foundation

The mechanism of long-short composite piled raft foundation was discussed. Assuming the relationship between shear stress and shear strain of the surrounding soil was elasto-plastic, shear displacement method was employed to establish the different explicit relational equations between the load and the displacement at the top of pile in either elastic or elasto-plastic period. Then Mylonakis ＆ Gazetas model was introduced to simulate the interaction between two piles or between piles and soil. Considering the effect of cushion, the flexible coefficients of interaction were provided, With the addition of a relevant program, the settlement calculation for long-short composite piled raft foundation was developed which could be used to account for the interaction of piles, soil and cushion. Finally, the calculation method was used to analyze an engineering example. The calculated value of settlement is 10.2 ram, which is close to the observed value 8.8 mm.

Rapid Springback Compensation for Age Forming Based on Quasi Newton Method

Iterative methods based on finite element simulation are effective approaches to design mold shape to compensate springback in sheet metal forming. However, convergence rate of iterative methods is difficult to improve greatly. To increase the springback compensate speed of designing age forming mold, process of calculating springback for a certain mold with finite element method is analyzed. Springback compensation is abstracted as finding a solution for a set of nonlinear functions and a springback compensation algorithm is presented on the basis of quasi Newton method. The accuracy of algorithm is verified by developing an ABAQUS secondary development program with MATLAB. Three rectangular integrated panels of dimensions 710 mmx750 mm integrated panels with intersected ribs of 10 mm are selected to perform case studies. The algorithm is used to compute mold contours for the panels with cylinder, sphere and saddle contours respectively and it takes 57%, 22% and 33% iterations as compared to that of displacement adjustment （DA） method. At the end of iterations, maximum deviations on the three panels are 0.618 4 mm, 0.624 1 mm and 0.342 0 mm that are smaller than the deviations determined by DA method （0.740 8 mm, 0.740 8 mm and 0.713 7 mm respectively）. In following experimental verification, mold contour for another integrated panel with 400 ram^380 mm size is designed by the algorithm. Then the panel is age formed in an autoclave and measured by a three dimensional digital measurement devise. Deviation between measuring results and the panel＇s design contour is less than 1 mm. Finally, the iterations with different mesh sizes （40 mm, 35 mm, 30 mm, 25 mm, 20 mm） in finite element models are compared and found no considerable difference. Another possible compensation method, Broyden-Fletcher-Shanmo method, is also presented based on the solving nonlinear fimctions idea. The Broyden-Fletcher-Shanmo method is employed to compute mold contour for the second panel. It only takes 50% iterations compared to that of DA. The proposed method can serve a faster mold contour compensation method for sheet metal forming.