A novel box-counting method for quantitative fractal analysis of threedimensional pore characteristics in sandstone

Fractal theory offers a powerful tool for the precise description and quantification of the complex pore structures in reservoir rocks,crucial for understanding the storage and migration characteristics of media within these rocks.Faced with the challenge of calculating the three-dimensional fractal dimensions of rock porosity,this study proposes an innovative computational process that directly calculates the three-dimensional fractal dimensions from a geometric perspective.By employing a composite denoising approach that integrates Fourier transform(FT)and wavelet transform(WT),coupled with multimodal pore extraction techniques such as threshold segmentation,top-hat transformation,and membrane enhancement,we successfully crafted accurate digital rock models.The improved box-counting method was then applied to analyze the voxel data of these digital rocks,accurately calculating the fractal dimensions of the rock pore distribution.Further numerical simulations of permeability experiments were conducted to explore the physical correlations between the rock pore fractal dimensions,porosity,and absolute permeability.The results reveal that rocks with higher fractal dimensions exhibit more complex pore connectivity pathways and a wider,more uneven pore distribution,suggesting that the ideal rock samples should possess lower fractal dimensions and higher effective porosity rates to achieve optimal fluid transmission properties.The methodology and conclusions of this study provide new tools and insights for the quantitative analysis of complex pores in rocks and contribute to the exploration of the fractal transport properties of media within rocks.

A symmetric substructuring method for analyzing the natural frequencies of conical origami structures

Conical origami structures are characterized by their substantial out-of-plane stiffness and energy-absorptioncapacity.Previous investigations have commonly focused on the static characteristics of these lightweight struc-tures.However,the efficient analysis of the natural vibrations of these structures is pivotal for designing conicalorigami structures with programmable stiffness and mass.In this paper,we propose a novel method to analyzethe natural vibrations of such structures by combining a symmetric substructuring method(SSM)and a gener-alized eigenvalue analysis.SSM exploits the inherent symmetry of the structure to decompose it into a finiteset of repetitive substructures.In doing so,we reduce the dimensions of matrices and improve computationalefficiency by adopting the stiffness and mass matrices of the substructures in the generalized eigenvalue analysis.Finite element simulations of pin-jointed models are used to validate the computational results of the proposedapproach.Moreover,the parametric analysis of the structures demonstrates the influences of the number of seg-ments along the circumference and the radius of the cone on the structural mass and natural frequencies of thestructures.Furthermore,we present a comparison between six-fold and four-fold conical origami structures anddiscuss the influence of various geometric parameters on their natural frequencies.This study provides a strategyfor efficiently analyzing the natural vibration of symmetric origami structures and has the potential to contributeto the efficient design and customization of origami metastructures with programmable stiffness.

Size-dependent effect on biaxial and shear nonlinear buckling analysis of nonlocal isotropic and orthotropic micro-plate based on surface stress and modified couple stress theories using differential quadrature method

The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.

Quantitative Analysis of Seeing with Height and Time at Muztagh-Ata Site Based on ERA5 Database

Seeing is an important index to evaluate the quality of an astronomical site.To estimate seeing at the Muztagh-Ata site with height and time quantitatively,the European Centre for Medium-Range Weather Forecasts reanalysis database(ERA5)is used.Seeing calculated from ERA5 is compared consistently with the Differential Image Motion Monitor seeing at the height of 12 m.Results show that seeing decays exponentially with height at the Muztagh-Ata site.Seeing decays the fastest in fall in 2021 and most slowly with height in summer.The seeing condition is better in fall than in summer.The median value of seeing at 12 m is 0.89 arcsec,the maximum value is1.21 arcsec in August and the minimum is 0.66 arcsec in October.The median value of seeing at 12 m is 0.72arcsec in the nighttime and 1.08 arcsec in the daytime.Seeing is a combination of annual and about biannual variations with the same phase as temperature and wind speed indicating that seeing variation with time is influenced by temperature and wind speed.The Richardson number Ri is used to analyze the atmospheric stability and the variations of seeing are consistent with Ri between layers.These quantitative results can provide an important reference for a telescopic observation strategy.

Calculation and Analysis of TVMS Considering Profile Shifts and Surface Wear Evolution Process of Spur Gear

