Application of GNSS-PPP on Dynamic Deformation Monitoring of Offshore Platforms

The real-time dynamic deformation monitoring of offshore platforms under environmental excitation is crucial to their safe operation.Although Global Navigation Satellite System-Precise Point Positioning(GNSS-PPP)has been considered for this purpose,its monitoring accuracy is relatively low.Moreover,the influence of background noise on the dynamic monitoring accuracy of GNSS-PPP remains unclear.Hence,it is imperative to further validate the feasibility of GNSS-PPP for deformation monitoring of offshore platforms.To address these concerns,vibration table tests with different amplitudes and frequencies are conducted.The results demonstrate that GNSS-PPP can effectively monitor horizontal vibration displacement as low as±30 mm,which is consistent with GNSS-RTK.Furthermore,the spectral characteristic of background noise in GNSS-PPP is similar to that of GNSS-RTK(Real Time Kinematic).Building on this observation,an improved Complete Ensemble Empirical Mode Decomposition with Adaptive Noise(CEEMDAN)has been proposed to de-noise the data and enhance the dynamic monitoring accuracy of GNSS-PPP.Field monitoring application research is also undertaken,successfully extracting and analyzing the dynamic deformation of an offshore platform structure under environmental excitation using GNSS-PPP monitoring in conjunction with improved CEEMDAN de-noising.By comparing the de-noised dynamic deformation trajectories of the offshore platform during different periods,it is observed that the platform exhibits reversible alternating vibration responses under environmental excitation,with more pronounced displacement deformation in the direction of load action.The research results confirm the feasibility and potential of GNSS-PPP for dynamic deformation monitoring of offshore platforms.

Flow characteristics and hot workability of a typical low-alloy high-strength steel during multi-pass deformation

Heavy components of low-alloy high-strength(LAHS) steels are generally formed by multi-pass forging. It is necessary to explore the flow characteristics and hot workability of LAHS steels during the multi-pass forging process, which is beneficial to the formulation of actual processing parameters. In the study, the multi-pass hot compression experiments of a typical LAHS steel are carried out at a wide range of deformation temperatures and strain rates. It is found that the work hardening rate of the experimental material depends on deformation parameters and deformation passes, which is ascribed to the impacts of static and dynamic softening behaviors. A new model is established to describe the flow characteristics at various deformation passes. Compared to the classical Arrhenius model and modified Zerilli and Armstrong model, the newly proposed model shows higher prediction accuracy with a confidence level of 0.98565. Furthermore, the connection between power dissipation efficiency(PDE) and deformation parameters is revealed by analyzing the microstructures. The PDE cannot be utilized to reflect the efficiency of energy dissipation for microstructure evolution during the entire deformation process, but only to assess the efficiency of energy dissipation for microstructure evolution in a specific deformation parameter state.As a result, an integrated processing map is proposed to better study the hot workability of the LAHS steel, which considers the effects of instability factor(IF), PDE, and distribution and size of grains. The optimized processing parameters for the multi-pass deformation process are the deformation parameters of 1223–1318 K and 0.01–0.08 s^(-1). Complete dynamic recrystallization occurs within the optimized processing parameters with an average grain size of 18.36–42.3 μm. This study will guide the optimization of the forging process of heavy components.

Liquefaction susceptibility and deformation characteristics of saturated coral sandy soils subjected to cyclic loadings-a critical review

Coral sandy soils widely exist in coral island reefs and seashores in tropical and subtropical regions.Due to the unique marine depositional environment of coral sandy soils,the engineering characteristics and responses of these soils subjected to monotonic and cyclic loadings have been a subject of intense interest among the geotechnical and earthquake engineering communities.This paper critically reviews the progress of experimental investigations on the undrained behavior of coral sandy soils under monotonic and cyclic loadings over the last three decades.The focus of coverage includes the contractive-dilative behavior,the pattern of excess pore-water pressure(EPWP)generation and the liquefaction mechanism and liquefaction resistance,the small-strain shear modulus and strain-dependent shear modulus and damping,the cyclic softening feature,and the anisotropic characteristics of undrained responses of saturated coral sandy soils.In particular,the advances made in the past decades are reviewed from the following aspects:(1)the characterization of factors that impact the mechanism and patterns of EPWP build-up;(2)the identification of liquefaction triggering in terms of the apparent viscosity and the average flow coefficient;(3)the establishment of the invariable form of strain-based,stress-based,or energy-based EPWP ratio formulas and the unique relationship between the new proxy of liquefaction resistance and the number of cycles required to reach liquefaction;(4)the establishment of the invariable form of the predictive formulas of small strain modulus and strain-dependent shear modulus;and(5)the investigation on the effects of stress-induced anisotropy on liquefaction susceptibility and dynamic deformation characteristics.Insights gained through the critical review of these advances in the past decades offer a perspective for future research to further resolve the fundamental issues concerning the liquefaction mechanism and responses of coral sandy sites subjected to cyclic loadings associated with seismic events in marine environments.

