Dynamical Modeling and Dynamic Characteristics Analysisof a Coaxial Dual-Rotor System

The dual-rotor structure serves as the primary source of vibration in aero-engines. Understanding itsdynamical model and analyzing dynamic characteristics, such as critical speed and unbalanced response, arecrucial for rotor system dynamics. Previous work introduced a coaxial dual-rotor-support scheme for aeroengines,and a physical model featuring a high-speed flexible inner rotor with a substantial length-to-diameter ratiowas designed. Then a finite element (FE) dynamic model based on the Timoshenko beam elements and rigid bodykinematics of the dual-rotor system is modeled, with the Newmark method and Newton–Raphson method used forthe numerical calculation to study the dynamic characteristics of the system. Three different simulation models,including beam-based FE (1D) model, solid-based FE (3D) model, and transfer matrix model, were designed tostudy the characteristics of mode and the critical speed characteristic of the dual-rotor system. The unbalancedresponse of the dual-rotor system was analyzed to study the influence of mass unbalance on the rotor system. Theeffect of different disk unbalance phases and different speed ratios on the dynamic characteristics of the dual-rotorsystem was investigated in detail. The experimental result shows that the beam-based FE model is effective andsuitable for studying the dual-rotor system.

Time-domain dynamic constitutive model suitable for mucky soil site seismic response

Soil nonlinear behavior displays noticeable effects on the site seismic response.This study proposes a new functional expression of the skeleton curve to replace the hyperbolic skeleton curve.By integrating shear modulus and combining the dynamic skeleton curve and the damping degradation coefficient,the constitutive equation of the logarithmic dynamic skeleton can be obtained,which considers the damping effect in a soil dynamics problem.Based on the finite difference method and the multi-transmitting boundary condition,a 1D site seismic response analysis program called Soilresp1D has been developed herein and used to analyze the time-domain seismic response in three types of sites.At the same time,this study also provides numerical simulation results based on the hyperbolic constitutive model and the equivalent linear method.The results verify the rationality of the new soil dynamic constitutive model.It can analyze the mucky soil site nonlinear seismic response,reflecting the deformation characteristics and damping effect of the silty soil.The hysteresis loop area is more extensive,and the residual strain is evident.

Dynamically triggered seismicity on a tectonic scale:A review

Earthquakes triggered by dynamic disturbances have been confirmed by numerous observations and experiments.In the past several decades,earthquake triggering has attracted increasing attention of scholars in relation to exploring the mechanism of earthquake triggering,earthquake prediction,and the desire to use the mechanism of earthquake triggering to reduce,prevent,or trigger earthquakes.Natural earthquakes and large‐scale explosions are the most common sources of dynamic disturbances that trigger earthquakes.In the past several decades,some models have been developed,including static,dynamic,quasi‐static,and other models.Some reviews have been published,but explosiontriggered seismicity was not included.In recent years,some new results on earthquake triggering have emerged.Therefore,this paper presents a new review to reflect the new results and include the content of explosion‐triggered earthquakes for the reference of scholars in this area.Instead of a complete review of the relevant literature,this paper primarily focuses on the main aspects of dynamic earthquake triggering on a tectonic scale and makes some suggestions on issues that need to be resolved in this area in the future.

Hybrid Dynamic Variables-Dependent Event-Triggered Fuzzy Model Predictive Control

This article focuses on dynamic event-triggered mechanism(DETM)-based model predictive control(MPC) for T-S fuzzy systems.A hybrid dynamic variables-dependent DETM is carefully devised,which includes a multiplicative dynamic variable and an additive dynamic variable.The addressed DETM-based fuzzy MPC issue is described as a “min-max” optimization problem(OP).To facilitate the co-design of the MPC controller and the weighting matrix of the DETM,an auxiliary OP is proposed based on a new Lyapunov function and a new robust positive invariant(RPI) set that contain the membership functions and the hybrid dynamic variables.A dynamic event-triggered fuzzy MPC algorithm is developed accordingly,whose recursive feasibility is analysed by employing the RPI set.With the designed controller,the involved fuzzy system is ensured to be asymptotically stable.Two examples show that the new DETM and DETM-based MPC algorithm have the advantages of reducing resource consumption while yielding the anticipated performance.

Adaptive Optimal Discrete-Time Output-Feedback Using an Internal Model Principle and Adaptive Dynamic Programming

