Pulses in ground motions identified through surface partial matching and their impact on seismic rocking consequence

In seismology and earthquake engineering,it is fundamental to identify and characterize the pulse-like features in pulse-type ground motions.To capture the pulses that dominate structural responses,this study establishes congruence and shift relationships between response spectrum surfaces.A similarity search between spectrum surfaces,supplemented with a similarity search in time series,has been applied to characterize the pulse-like features in pulse-type ground motions.The identified pulses are tested in predicting the rocking consequences of slender rectangular blocks under the original ground motions.Generally,the prediction is promising for the majority of the ground motions where the dominant pulse is correctly identified.

Failure mechanisms and destruction characteristics of cemented coal gangue backfill under compression effect of non-uniform load

Backfill mining is one of the most important technical means for controlling strata movement and reducing surface subsidence and environmental damage during exploitation of underground coal resources. Ensuring the stability of the backfill bodies is the primary prerequisite for maintaining the safety of the backfilling working face, and the loading characteristics of backfill are closely related to the deformation and subsidence of the roof. Elastic thin plate model was used to explore the non-uniform subsidence law of the roof, and then the non-uniform distribution characteristics of backfill bodies’ load were revealed. Through a self-developed non-uniform loading device combined with acoustic emission (AE) and digital image correlation (DIC) monitoring technology, the synergistic dynamic evolution law of the bearing capacity, apparent crack, and internal fracture of cemented coal gangue backfills (CCGBs) under loads with different degrees of non-uniformity was deeply explored. The results showed that: 1) The uniaxial compressive strength (UCS) of CCGB increased and then decreased with an increase in the degree of non-uniformity of load (DNL). About 40% of DNL was the inflection point of DNL-UCS curve and when DNL exceeded 40%, the strength decreased in a cliff-like manner;2) A positive correlation was observed between the AE ringing count and UCS during the loading process of the specimen, which was manifested by a higher AE ringing count of the high-strength specimen. 3) Shear cracks gradually increased and failure mode of specimens gradually changed from “X” type dominated by tension cracks to inverted “Y” type dominated by shear cracks with an increase in DNL, and the crack opening displacement at the peak stress decreased and then increased. The crack opening displacement at 40% of the DNL was the smallest. This was consistent with the judgment of crack size based on the AE b-value, i. e., it showed the typical characteristics of “small b-value-large crack and large b-value-small crack”. The research results are of significance for preventing the instability and failure of backfill.

Consolidation of high replacement ratio stone column-reinforced ground:Analytical solutions incorporating clogging effect

The utilization of stone columns has emerged as a popular ground improvement strategy,whereas the drainage performance can be adversely hampered by clogging effect.Despite the ample progress of calculation methods for the consolidation of stone column-improved ground,theoretical investigations into the clogging effect have not been thoroughly explored.Furthermore,it is imperative to involve the column consolidation deformation to mitigate computational error on the consolidation of composite ground with high replacement ratios.In this context,an analytical model accounting for the initial clogging and coupled time and depth-dependent clogging of stone columns is established.Then,the resulting governing equations and analytical solutions are obtained under a new flow continuity relationship to incorporate column consolidation deformation.The accuracy and reliability of the proposed model are illustrated by degradation analysis and case studies with good agreements.Subsequently,the computed results of the current study are juxtaposed against the existing models,and an in-depth assessment of the impacts of several crucial parameters on the consolidation behavior is conducted.The results reveal that ignoring column consolidation deformation leads to an overestimate of the consolidation rate,with maximum error reaching up to 16%as the replacement ratio increases.Furthermore,the initial clogging also has a significant influence on the consolidation performance.Additionally,the increment of depth and time-clogging factors a and b will induce a noticeable retardation of the consolidation process,particularly in the later stage.

Study of the pressure transient behavior of directional wells considering the effect of non-uniform flux distribution