Coseismic deformation and fault slip distribution of the 2023 M_(W)7.8 and M_(W)7.6 earthquakes in Türkiye

On February 6,2023,a devastating earthquake with a moment magnitude of M_(W)7.8 struck the town of Pazarcik in south-central Türkiye,followed by another powerful earthquake with a moment magnitude of M_(W)7.6 that struck the nearby city of Elbistan 9 h later.To study the characteristics of surface deformation caused by this event and the influence of fault rupture,this study calculated the static coseismic deformation of 56 stations and dynamic displacement waveforms of 15 stations using data from the Turkish national fixed global navigation satellite system(GNSS)network.A maximum static coseismic displacement of 0.38 m for the M_(W)7.8 Kahramanmaras earthquake was observed at station ANTE,36 km from the epicenter,and a maximum dynamic coseismic displacement of 4.4 m for the M_(W)7.6 Elbistan earthquake was observed at station EKZ1,5 km from the epicenter.The rupture-slip distributions of the two earthquakes were inverted using GNSS coseismic deformation as a constraint.The results showed that the Kahramanmaras earthquake rupture segment was distinct and exposed on the ground,resulting in significant rupture slip along the Amanos and Pazarcik fault segments of the East Anatolian Fault.The maximum slip in the Pazarcik fault segment was 10.7 m,and rupture occurred at depths of 0–15 km.In the Cardak fault region,the Elbistan earthquake caused significant ruptures at depths of 0–12 km,with the largest amount of slip reaching 11.6 m.The Coulomb stress change caused by the Kahramanmaras earthquake rupture along the Cardak fault segment was approximately 2 bars,and the area of increased Coulomb stress corresponded to the subsequent rupture region of the M_(W)7.6 earthquake.Thus,it is likely that the M_(W)7.8 earthquake triggered or promoted the M_(W)7.6 earthquake.Based on the cumulative stress impact of the M_(W)7.8 and M_(W)7.6 events,the southwestern segment of the East Anatolian Fault,specifically the Amanos fault segment,experienced a Coulomb rupture stress change exceeding 2 bars,warranting further attention to assess its future seismic hazard risk.

Deformation monitoring of long-span railway bridges based on SBAS-InSAR technology

The deformation monitoring of long-span railway bridges is significant to ensure the safety of human life and property.The interferometric synthetic aperture radar(In SAR)technology has the advantage of high accuracy in bridge deformation monitoring.This study monitored the deformation of the Ganjiang Super Bridge based on the small baseline subsets(SBAS)In SAR technology and Sentinel-1A data.We analyzed the deformation results combined with bridge structure,temperature,and riverbed sediment scouring.The results are as follows:(1)The Ganjiang Super Bridge area is stable overall,with deformation rates ranging from-15.6 mm/yr to 10.7 mm/yr(2)The settlement of the Ganjiang Super Bridge deck gradually increases from the bridge tower toward the main span,which conforms to the typical deformation pattern of a cable-stayed bridge.(3)The sediment scouring from the riverbed cause the serious settlement on the bridge’s east side compared with that on the west side.(4)The bridge deformation negatively correlates with temperature,with a faster settlement at a higher temperature and a slow rebound trend at a lower temperature.The study findings can provide scientific data support for the health monitoring of long-span railway bridges.

Violent collisions can reveal hexadecapole deformation of nuclei

A new observable in heavy ion collision experiments was identified to be sensitive to the hexadecapole deformation of the colliding nuclei.Such deformation is difficult to measure in traditional nuclear electric transition measurements,as it is often overwhelmed by the nuclear quadrupole deformation.This opens the door to gain new insight into nuclear structure with experiments that were designed to study hot and dense nuclear matter.

