Phonons are the quantum mechanical descriptions of vibrational modes that manifest themselves in many physical properties of condensed matter systems. As the size of electronic devices continues to decrease below mean...Phonons are the quantum mechanical descriptions of vibrational modes that manifest themselves in many physical properties of condensed matter systems. As the size of electronic devices continues to decrease below mean free paths of acoustic phonons, the engineering of phonon spectra at the nanoscale becomes an important topic. Phonon manipulation allows for active control and management of heat fow, enabling functions such as regulated heat transport. At the same time, phonon transmission, as a novel signal transmission method, holds great potential to revolutionize modern industry like microelectronics technology, and boasts wide-ranging applications. Unlike fermions such as electrons, polarity regulation is difficult to act on phonons as bosons, making the development of effective phonon modulation methods a daunting task.This work reviews the development of phonon engineering and strategies of phonon manipulation at different scales, reports the latest research progress of nanophononic devices such as thermal rectifiers, thermal transistors, thermal memories, and thermoelectric devices,and analyzes the phonon transport mechanisms involved. Lastly, we survey feasible perspectives and research directions of phonon engineering. Thermoelectric analogies, external field regulation, and acousto-optic co-optimization are expected to become future research hotspots.展开更多
Piezoelectric semiconductors(PSs)possess both semiconducting properties and piezoelectric coupling effects,making them optimal building blocks for semiconductor devices.PS fiber-like structures have wide applications ...Piezoelectric semiconductors(PSs)possess both semiconducting properties and piezoelectric coupling effects,making them optimal building blocks for semiconductor devices.PS fiber-like structures have wide applications in multi-functional semiconductor devices.In this paper,a one-dimensional(1D)theoretical model is established to describe the piezotronic responses of a PS fiber under gradient temperature changes.The theoretical model aims to explain the mechanism behind the resistance change caused by such gradient temperature changes.Numerical results demonstrate that a gradient temperature change significantly affects the physical fields within the PS fiber,and can induce changes in its surface resistance.It provides important theoretical guidance on the development of piezotronic devices that are sensitive to temperature effects.展开更多
In the practical environment,it is very common for the simultaneous occurrence of base excitation and crosswind.Scavenging the combined energy of vibration and wind with a single energy harvesting structure is fascina...In the practical environment,it is very common for the simultaneous occurrence of base excitation and crosswind.Scavenging the combined energy of vibration and wind with a single energy harvesting structure is fascinating.For this purpose,the effects of the wind speed and random excitation level are investigated with the stochastic averaging method(SAM)based on the energy envelope.The results of the analytical prediction are verified with the Monte-Carlo method(MCM).The numerical simulation shows that the introduction of wind can reduce the critical excitation level for triggering an inter-well jump and make a bi-stable energy harvester(BEH)realize the performance enhancement for a weak base excitation.However,as the strength of the wind increases to a particular level,the influence of the random base excitation on the dynamic responses is weakened,and the system exhibits a periodic galloping response.A comparison between a BEH and a linear energy harvester(LEH)indicates that the BEH demonstrates inferior performance for high-speed wind.Relevant experiments are conducted to investigate the validity of the theoretical prediction and numerical simulation.The experimental findings also show that strong random excitation is favorable for the BEH in the range of low wind speeds.However,as the speed of the incoming wind is up to a particular level,the disadvantage of the BEH becomes clear and evident.展开更多
Marine umbilical is one of the key equipment for subsea oil and gas exploitation,which is usually integrated by a great number of different functional components with multi-layers.The layout of these components direct...Marine umbilical is one of the key equipment for subsea oil and gas exploitation,which is usually integrated by a great number of different functional components with multi-layers.The layout of these components directly affects manufacturing,operation and storage performances of the umbilical.For the multi-layer cross-sectional layout design of the umbilical,a quantifiable multi-objective optimization model is established according to the operation and storage requirements.Considering the manufacturing factors,the multi-layering strategy based on contact point identification is introduced for a great number of functional components.Then,the GA-GLM global optimization algorithm is proposed combining the genetic algorithm and the generalized multiplier method,and the selection operator of the genetic algorithm is improved based on the steepest descent method.Genetic algorithm is used to find the optimal solution in the global space,which can converge from any initial layout to the feasible layout solution.The feasible layout solution is taken as the initial value of the generalized multiplier method for fast and accurate solution.Finally,taking umbilicals with a great number of components as examples,the results show that the cross-sectional performance of the umbilical obtained by optimization algorithm is better and the solution efficiency is higher.Meanwhile,the multi-layering strategy is effective and feasible.The design method proposed in this paper can quickly obtain the optimal multi-layer cross-sectional layout,which replaces the manual design,and provides useful reference and guidance for the umbilical industry.展开更多
Besides exhibiting excellent capabilities such as energy absorption,phase-transforming metamaterials offer a vast design space for achieving nonlinear constitutive relations.This is facilitated by switching between di...Besides exhibiting excellent capabilities such as energy absorption,phase-transforming metamaterials offer a vast design space for achieving nonlinear constitutive relations.This is facilitated by switching between different patterns under deformation.However,the related inverse design problem is quite challenging,due to the lack of appropriate mathematical formulation and the convergence issue in the post-buckling analysis of intermediate designs.In this work,periodic unit cells are explicitly described by the moving morphable voids method and effectively analyzed by eliminating the degrees of freedom in void regions.Furthermore,by exploring the Pareto frontiers between error and cost,an inverse design formulation is proposed for unit cells.This formulation aims to achieve a prescribed constitutive curve and is validated through numerical examples and experimental results.The design approach presented here can be extended to the inverse design of other types of mechanical metamaterials with prescribed nonlinear effective properties.展开更多
High spatiotemporal resolution brain electrical signals are critical for basic neuroscience research and high-precision focus diagnostic localization,as the spatial scale of some pathologic signals is at the submillim...High spatiotemporal resolution brain electrical signals are critical for basic neuroscience research and high-precision focus diagnostic localization,as the spatial scale of some pathologic signals is at the submillimeter or micrometer level.This entails connecting hundreds or thousands of electrode wires on a limited surface.This study reported a class of flexible,ultrathin,highdensity electrocorticogram(ECoG)electrode arrays.The challenge of a large number of wiring arrangements was overcome by a laminated structure design and processing technology improvement.The flexible,ultrathin,high-density ECoG electrode array was conformably attached to the cortex for reliable,high spatial resolution electrophysiologic recordings.The minimum spacing between electrodes was 15μm,comparable to the diameter of a single neuron.Eight hundred electrodes were prepared with an electrode density of 4444 mm^(-2).In focal epilepsy surgery,the flexible,high-density,laminated ECoG electrode array with 36 electrodes was applied to collect epileptic spike waves inrabbits,improving the positioning accuracy of epilepsy lesions from the centimeter to the submillimeter level.The flexible,high-density,laminated ECoG electrode array has potential clinical applications in intractable epilepsy and other neurologic diseases requiring high-precision electroencephalogram acquisition.展开更多
Entropy production in quasi-isentropic compression (QIC) is critically important for understanding the properties of materials under extremeconditions. However, the origin and accurate quantification of entropy in thi...Entropy production in quasi-isentropic compression (QIC) is critically important for understanding the properties of materials under extremeconditions. However, the origin and accurate quantification of entropy in this situation remain long-standing challenges. In this work, a framework is established for the quantification of entropy production and partition, and their relation to microstructural change in QIC. Cu50Zr50is taken as a model material, and its compression is simulated by molecular dynamics. On the basis of atomistic simulation-informed physicalproperties and free energy, the thermodynamic path is recovered, and the entropy production and its relation to microstructural change aresuccessfully quantified by the proposed framework. Contrary to intuition, entropy production during QIC of metallic glasses is relativelyinsensitive to the strain rate ˙γ when ˙γ ranges from 7.5 × 10^(8) to 2 × 10^(9)/s, which are values reachable in QIC experiments, with a magnitudeof the order of 10^(−2)kB/atom per GPa. However, when ˙γ is extremely high (>2 × 10^(9)/s), a notable increase in entropy production rate with˙γ is observed. The Taylor–Quinney factor is found to vary with strain but not with strain rate in the simulated regime. It is demonstrated thatentropy production is dominated by the configurational part, compared with the vibrational part. In the rate-insensitive regime, the increase inconfigurational entropy exhibits a linear relation to the Shannon-entropic quantification of microstructural change, and a stretched exponential relation to the Taylor–Quinney factor. The quantification of entropy is expected to provide thermodynamic insights into the fundamentalrelation between microstructure evolution and plastic dissipation.展开更多
Recently,numerous studies have demonstrated that the physics-informed neural network(PINN)can effectively and accurately resolve hyperelastic finite deformation problems.In this paper,a PINN framework for tackling hyp...Recently,numerous studies have demonstrated that the physics-informed neural network(PINN)can effectively and accurately resolve hyperelastic finite deformation problems.In this paper,a PINN framework for tackling hyperelastic-magnetic coupling problems is proposed.Since the solution space consists of two-phase domains,two separate networks are constructed to independently predict the solution for each phase region.In addition,a conscious point allocation strategy is incorporated to enhance the prediction precision of the PINN in regions characterized by sharp gradients.With the developed framework,the magnetic fields and deformation fields of magnetorheological elastomers(MREs)are solved under the control of hyperelastic-magnetic coupling equations.Illustrative examples are provided and contrasted with the reference results to validate the predictive accuracy of the proposed framework.Moreover,the advantages of the proposed framework in solving hyperelastic-magnetic coupling problems are validated,particularly in handling small data sets,as well as its ability in swiftly and precisely forecasting magnetostrictive motion.展开更多
The damage evolution process of non-penetrating cracks often causes some unexpected engineering disasters.Gypsum specimens containing non-penetrating crack(s)are used to study the damage evolution and characteristics ...The damage evolution process of non-penetrating cracks often causes some unexpected engineering disasters.Gypsum specimens containing non-penetrating crack(s)are used to study the damage evolution and characteristics under cyclic loading.The results show that under cyclic loading,the relationship between the number of non-penetrating crack(s)and the characteristic parameters(cyclic number,peak stress,peak strain,failure stress,and failure strain)of the pre-cracked specimens can be represented by a decreasing linear function.The damage evolution equation is fitted by calibrating the accumulative plastic strain for each cycle,and the damage constitutive equation is proposed by the concept of effective stress.Additionally,non-penetrating cracks are more likely to cause uneven stress distribution,damage accumulation,and local failure of specimen.The local failure can change the stress distribution and relieve the inhibition of non-penetrating crack extension and eventually cause a dramatic destruction of the specimen.Therefore,the evolution process caused by non-penetrating cracks can be regarded as one of the important reasons for inducing rockburst.These results are expected to improve the understanding of the process of spalling formation and rockburst and can be used to analyze the stability of rocks or rock structures.展开更多
The micro-ablation processes and morphological evolution of ablative craters on single-crystal magnesium under subpicosecond laser irradiation are investigated using molecular dynamics(MD) simulations and experiments....The micro-ablation processes and morphological evolution of ablative craters on single-crystal magnesium under subpicosecond laser irradiation are investigated using molecular dynamics(MD) simulations and experiments.The simulation results exhibit that the main failure mode of single-crystal Mg film irradiated by a low fluence and long pulse width laser is the ejection of surface atoms,which has laser-induced high stress.However,under high fluence and short pulse width laser irradiation,the main damage mechanism is nucleation fracture caused by stress wave reflection and superposition at the bottom of the film.In addition,Mg[0001] has higher pressure sensitivity and is more prone to ablation than Mg[0001].The evolution equation of crater depth is established using multi-pulse laser ablation simulation and verified by experiments.The results show that,under multiple pulsed laser irradiation,not only does the crater depth increase linearly with the pulse number,but also the quadratic term and constant term of the fitted crater profile curve increase linearly.展开更多
This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and ...This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and other surface of the PSC half-space are assumed to be electrically insulating.By the Hankel integral transformation,the problem is reduced to the Fredholm integral equation of the second kind.This equation is solved numerically to obtain the indentation behaviors of the PSC half-space,mainly including the indentation force-depth relation and the electric potential-depth relation.The results show that the effect of the semiconductor property on the indentation responses is limited within a certain range of variation of the steady carrier concentration.The dependence of indentation behavior on material properties is also analyzed by two different kinds of PSCs.Finite element simulations are conducted to verify the results calculated by the integral equation technique,and good agreement is demonstrated.展开更多
The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is nece...The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is necessary to accurately predict the shakedown domains of these materials. The static shakedown theorem, also known as Melan's theorem, is a fundamental method used to predict the shakedown domains of structures and materials. Within this method, a key aspect lies in the construction and application of an appropriate self-equilibrium stress field(SSF). In the structural shakedown analysis, the SSF is typically constructed by governing equations that satisfy no external force(NEF) boundary conditions. However, we discover that directly applying these governing equations is not suitable for the shakedown analysis of heterogeneous materials. Researchers must consider the requirements imposed by the Hill-Mandel condition for boundary conditions and the physical significance of representative volume elements(RVEs). This paper addresses this issue and demonstrates that the sizes of SSFs vary under different boundary conditions, such as uniform displacement boundary conditions(DBCs), uniform traction boundary conditions(TBCs), and periodic boundary conditions(PBCs). As a result, significant discrepancies arise in the predicted shakedown domain sizes of heterogeneous materials. Built on the demonstrated relationship between SSFs under different boundary conditions, this study explores the conservative relationships among different shakedown domains, and provides proof of the relationship between the elastic limit(EL) factors and the shakedown loading factors under the loading domain of two load vertices. By utilizing numerical examples, we highlight the conservatism present in certain results reported in the existing literature. Among the investigated boundary conditions, the obtained shakedown domain is the most conservative under TBCs.Conversely, utilizing PBCs to construct an SSF for the shakedown analysis leads to less conservative lower bounds, indicating that PBCs should be employed as the preferred boundary conditions for the shakedown analysis of heterogeneous materials.展开更多
Soft robot incarnates its unique advantages in deep-sea exploration,but grapples with high hydrostatic pressure’s unpredictable impact on its mechanical performances.In our previous work,a self-powered soft robot sho...Soft robot incarnates its unique advantages in deep-sea exploration,but grapples with high hydrostatic pressure’s unpredictable impact on its mechanical performances.In our previous work,a self-powered soft robot showed excellent work performance in the Mariana Trench at a depth of 11000 m,yet experienced notable degradation in deforming capability.Here,we propose a magnetic loading method for characterizing elastomer’s mechanical properties under extremely high hydrostatic pressure of up to 120 MPa.This method facilitates remote loading and enables in-situ observation,so that the dimensions and deformation at high hydrostatic pressure are obtained and used for calculations.The results reveal that the Young’s modulus of Polydimethylsiloxane(PDMS)monotonously increases with pressure.It is found that the relative increase in Young’s modulus is determined by its initial value,which is 8% for an initial Young’s modulus of 2200 kPa and 38% for 660 kPa.The relation between initial Young’s modulus and relevant increase can be fitted by an exponential function.The bulk modulus of PDMS is about 1.4 GPa at 20℃ and is barely affected by hydrostatic pressure.The method can quantify alterations in the mechanical properties of elastomers induced by hydrostatic pressure,and provide guidance for the design of soft robots which serve in extreme pressure environment.展开更多
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...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.展开更多
Hunting stability is an important performance criterion in railway vehicles.This study proposes an incorporation of a bio-inspired limb-like structure(LLS)-based nonlinear damping into the motor suspension system for ...Hunting stability is an important performance criterion in railway vehicles.This study proposes an incorporation of a bio-inspired limb-like structure(LLS)-based nonlinear damping into the motor suspension system for traction units to improve the nonlinear critical speed and hunting stability of high-speed trains(HSTs).Initially,a vibration transmission analysis is conducted on a HST vehicle and a metro vehicle that suffered from hunting motion to explore the effect of different motor suspension systems from on-track tests.Subsequently,a simplified lateral dynamics model of an HST bogie is established to investigate the influence of the motor suspension on the bogie hunting behavior.The bifurcation analysis is applied to optimize the motor suspension parameters for high critical speed.Then,the nonlinear damping of the bio-inspired LLS,which has a positive correlation with the relative displacement,can further improve the modal damping of hunting motion and nonlinear critical speed compared with the linear motor suspension system.Furthermore,a comprehensive numerical model of a high-speed train,considering all nonlinearities,is established to investigate the influence of different types of motor suspension.The simulation results are well consistent with the theoretical analysis.The benefits of employing nonlinear damping of the bio-inspired LLS into the motor suspension of HSTs to enhance bogie hunting stability are thoroughly validated.展开更多
Recently,a Schwarz crystal structure with curved grain boundaries(GBs)constrained by twin-boundary(TB)networks was discovered in nanocrystalline Cu through experiments and atomistic simulations.Nanocrystalline Cu with...Recently,a Schwarz crystal structure with curved grain boundaries(GBs)constrained by twin-boundary(TB)networks was discovered in nanocrystalline Cu through experiments and atomistic simulations.Nanocrystalline Cu with nanosized Schwarz crystals exhibited high strength and excellent thermal stability.However,the grainsize effect and associated deformation mechanisms of Schwarz nanocrystals remain unknown.Here,we performed large-scale atomistic simulations to investigate the deformation behaviors and grain-size effect of nanocrystalline Cu with Schwarz crystals.Our simulations showed that similar to regular nanocrystals,Schwarz nanocrystals exhibit a strengthening-softening transition with decreasing grain size.The critical grain size in Schwarz nanocrystals is smaller than that in regular nanocrystals,leading to a maximum strength higher than that of regular nanocrystals.Our simulations revealed that the softening in Schwarz nanocrystals mainly originates from TB migration(or detwinning)and annihilation of GBs,rather than GB-mediated processes(including GB migration,sliding and diffusion)dominating the softening in regular nanocrystals.Quantitative analyses of simulation data further showed that compared with those in regular nanocrystals,the GB-mediated processes in Schwarz nanocrystals are suppressed,which is related to the low volume fraction of amorphous-like GBs and constraints of TB networks.The smaller critical grain size arises from the suppression of GB-mediated processes.展开更多
In the context of global mean square error concerning the number of random variables in the representation,the Karhunen–Loève(KL)expansion is the optimal series expansion method for random field discretization.T...In the context of global mean square error concerning the number of random variables in the representation,the Karhunen–Loève(KL)expansion is the optimal series expansion method for random field discretization.The computational efficiency and accuracy of the KL expansion are contingent upon the accurate resolution of the Fredholm integral eigenvalue problem(IEVP).The paper proposes an interpolation method based on different interpolation basis functions such as moving least squares(MLS),least squares(LS),and finite element method(FEM)to solve the IEVP.Compared with the Galerkin method based on finite element or Legendre polynomials,the main advantage of the interpolation method is that,in the calculation of eigenvalues and eigenfunctions in one-dimensional random fields,the integral matrix containing covariance function only requires a single integral,which is less than a two-folded integral by the Galerkin method.The effectiveness and computational efficiency of the proposed interpolation method are verified through various one-dimensional examples.Furthermore,based on theKL expansion and polynomial chaos expansion,the stochastic analysis of two-dimensional regular and irregular domains is conducted,and the basis function of the extended finite element method(XFEM)is introduced as the interpolation basis function in two-dimensional irregular domains to solve the IEVP.展开更多
This study proposes a comprehensive,coupled thermomechanical model that replaces local spatial derivatives in classical differential thermomechanical equations with nonlocal integral forms derived from the peridynamic...This study proposes a comprehensive,coupled thermomechanical model that replaces local spatial derivatives in classical differential thermomechanical equations with nonlocal integral forms derived from the peridynamic differential operator(PDDO),eliminating the need for calibration procedures.The model employs a multi-rate explicit time integration scheme to handle varying time scales in multi-physics systems.Through simulations conducted on granite and ceramic materials,this model demonstrates its effectiveness.It successfully simulates thermal damage behavior in granite arising from incompatible mineral expansion and accurately calculates thermal crack propagation in ceramic slabs during quenching.To account for material heterogeneity,the model utilizes the Shuffle algorithm andWeibull distribution,yielding results that align with numerical simulations and experimental observations.This coupled thermomechanical model shows great promise for analyzing intricate thermomechanical phenomena in brittle materials.展开更多
Inspired by the excellent adhesion performances of setae structure from organisms,micro/nano-pillar array has become one of the paradigms for adhesive surfaces.The micropillar arrays are composed of the resin pillars ...Inspired by the excellent adhesion performances of setae structure from organisms,micro/nano-pillar array has become one of the paradigms for adhesive surfaces.The micropillar arrays are composed of the resin pillars for adhesion and the substrate with different elastic modulus for supporting.The stress singularity at the bi-material corner between the pillars and the substrate can induce the failure of the micropillar-substrate corner and further hinder the fabrication and application of micropillar arrays,yet the design for the stability of the micropillar array lacks systematical and quantitative guides.In this work,we develop a semi-analytical method to provide the full expressions for the stress distribution within the bi-material corner combining analytical derivations and numerical calculations.The predictions for the stress within the singularity field can be obtained based on the full expressions of the stress.The good agreement between the predictions and the FEM results demonstrates the high reliability of our method.By adopting the strain energy density factor approach,the stability of the pillar-substrate corner is assessed by predicting the failure at the corner.For the elastic mismatch between the pillar and substrate given in this paper,the stability can be improved by increasing the ratio of the shear modulus of the substrate to that of the micropillar.Our study provides accurate predictions for the stress distribution at the bi-material corner and can guide the optimization of material combinations of the pillars and the substrate for more stable bioinspired dry adhesives.展开更多
Due to technical limitations,existing vibration isolation and energy harvesting(VIEH)devices have poor performance at low frequency.This paper proposes a new multilink-spring mechanism(MLSM)that can be used to solve t...Due to technical limitations,existing vibration isolation and energy harvesting(VIEH)devices have poor performance at low frequency.This paper proposes a new multilink-spring mechanism(MLSM)that can be used to solve this problem.The VIEH performance of the MLSM under harmonic excitation and Gaussian white noise was analyzed.It was found that the MLSM has good vibration isolation performance for low-frequency isolation and the frequency band can be widened by adjusting parameters to achieve a higher energy harvesting power.By comparison with two special cases,the results show that the MLSM is basically the same as the other two oscillators in terms of vibration isolation but has better energy harvesting performance under multistable characteristics.The MLSM is expected to reduce the impact of vibration on high-precision sensitive equipment in some special sites such as subways and mines,and at the same time supply power to structural health monitoring devices.展开更多
基金supported by the National Natural Science Foundation of China Grant Nos. 52276072 and 51976096。
文摘Phonons are the quantum mechanical descriptions of vibrational modes that manifest themselves in many physical properties of condensed matter systems. As the size of electronic devices continues to decrease below mean free paths of acoustic phonons, the engineering of phonon spectra at the nanoscale becomes an important topic. Phonon manipulation allows for active control and management of heat fow, enabling functions such as regulated heat transport. At the same time, phonon transmission, as a novel signal transmission method, holds great potential to revolutionize modern industry like microelectronics technology, and boasts wide-ranging applications. Unlike fermions such as electrons, polarity regulation is difficult to act on phonons as bosons, making the development of effective phonon modulation methods a daunting task.This work reviews the development of phonon engineering and strategies of phonon manipulation at different scales, reports the latest research progress of nanophononic devices such as thermal rectifiers, thermal transistors, thermal memories, and thermoelectric devices,and analyzes the phonon transport mechanisms involved. Lastly, we survey feasible perspectives and research directions of phonon engineering. Thermoelectric analogies, external field regulation, and acousto-optic co-optimization are expected to become future research hotspots.
基金Project supported by the National Natural Science Foundation of China (Nos.12172326 and 11972319)the National Key Research and Development Program of China (No.2020YFA0711700)the Natural Science Foundation of Zhejiang Province of China (No.LR21A020002)。
文摘Piezoelectric semiconductors(PSs)possess both semiconducting properties and piezoelectric coupling effects,making them optimal building blocks for semiconductor devices.PS fiber-like structures have wide applications in multi-functional semiconductor devices.In this paper,a one-dimensional(1D)theoretical model is established to describe the piezotronic responses of a PS fiber under gradient temperature changes.The theoretical model aims to explain the mechanism behind the resistance change caused by such gradient temperature changes.Numerical results demonstrate that a gradient temperature change significantly affects the physical fields within the PS fiber,and can induce changes in its surface resistance.It provides important theoretical guidance on the development of piezotronic devices that are sensitive to temperature effects.
基金Project supported by the National Natural Science Foundation of China(Nos.12272355,1202520411902294)+1 种基金the Opening Foundation of Shanxi Provincial Key Laboratory for Advanced Manufacturing Technology of China(No.XJZZ202304)the Shanxi Provincial Graduate Innovation Project of China(No.2023KY629)。
文摘In the practical environment,it is very common for the simultaneous occurrence of base excitation and crosswind.Scavenging the combined energy of vibration and wind with a single energy harvesting structure is fascinating.For this purpose,the effects of the wind speed and random excitation level are investigated with the stochastic averaging method(SAM)based on the energy envelope.The results of the analytical prediction are verified with the Monte-Carlo method(MCM).The numerical simulation shows that the introduction of wind can reduce the critical excitation level for triggering an inter-well jump and make a bi-stable energy harvester(BEH)realize the performance enhancement for a weak base excitation.However,as the strength of the wind increases to a particular level,the influence of the random base excitation on the dynamic responses is weakened,and the system exhibits a periodic galloping response.A comparison between a BEH and a linear energy harvester(LEH)indicates that the BEH demonstrates inferior performance for high-speed wind.Relevant experiments are conducted to investigate the validity of the theoretical prediction and numerical simulation.The experimental findings also show that strong random excitation is favorable for the BEH in the range of low wind speeds.However,as the speed of the incoming wind is up to a particular level,the disadvantage of the BEH becomes clear and evident.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52001088,52271269,U1906233)the Natural Science Foundation of Heilongjiang Province(Grant No.LH2021E050)+2 种基金the State Key Laboratory of Ocean Engineering(Grant No.GKZD010084)Liaoning Province’s Xing Liao Talents Program(Grant No.XLYC2002108)Dalian City Supports Innovation and Entrepreneurship Projects for High-Level Talents(Grant No.2021RD16)。
文摘Marine umbilical is one of the key equipment for subsea oil and gas exploitation,which is usually integrated by a great number of different functional components with multi-layers.The layout of these components directly affects manufacturing,operation and storage performances of the umbilical.For the multi-layer cross-sectional layout design of the umbilical,a quantifiable multi-objective optimization model is established according to the operation and storage requirements.Considering the manufacturing factors,the multi-layering strategy based on contact point identification is introduced for a great number of functional components.Then,the GA-GLM global optimization algorithm is proposed combining the genetic algorithm and the generalized multiplier method,and the selection operator of the genetic algorithm is improved based on the steepest descent method.Genetic algorithm is used to find the optimal solution in the global space,which can converge from any initial layout to the feasible layout solution.The feasible layout solution is taken as the initial value of the generalized multiplier method for fast and accurate solution.Finally,taking umbilicals with a great number of components as examples,the results show that the cross-sectional performance of the umbilical obtained by optimization algorithm is better and the solution efficiency is higher.Meanwhile,the multi-layering strategy is effective and feasible.The design method proposed in this paper can quickly obtain the optimal multi-layer cross-sectional layout,which replaces the manual design,and provides useful reference and guidance for the umbilical industry.
基金supported by the National Natural Science Foun-dation of China(Grant Nos.12002073 and 12372122)the National Key Research and Development Plan of China(Grant No.2020YFB 1709401)+2 种基金the Science Technology Plan of Liaoning Province(Grant No.2023JH2/101600044)the Liaoning Revitalization Talents Pro-gram(Grant No.XLYC2001003)111 Project of China(Grant No.B14013).
文摘Besides exhibiting excellent capabilities such as energy absorption,phase-transforming metamaterials offer a vast design space for achieving nonlinear constitutive relations.This is facilitated by switching between different patterns under deformation.However,the related inverse design problem is quite challenging,due to the lack of appropriate mathematical formulation and the convergence issue in the post-buckling analysis of intermediate designs.In this work,periodic unit cells are explicitly described by the moving morphable voids method and effectively analyzed by eliminating the degrees of freedom in void regions.Furthermore,by exploring the Pareto frontiers between error and cost,an inverse design formulation is proposed for unit cells.This formulation aims to achieve a prescribed constitutive curve and is validated through numerical examples and experimental results.The design approach presented here can be extended to the inverse design of other types of mechanical metamaterials with prescribed nonlinear effective properties.
基金support of the National Natural Science Foundation of China(Nos.U20A6001,12002190,11972207,and 11921002)the Fundamental Research Funds for the Central Universities,China(No.SWUKQ22029)the Chongqing Natural Science Foundation of China(No.CSTB2022NSCQ-MSX1635).
文摘High spatiotemporal resolution brain electrical signals are critical for basic neuroscience research and high-precision focus diagnostic localization,as the spatial scale of some pathologic signals is at the submillimeter or micrometer level.This entails connecting hundreds or thousands of electrode wires on a limited surface.This study reported a class of flexible,ultrathin,highdensity electrocorticogram(ECoG)electrode arrays.The challenge of a large number of wiring arrangements was overcome by a laminated structure design and processing technology improvement.The flexible,ultrathin,high-density ECoG electrode array was conformably attached to the cortex for reliable,high spatial resolution electrophysiologic recordings.The minimum spacing between electrodes was 15μm,comparable to the diameter of a single neuron.Eight hundred electrodes were prepared with an electrode density of 4444 mm^(-2).In focal epilepsy surgery,the flexible,high-density,laminated ECoG electrode array with 36 electrodes was applied to collect epileptic spike waves inrabbits,improving the positioning accuracy of epilepsy lesions from the centimeter to the submillimeter level.The flexible,high-density,laminated ECoG electrode array has potential clinical applications in intractable epilepsy and other neurologic diseases requiring high-precision electroencephalogram acquisition.
基金supported by the NSAF under Grant No.U1830206,the National Key R&D Program of China under Grant No.2017YFA0403200the National Natural Science Foundation of China under Grant Nos.11874424 and 12104507the Science and Technology Innovation Program of Hunan Province under Grant No.2021RC4026.
文摘Entropy production in quasi-isentropic compression (QIC) is critically important for understanding the properties of materials under extremeconditions. However, the origin and accurate quantification of entropy in this situation remain long-standing challenges. In this work, a framework is established for the quantification of entropy production and partition, and their relation to microstructural change in QIC. Cu50Zr50is taken as a model material, and its compression is simulated by molecular dynamics. On the basis of atomistic simulation-informed physicalproperties and free energy, the thermodynamic path is recovered, and the entropy production and its relation to microstructural change aresuccessfully quantified by the proposed framework. Contrary to intuition, entropy production during QIC of metallic glasses is relativelyinsensitive to the strain rate ˙γ when ˙γ ranges from 7.5 × 10^(8) to 2 × 10^(9)/s, which are values reachable in QIC experiments, with a magnitudeof the order of 10^(−2)kB/atom per GPa. However, when ˙γ is extremely high (>2 × 10^(9)/s), a notable increase in entropy production rate with˙γ is observed. The Taylor–Quinney factor is found to vary with strain but not with strain rate in the simulated regime. It is demonstrated thatentropy production is dominated by the configurational part, compared with the vibrational part. In the rate-insensitive regime, the increase inconfigurational entropy exhibits a linear relation to the Shannon-entropic quantification of microstructural change, and a stretched exponential relation to the Taylor–Quinney factor. The quantification of entropy is expected to provide thermodynamic insights into the fundamentalrelation between microstructure evolution and plastic dissipation.
基金supported by the National Natural Science Foundation of China(Nos.12072105 and 11932006)。
文摘Recently,numerous studies have demonstrated that the physics-informed neural network(PINN)can effectively and accurately resolve hyperelastic finite deformation problems.In this paper,a PINN framework for tackling hyperelastic-magnetic coupling problems is proposed.Since the solution space consists of two-phase domains,two separate networks are constructed to independently predict the solution for each phase region.In addition,a conscious point allocation strategy is incorporated to enhance the prediction precision of the PINN in regions characterized by sharp gradients.With the developed framework,the magnetic fields and deformation fields of magnetorheological elastomers(MREs)are solved under the control of hyperelastic-magnetic coupling equations.Illustrative examples are provided and contrasted with the reference results to validate the predictive accuracy of the proposed framework.Moreover,the advantages of the proposed framework in solving hyperelastic-magnetic coupling problems are validated,particularly in handling small data sets,as well as its ability in swiftly and precisely forecasting magnetostrictive motion.
基金supported by the National Natural Science Foundation of China(Nos.52204092 and 52274203).
文摘The damage evolution process of non-penetrating cracks often causes some unexpected engineering disasters.Gypsum specimens containing non-penetrating crack(s)are used to study the damage evolution and characteristics under cyclic loading.The results show that under cyclic loading,the relationship between the number of non-penetrating crack(s)and the characteristic parameters(cyclic number,peak stress,peak strain,failure stress,and failure strain)of the pre-cracked specimens can be represented by a decreasing linear function.The damage evolution equation is fitted by calibrating the accumulative plastic strain for each cycle,and the damage constitutive equation is proposed by the concept of effective stress.Additionally,non-penetrating cracks are more likely to cause uneven stress distribution,damage accumulation,and local failure of specimen.The local failure can change the stress distribution and relieve the inhibition of non-penetrating crack extension and eventually cause a dramatic destruction of the specimen.Therefore,the evolution process caused by non-penetrating cracks can be regarded as one of the important reasons for inducing rockburst.These results are expected to improve the understanding of the process of spalling formation and rockburst and can be used to analyze the stability of rocks or rock structures.
文摘The micro-ablation processes and morphological evolution of ablative craters on single-crystal magnesium under subpicosecond laser irradiation are investigated using molecular dynamics(MD) simulations and experiments.The simulation results exhibit that the main failure mode of single-crystal Mg film irradiated by a low fluence and long pulse width laser is the ejection of surface atoms,which has laser-induced high stress.However,under high fluence and short pulse width laser irradiation,the main damage mechanism is nucleation fracture caused by stress wave reflection and superposition at the bottom of the film.In addition,Mg[0001] has higher pressure sensitivity and is more prone to ablation than Mg[0001].The evolution equation of crater depth is established using multi-pulse laser ablation simulation and verified by experiments.The results show that,under multiple pulsed laser irradiation,not only does the crater depth increase linearly with the pulse number,but also the quadratic term and constant term of the fitted crater profile curve increase linearly.
基金Project supported by the National Natural Science Foundation of China(Nos.12072209,U21A2043012192211)+1 种基金the Natural Science Foundation of Hebei Province of China(No.A2020210009)the S&T Program of Hebei Province of China(No.225676162GH)。
文摘This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and other surface of the PSC half-space are assumed to be electrically insulating.By the Hankel integral transformation,the problem is reduced to the Fredholm integral equation of the second kind.This equation is solved numerically to obtain the indentation behaviors of the PSC half-space,mainly including the indentation force-depth relation and the electric potential-depth relation.The results show that the effect of the semiconductor property on the indentation responses is limited within a certain range of variation of the steady carrier concentration.The dependence of indentation behavior on material properties is also analyzed by two different kinds of PSCs.Finite element simulations are conducted to verify the results calculated by the integral equation technique,and good agreement is demonstrated.
基金Project supported by the National Natural Science Foundation of China (Nos. 52075070 and12302254)the Dalian City Supports Innovation and Entrepreneurship Projects for High-Level Talents (No. 2021RD16)the Liaoning Revitalization Talents Program (No. XLYC2002108)。
文摘The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is necessary to accurately predict the shakedown domains of these materials. The static shakedown theorem, also known as Melan's theorem, is a fundamental method used to predict the shakedown domains of structures and materials. Within this method, a key aspect lies in the construction and application of an appropriate self-equilibrium stress field(SSF). In the structural shakedown analysis, the SSF is typically constructed by governing equations that satisfy no external force(NEF) boundary conditions. However, we discover that directly applying these governing equations is not suitable for the shakedown analysis of heterogeneous materials. Researchers must consider the requirements imposed by the Hill-Mandel condition for boundary conditions and the physical significance of representative volume elements(RVEs). This paper addresses this issue and demonstrates that the sizes of SSFs vary under different boundary conditions, such as uniform displacement boundary conditions(DBCs), uniform traction boundary conditions(TBCs), and periodic boundary conditions(PBCs). As a result, significant discrepancies arise in the predicted shakedown domain sizes of heterogeneous materials. Built on the demonstrated relationship between SSFs under different boundary conditions, this study explores the conservative relationships among different shakedown domains, and provides proof of the relationship between the elastic limit(EL) factors and the shakedown loading factors under the loading domain of two load vertices. By utilizing numerical examples, we highlight the conservatism present in certain results reported in the existing literature. Among the investigated boundary conditions, the obtained shakedown domain is the most conservative under TBCs.Conversely, utilizing PBCs to construct an SSF for the shakedown analysis leads to less conservative lower bounds, indicating that PBCs should be employed as the preferred boundary conditions for the shakedown analysis of heterogeneous materials.
基金supported in part by the National Natural Science Foundation of China(52205424)in part by National Natural Science Foundation of China(T2125009,92048302)+2 种基金in part by Laoshan laboratory(Grant No.LSKJ202205300)in part by‘Pioneer’R&D Program of Zhejiang(Grant No.2023C03007)in part by the Zhejiang Provincial Natural Science Foundation of China(LY23A020001).
文摘Soft robot incarnates its unique advantages in deep-sea exploration,but grapples with high hydrostatic pressure’s unpredictable impact on its mechanical performances.In our previous work,a self-powered soft robot showed excellent work performance in the Mariana Trench at a depth of 11000 m,yet experienced notable degradation in deforming capability.Here,we propose a magnetic loading method for characterizing elastomer’s mechanical properties under extremely high hydrostatic pressure of up to 120 MPa.This method facilitates remote loading and enables in-situ observation,so that the dimensions and deformation at high hydrostatic pressure are obtained and used for calculations.The results reveal that the Young’s modulus of Polydimethylsiloxane(PDMS)monotonously increases with pressure.It is found that the relative increase in Young’s modulus is determined by its initial value,which is 8% for an initial Young’s modulus of 2200 kPa and 38% for 660 kPa.The relation between initial Young’s modulus and relevant increase can be fitted by an exponential function.The bulk modulus of PDMS is about 1.4 GPa at 20℃ and is barely affected by hydrostatic pressure.The method can quantify alterations in the mechanical properties of elastomers induced by hydrostatic pressure,and provide guidance for the design of soft robots which serve in extreme pressure environment.
基金Project supported by the National Natural Science Foundation of China (Nos.12072119,12325201,and 52205594)the China National Postdoctoral Program for Innovative Talents (No.BX20220118)。
文摘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.
基金the National Natural Science Foundation of China (Nos. 52388102, 52072317 and U2268210)the State Key Laboratory of Rail Transit Vehicle System (No. 2024RVL-T12)
文摘Hunting stability is an important performance criterion in railway vehicles.This study proposes an incorporation of a bio-inspired limb-like structure(LLS)-based nonlinear damping into the motor suspension system for traction units to improve the nonlinear critical speed and hunting stability of high-speed trains(HSTs).Initially,a vibration transmission analysis is conducted on a HST vehicle and a metro vehicle that suffered from hunting motion to explore the effect of different motor suspension systems from on-track tests.Subsequently,a simplified lateral dynamics model of an HST bogie is established to investigate the influence of the motor suspension on the bogie hunting behavior.The bifurcation analysis is applied to optimize the motor suspension parameters for high critical speed.Then,the nonlinear damping of the bio-inspired LLS,which has a positive correlation with the relative displacement,can further improve the modal damping of hunting motion and nonlinear critical speed compared with the linear motor suspension system.Furthermore,a comprehensive numerical model of a high-speed train,considering all nonlinearities,is established to investigate the influence of different types of motor suspension.The simulation results are well consistent with the theoretical analysis.The benefits of employing nonlinear damping of the bio-inspired LLS into the motor suspension of HSTs to enhance bogie hunting stability are thoroughly validated.
基金the financial support from National Natural Science Foundation of China (Grants Nos.12325203,91963117,and 11921002)。
文摘Recently,a Schwarz crystal structure with curved grain boundaries(GBs)constrained by twin-boundary(TB)networks was discovered in nanocrystalline Cu through experiments and atomistic simulations.Nanocrystalline Cu with nanosized Schwarz crystals exhibited high strength and excellent thermal stability.However,the grainsize effect and associated deformation mechanisms of Schwarz nanocrystals remain unknown.Here,we performed large-scale atomistic simulations to investigate the deformation behaviors and grain-size effect of nanocrystalline Cu with Schwarz crystals.Our simulations showed that similar to regular nanocrystals,Schwarz nanocrystals exhibit a strengthening-softening transition with decreasing grain size.The critical grain size in Schwarz nanocrystals is smaller than that in regular nanocrystals,leading to a maximum strength higher than that of regular nanocrystals.Our simulations revealed that the softening in Schwarz nanocrystals mainly originates from TB migration(or detwinning)and annihilation of GBs,rather than GB-mediated processes(including GB migration,sliding and diffusion)dominating the softening in regular nanocrystals.Quantitative analyses of simulation data further showed that compared with those in regular nanocrystals,the GB-mediated processes in Schwarz nanocrystals are suppressed,which is related to the low volume fraction of amorphous-like GBs and constraints of TB networks.The smaller critical grain size arises from the suppression of GB-mediated processes.
基金The authors gratefully acknowledge the support provided by the Postgraduate Research&Practice Program of Jiangsu Province(Grant No.KYCX18_0526)the Fundamental Research Funds for the Central Universities(Grant No.2018B682X14)Guangdong Basic and Applied Basic Research Foundation(No.2021A1515110807).
文摘In the context of global mean square error concerning the number of random variables in the representation,the Karhunen–Loève(KL)expansion is the optimal series expansion method for random field discretization.The computational efficiency and accuracy of the KL expansion are contingent upon the accurate resolution of the Fredholm integral eigenvalue problem(IEVP).The paper proposes an interpolation method based on different interpolation basis functions such as moving least squares(MLS),least squares(LS),and finite element method(FEM)to solve the IEVP.Compared with the Galerkin method based on finite element or Legendre polynomials,the main advantage of the interpolation method is that,in the calculation of eigenvalues and eigenfunctions in one-dimensional random fields,the integral matrix containing covariance function only requires a single integral,which is less than a two-folded integral by the Galerkin method.The effectiveness and computational efficiency of the proposed interpolation method are verified through various one-dimensional examples.Furthermore,based on theKL expansion and polynomial chaos expansion,the stochastic analysis of two-dimensional regular and irregular domains is conducted,and the basis function of the extended finite element method(XFEM)is introduced as the interpolation basis function in two-dimensional irregular domains to solve the IEVP.
基金supported by the University Natural Science Foundation of Jiangsu Province(Grant No.23KJB130004)the National Natural Science Foundation of China(Grant Nos.11932006,U1934206,12172121,12002118).
文摘This study proposes a comprehensive,coupled thermomechanical model that replaces local spatial derivatives in classical differential thermomechanical equations with nonlocal integral forms derived from the peridynamic differential operator(PDDO),eliminating the need for calibration procedures.The model employs a multi-rate explicit time integration scheme to handle varying time scales in multi-physics systems.Through simulations conducted on granite and ceramic materials,this model demonstrates its effectiveness.It successfully simulates thermal damage behavior in granite arising from incompatible mineral expansion and accurately calculates thermal crack propagation in ceramic slabs during quenching.To account for material heterogeneity,the model utilizes the Shuffle algorithm andWeibull distribution,yielding results that align with numerical simulations and experimental observations.This coupled thermomechanical model shows great promise for analyzing intricate thermomechanical phenomena in brittle materials.
基金supported by the National Natural Science Foundation of China(Grant No.12272276)the Fundamental Research Funds for the Central Universities(Grant No.2042023kf0194)the Shenzhen Science and Technology Program(Grant No.JCYJ20220530140606013)。
文摘Inspired by the excellent adhesion performances of setae structure from organisms,micro/nano-pillar array has become one of the paradigms for adhesive surfaces.The micropillar arrays are composed of the resin pillars for adhesion and the substrate with different elastic modulus for supporting.The stress singularity at the bi-material corner between the pillars and the substrate can induce the failure of the micropillar-substrate corner and further hinder the fabrication and application of micropillar arrays,yet the design for the stability of the micropillar array lacks systematical and quantitative guides.In this work,we develop a semi-analytical method to provide the full expressions for the stress distribution within the bi-material corner combining analytical derivations and numerical calculations.The predictions for the stress within the singularity field can be obtained based on the full expressions of the stress.The good agreement between the predictions and the FEM results demonstrates the high reliability of our method.By adopting the strain energy density factor approach,the stability of the pillar-substrate corner is assessed by predicting the failure at the corner.For the elastic mismatch between the pillar and substrate given in this paper,the stability can be improved by increasing the ratio of the shear modulus of the substrate to that of the micropillar.Our study provides accurate predictions for the stress distribution at the bi-material corner and can guide the optimization of material combinations of the pillars and the substrate for more stable bioinspired dry adhesives.
基金Project supported by Guangdong Basic and Applied Basic Research Foundation(Grant Nos.2022A1515010967 and 2023A1515012821)the National Natural Science Foundation of China(Grant Nos.12002272 and 12272293)Opening Project of Applied Mechanics and Structure Safety Key Laboratory of Sichuan Province(Grant No.SZDKF-202101)。
文摘Due to technical limitations,existing vibration isolation and energy harvesting(VIEH)devices have poor performance at low frequency.This paper proposes a new multilink-spring mechanism(MLSM)that can be used to solve this problem.The VIEH performance of the MLSM under harmonic excitation and Gaussian white noise was analyzed.It was found that the MLSM has good vibration isolation performance for low-frequency isolation and the frequency band can be widened by adjusting parameters to achieve a higher energy harvesting power.By comparison with two special cases,the results show that the MLSM is basically the same as the other two oscillators in terms of vibration isolation but has better energy harvesting performance under multistable characteristics.The MLSM is expected to reduce the impact of vibration on high-precision sensitive equipment in some special sites such as subways and mines,and at the same time supply power to structural health monitoring devices.