In order to address the output feedback issue for linear discrete-time systems, this work suggests a brand-new adaptive dynamic programming(ADP) technique based on the internal model principle(IMP). The proposed method, termed as IMP-ADP, does not require complete state feedback-merely the measurement of input and output data. More specifically, based on the IMP, the output control problem can first be converted into a stabilization problem. We then design an observer to reproduce the full state of the system by measuring the inputs and outputs. Moreover, this technique includes both a policy iteration algorithm and a value iteration algorithm to determine the optimal feedback gain without using a dynamic system model. It is important that with this concept one does not need to solve the regulator equation. Finally, this control method was tested on an inverter system of grid-connected LCLs to demonstrate that the proposed method provides the desired performance in terms of both tracking and disturbance rejection.

Dynamic Modeling and Experimental Verification of an RPR Type Compliant Paralle Mechanism with Low Orders

Efficiency of calculating a dynamic response is an important point of the compliant mechanism for posture adjustment.Dynamic modeling with low orders of a 2R1T compliant parallel mechanism is studied in the paper.The mechanism with two out-of-plane rotational and one lifting degrees of freedom(DoFs)plays an important role in posture adjustment.Based on elastic beam theory,the stiffness matrix and mass matrix of the beam element are established where the moment of inertia is considered.To improve solving efficiency,a dynamic model with low orders of the mechanism is established based on a modified modal synthesis method.Firstly,each branch of the RPR type mechanism is divided into a substructure.Subsequently,a set of hypothetical modes of each substructure is obtained based on the C-B method.Finally,dynamic equation of the whole mechanism is established by the substructure assembly.A dynamic experiment is conducted to verify the dynamic characteristics of the compliant mechanism.

Dynamics Modeling and Parameter Identification for a Coupled-Drive Dual-Arm Nursing Robot

A dual-arm nursing robot can gently lift patients and transfer them between a bed and a wheelchair.With its lightweight design,high load-bearing capacity,and smooth surface,the coupled-drive joint is particularly well suited for these robots.However,the coupled nature of the joint disrupts the direct linear relationship between the input and output torques,posing challenges for dynamic modeling and practical applications.This study investigated the transmission mechanism of this joint and employed the Lagrangian method to construct a dynamic model of its internal dynamics.Building on this foundation,the Newton-Euler method was used to develop a dynamic model for the entire robotic arm.A continuously differentiable friction model was incorporated to reduce the vibrations caused by speed transitions to zero.An experimental method was designed to compensate for gravity,inertia,and modeling errors to identify the parameters of the friction model.This method establishes a mapping relationship between the friction force and motor current.In addition,a Fourier series-based excitation trajectory was developed to facilitate the identification of the dynamic model parameters of the robotic arm.Trajectory tracking experiments were conducted during the experimental validation phase,demonstrating the high accuracy of the dynamic model and the parameter identification method for the robotic arm.This study presents a dynamic modeling and parameter identification method for coupled-drive joint robotic arms,thereby establishing a foundation for motion control in humanoid nursing robots.

Phase diagram and quench dynamics of a periodically drivenHaldane model

We investigate a periodically driven Haldane model subjected to a two-stage driving scheme in the form of a step function.By using the Floquet theory,we obtain the topological phase diagram of the system.We also find that anomalous Floquet topological phases exist in the system.Focusing on examining the quench dynamics among topological phases,we analyze the site distribution of the 0-mode and p-mode edge states in long-period evolution after a quench.The results demonstrate that,under certain conditions,the site distribution of the 0-mode can be confined at the edge even in long-period evolution.Additionally,both the 0-mode and p-mode can recover and become confined at the edge in long-period evolution when the post-quench parameters(T,M_(2) /M_(1))in the phase diagram cross away from the phase boundary (M_(2)/ M_(1))=(6√3t2)/ M_(1)−1.Furthermore,we conclude that whether the edge state is confined at the edge in the long-period evolution after a quench depends on the similarity of the edge states before and after the quench.Our findings reveal some new characteristics of quench dynamics in a periodically driven system.

A semi-analytical pressure and rate transient analysis model for inner boundary and propped fractures exhibiting dynamic behavior under long-term production conditions