Effect of deformation parameters on the austenite dynamic recrystallization behavior of a eutectoid pearlite rail steel

Understandings of the effect of hot deformation parameters close to the practical production line on grain refinement are crucial for enhancing both the strength and toughness of future rail steels.In this work,the austenite dynamic recrystallization(DRX)behaviors of a eutectoid pearlite rail steel were studied using a thermo-mechanical simulator with hot deformation parameters frequently employed in rail production lines.The single-pass hot deformation results reveal that the prior austenite grain sizes(PAGSs)for samples with different deformation reductions decrease initially with an increase in deformation temperature.However,once the deformation temperature is beyond a certain threshold,the PAGSs start to increase.It can be attributed to the rise in DRX volume fraction and the increase of DRX grain with deformation temperature,respectively.Three-pass hot deformation results show that the accumulated strain generated in the first and second deformation passes can increase the extent of DRX.In the case of complete DRX,PAGS is predominantly determined by the deformation temperature of the final pass.It suggests a strategic approach during industrial production where part of the deformation reduction in low temperature range can be shifted to the medium temperature range to release rolling mill loads.

Study on the vertical deformations induced by terrestrial water storage changes in Huang-Huai-Hai river basin

Terrestrial water storage(TWs)variations are associated with water mass movements,which may cause the deformation displacements of the Global Navigation Satellite System(GNSS)stations.This study investigates the spatio-temporal Tws variations and addresses the relationship between deformation variations observed in the Huang-Huai-Hai River Basin(HHHRB)and local hydrological features.Results indicate that the vertical velocities at the GNSS stations induced by TWS changes are relatively small,and the impacts of the terrestrial water storage changes are mainly reflected in the changes of seasonal characteristics.Although there is a downward TWS trend from 2011 to 2022 in most HHHRB areas,velocities from the vertical displacements of both Gravity Recovery and Climate Experiment(GRACE)and GRACE Follow-On(GFO)and the GNSS reflect that the HHHRB is undergoing an uplift process,while the magnitude of the GRACE/GFO derived velocities is much smaller than that of the GNSS solutions.Common hydrological deformations estimated from GRACE/GFO and GNSS measurements reveal that the TWS-derived displacements can explain 54.5%of the GNSS seasonal variations,with the phases of terrestrial water storage advancing by about one month relative to GNss common signal phases.Moreover,the decrease of the groundwater storage in the HHHRB has been accelerating since 2008.After reaching its lowest level around mid-2020,it began to rise rapidly,which might be closely related to the implementation of the South-North Water Transfer Central Project.

Room and cryogenic deformation behavior of AZ61 and AZ61-xCaO(x=0.5,1 wt.%)alloy

This study investigates the influence of CaO(0.5,1(wt.%))alloying on the microstructural evolution,texture development and deformation behavior of AZ61 magnesium alloy.The uniaxial tension tests at room(RT)and cryogenic(CT,-150℃)temperature were performed to investigate the twinability and dislocation behavior and its consequent effect on flow stress,ductility and strain hardening rate.The results showed that the AZ61-1CaO exhibited superior strength/ductility synergy at RT with a yield strength(YS)of 223 MPa and a ductility of 23% as compared to AZ61(178 MPa,18.5%)and AZ61-0.5CaO(198 MPa,21%).Similar trend was witnessed for all the samples during CT deformation,where increase in the YS and decrease in ductility were observed.The Mtex tools based in-grain misorientation axis(IGMA)analysis of RT deformed samples revealed the higher activities of prismatic slip in AZ61-CaO,which led to superior ductility.Moreover,subsequent EBSD analysis of CT deformed samples showed the increased fraction of fine{10-12}tension twins and nucleation of multiple{10-12}twin variants caused by higher local stress concentration at the grain boundaries,which imposed the strengthening by twin-twin interaction.Lastly,the detailed investigations on strengthening contributors showed that the dislocation strengthening has the highest contribution towards strength in all samples.

Mechanical Behavior and Microstructure Evolution during Tensile Deformation of Twinning Induced Plasticity Steel Processed by Warm Forgings

