Isogeometric analysis (IGA) is known to showadvanced features compared to traditional finite element approaches.Using IGA one may accurately obtain the geometrically nonlinear bending behavior of plates with functiona...Isogeometric analysis (IGA) is known to showadvanced features compared to traditional finite element approaches.Using IGA one may accurately obtain the geometrically nonlinear bending behavior of plates with functionalgrading (FG). However, the procedure is usually complex and often is time-consuming. We thus put forward adeep learning method to model the geometrically nonlinear bending behavior of FG plates, bypassing the complexIGA simulation process. A long bidirectional short-term memory (BLSTM) recurrent neural network is trainedusing the load and gradient index as inputs and the displacement responses as outputs. The nonlinear relationshipbetween the outputs and the inputs is constructed usingmachine learning so that the displacements can be directlyestimated by the deep learning network. To provide enough training data, we use S-FSDT Von-Karman IGA andobtain the displacement responses for different loads and gradient indexes. Results show that the recognition erroris low, and demonstrate the feasibility of deep learning technique as a fast and accurate alternative to IGA formodeling the geometrically nonlinear bending behavior of FG plates.展开更多
Geometrical configurations play a crucial role in dual-atom catalysts(DACs)for electrocatalytic applications.Significant progress has been made to design DACs electrocatalysts with various geometri-cal configurations,...Geometrical configurations play a crucial role in dual-atom catalysts(DACs)for electrocatalytic applications.Significant progress has been made to design DACs electrocatalysts with various geometri-cal configurations,but in-depth understanding the relationship between geometrical configurations and metal-metal interaction mechanisms for designing targeted DACs is still required.In this review,the recent progress in engineering of geometrical configurations of DACs is systematically summarized.Based on the polarity of geometrical configuration,DACs can be classified into two different types that are homonuclear and heteronuclear DACs.Furthermore,with regard to the geometrical configurations of the active sites,homonuclear DACs are identified into adjacent and bridged configurations,and heteronuclear DACs can be classified into adjacent,bridged,and separated configurations.Subsequently,metal-metal interactions in DACs with different geometrical configurations are introduced.Additionally,the applications of DACs in different electrocatalytic reactions are discussed,including the oxygen reduction reaction(ORR),oxygen evolution reaction(OER),hydrogen evolution reaction(HER),and other catalysis.Finally,the future challenges and perspectives for advancements in DACs are high-lighted.This review aims to provide inspiration for the design of highly effcient DACs towards energy relatedapplications.展开更多
We present a class of preconditioners for the linear systems resulting from a finite element or discontinuous Galerkin discretizations of advection-dominated problems.These preconditioners are designed to treat the ca...We present a class of preconditioners for the linear systems resulting from a finite element or discontinuous Galerkin discretizations of advection-dominated problems.These preconditioners are designed to treat the case of geometrically localized stiffness,where the convergence rates of iterative methods are degraded in a localized subregion of the mesh.Slower convergence may be caused by a number of factors,including the mesh size,anisotropy,highly variable coefficients,and more challenging physics.The approach taken in this work is to correct well-known preconditioners such as the block Jacobi and the block incomplete LU(ILU)with an adaptive inner subregion iteration.The goal of these preconditioners is to reduce the number of costly global iterations by accelerating the convergence in the stiff region by iterating on the less expensive reduced problem.The tolerance for the inner iteration is adaptively chosen to minimize subregion-local work while guaranteeing global convergence rates.We present analysis showing that the convergence of these preconditioners,even when combined with an adaptively selected tolerance,is independent of discretization parameters(e.g.,the mesh size and diffusion coefficient)in the subregion.We demonstrate significant performance improvements over black-box preconditioners when applied to several model convection-diffusion problems.Finally,we present performance results of several variations of iterative subregion correction preconditioners applied to the Reynolds number 2.25×10^(6)fluid flow over the NACA 0012 airfoil,as well as massively separated flow at 30°angle of attack.展开更多
In situ tensile testing in a scanning electron microscope(SEM)in conjunction with high-resolution electron backscatter diffraction(HR-EBSD)under load was used to characterize the evolution of geometrically necessary d...In situ tensile testing in a scanning electron microscope(SEM)in conjunction with high-resolution electron backscatter diffraction(HR-EBSD)under load was used to characterize the evolution of geometrically necessary dislocation(GND)densities at individual grain boundaries as a function of applied strain in a polycrystalline Mg-4Al alloy.The increase in GND density was investigated at plastic strains of 0%,0.6%,2.2%,3.3% from the area including 76 grains and correlated with(i)geometric compatibility between slip systems across grain boundaries,and(ii)plastic incompatibility.We develop expressions for the grain boundary GND density evolution as a function of plastic strain and plastic incompatibility,from which uniaxial tensile stress-strain response of polycrystalline Mg-4Al are computed and compared with experimental measurement.The findings in this study contribute to understanding the mechanisms governing the strain hardening response of single-phase polycrystalline alloys and more reliable prediction of mechanical behaviors in diverse microstructures.展开更多
Modern additive manufacturing processes enable fabricating architected cellular materials of complex shape,which can be used for different purposes.Among them,lattice structures are increasingly used in applications r...Modern additive manufacturing processes enable fabricating architected cellular materials of complex shape,which can be used for different purposes.Among them,lattice structures are increasingly used in applications requiring a compromise among lightness and suited mechanical properties,like improved energy absorption capacity and specific stiffness-to-weight and strength-to-weight ratios.A dedicated modeling strategy to assess the energy absorption capacity of lattice structures under uni-axial compression loading is presented in this work.The numerical model is developed in a non-linear framework accounting for the strain rate effect on the mechanical responses of the lattice structure.Four geometries,i.e.,cubic body centered cell,octet cell,rhombic-dodecahedron and truncated cuboctahedron 2+,are investigated.Specifically,the influence of the relative density of the representative volume element of each geometry,the strain-rate dependency of the bulk material and of the presence of the manufacturing process-induced geometrical imperfections on the energy absorption capacity of the lattice structure is investigated.The main outcome of this study points out the importance of correctly integrating geometrical imperfections into the modeling strategy when shock absorption applications are aimed for.展开更多
Ion cyclotron resonance heating(ICRH)stands out as a widely utilized and cost-effective auxiliary method for plasma heating,bearing significant importance in achieving high-performance discharges in p-^(11)B plasmas.I...Ion cyclotron resonance heating(ICRH)stands out as a widely utilized and cost-effective auxiliary method for plasma heating,bearing significant importance in achieving high-performance discharges in p-^(11)B plasmas.In light of the specific context of p-^(11)B plasma in the EHL-2 device,we conducted a comprehensive scan of the fundamental physical parameters of the antenna using the full-wave simulation program TORIC.Our preliminary result indicated that for p-^(11)B plasma,optimal ion heating parameters include a frequency of 40 MHz,with a high toroidal mode number like N_(?)=28 to heat the majority H ions.In addition,we discussed the impact of concentration of minority ion species on ion cyclotron resonance heating when^(11)B serves as the heavy minority species.The significant difference in charge-to-mass ratio between boron and hydrogen ions results in a considerable distance between the hybrid resonance layer and the tow inverted cyclotron resonance layer,necessitating a quite low boron ion concentration to achieve effective minority heating.We also considered another method of direct heating of hydrogen ions in the presence of boron ion minority.It is found that at appropriate boron ion concentrations(X(^(11)B)~17%),the position of the hybrid resonance layer approaches that of the hydrogen ion cyclotron resonance layer,thereby altering the polarization at this position and significantly enhancing hydrogen ion fundamental absorption.展开更多
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.展开更多
The study of Electromagnetic Compatibility is essential to ensure the harmonious operation of electronic equipment in a shared environment. The basic principles of Electromagnetic Compatibility focus on the ability of...The study of Electromagnetic Compatibility is essential to ensure the harmonious operation of electronic equipment in a shared environment. The basic principles of Electromagnetic Compatibility focus on the ability of devices to withstand electromagnetic disturbances and not produce disturbances that could affect other systems. Imperceptible in most work situations, electromagnetic fields can, beyond certain thresholds, have effects on human health. The objective of the present article is focused on the modeling analysis of the influence of geometric parameters of industrial static converters radiated electromagnetic fields using Maxwell’s equations. To do this we used the analytical formalism for calculating the electromagnetic field emitted by a filiform conductor, to model the electromagnetic radiation of this device in the spatio-temporal domain. The interactions of electromagnetic waves with human bodies are complex and depend on several factors linked to the characteristics of the incident wave. To model these interactions, we implemented the physical laws of electromagnetic wave propagation based on Maxwell’s and bio-heat equations to obtain consistent results. These obtained models allowed us to evaluate the spatial profile of induced current and temperature of biological tissue during exposure to electromagnetic waves generated by this system. The simulation 2D results obtained from computer tools show that the temperature variation and current induced by the electromagnetic field can have a very significant influence on the life of biological tissue. The paper provides a comprehensive analysis using advanced mathematical models to evaluate the influence of electromagnetic fields. The findings have direct implications for workplace safety, potentially influencing standards and regulations concerning electromagnetic exposure in industrial settings.展开更多
Due to the importance of vibration effects on the functional accuracy of mechanical systems,this research aims to develop a precise model of a nonlinearly vibrating single-link mobile flexible manipulator.The manipula...Due to the importance of vibration effects on the functional accuracy of mechanical systems,this research aims to develop a precise model of a nonlinearly vibrating single-link mobile flexible manipulator.The manipulator consists of an elastic arm,a rotary motor,and a rigid carrier,and undergoes general in-plane rigid body motion along with elastic transverse deformation.To accurately model the elastic behavior,Timoshenko’s beam theory is used to describe the flexible arm,which accounts for rotary inertia and shear deformation effects.By applying Newton’s second law,the nonlinear governing equations of motion for the manipulator are derived as a coupled system of ordinary differential equations(ODEs)and partial differential equations(PDEs).Then,the assumed mode method(AMM)is used to solve this nonlinear system of governing equations with appropriate shape functions.The assumed modes can be obtained after solving the characteristic equation of a Timoshenko beam with clamped boundary conditions at one end and an attached mass/inertia at the other.In addition,the effect of the transverse vibration of the inextensible arm on its axial behavior is investigated.Despite the axial rigidity,the effect makes the rigid body dynamics invalid for the axial behavior of the arm.Finally,numerical simulations are conducted to evaluate the performance of the developed model,and the results are compared with those obtained by the finite element approach.The comparison confirms the validity of the proposed dynamic model for the system.According to the mentioned features,this model can be reliable for investigating the system’s vibrational behavior and implementing vibration control algorithms.展开更多
Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes ...Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.展开更多
A novel approach to compute the high frequency radar cross-section (RCS) of complex targets is described in this paper.From the three views or the sectional views of the target, target is geometrically modeled by non-...A novel approach to compute the high frequency radar cross-section (RCS) of complex targets is described in this paper.From the three views or the sectional views of the target, target is geometrically modeled by non-uniform rational B-spline (NURBS) parametric surfaces using the software CNFEOV developed by oneself which constructs NURBS representation of complex target from engineering orthographic views. RCS is obtained through PO, PTD, MEC and IBC techniques. When calculating RCS of the target, it is necessary to get the unit normal vector to surface illumi- nated by radar and the value Z which is the distance from the point on the surface to radar. ln this novel approach, the unit normal vector to the surface can be obtained either by the Phong rendering model, in which the color components (RGB) of every pixel on the image are equal to the coordinate components of the normal, or by the NURBS expressions. The value Z can be achieved by software or hardware Z-buffer. The effects of the size of image on the RCS of target are discussed and the correct method is recommended. The RCS of the perfect conducting sphere, cylinder and dihedral as well as the coated cylinder, as some examples, are computed. The accuracy of the method is verified by comparing the numerical results with those obtained by using other methods.展开更多
A novel method for modeling cellular materials is proposed based on MATLAB image processing and synchrotron X-ray computed tomography scan- ning to obtain an accurate calculation result of aluminum foam based on finit...A novel method for modeling cellular materials is proposed based on MATLAB image processing and synchrotron X-ray computed tomography scan- ning to obtain an accurate calculation result of aluminum foam based on finite element model. The maximum entropy algorithm is employed to obtain the bina- rization image, and the median filtering algorithm is used to reduce the noise after binarization. The external contour and internal pores boundary is extracted by the "edge" function in MATLAB, and the geometrical model is reconstructed. A two-step mesh algorithm is adopted to mesh the reconstructed geometrical model. Accordingly, the finite element model of aluminum foam is established by the proposed method based on reconstruction geometrical model. The compression behavior of aluminum foam is obtained at 25℃, 100℃, 200℃ by ABAQUS, and good agreements with experiments are achieved by applying the present recon- struction algorithm and modeling method.展开更多
This study presents a high-speed geometrically nonlinear flutter analysis calculation method based on the highprecision computational fluid dynamics/computational structural dynamics methods.In the proposed method,the...This study presents a high-speed geometrically nonlinear flutter analysis calculation method based on the highprecision computational fluid dynamics/computational structural dynamics methods.In the proposed method,the aerodynamic simulation was conducted based on computational fluid dynamics,and the structural model was established using the nonlinear finite element model and tangential stiffness matrix.First,the equilibrium position was obtained using the nonlinear static aeroelastic iteration.Second,the structural modal under a steady aerodynamic load was extracted.Finally,the generalized displacement time curve was obtained by coupling the unsteady aerodynamics and linearized structure motion equations.Moreover,if the flutter is not at a critical state,the incoming flow dynamic pressure needs to be changed,and the above steps must be repeated until the vibration amplitude are equal.Furthermore,the high-speed geometrically nonlinear flutter of the wing-body assemblymodel with a high-aspect ratio was investigated,and the correctness of the method was verified using high-speed wind tunnel experiments.The results showed that the geometric nonlinearity of the large deformation of the wing caused in-plane bending to become a key factor in flutter characteristics and significantly decreased the dynamic pressure and frequency of the nonlinear flutter compared to those of the linear flutter.展开更多
Auxetic materials are cellular materials with a unique property of negative Poisson’s ratio.The auxeticity and performance of these metamaterials utterly depend on the geometrical parameters and loading direction.For...Auxetic materials are cellular materials with a unique property of negative Poisson’s ratio.The auxeticity and performance of these metamaterials utterly depend on the geometrical parameters and loading direction.For the first time,the quasi-static uniaxial compression performance of fused filament fabricated re-entrant diamond auxetic metamaterial is evaluated in the x-direction(in-plane)and z-direction(out-of-plane).The most commonly used thermoplastic feedstock,Acrylonitrile butadiene styrene,is considered a material of choice.The effect of influential geometrical parameters of the re-entrant diamond structure and printing parameter is systematically studied using Taguchi’s design of experiments.Grey-based multi-objective optimisation technique has been adopted to arrive at the optimal structure.Efforts are made to improve the stiffness and strength of the structure with fibre reinforcements.Micro glass fibre reinforcements have enhanced specific strength and stiffness in both in-plane and out-ofplane directions.A sevenfold and thirteen times increase in specific strength and energy absorption is evident for glass fibre-reinforced structures in out-of-plane directions compared to in-plane ones.Proper tuning of geometrical parameters of the re-entrant diamond structure can result in a Poisson’s ratio of up to-3.49 when tested in the x-direction.The parametric study has illustrated the tailorability of the structure according to the application requirements.The statistical study has signified each considered parameter’s contribution to the compression performance characteristics of the 3D printed re-entrant diamond auxetic metamaterial.展开更多
The Kuqa fold-and-thrust belt exhibits apparent structural variation in the western and eastern zone.Two salt layer act as effective decollements and influence the varied deformation.In this study,detailed seismic int...The Kuqa fold-and-thrust belt exhibits apparent structural variation in the western and eastern zone.Two salt layer act as effective decollements and influence the varied deformation.In this study,detailed seismic interpretations and analog modeling are presented to construct the suprasalt and subsalt structures in the transfer zone of the middle Kuqa and investigate the influence of the two salt layers.The results reveal that the relationship of the two salt layers changes from separated to connected,and then overlapped toward the foreland in the transfer zone.Different structural models are formed in the suprasalt and subsalt units due to the interaction of the two salt layers.The imbricate thrust faults form two broom-like fault systems in the subsalt units.The suprasalt units develop detached folds terminating toward the east in the region near the orogenic belt.Whereas,two offset anticlines with different trends develop at the frontal edge of the lower salt layer and the trailing edge of the upper salt layer,respectively.According to exploration results in this region,the relationship between suprasalt and subsalt structures has an influence on hydrocarbon accumulation.We believe that the connected deformation contains high-risk plays while the decoupled deformation contains well-preserved plays.展开更多
The perturbed Riemann problem for a hyperbolic system of conservation laws arising in geometrical optics with three constant initial states is solved.By studying the interactions among of the delta-shock,vacuum,and co...The perturbed Riemann problem for a hyperbolic system of conservation laws arising in geometrical optics with three constant initial states is solved.By studying the interactions among of the delta-shock,vacuum,and contact discontinuity,fourteen kinds of structures of Riemann solutions are obtained.The compound wave solutions consisting of delta-shocks,vacuums,and contact discontinuities are found.The single and double closed vacuum cavitations develop in solutions.Furthermore,it is shown that the solutions of the Riemann problem for the geometrical optics system are stable under certain perturbation of the initial data.Finally,the numerical results completely coinciding with theoretical analysis are presented.展开更多
Although most pipes are restrained by retaining clips in aircraft,the influence of the clip parameters on the vibration of the fluid-conveying pipe has not been revealed.By considering the clip width,a new dynamic mod...Although most pipes are restrained by retaining clips in aircraft,the influence of the clip parameters on the vibration of the fluid-conveying pipe has not been revealed.By considering the clip width,a new dynamic model of a fluid-conveying pipe restrained by an intermediate clip is established in this paper.To demonstrate the necessity of the proposed model,a half pipe model is established by modeling the clip as one end.By comparing the two models,it is found that the half pipe model overestimates the critical velocity and may estimate the dynamical behavior of the pipe incorrectly.In addition,with the increase in the clip stiffness,the conversion processes of the first two modes of the pipe are shown.Furthermore,by ignoring the width of the clip,the effect of the flow velocity on the accuracy of a concentrated restraint clip model is presented.When the flow velocity is close to the critical velocity,the accuracy of the concentrated restraint clip model significantly reduces,especially when the width of the clip is large.In general,the contribution of this paper is to establish a dynamic model of the fluid-conveying pipe which can describe the influence of the clip parameters,and to demonstrate the necessity of this model.展开更多
Degradation and overstress failures occur in many electronic systems in which the operation load and environmental conditions are complex.The dependency of them called dependent competing failure process(DCFP),has bee...Degradation and overstress failures occur in many electronic systems in which the operation load and environmental conditions are complex.The dependency of them called dependent competing failure process(DCFP),has been widely studied.Electronic system may experience mutual effects of degradation and shocks,they are considered to be interdependent.Both the degradation and the shock processes will decrease the limit of system and cause cumulative effect.Finally,the competition of hard and soft failure will cause the system failure.Based on the failure mechanism accumulation theory,this paper constructs the shock-degradation acceleration and the threshold descent model,and a system reliability model established by using these two models.The mutually DCFP effect of electronic system interaction has been decomposed into physical correlation of failure,including acceleration,accumulation and competition.As a case,a reliability of electronic system in aeronautical system has been analyzed with the proposed method.The method proposed is based on failure physical evaluation,and could provide important reference for quantitative evaluation and design improvement of the newly designed system in case of data deficiency.展开更多
基金the National Natural Science Foundation of China(NSFC)under Grant Nos.12272124 and 11972146.
文摘Isogeometric analysis (IGA) is known to showadvanced features compared to traditional finite element approaches.Using IGA one may accurately obtain the geometrically nonlinear bending behavior of plates with functionalgrading (FG). However, the procedure is usually complex and often is time-consuming. We thus put forward adeep learning method to model the geometrically nonlinear bending behavior of FG plates, bypassing the complexIGA simulation process. A long bidirectional short-term memory (BLSTM) recurrent neural network is trainedusing the load and gradient index as inputs and the displacement responses as outputs. The nonlinear relationshipbetween the outputs and the inputs is constructed usingmachine learning so that the displacements can be directlyestimated by the deep learning network. To provide enough training data, we use S-FSDT Von-Karman IGA andobtain the displacement responses for different loads and gradient indexes. Results show that the recognition erroris low, and demonstrate the feasibility of deep learning technique as a fast and accurate alternative to IGA formodeling the geometrically nonlinear bending behavior of FG plates.
基金supported by the Natural Science Foundation of China (22179062,52125202,and U2004209)the Natural Science Foundation of Jiangsu Province (BK20230035)+1 种基金the Fundamental Research Funds for the Central Universities (30922010303)the Intergovernmental Cooperation Projects in the National Key Research and Development Plan of the Ministry of Science and Technology of PRC (2022YFE0196800)
文摘Geometrical configurations play a crucial role in dual-atom catalysts(DACs)for electrocatalytic applications.Significant progress has been made to design DACs electrocatalysts with various geometri-cal configurations,but in-depth understanding the relationship between geometrical configurations and metal-metal interaction mechanisms for designing targeted DACs is still required.In this review,the recent progress in engineering of geometrical configurations of DACs is systematically summarized.Based on the polarity of geometrical configuration,DACs can be classified into two different types that are homonuclear and heteronuclear DACs.Furthermore,with regard to the geometrical configurations of the active sites,homonuclear DACs are identified into adjacent and bridged configurations,and heteronuclear DACs can be classified into adjacent,bridged,and separated configurations.Subsequently,metal-metal interactions in DACs with different geometrical configurations are introduced.Additionally,the applications of DACs in different electrocatalytic reactions are discussed,including the oxygen reduction reaction(ORR),oxygen evolution reaction(OER),hydrogen evolution reaction(HER),and other catalysis.Finally,the future challenges and perspectives for advancements in DACs are high-lighted.This review aims to provide inspiration for the design of highly effcient DACs towards energy relatedapplications.
文摘We present a class of preconditioners for the linear systems resulting from a finite element or discontinuous Galerkin discretizations of advection-dominated problems.These preconditioners are designed to treat the case of geometrically localized stiffness,where the convergence rates of iterative methods are degraded in a localized subregion of the mesh.Slower convergence may be caused by a number of factors,including the mesh size,anisotropy,highly variable coefficients,and more challenging physics.The approach taken in this work is to correct well-known preconditioners such as the block Jacobi and the block incomplete LU(ILU)with an adaptive inner subregion iteration.The goal of these preconditioners is to reduce the number of costly global iterations by accelerating the convergence in the stiff region by iterating on the less expensive reduced problem.The tolerance for the inner iteration is adaptively chosen to minimize subregion-local work while guaranteeing global convergence rates.We present analysis showing that the convergence of these preconditioners,even when combined with an adaptively selected tolerance,is independent of discretization parameters(e.g.,the mesh size and diffusion coefficient)in the subregion.We demonstrate significant performance improvements over black-box preconditioners when applied to several model convection-diffusion problems.Finally,we present performance results of several variations of iterative subregion correction preconditioners applied to the Reynolds number 2.25×10^(6)fluid flow over the NACA 0012 airfoil,as well as massively separated flow at 30°angle of attack.
基金supported by the U.S.Department of Energy,Office of Basic Energy Sciences,Division of Materials Sciences and Engineering under Award#DE-SC0008637 as part of the Center for PRedictive Integrated Materials Science(PRISMS Center)at the University of Michigan。
文摘In situ tensile testing in a scanning electron microscope(SEM)in conjunction with high-resolution electron backscatter diffraction(HR-EBSD)under load was used to characterize the evolution of geometrically necessary dislocation(GND)densities at individual grain boundaries as a function of applied strain in a polycrystalline Mg-4Al alloy.The increase in GND density was investigated at plastic strains of 0%,0.6%,2.2%,3.3% from the area including 76 grains and correlated with(i)geometric compatibility between slip systems across grain boundaries,and(ii)plastic incompatibility.We develop expressions for the grain boundary GND density evolution as a function of plastic strain and plastic incompatibility,from which uniaxial tensile stress-strain response of polycrystalline Mg-4Al are computed and compared with experimental measurement.The findings in this study contribute to understanding the mechanisms governing the strain hardening response of single-phase polycrystalline alloys and more reliable prediction of mechanical behaviors in diverse microstructures.
文摘Modern additive manufacturing processes enable fabricating architected cellular materials of complex shape,which can be used for different purposes.Among them,lattice structures are increasingly used in applications requiring a compromise among lightness and suited mechanical properties,like improved energy absorption capacity and specific stiffness-to-weight and strength-to-weight ratios.A dedicated modeling strategy to assess the energy absorption capacity of lattice structures under uni-axial compression loading is presented in this work.The numerical model is developed in a non-linear framework accounting for the strain rate effect on the mechanical responses of the lattice structure.Four geometries,i.e.,cubic body centered cell,octet cell,rhombic-dodecahedron and truncated cuboctahedron 2+,are investigated.Specifically,the influence of the relative density of the representative volume element of each geometry,the strain-rate dependency of the bulk material and of the presence of the manufacturing process-induced geometrical imperfections on the energy absorption capacity of the lattice structure is investigated.The main outcome of this study points out the importance of correctly integrating geometrical imperfections into the modeling strategy when shock absorption applications are aimed for.
基金supported by Shenzhen Municipal Collaborative Innovation Technology Program-International Science and Technology(S&T)Cooperation Project(No.GJHZ20220913142609017)Shenzhen Science and Technology Innovation Commission Key Technical Project(No.JSGG20210713091539014)+1 种基金Ling Chuang Research Project of China National Nuclear Corporationthe“Fourteen Five-Year Plan”Basic Technological Research Project(No.JSZL2022XXXX001)。
文摘Ion cyclotron resonance heating(ICRH)stands out as a widely utilized and cost-effective auxiliary method for plasma heating,bearing significant importance in achieving high-performance discharges in p-^(11)B plasmas.In light of the specific context of p-^(11)B plasma in the EHL-2 device,we conducted a comprehensive scan of the fundamental physical parameters of the antenna using the full-wave simulation program TORIC.Our preliminary result indicated that for p-^(11)B plasma,optimal ion heating parameters include a frequency of 40 MHz,with a high toroidal mode number like N_(?)=28 to heat the majority H ions.In addition,we discussed the impact of concentration of minority ion species on ion cyclotron resonance heating when^(11)B serves as the heavy minority species.The significant difference in charge-to-mass ratio between boron and hydrogen ions results in a considerable distance between the hybrid resonance layer and the tow inverted cyclotron resonance layer,necessitating a quite low boron ion concentration to achieve effective minority heating.We also considered another method of direct heating of hydrogen ions in the presence of boron ion minority.It is found that at appropriate boron ion concentrations(X(^(11)B)~17%),the position of the hybrid resonance layer approaches that of the hydrogen ion cyclotron resonance layer,thereby altering the polarization at this position and significantly enhancing hydrogen ion fundamental absorption.
基金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 study of Electromagnetic Compatibility is essential to ensure the harmonious operation of electronic equipment in a shared environment. The basic principles of Electromagnetic Compatibility focus on the ability of devices to withstand electromagnetic disturbances and not produce disturbances that could affect other systems. Imperceptible in most work situations, electromagnetic fields can, beyond certain thresholds, have effects on human health. The objective of the present article is focused on the modeling analysis of the influence of geometric parameters of industrial static converters radiated electromagnetic fields using Maxwell’s equations. To do this we used the analytical formalism for calculating the electromagnetic field emitted by a filiform conductor, to model the electromagnetic radiation of this device in the spatio-temporal domain. The interactions of electromagnetic waves with human bodies are complex and depend on several factors linked to the characteristics of the incident wave. To model these interactions, we implemented the physical laws of electromagnetic wave propagation based on Maxwell’s and bio-heat equations to obtain consistent results. These obtained models allowed us to evaluate the spatial profile of induced current and temperature of biological tissue during exposure to electromagnetic waves generated by this system. The simulation 2D results obtained from computer tools show that the temperature variation and current induced by the electromagnetic field can have a very significant influence on the life of biological tissue. The paper provides a comprehensive analysis using advanced mathematical models to evaluate the influence of electromagnetic fields. The findings have direct implications for workplace safety, potentially influencing standards and regulations concerning electromagnetic exposure in industrial settings.
文摘Due to the importance of vibration effects on the functional accuracy of mechanical systems,this research aims to develop a precise model of a nonlinearly vibrating single-link mobile flexible manipulator.The manipulator consists of an elastic arm,a rotary motor,and a rigid carrier,and undergoes general in-plane rigid body motion along with elastic transverse deformation.To accurately model the elastic behavior,Timoshenko’s beam theory is used to describe the flexible arm,which accounts for rotary inertia and shear deformation effects.By applying Newton’s second law,the nonlinear governing equations of motion for the manipulator are derived as a coupled system of ordinary differential equations(ODEs)and partial differential equations(PDEs).Then,the assumed mode method(AMM)is used to solve this nonlinear system of governing equations with appropriate shape functions.The assumed modes can be obtained after solving the characteristic equation of a Timoshenko beam with clamped boundary conditions at one end and an attached mass/inertia at the other.In addition,the effect of the transverse vibration of the inextensible arm on its axial behavior is investigated.Despite the axial rigidity,the effect makes the rigid body dynamics invalid for the axial behavior of the arm.Finally,numerical simulations are conducted to evaluate the performance of the developed model,and the results are compared with those obtained by the finite element approach.The comparison confirms the validity of the proposed dynamic model for the system.According to the mentioned features,this model can be reliable for investigating the system’s vibrational behavior and implementing vibration control algorithms.
基金funded by a grant from Natural Sciences and Engineering Research Council of Canada (NSERC)the authors would like to acknowledge the Niobec mine (Saint-Honoré, QuébecQuébec)
文摘Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.
文摘A novel approach to compute the high frequency radar cross-section (RCS) of complex targets is described in this paper.From the three views or the sectional views of the target, target is geometrically modeled by non-uniform rational B-spline (NURBS) parametric surfaces using the software CNFEOV developed by oneself which constructs NURBS representation of complex target from engineering orthographic views. RCS is obtained through PO, PTD, MEC and IBC techniques. When calculating RCS of the target, it is necessary to get the unit normal vector to surface illumi- nated by radar and the value Z which is the distance from the point on the surface to radar. ln this novel approach, the unit normal vector to the surface can be obtained either by the Phong rendering model, in which the color components (RGB) of every pixel on the image are equal to the coordinate components of the normal, or by the NURBS expressions. The value Z can be achieved by software or hardware Z-buffer. The effects of the size of image on the RCS of target are discussed and the correct method is recommended. The RCS of the perfect conducting sphere, cylinder and dihedral as well as the coated cylinder, as some examples, are computed. The accuracy of the method is verified by comparing the numerical results with those obtained by using other methods.
基金supported by the National Natural Science Foundations of China(11202007,11232001,and 91216301)
文摘A novel method for modeling cellular materials is proposed based on MATLAB image processing and synchrotron X-ray computed tomography scan- ning to obtain an accurate calculation result of aluminum foam based on finite element model. The maximum entropy algorithm is employed to obtain the bina- rization image, and the median filtering algorithm is used to reduce the noise after binarization. The external contour and internal pores boundary is extracted by the "edge" function in MATLAB, and the geometrical model is reconstructed. A two-step mesh algorithm is adopted to mesh the reconstructed geometrical model. Accordingly, the finite element model of aluminum foam is established by the proposed method based on reconstruction geometrical model. The compression behavior of aluminum foam is obtained at 25℃, 100℃, 200℃ by ABAQUS, and good agreements with experiments are achieved by applying the present recon- struction algorithm and modeling method.
文摘This study presents a high-speed geometrically nonlinear flutter analysis calculation method based on the highprecision computational fluid dynamics/computational structural dynamics methods.In the proposed method,the aerodynamic simulation was conducted based on computational fluid dynamics,and the structural model was established using the nonlinear finite element model and tangential stiffness matrix.First,the equilibrium position was obtained using the nonlinear static aeroelastic iteration.Second,the structural modal under a steady aerodynamic load was extracted.Finally,the generalized displacement time curve was obtained by coupling the unsteady aerodynamics and linearized structure motion equations.Moreover,if the flutter is not at a critical state,the incoming flow dynamic pressure needs to be changed,and the above steps must be repeated until the vibration amplitude are equal.Furthermore,the high-speed geometrically nonlinear flutter of the wing-body assemblymodel with a high-aspect ratio was investigated,and the correctness of the method was verified using high-speed wind tunnel experiments.The results showed that the geometric nonlinearity of the large deformation of the wing caused in-plane bending to become a key factor in flutter characteristics and significantly decreased the dynamic pressure and frequency of the nonlinear flutter compared to those of the linear flutter.
文摘Auxetic materials are cellular materials with a unique property of negative Poisson’s ratio.The auxeticity and performance of these metamaterials utterly depend on the geometrical parameters and loading direction.For the first time,the quasi-static uniaxial compression performance of fused filament fabricated re-entrant diamond auxetic metamaterial is evaluated in the x-direction(in-plane)and z-direction(out-of-plane).The most commonly used thermoplastic feedstock,Acrylonitrile butadiene styrene,is considered a material of choice.The effect of influential geometrical parameters of the re-entrant diamond structure and printing parameter is systematically studied using Taguchi’s design of experiments.Grey-based multi-objective optimisation technique has been adopted to arrive at the optimal structure.Efforts are made to improve the stiffness and strength of the structure with fibre reinforcements.Micro glass fibre reinforcements have enhanced specific strength and stiffness in both in-plane and out-ofplane directions.A sevenfold and thirteen times increase in specific strength and energy absorption is evident for glass fibre-reinforced structures in out-of-plane directions compared to in-plane ones.Proper tuning of geometrical parameters of the re-entrant diamond structure can result in a Poisson’s ratio of up to-3.49 when tested in the x-direction.The parametric study has illustrated the tailorability of the structure according to the application requirements.The statistical study has signified each considered parameter’s contribution to the compression performance characteristics of the 3D printed re-entrant diamond auxetic metamaterial.
基金supported by the National Natural Science Foundation of China(Grant Nos.41572187,41972219,41927802 and 42072320)the China Postdoctoral Science Foundation(Grant No.2020M671432)。
文摘The Kuqa fold-and-thrust belt exhibits apparent structural variation in the western and eastern zone.Two salt layer act as effective decollements and influence the varied deformation.In this study,detailed seismic interpretations and analog modeling are presented to construct the suprasalt and subsalt structures in the transfer zone of the middle Kuqa and investigate the influence of the two salt layers.The results reveal that the relationship of the two salt layers changes from separated to connected,and then overlapped toward the foreland in the transfer zone.Different structural models are formed in the suprasalt and subsalt units due to the interaction of the two salt layers.The imbricate thrust faults form two broom-like fault systems in the subsalt units.The suprasalt units develop detached folds terminating toward the east in the region near the orogenic belt.Whereas,two offset anticlines with different trends develop at the frontal edge of the lower salt layer and the trailing edge of the upper salt layer,respectively.According to exploration results in this region,the relationship between suprasalt and subsalt structures has an influence on hydrocarbon accumulation.We believe that the connected deformation contains high-risk plays while the decoupled deformation contains well-preserved plays.
基金supported by the National Natural Science Foundation of China(12061084)the Natural Science Foundation of Yunnan Province(2019FY003007).
文摘The perturbed Riemann problem for a hyperbolic system of conservation laws arising in geometrical optics with three constant initial states is solved.By studying the interactions among of the delta-shock,vacuum,and contact discontinuity,fourteen kinds of structures of Riemann solutions are obtained.The compound wave solutions consisting of delta-shocks,vacuums,and contact discontinuities are found.The single and double closed vacuum cavitations develop in solutions.Furthermore,it is shown that the solutions of the Riemann problem for the geometrical optics system are stable under certain perturbation of the initial data.Finally,the numerical results completely coinciding with theoretical analysis are presented.
基金supported by the National Science Funds for Distinguished Young Scholars of China(No.12025204)the Shanghai Municipal Education Commission of China(No.2019-01-07-00-09-E00018)。
文摘Although most pipes are restrained by retaining clips in aircraft,the influence of the clip parameters on the vibration of the fluid-conveying pipe has not been revealed.By considering the clip width,a new dynamic model of a fluid-conveying pipe restrained by an intermediate clip is established in this paper.To demonstrate the necessity of the proposed model,a half pipe model is established by modeling the clip as one end.By comparing the two models,it is found that the half pipe model overestimates the critical velocity and may estimate the dynamical behavior of the pipe incorrectly.In addition,with the increase in the clip stiffness,the conversion processes of the first two modes of the pipe are shown.Furthermore,by ignoring the width of the clip,the effect of the flow velocity on the accuracy of a concentrated restraint clip model is presented.When the flow velocity is close to the critical velocity,the accuracy of the concentrated restraint clip model significantly reduces,especially when the width of the clip is large.In general,the contribution of this paper is to establish a dynamic model of the fluid-conveying pipe which can describe the influence of the clip parameters,and to demonstrate the necessity of this model.
基金supported by the National Natural Science Foundation of China(61503014,62073009)。
文摘Degradation and overstress failures occur in many electronic systems in which the operation load and environmental conditions are complex.The dependency of them called dependent competing failure process(DCFP),has been widely studied.Electronic system may experience mutual effects of degradation and shocks,they are considered to be interdependent.Both the degradation and the shock processes will decrease the limit of system and cause cumulative effect.Finally,the competition of hard and soft failure will cause the system failure.Based on the failure mechanism accumulation theory,this paper constructs the shock-degradation acceleration and the threshold descent model,and a system reliability model established by using these two models.The mutually DCFP effect of electronic system interaction has been decomposed into physical correlation of failure,including acceleration,accumulation and competition.As a case,a reliability of electronic system in aeronautical system has been analyzed with the proposed method.The method proposed is based on failure physical evaluation,and could provide important reference for quantitative evaluation and design improvement of the newly designed system in case of data deficiency.