Profile shift is a highly effective technique for optimizing the performance of spur gear transmission systems.However,tooth surface wear is inevitable during gear meshing due to inadequate lubrication and long-term operation.Both profile shift and tooth surface wear(TSW)can impact the meshing characteristics by altering the involute tooth profile.In this study,a tooth stiffness model of spur gears that incorporates profile shift,TSW,tooth deformation,tooth contact deformation,fillet-foundation deformation,and gear body structure coupling is established.This model efficiently and accurately determines the time-varying mesh stiffness(TVMS).Additionally,an improved wear depth prediction method for spur gears is developed,which takes into consideration the mutually prime teeth numbers and more accurately reflects actual gear meshing conditions.Results show that consideration of the mutual prime of teeth numbers will have a certain impact on the TSW process.Furthermore,the finite element method(FEM)is employed to accurately verify the values of TVMS and load sharing ratio(LSR)of profile-shifted gears and worn gears.This study quantitatively analyzes the effect of profile shift on the surface wear process,which suggests that gear profile shift can partially alleviate the negative effects of TSW.The contribution of this study provides valuable insights into the design and maintenance of spur gear systems.

Nonlinear dynamic modeling of planar moving Timoshenko beam considering non-rigid non-elastic axial effects

Due to the importance of vibration effects on the functional accuracy of mechanical systems,this research aims to develop a precise model of a nonlinearly vibrating single-link mobile flexible manipulator.The manipulator consists of an elastic arm,a rotary motor,and a rigid carrier,and undergoes general in-plane rigid body motion along with elastic transverse deformation.To accurately model the elastic behavior,Timoshenko’s beam theory is used to describe the flexible arm,which accounts for rotary inertia and shear deformation effects.By applying Newton’s second law,the nonlinear governing equations of motion for the manipulator are derived as a coupled system of ordinary differential equations(ODEs)and partial differential equations(PDEs).Then,the assumed mode method(AMM)is used to solve this nonlinear system of governing equations with appropriate shape functions.The assumed modes can be obtained after solving the characteristic equation of a Timoshenko beam with clamped boundary conditions at one end and an attached mass/inertia at the other.In addition,the effect of the transverse vibration of the inextensible arm on its axial behavior is investigated.Despite the axial rigidity,the effect makes the rigid body dynamics invalid for the axial behavior of the arm.Finally,numerical simulations are conducted to evaluate the performance of the developed model,and the results are compared with those obtained by the finite element approach.The comparison confirms the validity of the proposed dynamic model for the system.According to the mentioned features,this model can be reliable for investigating the system’s vibrational behavior and implementing vibration control algorithms.

Algorithms for Empirical Equations in Terms of the Cosmic Microwave Background Temperature

Previously, we presented several empirical equations using the cosmic microwave background (CMB) temperature. Next, we propose an empirical equation for the fine-structure constant. Considering the compatibility among these empirical equations, the CMB temperature (Tc) and gravitational constant (G) were calculated to be 2.726312 K and 6.673778 × 10−11 m3∙kg−1∙s−2, respectively. Every equation can be explained numerically in terms of the Compton length of an electron (λe), the Compton length of a proton (λp) and α. Furthermore, every equation can also be explained in terms of the Avogadro number and the number of electrons at 1 C. We show that every equation can be described in terms of the Planck constant. Then, the ratio of the gravitational force to the electric force can be uniquely determined with the assumption of minimum mass. In this report, we describe the algorithms used to explain these equations in detail. Thus, there are no dimension mismatch problems.

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

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

Effective stress analysis method of seismic response for high tailings dam

Based on the analysis method for tailings dam in upstream raising method presently used in metallurgy and nonferrous metals tailings depository in the world, an effective stress analysis method of seismic response for high tailings dam was developed according to the results of engineering geological exploration, static and dynamic test and stability analysis on Baizhishan tailing dam 113.5 m high. The law of generation, diffusion and dissipation of seismic pore water pressure during and after earthquake was investigated, and the results of tailings dam’s acceleration, seismic dynamic stress and pore water pressure were obtained. The results show that the seismic stability and liquefaction resistance of high tailings dam are strengthened remarkably, and the scope and depth of liquefaction area at the top of dam are reduced greatly. The interior stress is compressive stress, the stress level of every element is less than 1.0 and the safety coefficient of every element is greater than 1.0. The safety coefficient against liquefaction of every element of tailing dam is greater than 1.5 according to the effective stress analysis of seismic response by finite element method. The calculated results prove that liquefaction is the main reason of seismic failure of high tailing dams, and the effect of seismic inertia forces on high tailing dams’ stability during earthquake is secondary reason.

3D stability analysis method of concave slope based on the Bishop method