The mechanical behavior and microstructural evolution of an Fe-30Mn-3Al-3Si twinninginduced plasticity(TWIP)steel processed using warm forging was investigated.It is found that steel processed via warm forging improves comprehensive mechanical properties compared to the TWIP steel processed via cold rolling,with a high tensile strength(R_(m))of 793 MPa,a yield strength(R_(P))of 682 MPa,an extremely large R_(P)/R_(m)ratio as high as 0.86 as well as an excellent elongation rate of 46.8%.The microstructure observation demonstrates that steel processed by warm forging consists of large and elongated grains together with fine,equiaxed grains.Complicated micro-defect configurations were also observed within the steel,including dense dislocation networks and a few coarse deformation twins.As the plastic deformation proceeds,the densities of dislocations and deformation twins significantly increase.Moreover,a great number of slip lines could be observed in the elongated grains.These findings reveal that a much more dramatic interaction between microstructural defect and dislocations glide takes place in the forging sample,wherein the fine and equiaxed grains propagated dislocations more rapidly,together with initial defect configurations,are responsible for enhanced strength properties.Meanwhile,larger,elongated grains with more prevalently activated deformation twins result in high plasticity.

Unraveling engineering disturbance effects on deformation in red-bed mudstone railway cuttings:incorporating crack-facilitated moisture diffusion

Red-bed mudstone, prevalent in southwest China, poses a formidable challenge due to its hydrophilic clay minerals, resulting in expansion, deformation, and cracking upon exposure to moisture. This study addresses uplift deformation disasters in high-speed railways by employing a moisture diffusion-deformation-fracture coupling model based on the finite-discrete element method(FDEM). The model integrates the influence of cracks on moisture diffusion. The investigation into various excavation depths reveals a direct correlation between depth and the formation of tensile cracks at the bottom of the railway cutting. These cracks expedite moisture migration, significantly impacting the temporal and spatial evolution of the moisture field. Additionally, crack expansion dominates hygroscopic deformation, with the lateral coordinate of the crack zone determining peak vertical displacement. Furthermore, key factors influencing deformation in railway cuttings, including the swelling factor and initial moisture content at the bottom of the cutting, are explored. The number of tensile and shear cracks increases with greater excavation depth, particularly concerning shear cracks. Higher swelling factors and initial moisture contents result in an increased total number of cracks, predominantly shear cracks. Numerical calculations provide valuable insights, offering a scientific foundation and directional guidance for the precise prevention, control, prediction, and comprehensive treatment of mudstone-related issues in high-speed railways.

Kinematic deformation and intensity assessment of the 2021 Maduo M_(S)7.4 earthquake in Qinghai revealed by high-frequency GNSS

Rapid acquisition of the kinematic deformation field and seismic intensity distribution of large earthquakes is crucial for postseismic emergency rescue,disaster assessment,and future seismic risk research.The advancement of GNSS observation and data processing makes it play an important role in this field,especially the high-frequency GNSS.We used the differential positioning method to calculate the 1 HZ GNSS data from 98 sites within 1000 km of the M_(S)7.4 Maduo earthquake epicenter.The kinematic deformation field and the distribution of the seismic intensity by using the peak ground velocity derived from displacement waveforms were obtained.The results show that:1)Horizontal coseismic response deformation levels ranging from 25 mm to 301 mm can be observed within a 1000 km radius from the epicenter.Coseismic response deformation on the east and west sides shows bilateral asymmetry,which markedly differs from the symmetry presented by surface rupture.2)The seismic intensity obtained through high-frequency GNSS and field investigations exhibits good consistency of the scope and orientation in the high seismic intensity area,although the former is generally slightly smaller than the latter.3)There may exist obstacles on the eastern side of the seismogenic fault.The Maduo earthquake induced a certain tectonic stress loading effect on the western Kunlun Pass-Jiangcuo fault(KPJF)and Maqin-Maqu segment,resulting in higher seismic risk in the future.

Operando measurement of lattice deformation profiles of synchrotron radiation monochromator