During the production,the fluid in the vicinity of the directional well enters the wellbore with different rates,leading to non-uniform flux distribution along the directional well.However,in all existing studies,it is oversimplified to a uniform flux distribution,which can result in inaccurate results for field applications.Therefore,this paper proposes a semi-analytical model of a directional well based on the assumption of non-uniform flux distribution.Specifically,the direction well is discretized into a carefully chosen series of linear sources,such that the complex well trajectory can be captured and the nonuniform flux distribution along the wellbore can be considered to model the three-dimensional flow behavior.By using the finite difference method,we can obtain the numerical solutions of the transient flow within the wellbore.With the aid of Green's function method,we can obtain the analytical solutions of the transient flow from the matrix to the wellbore.The complete flow behavior of a directional well is perfectly represented by coupling the above two types of transient flow.Subsequently,on the basis of the proposed model,we conduct a comprehensive analysis of the pressure transient behavior of a directional well.The computation results show that the flux variation along the direction well has a significant effect on pressure responses.In addition,the directional well in an infinite reservoir may exhibit the following flow regimes:wellbore afterflow,transition flow,inclined radial flow,elliptical flow,horizontal linear flow,and horizontal radial flow.The horizontal linear flow can be observed only if the formation thickness is much smaller than the well length.Furthermore,a dip region that appears on the pressure derivative curve indicates the three-dimensional flow behavior near the wellbore.

Effect of non-uniform swelling on coal multiphysics during gas injection: The triangle approach

In current dual porosity/permeability models,there exists a fundamental assumption that the adsorption-induced swelling is distributed uniformly within the representative elementary volume (REV),irrespective of its internal structures and transient processes.However,both internal structures and transient processes can lead to the non-uniform swelling.In this study,we hypothesize that the non-uniform swelling is responsible for why coal permeability in experimental measurements is not only controlled by the effective stress but also is affected by the adsorption-induced swelling.We propose a concept of the swelling triangle composed of swelling paths to characterize the evolution of the non-uniform swelling and serve as a core link in coupled multiphysics.A swelling path is determined by a dimensionless volumetric ratio and a dimensionless swelling ratio.Different swelling paths have the same start and end point,and each swelling path represents a unique swelling case.The swelling path as the diagonal of the triangle represents the case of the uniform swelling while that as the two perpendicular boundaries represents the case of the localized swelling.The paths of all intermediate cases populate inside the triangle.The corresponding relations between the swelling path and the response of coal multiphysics are established by a non-uniform swelling coefficient.We define this method as the triangle approach and corresponding models as swelling path-based ones.The proposed concept and models are verified against a long-term experimental measurement of permeability and strains under constant effective stress.Our results demonstrate that during gas injection,coal multiphysics responses have a close dependence on the swelling path,and that in both future experiments and field predictions,this dependence must be considered.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

An improved non-uniform fast Fourier transform method for radio imaging of coronal mass ejections

Radioheliographs can obtain solar images at high temporal and spatial resolution,with a high dynamic range.These are among the most important instruments for studying solar radio bursts,understanding solar eruption events,and conducting space weather forecasting.This study aims to explore the effective use of radioheliographs for solar observations,specifically for imaging coronal mass ejections(CME),to track their evolution and provide space weather warnings.We have developed an imaging simulation program based on the principle of aperture synthesis imaging,covering the entire data processing flow from antenna configuration to dirty map generation.For grid processing,we propose an improved non-uniform fast Fourier transform(NUFFT)method to provide superior image quality.Using simulated imaging of radio coronal mass ejections,we provide practical recommendations for the performance of radioheliographs.This study provides important support for the validation and calibration of radioheliograph data processing,and is expected to profoundly enhance our understanding of solar activities.

An efficient physics-guided Bayesian framework for predicting ground settlement profile during excavations in clay

Recently,the application of Bayesian updating to predict excavation-induced deformation has proven successful and improved prediction accuracy significantly.However,updating the ground settlement profile,which is crucial for determining potential damage to nearby infrastructures,has received limited attention.To address this,this paper proposes a physics-guided simplified model combined with a Bayesian updating framework to accurately predict the ground settlement profile.The advantage of this model is that it eliminates the need for complex finite element modeling and makes the updating framework user-friendly.Furthermore,the model is physically interpretable,which can provide valuable references for construction adjustments.The effectiveness of the proposed method is demonstrated through two field case studies,showing that it can yield satisfactory predictions for the settlement profile.