In order to study the stability control mechanism of a concave slope with circular landslide, and remove the influence of differences in shape on slope stability, the limit analysis method of a simplified Bishop method was employed. The sliding body was divided into strips in a three-dimensional model, and the lateral earth pressure was put into mechanical analysis and the three-dimensional stability analysis methods applicable for circular sliding in concave slope were deduced. Based on geometric structure and the geological parameters of a concave slope, the influence rule of curvature radius and the top and bottom arch height on the concave slope stability were analyzed. The results show that the stability coefficient decreases after growth, first in the transition stage of slope shape from flat to concave, and it has been confirmed that there is a best size to make the slope stability factor reach a maximum. By contrast with average slope, the stability of a concave slope features a smaller range of ascension with slope height increase, which indicates that the enhancing effect of a concave slope is apparent only with lower slope heights.

Application of the homotopy analysis method to nonlinear characteristics of a piezoelectric semiconductor fiber

Based on the nonlinear constitutive equation,a piezoelectric semiconductor(PSC)fiber under axial loads and Ohmic contact boundary conditions is investigated.The analytical solutions of electromechanical fields are derived by the homotopy analysis method(HAM),indicating that the HAM is efficient for the nonlinear analysis of PSC fibers,along with a rapid rate of convergence.Furthermore,the nonlinear characteristics of electromechanical fields are discussed through numerical results.It is shown that the asymmetrical distribution of electromechanical fields is obvious under a symmetrical load,and the piezoelectric effect is weakened by an applied electric field.With the increase in the initial carrier concentration,the electric potential decreases,and owing to the screen-ing effect of electrons,the distribution of electromechanical fields tends to be symmetrical.

Comparison of Calibration Curve Method and Partial Least Square Method in the Laser Induced Breakdown Spectroscopy Quantitative Analysis

The Laser Induced Breakdown Spectroscopy (LIBS) is a fast, non-contact, no sample preparation analytic technology;it is very suitable for on-line analysis of alloy composition. In the copper smelting industry, analysis and control of the copper alloy concentration affect the quality of the products greatly, so LIBS is an efficient quantitative analysis tech- nology in the copper smelting industry. But for the lead brass, the components of Pb, Al and Ni elements are very low and the atomic emission lines are easily submerged under copper complex characteristic spectral lines because of the matrix effects. So it is difficult to get the online quantitative result of these important elements. In this paper, both the partial least squares (PLS) method and the calibration curve (CC) method are used to quantitatively analyze the laser induced breakdown spectroscopy data which is obtained from the standard lead brass alloy samples. Both the major and trace elements were quantitatively analyzed. By comparing the two results of the different calibration method, some useful results were obtained: both for major and trace elements, the PLS method was better than the CC method in quantitative analysis. And the regression coefficient of PLS method is compared with the original spectral data with background interference to explain the advantage of the PLS method in the LIBS quantitative analysis. Results proved that the PLS method used in laser induced breakdown spectroscopy was suitable for simultaneous quantitative analysis of different content elements in copper smelting industry.

Perturbation Analysis of the Attached Water Mass Effect on Offshore Platform Dynamic Response

Theoretically speaking, it is impossible to make the differential equation of motion uncoupled for the natural modes of a system in consideration of the attached water. The hydro-elastic structure is equal to the non-proportional damping system. In this paper a perturbation analysis method is put forward. The structure motion equation is strictly solved mathematically, and the non-proportional damping problem is transformed into a series of proportional damping ones in the superposition form. The paper also presents the calculation formula of the dynamic response of the structure being subjected to harmonic and arbitrary load. The convergence of the proposed method is also studied in this paper, and the corresponding convergence conditions are given. Finally, the proposed method is used to analyze the displacement response of a real offshore platform. The calculation results show that this method has the characteristics of high accuracy and fast convergence.

BRIDGE-TOUGHENING ANALYSIS OF THE FIBER REINFORCED COMPOSITE CONTAINING INTERPHASE EFFECT

A detailed fracture mechanics analysis of bridge-toughening in a fiber reinforced composite is presented in this paper. The integral equation governing bridge-toughening as well as crack opening displacement (COD) for the composite with interfacial layer is derived from the Castigliano's theorem and interface shear-lag model. A numerical result of the COD equation is obtained using the iteration solution of the second Fredholm integral equation. In order to investigate the effect of various parameters on the toughening, an approximate analytical solution of the equation is present and its error analysis is performed, which demonstrates the approximate solution to be appropriate. A parametric study of the influence of the crack length, interfacial shear modules, thickness of the interphase, fiber radius, fiber volume fraction and properties of materials on composite toughening is therefore carried out. The results are useful for experimental demonstration and toughening design including the fabrication process of the composite.

THE RATIONALISM THEORY AND ITS FINITE ELEMENT ANALYSIS METHOD OF SHELL STRUCTURES

In this paper, a kind of rationalism theory of shell is established which is of different mechanic characters in tension and in compression, and the finite element numerical analysis method is also described.