The loss of hydrocarbon production caused by the dynamic behavior of the inner boundary and propped fractures under long-term production conditions has been widely reported in recent studies.However,the quantitative relationships for the variations of the inner boundary and propped fractures have not been determined and incorporated in the semi-analytical models for the pressure and rate transient analysis.This work focuses on describing the variations of the inner boundary and propped fractures and capturing the typical characteristics from the pressure transient curves.A generalized semi-analytical model was developed to characterize the dynamic behavior of the inner boundary and propped fractures under long-term production conditions.The pressure-dependent length shrinkage coefficients,which quantify the length changes of the inner zone and propped fractures,are modified and incorporated into this multi-zone semi-analytical model.With simultaneous numerical iterations and numerical inversions in Laplace and real-time space,the transient solutions to pressure and rate behavior are determined in just a few seconds.The dynamic behavior of the inner boundary and propped fractures on transient pressure curves is divided into five periods:fracture bilinear flow(FR1),dynamic PFs flow(FR2),inner-area linear flow(FR3),dynamic inner boundary flow(FR4),and outer-area dominated linear flow(FR5).The early hump during FR2 period and a positive upward shift during FR4period are captured on the log-log pressure transient curves,reflecting the dynamic behavior of the inner boundary and propped fractures during the long-term production period.The transient pressure behavior will exhibit greater positive upward trend and the flow rate will be lower with the shrinkage of the inner boundary.The pressure derivative curve will be upward earlier as the inner boundary shrinks more rapidly.The lower permeability caused by the closure of un-propped fractures in the inner zone results in greater upward in pressure derivative curves.If the permeability loss for the dynamic behavior of the inner boundary caused by the closure of un-propped fractures is neglected,the flow rate will be overestimated in the later production period.

Characterization and identification towards dynamic-based electrical modeling of lithium-ion batteries

Lithium-ion batteries are widely recognized as a crucial enabling technology for the advancement of electric vehicles and energy storage systems in the grid.The design of battery state estimation and control algorithms in battery management systems is usually based on battery models,which interpret crucial battery dynamics through the utilization of mathematical functions.Therefore,the investigation of battery dynamics with the purpose of battery system identification has garnered considerable attention in the realm of battery research.Characterization methods in terms of linear and nonlinear response of lithium-ion batteries have emerged as a prominent area of study in this field.This review has undertaken an analysis and discussion of characterization methods,with a particular focus on the motivation of battery system identification.Specifically,this work encompasses the incorporation of frequency domain nonlinear characterization methods and dynamics-based battery electrical models.The aim of this study is to establish a connection between the characterization and identification of battery systems for researchers and engineers specialized in the field of batteries,with the intention of promoting the advancement of efficient battery technology for real-world applications.

A Neural-network-based Alternative Scheme to Include Nonhydrostatic Processes in an Atmospheric Dynamical Core

Here,a nonhydrostatic alternative scheme(NAS)is proposed for the grey zone where the nonhydrostatic impact on the atmosphere is evident but not large enough to justify the necessity to include an implicit nonhydrostatic solver in an atmospheric dynamical core.The NAS is designed to replace this solver,which can be incorporated into any hydrostatic models so that existing well-developed hydrostatic models can effectively serve for a longer time.Recent advances in machine learning(ML)provide a potential tool for capturing the main complicated nonlinear-nonhydrostatic relationship.In this study,an ML approach called a neural network(NN)was adopted to select leading input features and develop the NAS.The NNs were trained and evaluated with 12-day simulation results of dry baroclinic-wave tests by the Weather Research and Forecasting(WRF)model.The forward time difference of the nonhydrostatic tendency was used as the target variable,and the five selected features were the nonhydrostatic tendency at the last time step,and four hydrostatic variables at the current step including geopotential height,pressure in two different forms,and potential temperature,respectively.Finally,a practical NAS was developed with these features and trained layer by layer at a 20-km horizontal resolution,which can accurately reproduce the temporal variation and vertical distribution of the nonhydrostatic tendency.Corrected by the NN-based NAS,the improved hydrostatic solver at different horizontal resolutions can run stably for at least one month and effectively reduce most of the nonhydrostatic errors in terms of system bias,anomaly root-mean-square error,and the error of the wave spatial pattern,which proves the feasibility and superiority of this scheme.

Dynamic response of mountain tunnel,bridge,and embankment along the Sichuan-Tibet transportation corridor to active fault based on model tests