Thermal performance of cast-in-place piles with artificial ground freezing in permafrost regions

During the construction of cast-in-place piles in warm permafrost,the heat carried by concrete and the cement hydration reaction can cause strong thermal disturbance to the surrounding permafrost.Since the bearing capacity of the pile is quite small before the full freeze-back,the quick refreezing of the native soils surrounding the cast-in-place pile has become the focus of the infrastructure construction in permafrost.To solve this problem,this paper innovatively puts forward the application of the artificial ground freezing(AGF)method at the end of the curing period of cast-in-place piles in permafrost.A field test on the AGF was conducted at the Beiluhe Observation and Research Station of Frozen Soil Engineering and Environment(34°51.2'N,92°56.4'E)in the Qinghai Tibet Plateau(QTP),and then a 3-D numerical model was established to investigate the thermal performance of piles using AGF under different engineering conditions.Additionally,the long-term thermal performance of piles after the completion of AGF under different conditions was estimated.Field experiment results demonstrate that AGF is an effective method to reduce the refreezing time of the soil surrounding the piles constructed in permafrost terrain,with the ability to reduce the pile-soil interface temperatures to below the natural ground temperature within 3 days.Numerical results further prove that AGF still has a good cooling effect even under unfavorable engineering conditions such as high pouring temperature,large pile diameter,and large pile length.Consequently,the application of this method is meaningful to save the subsequent latency time and solve the problem of thermal disturbance in pile construction in permafrost.The research results are highly relevant for the spread of AGF technology and the rapid building of pile foundations in permafrost.

Ground threat prediction-based path planning of unmanned autonomous helicopter using hybrid enhanced artificial bee colony algorithm

Unmanned autonomous helicopter(UAH)path planning problem is an important component of the UAH mission planning system.Aiming to reduce the influence of non-complete ground threat information on UAH path planning,a ground threat prediction-based path planning method is proposed based on artificial bee colony(ABC)algorithm by collaborative thinking strategy.Firstly,a dynamic threat distribution probability model is developed based on the characteristics of typical ground threats.The dynamic no-fly zone of the UAH is simulated and established by calculating the distribution probability of ground threats in real time.Then,a dynamic path planning method for UAH is designed in complex environment based on the real-time prediction of ground threats.By adding the collision warning mechanism to the path planning model,the flight path could be dynamically adjusted according to changing no-fly zones.Furthermore,a hybrid enhanced ABC algorithm is proposed based on collaborative thinking strategy.The proposed algorithm applies the leader-member thinking mechanism to guide the direction of population evolution,and reduces the negative impact of local optimal solutions caused by collaborative learning update strategy,which makes the optimization performance of ABC algorithm more controllable and efficient.Finally,simulation results verify the feasibility and effectiveness of the proposed ground threat prediction path planning method.

Ground settlement and tunnel response due to twin-curved shield tunnelling in soft ground with small clear distance

Twin curved tunnels are often encountered in shield tunnelling,where significant complexities in densely exploited underground space are observed.In this study,the ground settlement and tunnel deformation due to twin-curved shield tunnelling in soft ground were investigated using numerical simulation and field monitoring.Different curvature radii of twin curved tunnels and subsequent effects of tunnel construction were considered to reveal the tunnelling effect on ground surface settlement and tunnel deformation.The results show that the settlement trough yields one offset towards inside of curved shield tunnelling.The location of settlement trough and maximum settlement were affected by curvature radius but except for the shape and width of settlement trough.Adjacent parallel twin-curved shield tunnelling could increase the offset of existing settlement trough and maximum settlement.Then,an empirical prediction of surface settlement trough due to twin-curved shield tunnelling with same tunnel diameters in soft clay was proposed,which was applicable to curvature radius less than 800 m.Finally,a minimum radius of 600 m of curvature tunnel was proposed in terms of allowable convergence deformation of tunnel.The result could provide guidance on safety evaluation for twin curved shield tunnelling construction.