This study presents a new method for characterizing the thermal lattice deformation of a monochromator with high precision under service conditions and first reports the operando measurements of nanoscale thermal lattice deformation on a double-crystal monochromator at different incident powers.The nanoscale thermal lattice deformation of the monochromator first crystal was obtained by analyzing the intensity of the distorted DuMond diagrams.DuMond diagrams of the 333 diffraction index,sensitive to lattice deformation,were obtained directly using a 2D detector and an analyzer crystal orthogonal to the monochromator.With increasing incident power and power density,the maximum height of the lattice deformation increased from 3.2 to 18.5 nm,and the deformation coefficient of the maximum height increased from 1.1 to 3.2 nm/W.The maximum relative standard deviation was 4.2%,and the maximum standard deviation was 0.1 nm.Based on the measured thermal deformations,the flux saturation phenomenon and critical point for the linear operation of the monochromator were predicted with increasing incident power.This study provides a simple solution to the problem of the lower precision of synchrotron radiation monochromator characterizations compared to simulations.

Finite Deformation, Finite Strain Nonlinear Dynamics and Dynamic Bifurcation in TVE Solids with Rheology

This paper presents a mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and ordered rate constitutive theories in Lagrangian description derived using entropy inequality and the representation theorem for thermoviscoelastic solids (TVES) with rheology. The CBL and the constitutive theories take into account finite deformation and finite strain deformation physics and are based on contravariant deviatoric second Piola-Kirchhoff stress tensor and its work conjugate covariant Green’s strain tensor and their material derivatives of up to order m and n respectively. All published works on nonlinear dynamics of TVES with rheology are mostly based on phenomenological mathematical models. In rare instances, some aspects of CBL are used but are incorrectly altered to obtain mass, stiffness and damping matrices using space-time decoupled approaches. In the work presented in this paper, we show that this is not possible using CBL of CCM for TVES with rheology. Thus, the mathematical models used currently in the published works are not the correct description of the physics of nonlinear dynamics of TVES with rheology. The mathematical model used in the present work is strictly based on the CBL of CCM and is thermodynamically and mathematically consistent and the space-time coupled finite element methodology used in this work is unconditionally stable and provides solutions with desired accuracy and is ideally suited for nonlinear dynamics of TVES with memory. The work in this paper is the first presentation of a mathematical model strictly based on CBL of CCM and the solution of the mathematical model is obtained using unconditionally stable space-time coupled computational methodology that provides control over the errors in the evolution. Both space-time coupled and space-time decoupled finite element formulations are considered for obtaining solutions of the IVPs described by the mathematical model and are presented in the paper. Factors or the physics influencing dynamic response and dynamic bifurcation for TVES with rheology are identified and are also demonstrated through model problem studies. A simple model problem consisting of a rod (1D) of TVES material with memory fixed at one end and subjected to harmonic excitation at the other end is considered to study nonlinear dynamics of TVES with rheology, frequency response as well as dynamic bifurcation phenomenon.

Impact of initial fluctuations and nuclear deformations in isobar collisions

Relativistic isobar^(96)_(44)Ru+^(96)_(44)Ru and^(96)_(40)Zr+^(96)_(40)Zrcollisions have revealed intricate differences in their nuclear size and shape,inspiring unconventional studies of nuclear structure using relativistic heavy ion collisions.In this study,we investigate the relative differences in the mean multiplicityR_(<Nch>)and the secondR_(ε2)and third-order eccentricityR_(ε3)between isobar collisions using initial state Glauber models.It is found that initial fluctuations and nuclear deformations have negligible effects on R_(<Nch>)in most central collisions,while both are important for the R_(ε2)and R_(ε3),the degree of which is sensitive to the underlying nucleonic or sub-nucleonic degree of freedom.These features,compared to real data,may probe the particle production mechanism and the physics underlying nuclear structure.

Improved spatio-temporal alignment measurement method for hull deformation

In this paper,an improved spatio-temporal alignment measurement method is presented to address the inertial matching measurement of hull deformation under the coexistence of time delay and large misalignment angle.Large misalignment angle and time delay often occur simultaneously and bring great challenges to the accurate measurement of hull deformation in space and time.The proposed method utilizes coarse alignment with large misalignment angle and time delay estimation of inertial measurement unit modeling to establish a brand-new spatiotemporal aligned hull deformation measurement model.In addition,two-step loop control is designed to ensure the accurate description of dynamic deformation angle and static deformation angle by the time-space alignment method of hull deformation.The experiments illustrate that the proposed method can effectively measure the hull deformation angle when time delay and large misalignment angle coexist.

Research on casing deformation prevention technology based on cementing slurry system optimization