The Sichuan-Tibet transportation corridor is prone to numerous active faults and frequent strong earthquakes.While extensive studies have individually explored the effect of active faults and strong earthquakes on different engineering structures,their combined effect remains unclear.This research employed multiple physical model tests to investigate the dynamic response of various engineering structures,including tunnels,bridges,and embankments,under the simultaneous influence of cumulative earthquakes and stick-slip misalignment of an active fault.The prototype selected for this study was the Kanding No.2 tunnel,which crosses the Yunongxi fault zone within the Sichuan-Tibet transportation corridor.The results demonstrated that the tunnel,bridge,and embankment exhibited amplification in response to the input seismic wave,with the amplification effect gradually decreasing as the input peak ground acceleration(PGA)increased.The PGAs of different engineering structures were weakened by the fault rupture zone.Nevertheless,the misalignment of the active fault may decrease the overall stiffness of the engineering structure,leading to more severe damage,with a small contribution from seismic vibration.Additionally,the seismic vibration effect might be enlarged with the height of the engineering structure,and the tunnel is supposed to have a smaller PGA and lower dynamic earth pressure compared to bridges and embankments in strong earthquake zones crossing active faults.The findings contribute valuable insights for evaluating the dynamic response of various engineering structures crossing an active fault and provide an experimental reference for secure engineering design in the challenging conditions of the Sichuan-Tibet transportation corridor.

Dynamics modeling and optimal control for multi-information diffusion in Social Internet of Things

As an ingenious convergence between the Internet of Things and social networks,the Social Internet of Things(SIoT)can provide effective and intelligent information services and has become one of the main platforms for people to spread and share information.Nevertheless,SIoT is characterized by high openness and autonomy,multiple kinds of information can spread rapidly,freely and cooperatively in SIoT,which makes it challenging to accurately reveal the characteristics of the information diffusion process and effectively control its diffusion.To this end,with the aim of exploring multi-information cooperative diffusion processes in SIoT,we first develop a dynamics model for multi-information cooperative diffusion based on the system dynamics theory in this paper.Subsequently,the characteristics and laws of the dynamical evolution process of multi-information cooperative diffusion are theoretically investigated,and the diffusion trend is predicted.On this basis,to further control the multi-information cooperative diffusion process efficiently,we propose two control strategies for information diffusion with control objectives,develop an optimal control system for the multi-information cooperative diffusion process,and propose the corresponding optimal control method.The optimal solution distribution of the control strategy satisfying the control system constraints and the control budget constraints is solved using the optimal control theory.Finally,extensive simulation experiments based on real dataset from Twitter validate the correctness and effectiveness of the proposed model,strategy and method.

Influence of Various Earth-Retaining Walls on the Dynamic Response Comparison Based on 3D Modeling

The present work aims to assess earthquake-induced earth-retaining(ER)wall displacement.This study is on the dynamics analysis of various earth-retaining wall designs in hollow precast concrete panels,reinforcement concrete facing panels,and gravity-type earth-retaining walls.The finite element(FE)simulations utilized a 3D plane strain condition to model full-scale ER walls and numerous nonlinear dynamics analyses.The seismic performance of differentmodels,which includes reinforcement concrete panels and gravity-type and hollowprecast concrete ER walls,was simulated and examined using the FE approach.It also displays comparative studies such as stress distribution,deflection of the wall,acceleration across the wall height,lateral wall displacement,lateral wall pressure,and backfill plastic strain.Three components of the created ER walls were found throughout this research procedure.One is a granular reinforcement backfill,while the other is a wall-facing panel and base foundation.The dynamic response effects of varied earth-retaining walls have also been studied.It was discovered that the facing panel of the model significantly impacts the earthquake-induced displacement of ER walls.The proposed analytical model’s validity has been evaluated and compared with the reinforcement concrete facing panels,gravity-type ER wall,scientifically available data,and American Association of State Highway and Transportation Officials(AASHTO)guidelines results based on FE simulation.The results of the observations indicate that the hollow prefabricated concrete ER wall is the most feasible option due to its lower displacement and high-stress distribution compared to the two types.The methodology and results of this study establish standards for future analogous investigations and professionals,particularly in light of the increasing computational capabilities of desktop computers.

Computational Fluid Dynamics Approach for Predicting Pipeline Response to Various Blast Scenarios: A Numerical Modeling Study