Artificial ground freezing of underground mines in cold regions using thermosyphons with air insulation

Current practice of underground artificial ground freezing(AGF)typically involves huge refrigeration systems of large economic and environmental costs.In this study,a novel AGF technique is proposed deploying available cold wind in cold regions.This is achieved by a static heat transfer device called thermosyphon equipped with an air insulation layer.A refrigeration unit can be optionally integrated to meet additional cooling requirements.The introduction of air insulation isolates the thermosyphon from ground zones where freezing is not needed,resulting in:(1)steering the cooling resources(cold wind or refrigeration)towards zones of interest;and(2)minimizing refrigeration load.This design is demonstrated using well-validated mathematical models from our previous work based on two-phase enthalpy method of the ground coupled with a thermal resistance network for the thermosyphon.Two Canadian mines are considered:the Cigar Lake Mine and the Giant Mine.The results show that our proposed design can speed the freezing time by 30%at the Giant Mine and by two months at the Cigar Lake Mine.Further,a cooling load of 2.4 GWh can be saved at the Cigar Lake Mine.Overall,this study provides mining practitioners with sustainable solutions of underground AGF.

Neural network study of the nuclear ground-state spin distribution within a random interaction ensemble

The distribution of the nuclear ground-state spin in a two-body random ensemble(TBRE)was studied using a general classification neural network(NN)model with two-body interaction matrix elements as input features and the corresponding ground-state spins as labels or output predictions.The quantum many-body system problem exceeds the capability of our optimized NNs in terms of accurately predicting the ground-state spin of each sample within the TBRE.However,our NN model effectively captured the statistical properties of the ground-state spin because it learned the empirical regularity of the ground-state spin distribution in TBRE,as discovered by physicists.

Crack mechanism of ground fissures in loess layer of Fenwei Basin, China

The Fenwei Basin, covered by loess, experiences severe ground fissure disasters. These disasters disrupt the continuity of the loess and pose significant threats to engineering construction safety along transportation routes. Nevertheless, the crack characteristics and the influence zone of ground fissures in the loess layer remain inadequately investigated. To effectively prevent and control ground fissure disasters, physical model tests and the PFC(particle flow code) numerical simulation method are used to investigate the crack mechanism of buried ground fissures in the loess layer. The results show that there are two main cracks in the layer profile, which have a Y-shape morphology. As the dip angle of the preset cracks increased from 60° to 90°, the main deformation zone at the surface gradually shifted towards the footwall. The process of crack propagation from depth to surface is divided into five stages. Additionally, the results confirm the accuracy of the width of the rupture zone d2in the footwall calculated by the cantilever beam theory. These findings can offer theoretical guidance for determining the avoidance distance of ground fissures in loess regions, as well as for implementing disaster prevention and corresponding control measures for various stages of buried ground fissure propagation.

Seismic response and correlation analysis of a pile-supported wharf to near-fault pulse-like ground motions

Earthquake investigations have shown that near-fault pulse-like(NF-P)ground motions have unique characteristics compared to near-fault non-pulse-like(NF-NP)and far-field(FF)ground motions.It is necessary to study the seismic response of pile-supported wharf(PSW)structures under NF-P ground motions.In this study,a three-dimensional finite element numerical model is created to simulate a PSW.By imparting three types of ground motion,the engineering demand parameters(EDPs)of PSW under NF-P ground motions were analyzed and compared,in which EDPs are the maximum displacement and bending moment of the piles.Twenty intensity measures(IMs)were selected to characterize the properties of ground motions.The correlation between IMs and EDPs was explored.The results show that the piles present larger displacement and bending moment under NF-P ground motions compared to NF-NP and FF ground motions.None of the IMs have a high correlation with EDPs under NF-P ground motions,and these IMs are more applicable to FF ground motions.The correlation coefficients between EDPs and IMs under three types of ground motion were obtained,which will provide a valuable reference for the seismic design of PSWs.

Structural failure analysis with CMS-based ground motion selection using innovative cost function and weight factors