The casing deformation prevention technology based on the optimization of cement slurry is proposed to reduce the casing deformation of shale oil and gas wells during hydraulic fracturing. In this paper, the fracture mechanism of hollow particles in cement sheath was firstly analyzed by discrete element method, and the effect of hollow particles in cement on casing deformation was investigated by laboratory experiment method. Finally, field test was carried out to verify the improvement effect of the casing deformation based on cement slurry modification. The results show that the formation displacement can be absorbed effectively by hollow particles inside the cement transferring the excessive deformation away from casing. The particles in the uncemented state provide deformation space during formation slipping. The casing with diameter of 139.7 mm could be passed through by bridge plug with the diameter of 99 mm when the mass ratio of particle/cement reaches 1:4. According to the field test feedback, the method based on optimization of cement slurry can effectively reduce the risk of casing deformation, and the recommended range of hollow microbeads content in the cement slurry is between 15% and 25%.

Swelling pressure evolution characterization of strong expansive soil considering the influence of reserved expansion deformation

Numerous engineering cases have demonstrated that the expansive soil channel slope remains susceptible to damage with the implementation of a rigid or closed protective structure. It is common for the protective structure to experience bulging failure due to excessive swelling pressure. To investigate the swelling pressure properties of expansive soil, the constant volume test was employed to study the influence of water content and reserved expansion deformation on the characteristics of swelling pressure in strong expansive soils, and also to explore the evolution mechanism of the swelling pressure. The findings demonstrate that the swelling pressure-time curve can be classified into swelling pressure-time softening and swelling pressure-time stability type. The swelling pressuretime curve of the specimen with low water content is the swelling pressure-time softening type, and the softening level will be weakened with increasing reserved expansion deformation. Besides, the maximum swelling pressure Psmax decreases with increasing water content and reserved expansion deformation, especially for expansion ratio η from 24% to 37%. The reserved deformation has little effect on reducing Psmax when it is beyond 7% of the expansion rate. The specimen with low water content has a more homogeneous structure due to the significant expansion-filling effect, and the fracture and reorganization of the aggregates in the specimens with low water content cause the swelling pressure-time softening behavior. In addition, the proposed swelling pressure-time curve prediction model has a good prediction on the test results. If necessary, a deformation space of about 7% expansion rate is recommended to be reserved in the engineering to reduce the swelling pressure except for keeping a stable water content.

Deformation Characteristics of Hydrate-Bearing Sediments

The safe and efficient development of natural gas hydrate requires a deep understanding of the deformation behaviors of reservoirs.In this study,a series of triaxial shearing tests are carried out to investigate the deformation properties of hydrate-bearing sediments.Variations of volumetric and lateral strains versus hydrate saturation are analyzed comprehensively.Results indicate that the sediments with high hydrate saturation show dilative behaviors,which lead to strain-softening characteristics during shearing.The volumetric strain curves have a tendency to transform gradually from dilatation to compression with the increase in effective confining pressure.An easy prediction model is proposed to describe the relationship between volumetric and axial strains.The model coefficientβis the key dominating factor for the shape of volumetric strain curves and can be determined by the hydrate saturation and stress state.Moreover,a modified model is established for the calculation of lateral strain.The corresponding determination method is provided for the easy estimation of model coefficients for medium sand sediments containing hydrate.This study provides a theoretical and experimental reference for deformation estimation in natural gas hydrate development.

An improved strain-softening constitutive model of granite considering the effect of crack deformation

This paper presents an improved strain-softening constitutive model considering the effect of crack deformation based on the triaxial cyclic loading and unloading test results.The improved model assumes that total strain is a combination of plastic,elastic,and crack strains.The constitutive relationship between the crack strain and the stress was further derived.The evolutions of mechanical parameters,i.e.strength parameters,dilation angle,unloading elastic modulus,and deformation parameters of crack,with the plastic strain and confining pressure were studied.With the increase in plastic strain,the cohesion,friction angle,dilation angle,and crack Poisson's ratio initially increase and subsequently decrease,and the unloading elastic modulus and the crack elastic modulus nonlinearly decrease.The increasing confining pressure enhances the strength and unloading elastic modulus,and decreases the dilation angle and Poisson's ratio of the crack.The theoretical triaxial compressive stress-strain curves were compared with the experimental results,and they present a good agreement with each other.The improved constitutive model can well reflect the nonlinear mechanical behavior of granite.