The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critica...The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.展开更多
Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation o...Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.展开更多
The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(R...The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.展开更多
Oxygenated volatile organic compounds(OVOCs) are key intermediates in the atmospheric photooxidation process. To further study the primary and secondary sources of OVOCs,their ambient levels were monitored using a pro...Oxygenated volatile organic compounds(OVOCs) are key intermediates in the atmospheric photooxidation process. To further study the primary and secondary sources of OVOCs,their ambient levels were monitored using a proton-transfer reaction mass spectrometer(PTR-MS) at an urban site in the Pearl River Delta of China. Continuous monitoring campaigns were conducted in the spring, summer, fall, and winter of 2016. Among the six types of OVOC species, the mean concentrations of methanol were the highest in each season(up to 13–20 ppbv), followed by those of acetone, acetaldehyde and acetic acid(approximately 2–4 ppbv), while those of formic acid and methyl ethyl ketone(MEK) were the lowest(approximately 1–2 ppbv). As observed from a diurnal variation chart, the OVOCs observed in Shenzhen may have been affected by numerous factors such as their primary and secondary sources and photochemical consumption. The photochemical age-based parameterization method was used to apportion the sources of ambient OVOCs. Methanol had significant anthropogenic primary sources but negligible anthropogenic secondary sources during all of the seasons. Acetone, MEK and acetic acid were mostly attributed to anthropogenic primary sources during each season with smaller contributions from anthropogenic secondary sources. Acetaldehyde had similar contributions from both anthropogenic secondary and anthropogenic primary sources throughout the year.Meanwhile, anthropogenic primary sources contributed the most to formic acid.展开更多
To effectively estimate the unknown aerodynamic parameters from the aircraft’s flight data,this paper proposes a novel aerodynamic parameter estimation method incorporating a stacked Long Short-Term Memory(LSTM)netwo...To effectively estimate the unknown aerodynamic parameters from the aircraft’s flight data,this paper proposes a novel aerodynamic parameter estimation method incorporating a stacked Long Short-Term Memory(LSTM)network model and the Levenberg-Marquardt(LM)method.The stacked LSTM network model was designed to realize the aircraft dynamics modeling by utilizing a frame of nonlinear functional mapping based entirely on the measured input-output data of the aircraft system without requiring explicit postulation of the dynamics.The LM method combines the already-trained LSTM network model to optimize the unknown aerodynamic parameters.The proposed method is applied by using the real flight data,generated by ATTAS aircraft and a bio-inspired morphing Unmanned Aerial Vehicle(UAV).The investigation reveals that for the two different flight data,the designed stacked LSTM network structure can maintain the efficacy of the network prediction capability only by appropriately adjusting the dropout rates of its hidden layers without changing other network parameters(i.e.,the initial weights,initial biases,number of hidden cells,time-steps,learning rate,and number of training iterations).Besides,the proposed method’s effectiveness and potential are demonstrated by comparing the estimated results of the ATTAS aircraft or the bio-inspired morphing UAV with the corresponding reference values or wind-tunnel results.展开更多
The centroid coordinate serves as a critical control parameter in motion systems,including aircraft,missiles,rockets,and drones,directly influencing their motion dynamics and control performance.Traditional methods fo...The centroid coordinate serves as a critical control parameter in motion systems,including aircraft,missiles,rockets,and drones,directly influencing their motion dynamics and control performance.Traditional methods for centroid measurement often necessitate custom equipment and specialized positioning devices,leading to high costs and limited accuracy.Here,we present a centroid measurement method that integrates 3D scanning technology,enabling accurate measurement of centroid across various types of objects without the need for specialized positioning fixtures.A theoretical framework for centroid measurement was established,which combined the principle of the multi-point weighing method with 3D scanning technology.The measurement accuracy was evaluated using a designed standard component.Experimental results demonstrate that the discrepancies between the theoretical and the measured centroid of a standard component with various materials and complex shapes in the X,Y,and Z directions are 0.003 mm,0.009 mm,and 0.105 mm,respectively,yielding a spatial deviation of 0.106 mm.Qualitative verification was conducted through experimental validation of three distinct types.They confirmed the reliability of the proposed method,which allowed for accurate centroid measurements of various products without requiring positioning fixtures.This advancement significantly broadened the applicability and scope of centroid measurement devices,offering new theoretical insights and methodologies for the measurement of complex parts and systems.展开更多
The dramatic rise in the number of people living in cities has made many environmental and social problems worse.The search for a productive method for disposing of solid waste is the most notable of these problems.Ma...The dramatic rise in the number of people living in cities has made many environmental and social problems worse.The search for a productive method for disposing of solid waste is the most notable of these problems.Many scholars have referred to it as a fuzzy multi-attribute or multi-criteria decision-making problem using various fuzzy set-like approaches because of the inclusion of criteria and anticipated ambiguity.The goal of the current study is to use an innovative methodology to address the expected uncertainties in the problem of solid waste site selection.The characteristics(or sub-attributes)that decision-makers select and the degree of approximation they accept for various options can both be indicators of these uncertainties.To tackle these problems,a novel mathematical structure known as the fuzzy parameterized possibility single valued neutrosophic hypersoft expert set(ρˆ-set),which is initially described,is integrated with a modified version of Sanchez’s method.Following this,an intelligent algorithm is suggested.The steps of the suggested algorithm are explained with an example that explains itself.The compatibility of solid waste management sites and systems is discussed,and rankings are established along with detailed justifications for their viability.This study’s strengths lie in its application of fuzzy parameterization and possibility grading to effectively handle the uncertainties embodied in the parameters’nature and alternative approximations,respectively.It uses specific mathematical formulations to compute the fuzzy parameterized degrees and possibility grades that are missing from the prior literature.It is simpler for the decisionmakers to look at each option separately because the decision is uncertain.Comparing the computed results,it is discovered that they are consistent and dependable because of their preferred properties.展开更多
The investigation of leachate-contaminated clay(LCC)is essential for landfill engineering assessment and achievement of sustainable development goals.Several static and dynamic laboratory tests,including unconfined co...The investigation of leachate-contaminated clay(LCC)is essential for landfill engineering assessment and achievement of sustainable development goals.Several static and dynamic laboratory tests,including unconfined compressive strength(UCS),California bearing ratio(CBR),and cyclic simple shear,are conducted.Cyclic simple shear experiments on LCCs were performed to evaluate the damping and shear modulus.The investigated factors are vertical load(VL),leachate content(LC),frequency(F),and shear strain(ShS)for LCC.Forensic-based investigation optimization(FBIO)and equilibrium optimizer algorithm(EOA)were utilized in addition to multiple types of ensemble models,including adaptive boosting(ADB),gradient boosting regression tree(GBRT),extreme gradient boosting(XGB) and random forest(RF).The comparison of the methods showed that GBRT-FBIO and XGB-EOA models outperformed other models for shear modulus and damping of LCC.The p-value less than 0.0001 shows the significance of the used models in the response surface methodology(RSM)method.展开更多
A novel porous shock absorption layer is put forward in this study, and the shock absorption performance of the porous shock absorption layer is evaluated based on three-dimensional pseudo-static analysis. The modifie...A novel porous shock absorption layer is put forward in this study, and the shock absorption performance of the porous shock absorption layer is evaluated based on three-dimensional pseudo-static analysis. The modified reaction acceleration method is adopted and validated in the three-dimensional model. Seven ground motions are selected and the peak ground acceleration is adjusted to 0.2 g, 0.4 g and 0.6 g. The impact of the void ratio and thickness of the porous shock absorption layer is studied, while the surrounding rock grade and tunnel depth are also investigated. The numerical results show that the porous shock absorption layer has good shock absorption performance and can effectively reduce the maximum internal force of the secondary lining, but it cannot reduce the maximum horizontal relative displacement of the secondary lining. The circumferential rubber strip in the porous shock absorption layer will reduce shock absorption performance. The results of parameter analysis indicate that the shock absorption performance of the porous shock absorption layer increases with the increase of the void ratio and thickness, and it has good shock absorption performance under different surrounding rock grades and tunnel depths.展开更多
Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods rema...Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.展开更多
For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for s...For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for solving the dynamical equations of the constrained Birkhoffian system is provided.First the differential equations of motion for the constrained Birkhoffian system as well as for the corresponding free Birkhoffian system are established.Secondly a system of auxiliary equations is constructed and the general solution of the equations is found.Finally by varying the parameters and utilizing the properties of the generalized canonical transformation of the Birkhoffian system the solution of the problem can be obtained.The proposed method reveals the inherent relationship between the solution of a free Birkhoffian system and that of a constrained Birkhoffian system. The research results are of universal significance which can be further used in a variety of constrained mechanical systems such as non-conservative systems and nonholonomic systems etc.展开更多
In the traditional strength reduction method,the cohesion and the friction angle adopt the same reduction parameter,resulting in equivalent proportional reduction.This method does not consider the different effects of...In the traditional strength reduction method,the cohesion and the friction angle adopt the same reduction parameter,resulting in equivalent proportional reduction.This method does not consider the different effects of the cohesion and friction angle on the stability of the same slope and is defective to some extent.Regarding this defect,a strength reduction method based on double reduction parameters,which adopts different reduction parameters,is proposed.The core of the double-parameter reduction method is the matching reduction principle of the slope with different angles.This principle is represented by the ratio of the reduction parameter of the cohesion to that of the friction angle,described as η.With the increase in the slopeangle,ηincreases; in particular,when the slope angle is 45°,tηis 1.0.Through the matching reduction principle,different safety margin factors can be calculated for the cohesion and friction angle.In combination with these two safety margin factors,a formula for calculating the overall safety factor of the slope is proposed,reflecting the different contributions of the cohesion and friction angle to the slope stability.Finally,it is shown that the strength reduction method based on double reduction parameters acquires a larger safety factor than the classic limit equilibrium method,but the calculation results are very close to those obtained by the limit equilibrium method.展开更多
The formation and development of typhoons are closely related to the disturbed low vortexes at the planetary boundary layer(PBL). The effects of five PBL parameterization schemes(PBL schemes hereinafter) on the trajec...The formation and development of typhoons are closely related to the disturbed low vortexes at the planetary boundary layer(PBL). The effects of five PBL parameterization schemes(PBL schemes hereinafter) on the trajectory,intensity, and distribution of physical quantities are studied using the mesoscale WRF model on Super Typhoon Sanba(2012) during its initial stage. Results show that the five PBL schemes exhibit significant different effects on the simulated intensity and path. The results simulated by QNSE and ACM2 without the Bogus method are close to the best track data in the numerical experiments. When the Bogus method is adopted, the simulated trajectories improve significantly because the initial field is close to the true data. Among the five PBL schemes, QNSE and ACM2 with the Bogus method present improved simulated path and intensity compared with the three other schemes. This finding indicates that the two schemes deal with the initial PBL process satisfactorily, especially in the formation and development of disturbed low vortexes. The differences in the treatment methods of the five PBL schemes affect the surface layer physical quantities and the middle and upper atmospheres during the middle to late periods of the typhoon.Although QNSE and ACM2 present better simulation results than other schemes, they exhibit a few differences in the internal structure of the typhoon. The results simulated by MYJ are worse, and this method may be unsuitable for studying the formation and development of typhoons.展开更多
The on-orbit parameter identification of a space structure can be used for the modification of a system dynamics model and controller coefficients. This study focuses on the estimation of a system state-space model fo...The on-orbit parameter identification of a space structure can be used for the modification of a system dynamics model and controller coefficients. This study focuses on the estimation of a system state-space model for a two-link space manipulator in the procedure of capturing an unknown object, and a recursive tracking approach based on the recursive predictor-based subspace identification(RPBSID) algorithm is proposed to identify the manipulator payload mass parameter. Structural rigid motion and elastic vibration are separated, and the dynamics model of the space manipulator is linearized at an arbitrary working point(i.e., a certain manipulator configuration).The state-space model is determined by using the RPBSID algorithm and matrix transformation. In addition, utilizing the identified system state-space model, the manipulator payload mass parameter is estimated by extracting the corresponding block matrix. In numerical simulations, the presented parameter identification method is implemented and compared with the classical algebraic algorithm and the recursive least squares method for different payload masses and manipulator configurations. Numerical results illustrate that the system state-space model and payload mass parameter of the two-link flexible space manipulator are effectively identified by the recursive subspace tracking method.展开更多
In order to obtain accurate probability integration method(PIM) parameters for surface movement of multi-panel mining, a genetic algorithm(GA) was used to optimize the parameters. As the measured surface movement is a...In order to obtain accurate probability integration method(PIM) parameters for surface movement of multi-panel mining, a genetic algorithm(GA) was used to optimize the parameters. As the measured surface movement is affected by more than one mining panel, traditional PIM parameter inversion model is difficult to ensure the reliability of the results due to the complexity of rock movement. With crossover,mutation and selection operators, GA can perform a global optimization search and has high computation efficiency. Compared with the pattern search algorithm, the fitness function can avoid falling into local minima traps. GA reduces the risk of local minima traps which improves the accuracy and reliability with the mutation mechanism. Application at Xuehu colliery shows that GA can be used to inverse the PIM parameters for multi-panel surface movement observation, and reliable results can be obtained. The research provides a new way for back-analysis of PIM parameters for mining subsidence under complex conditions.展开更多
An estimation approach using least squares method was presented for identificationof model parameters of pressure control in shield tunneling.The state equation ofthe pressure control system for shield tunneling was a...An estimation approach using least squares method was presented for identificationof model parameters of pressure control in shield tunneling.The state equation ofthe pressure control system for shield tunneling was analytically derived based on themass equilibrium principle that the entry mass of the pressure chamber from cutting headwas equal to excluding mass from the screw conveyor.The randomly observed noise wasnumerically simulated and mixed to simulated observation values of system responses.The numerical simulation shows that the state equation of the pressure control system forshield tunneling is reasonable and the proposed estimation approach is effective even ifthe random observation noise exists.The robustness of the controlling procedure is validatedby numerical simulation results.展开更多
The parameter X of the Muskingum method is a physical parameter that reflects the flood peak attenuation and hydrograph shape flattening of a diffusion wave in motion. In this paper, the historic process that hydrolog...The parameter X of the Muskingum method is a physical parameter that reflects the flood peak attenuation and hydrograph shape flattening of a diffusion wave in motion. In this paper, the historic process that hydrologists have undergone to find a physical explanation of this parameter is briefly discussed. Based on the fact that the Muskingum method is the second-order accuracy difference solution to the diffusion wave equation, its numerical stability condition is analyzed, and a conclusion is drawn： X ≤ 0.5 is the uniform condition satisfying the demands for its physical meaning and numerical stability. It is also pointed out that the methods that regard the sum of squares of differences between the calculated and observed discharges or stages as the objective function and the routing coefficients C0, C1 and C2 of the Muskingum method as the optimization parameters cannot guarantee the physical meaning of X.展开更多
This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to impro...This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.展开更多
Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate...Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate regularized matrix, and a singular decomposition method was used to generate regularization parameters. Numerical test results suggest that the regularized ambiguity float solution is more stable and reliable than the least-squares float solution. The mean square error matrix of the new method possesses a lower correlation than the variancecovariance matrix of the least-squares estimation. The size of the ambiguity search space is reduced and the search efficiency is improved. The success rate of the integer ambiguity searching process is improved significantly when the ambiguity resolution by using constraint equation method is used to determine the correct ambiguity integervector. The ambiguity resolution by using constraint equation method requires an initial input of the ambiguity float solution candidates which are obtained from the LAMBDA method in the new method. In addition, the observation time required to fix reliable integer ambiguities can he significantly reduced.展开更多
The accurate material physical properties, initial and boundary conditions are indispensable to the numerical simulation in the casting process, and they are related to the simulation accuracy directly. The inverse he...The accurate material physical properties, initial and boundary conditions are indispensable to the numerical simulation in the casting process, and they are related to the simulation accuracy directly. The inverse heat conduction method can be used to identify the mentioned above parameters based on the temperature measurement data. This paper presented a new inverse method according to Tikhonov regularization theory. A regularization functional was established and the regularization parameter was deduced, the Newton-Raphson iteration method was used to solve the equations. One detailed case was solved to identify the thermal conductivity and specific heat of sand mold and interfacial heat transfer coefficient (IHTC) at the meantime. This indicates that the regularization method is very efficient in decreasing the sensitivity to the temperature measurement data, overcoming the ill-posedness of the inverse heat conduction problem (IHCP) and improving the stability and accuracy of the results. As a general inverse method, it can be used to identify not only the material physical properties but also the initial and boundary conditions' parameters.展开更多
基金National Key R&D Program of China(No.2017YFB1304000)Fundamental Research Funds for the Central Universities,China(No.2232023G-05-1)。
文摘The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.
基金supported by the National Natural Science Foundation of China(Grant Nos.12372200 and 12072242).
文摘Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.
基金supported by the National Key Research and Development Program of China(No.2023YFC2907600)the National Natural Science Foundation of China(Nos.42077267,42277174 and 52074164)+2 种基金the Natural Science Foundation of Shandong Province,China(No.ZR2020JQ23)the Opening Project of State Key Laboratory of Explosion Science and Technology,Beijing Institute of Technology(No.KFJJ21-02Z)the Fundamental Research Funds for the Central Universities,China(No.2022JCCXSB03).
文摘The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.
基金supported by the Ministry of Science and Technology of China (Nos.2017YFC0210004,2014BAC21B01)the Science and Technology Plan of Shenzhen Municipality (Nos.JCYJ20170412150626172,JCYJ20160122105855253)
文摘Oxygenated volatile organic compounds(OVOCs) are key intermediates in the atmospheric photooxidation process. To further study the primary and secondary sources of OVOCs,their ambient levels were monitored using a proton-transfer reaction mass spectrometer(PTR-MS) at an urban site in the Pearl River Delta of China. Continuous monitoring campaigns were conducted in the spring, summer, fall, and winter of 2016. Among the six types of OVOC species, the mean concentrations of methanol were the highest in each season(up to 13–20 ppbv), followed by those of acetone, acetaldehyde and acetic acid(approximately 2–4 ppbv), while those of formic acid and methyl ethyl ketone(MEK) were the lowest(approximately 1–2 ppbv). As observed from a diurnal variation chart, the OVOCs observed in Shenzhen may have been affected by numerous factors such as their primary and secondary sources and photochemical consumption. The photochemical age-based parameterization method was used to apportion the sources of ambient OVOCs. Methanol had significant anthropogenic primary sources but negligible anthropogenic secondary sources during all of the seasons. Acetone, MEK and acetic acid were mostly attributed to anthropogenic primary sources during each season with smaller contributions from anthropogenic secondary sources. Acetaldehyde had similar contributions from both anthropogenic secondary and anthropogenic primary sources throughout the year.Meanwhile, anthropogenic primary sources contributed the most to formic acid.
基金co-supported by the National Natural Science Foundation of China(No.52192633)the Natural Science Foundation of Shaanxi Province,China(No.2022JC-03)the Fundamental Research Funds for the Central Universities,China(No.XJSJ23164)。
文摘To effectively estimate the unknown aerodynamic parameters from the aircraft’s flight data,this paper proposes a novel aerodynamic parameter estimation method incorporating a stacked Long Short-Term Memory(LSTM)network model and the Levenberg-Marquardt(LM)method.The stacked LSTM network model was designed to realize the aircraft dynamics modeling by utilizing a frame of nonlinear functional mapping based entirely on the measured input-output data of the aircraft system without requiring explicit postulation of the dynamics.The LM method combines the already-trained LSTM network model to optimize the unknown aerodynamic parameters.The proposed method is applied by using the real flight data,generated by ATTAS aircraft and a bio-inspired morphing Unmanned Aerial Vehicle(UAV).The investigation reveals that for the two different flight data,the designed stacked LSTM network structure can maintain the efficacy of the network prediction capability only by appropriately adjusting the dropout rates of its hidden layers without changing other network parameters(i.e.,the initial weights,initial biases,number of hidden cells,time-steps,learning rate,and number of training iterations).Besides,the proposed method’s effectiveness and potential are demonstrated by comparing the estimated results of the ATTAS aircraft or the bio-inspired morphing UAV with the corresponding reference values or wind-tunnel results.
基金supported by National Natural Science Foundation of China(No.52176122).
文摘The centroid coordinate serves as a critical control parameter in motion systems,including aircraft,missiles,rockets,and drones,directly influencing their motion dynamics and control performance.Traditional methods for centroid measurement often necessitate custom equipment and specialized positioning devices,leading to high costs and limited accuracy.Here,we present a centroid measurement method that integrates 3D scanning technology,enabling accurate measurement of centroid across various types of objects without the need for specialized positioning fixtures.A theoretical framework for centroid measurement was established,which combined the principle of the multi-point weighing method with 3D scanning technology.The measurement accuracy was evaluated using a designed standard component.Experimental results demonstrate that the discrepancies between the theoretical and the measured centroid of a standard component with various materials and complex shapes in the X,Y,and Z directions are 0.003 mm,0.009 mm,and 0.105 mm,respectively,yielding a spatial deviation of 0.106 mm.Qualitative verification was conducted through experimental validation of three distinct types.They confirmed the reliability of the proposed method,which allowed for accurate centroid measurements of various products without requiring positioning fixtures.This advancement significantly broadened the applicability and scope of centroid measurement devices,offering new theoretical insights and methodologies for the measurement of complex parts and systems.
文摘The dramatic rise in the number of people living in cities has made many environmental and social problems worse.The search for a productive method for disposing of solid waste is the most notable of these problems.Many scholars have referred to it as a fuzzy multi-attribute or multi-criteria decision-making problem using various fuzzy set-like approaches because of the inclusion of criteria and anticipated ambiguity.The goal of the current study is to use an innovative methodology to address the expected uncertainties in the problem of solid waste site selection.The characteristics(or sub-attributes)that decision-makers select and the degree of approximation they accept for various options can both be indicators of these uncertainties.To tackle these problems,a novel mathematical structure known as the fuzzy parameterized possibility single valued neutrosophic hypersoft expert set(ρˆ-set),which is initially described,is integrated with a modified version of Sanchez’s method.Following this,an intelligent algorithm is suggested.The steps of the suggested algorithm are explained with an example that explains itself.The compatibility of solid waste management sites and systems is discussed,and rankings are established along with detailed justifications for their viability.This study’s strengths lie in its application of fuzzy parameterization and possibility grading to effectively handle the uncertainties embodied in the parameters’nature and alternative approximations,respectively.It uses specific mathematical formulations to compute the fuzzy parameterized degrees and possibility grades that are missing from the prior literature.It is simpler for the decisionmakers to look at each option separately because the decision is uncertain.Comparing the computed results,it is discovered that they are consistent and dependable because of their preferred properties.
文摘The investigation of leachate-contaminated clay(LCC)is essential for landfill engineering assessment and achievement of sustainable development goals.Several static and dynamic laboratory tests,including unconfined compressive strength(UCS),California bearing ratio(CBR),and cyclic simple shear,are conducted.Cyclic simple shear experiments on LCCs were performed to evaluate the damping and shear modulus.The investigated factors are vertical load(VL),leachate content(LC),frequency(F),and shear strain(ShS)for LCC.Forensic-based investigation optimization(FBIO)and equilibrium optimizer algorithm(EOA)were utilized in addition to multiple types of ensemble models,including adaptive boosting(ADB),gradient boosting regression tree(GBRT),extreme gradient boosting(XGB) and random forest(RF).The comparison of the methods showed that GBRT-FBIO and XGB-EOA models outperformed other models for shear modulus and damping of LCC.The p-value less than 0.0001 shows the significance of the used models in the response surface methodology(RSM)method.
基金Science and Technology Plan Project of Xizang Autonomous Region,China under Grant No.XZ202501YD0007。
文摘A novel porous shock absorption layer is put forward in this study, and the shock absorption performance of the porous shock absorption layer is evaluated based on three-dimensional pseudo-static analysis. The modified reaction acceleration method is adopted and validated in the three-dimensional model. Seven ground motions are selected and the peak ground acceleration is adjusted to 0.2 g, 0.4 g and 0.6 g. The impact of the void ratio and thickness of the porous shock absorption layer is studied, while the surrounding rock grade and tunnel depth are also investigated. The numerical results show that the porous shock absorption layer has good shock absorption performance and can effectively reduce the maximum internal force of the secondary lining, but it cannot reduce the maximum horizontal relative displacement of the secondary lining. The circumferential rubber strip in the porous shock absorption layer will reduce shock absorption performance. The results of parameter analysis indicate that the shock absorption performance of the porous shock absorption layer increases with the increase of the void ratio and thickness, and it has good shock absorption performance under different surrounding rock grades and tunnel depths.
基金sponsorship of the National Natural Science Foundation of China(42430809,42274157,42030103,42404132)the Fund of State Key Laboratory of Deep Oil and Gas,China University of Petroleum(East China)(SKLDOG2024-ZYTS-02)+5 种基金the Postdoctoral Fellowship Program of CPSF(GZB20240850)the Postdoctoral Project of Qingdao(QDBSH20240102082)the Fundamental Research Funds for the Central Universities(24CX07004A,24CX06036A)the CNPC Innovation Fund(2024DQ02-0505,2024DQ02-0136)the Innovation fund project for graduate student of China University of Petroleum(East China)the Fundamental Research Funds for the Central Universities(24CX04002A).
文摘Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.
基金The National Natural Science Foundation of China(No.10972151,11272227)
文摘For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for solving the dynamical equations of the constrained Birkhoffian system is provided.First the differential equations of motion for the constrained Birkhoffian system as well as for the corresponding free Birkhoffian system are established.Secondly a system of auxiliary equations is constructed and the general solution of the equations is found.Finally by varying the parameters and utilizing the properties of the generalized canonical transformation of the Birkhoffian system the solution of the problem can be obtained.The proposed method reveals the inherent relationship between the solution of a free Birkhoffian system and that of a constrained Birkhoffian system. The research results are of universal significance which can be further used in a variety of constrained mechanical systems such as non-conservative systems and nonholonomic systems etc.
基金Project(KZCX2-YW-T12)supported by the Chinese Academy of Science,China
文摘In the traditional strength reduction method,the cohesion and the friction angle adopt the same reduction parameter,resulting in equivalent proportional reduction.This method does not consider the different effects of the cohesion and friction angle on the stability of the same slope and is defective to some extent.Regarding this defect,a strength reduction method based on double reduction parameters,which adopts different reduction parameters,is proposed.The core of the double-parameter reduction method is the matching reduction principle of the slope with different angles.This principle is represented by the ratio of the reduction parameter of the cohesion to that of the friction angle,described as η.With the increase in the slopeangle,ηincreases; in particular,when the slope angle is 45°,tηis 1.0.Through the matching reduction principle,different safety margin factors can be calculated for the cohesion and friction angle.In combination with these two safety margin factors,a formula for calculating the overall safety factor of the slope is proposed,reflecting the different contributions of the cohesion and friction angle to the slope stability.Finally,it is shown that the strength reduction method based on double reduction parameters acquires a larger safety factor than the classic limit equilibrium method,but the calculation results are very close to those obtained by the limit equilibrium method.
基金National Natural Science Foundation of China(41375033,41225018)National Basic Research Program of China(973 Program)(2013CB430100)Fundamental Research Funds for Central Universities of Lanzhou University(LZUJBKY-2013-K16)
文摘The formation and development of typhoons are closely related to the disturbed low vortexes at the planetary boundary layer(PBL). The effects of five PBL parameterization schemes(PBL schemes hereinafter) on the trajectory,intensity, and distribution of physical quantities are studied using the mesoscale WRF model on Super Typhoon Sanba(2012) during its initial stage. Results show that the five PBL schemes exhibit significant different effects on the simulated intensity and path. The results simulated by QNSE and ACM2 without the Bogus method are close to the best track data in the numerical experiments. When the Bogus method is adopted, the simulated trajectories improve significantly because the initial field is close to the true data. Among the five PBL schemes, QNSE and ACM2 with the Bogus method present improved simulated path and intensity compared with the three other schemes. This finding indicates that the two schemes deal with the initial PBL process satisfactorily, especially in the formation and development of disturbed low vortexes. The differences in the treatment methods of the five PBL schemes affect the surface layer physical quantities and the middle and upper atmospheres during the middle to late periods of the typhoon.Although QNSE and ACM2 present better simulation results than other schemes, they exhibit a few differences in the internal structure of the typhoon. The results simulated by MYJ are worse, and this method may be unsuitable for studying the formation and development of typhoons.
基金funded by the National Natural Science Foundation of China (Nos. 11572069 and 51775541)the China Postdoctoral Science Foundation (No. 2016M601354)
文摘The on-orbit parameter identification of a space structure can be used for the modification of a system dynamics model and controller coefficients. This study focuses on the estimation of a system state-space model for a two-link space manipulator in the procedure of capturing an unknown object, and a recursive tracking approach based on the recursive predictor-based subspace identification(RPBSID) algorithm is proposed to identify the manipulator payload mass parameter. Structural rigid motion and elastic vibration are separated, and the dynamics model of the space manipulator is linearized at an arbitrary working point(i.e., a certain manipulator configuration).The state-space model is determined by using the RPBSID algorithm and matrix transformation. In addition, utilizing the identified system state-space model, the manipulator payload mass parameter is estimated by extracting the corresponding block matrix. In numerical simulations, the presented parameter identification method is implemented and compared with the classical algebraic algorithm and the recursive least squares method for different payload masses and manipulator configurations. Numerical results illustrate that the system state-space model and payload mass parameter of the two-link flexible space manipulator are effectively identified by the recursive subspace tracking method.
基金provided by the National Natural Science Foundation of China(No.51404272)the Hunan Province Key Laboratory of Coal Resources Clean-Utilization and Mine Environment Protection(No.E21224)
文摘In order to obtain accurate probability integration method(PIM) parameters for surface movement of multi-panel mining, a genetic algorithm(GA) was used to optimize the parameters. As the measured surface movement is affected by more than one mining panel, traditional PIM parameter inversion model is difficult to ensure the reliability of the results due to the complexity of rock movement. With crossover,mutation and selection operators, GA can perform a global optimization search and has high computation efficiency. Compared with the pattern search algorithm, the fitness function can avoid falling into local minima traps. GA reduces the risk of local minima traps which improves the accuracy and reliability with the mutation mechanism. Application at Xuehu colliery shows that GA can be used to inverse the PIM parameters for multi-panel surface movement observation, and reliable results can be obtained. The research provides a new way for back-analysis of PIM parameters for mining subsidence under complex conditions.
基金Supported by the National Basic Research Program of China(2007CB714006)the National Natural Science Foundation of China(90815023)
文摘An estimation approach using least squares method was presented for identificationof model parameters of pressure control in shield tunneling.The state equation ofthe pressure control system for shield tunneling was analytically derived based on themass equilibrium principle that the entry mass of the pressure chamber from cutting headwas equal to excluding mass from the screw conveyor.The randomly observed noise wasnumerically simulated and mixed to simulated observation values of system responses.The numerical simulation shows that the state equation of the pressure control system forshield tunneling is reasonable and the proposed estimation approach is effective even ifthe random observation noise exists.The robustness of the controlling procedure is validatedby numerical simulation results.
基金supported by the Scientific and Technological Basic Research Grant of the Ministry of Science and Technology of China (Grant No. 2007FY140900)the Public Welfare Industry Special Fund Project of the Ministry of Water Resources of China (Grant No. 200801033)
文摘The parameter X of the Muskingum method is a physical parameter that reflects the flood peak attenuation and hydrograph shape flattening of a diffusion wave in motion. In this paper, the historic process that hydrologists have undergone to find a physical explanation of this parameter is briefly discussed. Based on the fact that the Muskingum method is the second-order accuracy difference solution to the diffusion wave equation, its numerical stability condition is analyzed, and a conclusion is drawn： X ≤ 0.5 is the uniform condition satisfying the demands for its physical meaning and numerical stability. It is also pointed out that the methods that regard the sum of squares of differences between the calculated and observed discharges or stages as the objective function and the routing coefficients C0, C1 and C2 of the Muskingum method as the optimization parameters cannot guarantee the physical meaning of X.
基金supported by the National Natural Science Foundation of China(No.51379039)the Excellent Young Scientists Fund(No.51222904)
文摘This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.
文摘Based on the structural characteristics of the double-differenced normal equation, a new method was proposed to resolve the ambiguity float solution through a selection of parameter weights to construct an appropriate regularized matrix, and a singular decomposition method was used to generate regularization parameters. Numerical test results suggest that the regularized ambiguity float solution is more stable and reliable than the least-squares float solution. The mean square error matrix of the new method possesses a lower correlation than the variancecovariance matrix of the least-squares estimation. The size of the ambiguity search space is reduced and the search efficiency is improved. The success rate of the integer ambiguity searching process is improved significantly when the ambiguity resolution by using constraint equation method is used to determine the correct ambiguity integervector. The ambiguity resolution by using constraint equation method requires an initial input of the ambiguity float solution candidates which are obtained from the LAMBDA method in the new method. In addition, the observation time required to fix reliable integer ambiguities can he significantly reduced.
文摘The accurate material physical properties, initial and boundary conditions are indispensable to the numerical simulation in the casting process, and they are related to the simulation accuracy directly. The inverse heat conduction method can be used to identify the mentioned above parameters based on the temperature measurement data. This paper presented a new inverse method according to Tikhonov regularization theory. A regularization functional was established and the regularization parameter was deduced, the Newton-Raphson iteration method was used to solve the equations. One detailed case was solved to identify the thermal conductivity and specific heat of sand mold and interfacial heat transfer coefficient (IHTC) at the meantime. This indicates that the regularization method is very efficient in decreasing the sensitivity to the temperature measurement data, overcoming the ill-posedness of the inverse heat conduction problem (IHCP) and improving the stability and accuracy of the results. As a general inverse method, it can be used to identify not only the material physical properties but also the initial and boundary conditions' parameters.