The selection and scaling of ground motion records is considered a primary and essential task in performing structural analysis and design.Conventional methods involve using ground motion models and a conditional spectrum to select ground motion records based on the target spectrum.This research demonstrates the influence of adopting different weighted factors for various period ranges during matching selected ground motions with the target hazard spectrum.The event data from the Next Generation Attenuation West 2(NGA-West 2)database is used as the basis for ground motion selection,and hazard de-aggregation is conducted to estimate the event parameters of interest,which are then used to construct the target intensity measure(IM).The target IMs are then used to select ground motion records with different weighted vector-valued objective functions.The weights are altered to account for the relative importance of IM in accordance with the structural analysis application of steel moment resisting frame(SMRF)buildings.Instead of an ordinary objective function for the matching spectrum,a novel model is introduced and compared with the conventional cost function.The results indicate that when applying the new cost function for ground motion selection,it places higher demands on structures compared to the conventional cost function.Moreover,submitting more weights to the first-mode period of structures increases engineering demand parameters.Findings demonstrate that weight factors allocated to different period ranges can successfully account for period elongation and higher mode effects.

Characteristics of Lightning Activity in Southeast China and its Relation to the Atmospheric Background

Based on the lightning observation data from the Fengyun-4A(FY-4A)Lightning Mapping Imager(FY-4A/LMI)and the Lightning Imaging Sensor(LIS)on the International Space Station(ISS),we extract the“event”type data as the lightning detection results.These observations are then compared with the cloud-to-ground(CG)lightning observation data from the China Meteorological Administration.This study focuses on the characteristics of lightning activity in Southeast China,primarily in Jiangxi Province and its adjacent areas,from April to September,2017–2022.In addition,with the fifth-generation European Centre for Medium-Range Weather Forecasts reanalysis data,we further delved into the potential factors influencing the distribution and variations in lightning activity and their primary related factors.Our findings indicate that the lightning frequency and density of the FY-4A/LMI,ISS-LIS and CG data are higher in southern and central Jiangxi,central Fujian Province,and western and central Guangdong Province,while they tend to be lower in eastern Hunan Province.In general,the high-value areas of lightning density for the FY-4A/LMI are located in inland mountainous areas.The lower the latitude is,the higher the CG lightning density is.High-value areas of the CG lightning density are more likely to be located in eastern Fujian and southeastern Zhejiang Province.However,the high-value areas of lightning density for the ISS-LIS are more dispersed,with a scattered distribution in inland mountainous areas and along the coast of eastern Fujian.Thus,the mountainous terrain is closely related to the high-value areas of the lightning density.The locations of the high-value areas of the lightning density for the FY-4A/LMI correspond well with those for the CG observations,and the seasonal variations are also consistent.In contrast,the distribution of the high-value areas of the lightning density for the ISS-LIS is more dispersed.The positions of the peak frequency of the FY-4A/LMI lightning and CG lightning contrast with local altitudes,primarily located at lower altitudes or near mountainsides.K-index and convective available potential energy(CAPE)can better reflect the local boundary layer conditions,where the lightning density is higher and lightning seasonal variations are apparent.There are strong correlations in the annual variations between the dew-point temperature(Td)and CG lightning frequency,and the monthly variations of the dew-point temperature and CAPE are also strongly correlated with monthly variations of CG lightning,while they are weakly correlated with the lightning frequency for the FY-4A/LMI and ISS-LIS.This result reflects that the CAPE shows a remarkable effect on the CG lightning frequency during seasonal transitions.

Physical and numerical investigations of target stratum selection for ground hydraulic fracturing of multiple hard roofs

Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based on engineering properties to simulate the gradual collapse of the roof during longwall top coal caving(LTCC).A numerical model is established using the material point method(MPM)and the strain-softening damage constitutive model according to the structure of the physical model.Numerical simulations are conducted to analyze the LTCC process under different hard roofs for ground hydraulic fracturing.The results show that ground hydraulic fracturing releases the energy and stress of the target stratum,resulting in a substantial lag in the fracturing of the overburden before collapse occurs in the hydraulic fracturing stratum.Ground hydraulic fracturing of a low hard roof reduces the lag effect of hydraulic fractures,dissipates the energy consumed by the fracture of the hard roof,and reduces the abutment stress.Therefore,it is advisable to prioritize the selection of the lower hard roof as the target stratum.