The Method for Optimum Estimation of COVID-19 Variant Type Virus Infection Status Analysis by the Multivariate Analysis Considering the Environmental Variability Impact in Japan

Currently, the estimated value of the effective reproduction number (ERN), which is an index for grasping the COVID-19 infection status, is used for important planning and evaluation of infection prevention measures. Since ERN in the Sequential SIR model fluctuates in multiple dimensions due to changes in the surrounding environment, it is difficult to set the appropriate accuracy of the uncertainty region of the estimated data. The challenge in this study is to build a mathematical model of infectious disease according to the characteristics and data characteristics of the infectious disease and select an appropriate estimation method. Highly accurate quantitative research that analyzes the validity of “how infectious diseases prevail” from an academic point of view is the key to prediction and estimation in appropriate infection situation analysis. In this study, we adopted a statistical multivariate analysis method (T method) that enables evaluation and prediction of important factors related to ERN estimation and analysis of phenomena that change in real time (time series analysis). It was clarified that it is possible to estimate with higher accuracy by applying the T method to the estimated value of ERN by the current SIR mathematical model.

The Effects of the Structure--Foundation Interaction in the Structural Response of a TBM Gallery

The effects of the interaction between a gallery lining and surrounding ground are evaluated. The gallery is circular with constant thickness and surrounding ground as geotechnical characteristics is soft one. Two types of ground are successively considered with bulk modulus K = 9,260 kPa and K = 4,630 kPa. The analysis is carried out for a gallery with 3.80 m inner diameter and 25 cm thickness. This under construction is part of the irrigation system of the ＂Los Rios＂ county （Ecuador）. Comprehensive computations in various hypotheses pointed out the significant effects of the interaction. The sectional stresses （M, N） in the gallery lining embedded in soft ground depend mainly of the gallery depth versus ground surface. Contrary, the mechanical characteristics of the surrounding ground resulted to have little influence.

Theoretical prediction of effective elastic constants for new intermetallic compound porous material

Based on microstructure analysis of the new Ti-A1 intermetallic compound porous material, a micromechanics model of heterogeneous Plateau porous structure was established and calculation formulas of elastic constants （including effective elastic modulus, effective shear elastic modulus and effective Poisson ratio） were derived by the energy method for this porous material. Calculation results show that both the effective elastic modulus and effective shear elastic modulus increase with the increase of the relative density while the effective Poisson ratio decreases. Compared with the currently-existing hexagonal honeycomb model and micromechanics model of composite materials, the micromechanics model of heterogeneous Plateau porous structure in this study is more suitable for characterizing the medium-density porous material and more accurate for predicting the effective elastic constants of the medium-density porous material. Moreover, the obtained explicit expressions of the effective elastic constants in term of the relative density rather than the microstructural parameters for the uniform and regular Plateau porous structure are more convenient to engineering application.

Predicting the distribution of reservoirs by applying the method of seismic data structure characteristics： Example from the eighth zone in Tahe Oilfield

Seismic data structure characteristics means the waveform character arranged in the time sequence at discrete data points in each 2-D or 3-D seismic trace. Hydrocarbon prediction using seismic data structure characteristics is a new reservoir prediction technique. When the main pay interval is in carbonate fracture and fissure-cavern type reservoirs with very strong inhomogeneity, there are some difficulties with hydrocarbon prediction. Because of the special geological conditions of the eighth zone in the Tahe oil field, we apply seismic data structure characteristics to hydrocarbon prediction for the Ordovician reservoir in this zone. We divide the area oil zone into favorable and unfavorable blocks. Eighteen well locations were proposed in the favorable oil block, drilled, and recovered higher output of oil and gas.

Probabilistic model and analysis of coupled train-ballasted track-subgrade system with uncertain structural parameters

Random dynamic responses caused by the uncertainty of structural parameters of the coupled train-ballasted track-subgrade system under train loading can pose safety concerns to the train operation.This paper introduced a computational model for analyzing probabilistic dynamic responses of three-dimensional(3D)coupled train-ballasted track-subgrade system(TBTSS),where the coupling effects of uncertain rail irregularities,stiffness and damping properties of ballast and subgrade layers were simultaneously considered.The number theoretical method(NTM)was employed to design discrete points for the multi-dimensional stochastic parameters.The time-histories of stochastic dynamic vibrations of the TBSS with systematically uncertain structural parameters were calculated accurately and efficiently by employing the probability density evolution method(PDEM).The model-predicted results were consistent with those by the Monte Carlo simulation method.A sensitivity study was performed to assess the relative importance of those uncertain structural parameters,based on which a case study was presented to explore the stochastic probability evolution mechanism of such train-ballasted track-subgrade system.