Recent industrial explosions globally have intensified the focus in mechanical engineering on designing infras-tructure systems and networks capable of withstanding blast loading.Initially centered on high-profile facilities such as embassies and petrochemical plants,this concern now extends to a wider array of infrastructures and facilities.Engineers and scholars increasingly prioritize structural safety against explosions,particularly to prevent disproportionate collapse and damage to nearby structures.Urbanization has further amplified the reliance on oil and gas pipelines,making them vital for urban life and prime targets for terrorist activities.Consequently,there is a growing imperative for computational engineering solutions to tackle blast loading on pipelines and mitigate associated risks to avert disasters.In this study,an empty pipe model was successfully validated under contact blast conditions using Abaqus software,a powerful tool in mechanical engineering for simulating blast effects on buried pipelines.Employing a Eulerian-Lagrangian computational fluid dynamics approach,the investigation extended to above-surface and below-surface blasts at standoff distances of 25 and 50 mm.Material descriptions in the numerical model relied on Abaqus’default mechanical models.Comparative analysis revealed varying pipe performance,with deformation decreasing as explosion-to-pipe distance increased.The explosion’s location relative to the pipe surface notably influenced deformation levels,a key finding highlighted in the study.Moreover,quantitative findings indicated varying ratios of plastic dissipation energy(PDE)for different blast scenarios compared to the contact blast(P0).Specifically,P1(25 mm subsurface blast)and P2(50 mm subsurface blast)showed approximately 24.07%and 14.77%of P0’s PDE,respectively,while P3(25 mm above-surface blast)and P4(50 mm above-surface blast)exhibited lower PDE values,accounting for about 18.08%and 9.67%of P0’s PDE,respectively.Utilising energy-absorbing materials such as thin coatings of ultra-high-strength concrete,metallic foams,carbon fiber-reinforced polymer wraps,and others on the pipeline to effectively mitigate blast damage is recommended.This research contributes to the advancement of mechanical engineering by providing insights and solutions crucial for enhancing the resilience and safety of underground pipelines in the face of blast events.

Dynamic optimal allocation of energy storage systems integrated within photovoltaic based on a dual timescale dynamics model

Energy storage systems(ESSs)operate as independent market participants and collaborate with photovoltaic(PV)generation units to enhance the flexible power supply capabilities of PV units.However,the dynamic variations in the profitability of ESSs in the electricity market are yet to be fully understood.This study introduces a dual-timescale dynamics model that integrates a spot market clearing(SMC)model into a system dynamics(SD)model to investigate the profit-aware capacity growth of ESSs and compares the profitability of independent energy storage systems(IESSs)with that of an ESS integrated within a PV(PV-ESS).Furthermore,this study aims to ascertain the optimal allocation of the PV-ESS.First,SD and SMC models were set up.Second,the SMC model simulated on an hourly timescale was incorporated into the SD model as a subsystem,a dual-timescale model was constructed.Finally,a development simulation and profitability analysis was conducted from 2022 to 2040 to reveal the dynamic optimal range of PV-ESS allocation.Additionally,negative electricity prices were considered during clearing processes.The simulation results revealed differences in profitability and capacity growth between IESS and PV-ESS,helping grid investors and policymakers to determine the boundaries of ESSs and dynamic optimal allocation of PV-ESSs.

A novel approach to the dynamic response analysis of Euler-Bernoulli beams resting on a Winkler soil model and subjected to impact loads

This work presents a novel approach to the dynamic response analysis of a Euler-Bernoulli beam resting on a Winkler soil model and subjected to an impact loading.The approach considers that damping has much less importance in controlling the maximum response to impulsive loadings because the maximum response is reached in a very short time,before the damping forces can dissipate a significant portion of the energy input into the system.The development of two sine series solutions,relating to different types of impulsive loadings,one involving a single concentrated force and the other a distributed line load,are presented.This study revealed that when a simply supported Euler-Bernoulli beam,resting on a Winkler soil model,is subject to an impact load,the resulting vertical displacements,bending moments and shear forces produced along the span of the beam are considerably affected.In particular,the quantification of this effect is best observed,relative to the corresponding static solution,via an amplification factor.The computed impact amplification factors,for the sub-grade moduli used in this study,were in magnitude greater than 2,thus confirming the multiple-degree-of-freedom nature of the problem.

Improving the accuracy of precipitation estimates in a typical inland arid area of China using a dynamic Bayesian model averaging approach