Comparative Analysis of Energy Characteristics of Two Southwest Vortices in Sichuan Under Similar Circulation Backgrounds

Based on ERA5 reanalysis data,the present study analyzed the thermal energy development mechanism and kinetic energy conversion characteristics of two extreme rainstorm processes in relation to the shallow southwest vortex in the warm-sector during a“rain-generated vortex”process and the deep southwest vortex in a“vortex-generated rain”process.The findings were as follows:(1)During the extreme rainstorm on August 11,2020(hereinafter referred to as the“8·11”process),intense surface heating and a high-energy unstable environment were observed.The mesoscale convergence system triggered convection to produce heavy rainfall,and the release of latent condensation heat generated by the rainfall promoted the formation of a southwest vortex.The significant increase(decrease)in atmospheric diabatic heating and kinetic energy preceded the increase(decrease)in vorticity.By contrast,the extreme rainstorm on August 16,2020(hereinafter referred to as the“8·16”process)involved the generation of southwest vortex in a low-energy and highhumidity environment.The dynamic uplift of the southwest vortex triggered rainfall,and the release of condensation latent heat from rainfall further strengthened the development of the southwest vortex.The significant increase(decrease)in atmospheric diabatic heating and kinetic energy exhibited a delayed progression compared to the increase(decrease)in vorticity.(2)The heating effect around the southwest vortex region was non-uniform,and the heating intensity varied in different stages.In the“8·11”process,the heating effect was the strongest in the initial stage,but weakened during the vortex's development.On the contrary,the heating effect was initially weak in the“8·16”process,and intensified during the development stage.(3)The available potential energy of the“8·11”process significantly increased in kinetic energy converted from rotational and divergent winds through baroclinic action,and the divergent wind energy continued to convert into rotational wind energy.By contrast,the“8·16”process involved the conversion of rotational wind energy into divergent wind energy,which in turn converted kinetic energy back into available potential energy,thereby impeding the further development and maintenance of the southwest vortex.

The deterministic condition for the ground reaction force acting point on the combined knee valgus and tibial internal rotation moments in early phase of cutting maneuvers in female athletes

Background:Combined knee valgus and tibial internal rotation(VL+IR)moments have been shown to stress the anterior cruciate ligament(ACL)in several in vitro cadaveric studies.To utilize this knowledge for non-contact ACL injury prevention in sports,it is necessary to elucidate how the ground reaction force(GRF)acting point(center of pressure(CoP))in the stance foot produces combined knee VL+IR moments in risky maneuvers,such as cuttings.However,the effects of the GRF acting point on the development of the combined knee VL+IR moment in cutting are still unknown.Methods:We first established the deterministic mechanical condition that the CoP position relative to the tibial rotational axis differentiates the GRF vector’s directional probability for developing the combined knee VL+IR moment,and theoretically predicted that when the CoP is posterior to the tibial rotational axis,the GRF vector is more likely to produce the combined knee VL+IR moment than when the CoP is anterior to the tibial rotational axis.Then,we tested a stochastic aspect of our theory in a lab-controlled in vivo experiment.Fourteen females performed 60˚cutting under forefoot/rearfoot strike conditions(10 trials each).The positions of lower limb markers and GRF data were measured,and the knee moment due to GRF vector was calculated.The trials were divided into anterior-and posterior-CoP groups depending on the CoP position relative to the tibial rotational axis at each 10 ms interval from 0 to 100 ms after foot strike,and the occurrence rate of the combined knee VL+IR moment was compared between trial groups.Results:The posterior-CoP group showed significantly higher occurrence rates of the combined knee VL+IR moment(maximum of 82.8%)at every time point than those of the anterior-CoP trials,as theoretically predicted by the deterministic mechanical condition.Conclusion:The rearfoot strikes inducing the posterior CoP should be avoided to reduce the risk of non-contact ACL injury associated with the combined knee VL+IR stress.