Xinjiang Uygur Autonomous Region is a typical inland arid area in China with a sparse and uneven distribution of meteorological stations,limited access to precipitation data,and significant water scarcity.Evaluating and integrating precipitation datasets from different sources to accurately characterize precipitation patterns has become a challenge to provide more accurate and alternative precipitation information for the region,which can even improve the performance of hydrological modelling.This study evaluated the applicability of widely used five satellite-based precipitation products(Climate Hazards Group InfraRed Precipitation with Station(CHIRPS),China Meteorological Forcing Dataset(CMFD),Climate Prediction Center morphing method(CMORPH),Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Climate Data Record(PERSIANN-CDR),and Tropical Rainfall Measuring Mission Multi-satellite Precipitation Analysis(TMPA))and a reanalysis precipitation dataset(ECMWF Reanalysis v5-Land Dataset(ERA5-Land))in Xinjiang using ground-based observational precipitation data from a limited number of meteorological stations.Based on this assessment,we proposed a framework that integrated different precipitation datasets with varying spatial resolutions using a dynamic Bayesian model averaging(DBMA)approach,the expectation-maximization method,and the ordinary Kriging interpolation method.The daily precipitation data merged using the DBMA approach exhibited distinct spatiotemporal variability,with an outstanding performance,as indicated by low root mean square error(RMSE=1.40 mm/d)and high Person's correlation coefficient(CC=0.67).Compared with the traditional simple model averaging(SMA)and individual product data,although the DBMA-fused precipitation data were slightly lower than the best precipitation product(CMFD),the overall performance of DBMA was more robust.The error analysis between DBMA-fused precipitation dataset and the more advanced Integrated Multi-satellite Retrievals for Global Precipitation Measurement Final(IMERG-F)precipitation product,as well as hydrological simulations in the Ebinur Lake Basin,further demonstrated the superior performance of DBMA-fused precipitation dataset in the entire Xinjiang region.The proposed framework for solving the fusion problem of multi-source precipitation data with different spatial resolutions is feasible for application in inland arid areas,and aids in obtaining more accurate regional hydrological information and improving regional water resources management capabilities and meteorological research in these regions.

Nonlinear dynamics of a circular curved cantilevered pipe conveying pulsating fluid based on the geometrically exact model

Due to the novel applications of flexible pipes conveying fluid in the field of soft robotics and biomedicine,the investigations on the mechanical responses of the pipes have attracted considerable attention.The fluid-structure interaction(FSI)between the pipe with a curved shape and the time-varying internal fluid flow brings a great challenge to the revelation of the dynamical behaviors of flexible pipes,especially when the pipe is highly flexible and usually undergoes large deformations.In this work,the geometrically exact model(GEM)for a curved cantilevered pipe conveying pulsating fluid is developed based on the extended Hamilton's principle.The stability of the curved pipe with three different subtended angles is examined with the consideration of steady fluid flow.Specific attention is concentrated on the large-deformation resonance of circular pipes conveying pulsating fluid,which is often encountered in practical engineering.By constructing bifurcation diagrams,oscillating shapes,phase portraits,time traces,and Poincarémaps,the dynamic responses of the curved pipe under various system parameters are revealed.The mean flow velocity of the pulsating fluid is chosen to be either subcritical or supercritical.The numerical results show that the curved pipe conveying pulsating fluid can exhibit rich dynamical behaviors,including periodic and quasi-periodic motions.It is also found that the preferred instability type of a cantilevered curved pipe conveying steady fluid is mainly in the flutter of the second mode.For a moderate value of the mass ratio,however,a third-mode flutter may occur,which is quite different from that of a straight pipe system.

Anisotropic dynamic permeability model for porous media

Based on the tortuous capillary network model,the relationship between anisotropic permeability and rock normal strain,namely the anisotropic dynamic permeability model(ADPM),was derived and established.The model was verified using pore-scale flow simulation.The uniaxial strain process was calculated and the main factors affecting permeability changes in different directions in the deformation process were analyzed.In the process of uniaxial strain during the exploitation of layered oil and gas reservoirs,the effect of effective surface porosity on the permeability in all directions is consistent.With the decrease of effective surface porosity,the sensitivity of permeability to strain increases.The sensitivity of the permeability perpendicular to the direction of compression to the strain decreases with the increase of the tortuosity,while the sensitivity of the permeability in the direction of compression to the strain increases with the increase of the tortuosity.For layered reservoirs with the same initial tortuosity in all directions,the tortuosity plays a decisive role in the relative relationship between the variations of permeability in all directions during pressure drop.When the tortuosity is less than 1.6,the decrease rate of horizontal permeability is higher than that of vertical permeability,while the opposite is true when the tortuosity is greater than 1.6.This phenomenon cannot be represented by traditional dynamic permeability model.After the verification by experimental data of pore-scale simulation,the new model has high fitting accuracy and can effectively characterize the effects of deformation in different directions on the permeability in all directions.