The macroscopic mechanical properties of rocks are significantly influenced by their microstructure.As a material bonded by mineral grains,the grain morphology of crystalline rock is the primary factor influencing the...The macroscopic mechanical properties of rocks are significantly influenced by their microstructure.As a material bonded by mineral grains,the grain morphology of crystalline rock is the primary factor influencing the strength.However,most strength criteria neglect the strength variations caused by different grain characteristics in rocks.Furthermore,the traditional linear criteria tend to overestimate tensile strength and exhibit apex singularity.To address these shortcomings,a piecewise strength criterion that considers the grain size effect has been proposed.A part of an ellipse was employed to construct the envelope of the tensive-shear region on the meridian plane,to accurately reproduce the low tensile-compressive strength ratio.Based on the analysis of experimental data,both linear and exponential modification functions that account for grain size effects were integrated into the proposed criterion.The corresponding finite element algorithm has been implemented.The accuracy and applicability of the proposed criterion were validated by comparing with the experimental data.展开更多
Permanent magnet linear synchronous motor(PMLSM)has the advantages of high thrust density and good control accuracy,which can be applied in high-power and high-speed occasions.In this paper,the analytical models are e...Permanent magnet linear synchronous motor(PMLSM)has the advantages of high thrust density and good control accuracy,which can be applied in high-power and high-speed occasions.In this paper,the analytical models are established to obtain the electromagnetic performance for the PMLSMs with dual secondaries and dual primaries.The air-gap flux density and the electromagnetic thrust are also obtained by the finite element model to verify theoretical analysis.Besides,an improved structure is also put forward in order to suppress the thrust fluctuation of the PMLSM.Finally,the advantages and disadvantages of two PMLSMs topologies are listed.These analyses would provide a guide for the design of PMLSMs applied in high-power and high-speed occasions.展开更多
Due to the large rotor eddy current loss and low thermal conductivity of carbon fiber sleeve,the high temperature usually occurs in high speed permanent magnet machines(HSPMMs)at the rated operation condition,resultin...Due to the large rotor eddy current loss and low thermal conductivity of carbon fiber sleeve,the high temperature usually occurs in high speed permanent magnet machines(HSPMMs)at the rated operation condition,resulting in irreversible demagnetization of the permanent magnet(PM).To obtain low rotor temperature,two novel rotor structures with low rotor eddy current loss are proposed in this paper.With the output torque and air gap flux density unchanged,the performance of HSPMMs with the two proposed rotor structures are analyzed based on finite element algorithm(FEA),including eddy current loss and temperature.Finally,the appropriate parameters of the proposed rotor structures are selected,and the electromagnetic(EM)performance,rotor stress and temperature are compared with those of the conventional rotor structure.Index Terms-Eddy current loss,finite element algorithm(FEA),electromagnetic(EM)performance,high speed permanent magnet machines(HSPMMs).展开更多
For the low-order finite element pair P1P1,based on full domain partition technique,a parallel pressure projection stabilized finite element algorithm for the Stokes equation with nonlinear slip boundary conditions is...For the low-order finite element pair P1P1,based on full domain partition technique,a parallel pressure projection stabilized finite element algorithm for the Stokes equation with nonlinear slip boundary conditions is designed and analyzed.From the definition of the subdifferential,the variational formulation of this equation is the variational inequality problem of the second kind.Each subproblem is a global problem on the composite grid,which is easy to program and implement.The optimal error estimates of the approximate solutions are obtained by theoretical analysis since the appropriate stabilization parameter is chosen.Finally,some numerical results are given to demonstrate the hight efficiency of the parallel stabilized finite element algorithm.展开更多
With advantages of strong drive capability,nested-loop secondary linear machine(NLS-LM)has great potentiality in linear metro.For its secondary structure with multiple loops,it is difficult to calculate the electromag...With advantages of strong drive capability,nested-loop secondary linear machine(NLS-LM)has great potentiality in linear metro.For its secondary structure with multiple loops,it is difficult to calculate the electromagnetic thrust of NLS-LM reasonably.Hence,in this paper,one thrust calculation method is proposed considering variable loop inductance and transient loop current.Firstly,to establish the secondary winding function,the modeling domain is confined to a limited range,and the equivalent loop span is employed by analyzing the coupling relationship between primary and secondary.Then,in order to obtain the secondary flux density,the transient secondary current is solved based on the loop impedance and induced voltage.Finally,the electromagnetic thrust can be calculated reasonably by the given primary current sheet and the calculated secondary flux density.Comprehensive simulations and experiments have demonstrated the effectiveness of the proposed method.展开更多
Uncontrolled coal fires are natural disasters that may cause mineral loss and environmental damage.The traditional loop source transient electromagnetic method can effectively detect the low-resistivity region of coal...Uncontrolled coal fires are natural disasters that may cause mineral loss and environmental damage.The traditional loop source transient electromagnetic method can effectively detect the low-resistivity region of coal fires,but its detection efficiency is not so good for high-resistivity regions.In view of this limitation,a technique based on electrical source transient electromagnetics is proposed in this paper to detect high-resistivity regions in the spontaneous combustion process of coal.Considering the complex geometry of the coal fire area,an unstructured tetrahedral grid is used in this study to realize the spatial discretization of the model,and solve the electromagnetic field based on a vector finite element algorithm.Numerical analysis is used to investigate methods for detecting coal fires and the characteristics of effective anomalies are further examined to provide guidance for practical detection.展开更多
Nested-loop secondary linear doubly-fed machine(NLS-LDFM) is a novel linear machine evolved from rotary brushless doubly-fed induction machine, which has a good application prospect in linear metro. In order to analyz...Nested-loop secondary linear doubly-fed machine(NLS-LDFM) is a novel linear machine evolved from rotary brushless doubly-fed induction machine, which has a good application prospect in linear metro. In order to analyze the performance of NLS-LDFM, the mechanism and action rules of end effects are investigated in this paper. Firstly, the mechanism of static and dynamic end effects is analyzed in aspect of direct coupling, winding asymmetry and transient secondary current. Furthermore, based on the winding theory for short primary linear machines, the machine parameters are established qualitatively considering pulsating magnetic field of NLS-LDFM. Finally, the NLS-LDFM performance analysis is supplemented by the finite element algorithm(FEA) simulation and experiments under different operating conditions.展开更多
An efficient and accurate solution algorithm was proposed for 1-D unsteady flow problems widely existing in hydraulic engineering. Based on the split-characteristic finite element method, the numerical model with the ...An efficient and accurate solution algorithm was proposed for 1-D unsteady flow problems widely existing in hydraulic engineering. Based on the split-characteristic finite element method, the numerical model with the Saint-Venant equations of 1-D unsteady flows was established. The assembled f＇mite element equations were solved with the tri-diagonal matrix algorithm. In the semi-implicit and explicit scheme, the critical time step of the method was dependent on the space step and flow velocity, not on the wave celerity. The method was used to eliminate the restriction due to the wave celerity for the computational analysis of unsteady open-channel flows. The model was verified by the experimental data and theoretical solution and also applied to the simulation of the flow in practical river networks. It shows that the numerical method has high efficiency and accuracy and can be used to simulate 1-D steady flows, and unsteady flows with shock waves or flood waves. Compared with other numerical methods, the algorithm of this method is simpler with higher accuracy, less dissipation, higher computation efficiency and less computer storage.展开更多
Based on the finite element method(FEM), some iterative methods related to different Reynolds numbers are designed and analyzed for solving the 2D/3D stationary incompressible magnetohydrodynamics(MHD) numerically. Tw...Based on the finite element method(FEM), some iterative methods related to different Reynolds numbers are designed and analyzed for solving the 2D/3D stationary incompressible magnetohydrodynamics(MHD) numerically. Two-level finite element iterative methods, consisting of the classical m-iteration methods on a coarse grid and corrections on a fine grid, are designed to solve the system at low Reynolds numbers under the strong uniqueness condition. One-level Oseen-type iterative method is investigated on a fine mesh at high Reynolds numbers under the weak uniqueness condition. Furthermore, the uniform stability and convergence of these methods with respect to equation parameters R_e, R_m, S_c, mesh sizes h, H and iterative step m are provided. Finally, the efficiency of the proposed methods is confirmed by numerical investigations.展开更多
Two algorithms of computing stress increment by using the elasto-plasticity constitutive model are firstly formulated, which are the Euler integration method and the radial return method. Hill＇48 anisotropic yield cr...Two algorithms of computing stress increment by using the elasto-plasticity constitutive model are firstly formulated, which are the Euler integration method and the radial return method. Hill＇48 anisotropic yield criterion is used. The Euler integration method can not obtain more accurate computation of the stress increment as the radial return method unless enough subintervals are taken,by which the Euler integration method will take excessive computing time. Without docreasing any accuracy, the radial return method can save much time. Finally, a square cup deep drawing from NUMISHEET＇93 benchmarks is simulated with a self-developed code SheetForrn in order to investigate the accuracy and efficiency of the radial return method.展开更多
A mesh-less Refined Integral Algorithm （RIA） of Boundary Element Method （BEM） is proposed to accurately solve the Helmholtz Integral Equation （HIE）.The convergence behavior and the practicability of the method a...A mesh-less Refined Integral Algorithm （RIA） of Boundary Element Method （BEM） is proposed to accurately solve the Helmholtz Integral Equation （HIE）.The convergence behavior and the practicability of the method are validated.Computational Fluid Dynamics （CFD）,Finite Element Method （FEM） and RIA are used to predict the propeller excited underwater noise of the submarine hull structure.Firstly the propeller and submarine＇s flows are independently validated,then the self propulsion of the ＂submarine＋propeller＂ system is simulated via CFD and the balanced point of the system is determined as well as the self propulsion factors.Secondly,the transient response of the ＂submarine＋ propeller＂ system is analyzed at the balanced point,and the propeller thrust and torque excitations are calculated.Thirdly the thrust and the torque excitations of the propeller are loaded on the submarine,respectively,to calculate the acoustic response,and the sound power and the main peak frequencies are obtained.Results show that：（1） the thrust mainly excites the submarine axial mode and the high frequency area appears at the two conical-type ends,while the torque mainly excites the circumferential mode and the high frequency area appears at the broadside of the cylindrical section,but with rather smaller sound power and radiation efficiency than the former,（2） the main sound source appears at BPF and 2BPF and comes from the harmonic propeller excitations.So,the main attention should be paid on the thrust excitation control for the sound reduction of the propeller excited submarine structure.展开更多
基金Project(2023YFC2907403)supported by the National Key R&D Program of ChinaProject(52074021)supported by the National Natural Science Foundation of China+1 种基金Project(2242045)supported by Beijing Natural Science Foundation,ChinaProject(ZD202216)supported by the Beijing Association of Higher Education,China。
文摘The macroscopic mechanical properties of rocks are significantly influenced by their microstructure.As a material bonded by mineral grains,the grain morphology of crystalline rock is the primary factor influencing the strength.However,most strength criteria neglect the strength variations caused by different grain characteristics in rocks.Furthermore,the traditional linear criteria tend to overestimate tensile strength and exhibit apex singularity.To address these shortcomings,a piecewise strength criterion that considers the grain size effect has been proposed.A part of an ellipse was employed to construct the envelope of the tensive-shear region on the meridian plane,to accurately reproduce the low tensile-compressive strength ratio.Based on the analysis of experimental data,both linear and exponential modification functions that account for grain size effects were integrated into the proposed criterion.The corresponding finite element algorithm has been implemented.The accuracy and applicability of the proposed criterion were validated by comparing with the experimental data.
基金This work has been partly supported in part by National Natural Science Foundation of China under Grants 51877093 and 51707079in part by the National Key Research and Development Program of China under Grant YS2018YFGH000200,in part by the Key Technical Innovation Program of Hubei Province under Grant 2019AAA026.
文摘Permanent magnet linear synchronous motor(PMLSM)has the advantages of high thrust density and good control accuracy,which can be applied in high-power and high-speed occasions.In this paper,the analytical models are established to obtain the electromagnetic performance for the PMLSMs with dual secondaries and dual primaries.The air-gap flux density and the electromagnetic thrust are also obtained by the finite element model to verify theoretical analysis.Besides,an improved structure is also put forward in order to suppress the thrust fluctuation of the PMLSM.Finally,the advantages and disadvantages of two PMLSMs topologies are listed.These analyses would provide a guide for the design of PMLSMs applied in high-power and high-speed occasions.
基金This work has been partly supported by the National Natural Science Foundation of China(NSFC 51877093,51807075)National Key Research and Development Program of China(Project ID:2018YFE0100200)+2 种基金Fundamental Research Funds for the Central Universities(2019kfyXMBZ031)Project funded by China Postdoctoral Science Foundation(2019M652640)Natural Science Basic Research Plan in Shaanxi Province of China under Grant 2018JQ5009 and Scientific Research Program Funded by Shaanxi Provincial Education Department of China under Grant No.18JK0398。
文摘Due to the large rotor eddy current loss and low thermal conductivity of carbon fiber sleeve,the high temperature usually occurs in high speed permanent magnet machines(HSPMMs)at the rated operation condition,resulting in irreversible demagnetization of the permanent magnet(PM).To obtain low rotor temperature,two novel rotor structures with low rotor eddy current loss are proposed in this paper.With the output torque and air gap flux density unchanged,the performance of HSPMMs with the two proposed rotor structures are analyzed based on finite element algorithm(FEA),including eddy current loss and temperature.Finally,the appropriate parameters of the proposed rotor structures are selected,and the electromagnetic(EM)performance,rotor stress and temperature are compared with those of the conventional rotor structure.Index Terms-Eddy current loss,finite element algorithm(FEA),electromagnetic(EM)performance,high speed permanent magnet machines(HSPMMs).
基金supported by the Natural Science Foundation of China(No.11361016)the Basic and Frontier Explore Program of Chongqing Municipality,China(No.cstc2018jcyjAX0305)Funds for the Central Universities(No.XDJK2018B032).
文摘For the low-order finite element pair P1P1,based on full domain partition technique,a parallel pressure projection stabilized finite element algorithm for the Stokes equation with nonlinear slip boundary conditions is designed and analyzed.From the definition of the subdifferential,the variational formulation of this equation is the variational inequality problem of the second kind.Each subproblem is a global problem on the composite grid,which is easy to program and implement.The optimal error estimates of the approximate solutions are obtained by theoretical analysis since the appropriate stabilization parameter is chosen.Finally,some numerical results are given to demonstrate the hight efficiency of the parallel stabilized finite element algorithm.
基金supported in part by the National Natural Science Foundation of China under Grants 52277050the Shenzhen International Collaboration under Grant GJHZ20210705142539007。
文摘With advantages of strong drive capability,nested-loop secondary linear machine(NLS-LM)has great potentiality in linear metro.For its secondary structure with multiple loops,it is difficult to calculate the electromagnetic thrust of NLS-LM reasonably.Hence,in this paper,one thrust calculation method is proposed considering variable loop inductance and transient loop current.Firstly,to establish the secondary winding function,the modeling domain is confined to a limited range,and the equivalent loop span is employed by analyzing the coupling relationship between primary and secondary.Then,in order to obtain the secondary flux density,the transient secondary current is solved based on the loop impedance and induced voltage.Finally,the electromagnetic thrust can be calculated reasonably by the given primary current sheet and the calculated secondary flux density.Comprehensive simulations and experiments have demonstrated the effectiveness of the proposed method.
基金Supported by National Key R&D Program of China(No.2018YFC0807900).
文摘Uncontrolled coal fires are natural disasters that may cause mineral loss and environmental damage.The traditional loop source transient electromagnetic method can effectively detect the low-resistivity region of coal fires,but its detection efficiency is not so good for high-resistivity regions.In view of this limitation,a technique based on electrical source transient electromagnetics is proposed in this paper to detect high-resistivity regions in the spontaneous combustion process of coal.Considering the complex geometry of the coal fire area,an unstructured tetrahedral grid is used in this study to realize the spatial discretization of the model,and solve the electromagnetic field based on a vector finite element algorithm.Numerical analysis is used to investigate methods for detecting coal fires and the characteristics of effective anomalies are further examined to provide guidance for practical detection.
基金supported in part by the National Natural Science Foundations of China under Grants 52277050 and 51877093the fund from Science,Technology,Shenzhen International Collaboration under Grant GJHZ20210705142539007+1 种基金the Key Research and Development Program of Sichuan Province under Grant 2021YFG0081the fund from Science,Technology and Innovation Commission of Shenzhen Municipality under Grant JCYJ20190809101205546。
文摘Nested-loop secondary linear doubly-fed machine(NLS-LDFM) is a novel linear machine evolved from rotary brushless doubly-fed induction machine, which has a good application prospect in linear metro. In order to analyze the performance of NLS-LDFM, the mechanism and action rules of end effects are investigated in this paper. Firstly, the mechanism of static and dynamic end effects is analyzed in aspect of direct coupling, winding asymmetry and transient secondary current. Furthermore, based on the winding theory for short primary linear machines, the machine parameters are established qualitatively considering pulsating magnetic field of NLS-LDFM. Finally, the NLS-LDFM performance analysis is supplemented by the finite element algorithm(FEA) simulation and experiments under different operating conditions.
基金Project supported by the National Nature Science Foundation of China (Grant No.50479068) the Program for New Century Excellent Talents in Universities (Grant No. NCET-04-0494).
文摘An efficient and accurate solution algorithm was proposed for 1-D unsteady flow problems widely existing in hydraulic engineering. Based on the split-characteristic finite element method, the numerical model with the Saint-Venant equations of 1-D unsteady flows was established. The assembled f＇mite element equations were solved with the tri-diagonal matrix algorithm. In the semi-implicit and explicit scheme, the critical time step of the method was dependent on the space step and flow velocity, not on the wave celerity. The method was used to eliminate the restriction due to the wave celerity for the computational analysis of unsteady open-channel flows. The model was verified by the experimental data and theoretical solution and also applied to the simulation of the flow in practical river networks. It shows that the numerical method has high efficiency and accuracy and can be used to simulate 1-D steady flows, and unsteady flows with shock waves or flood waves. Compared with other numerical methods, the algorithm of this method is simpler with higher accuracy, less dissipation, higher computation efficiency and less computer storage.
基金National Natural Science Foundation of China (Grant Nos. 11271298 and 11362021)
文摘Based on the finite element method(FEM), some iterative methods related to different Reynolds numbers are designed and analyzed for solving the 2D/3D stationary incompressible magnetohydrodynamics(MHD) numerically. Two-level finite element iterative methods, consisting of the classical m-iteration methods on a coarse grid and corrections on a fine grid, are designed to solve the system at low Reynolds numbers under the strong uniqueness condition. One-level Oseen-type iterative method is investigated on a fine mesh at high Reynolds numbers under the weak uniqueness condition. Furthermore, the uniform stability and convergence of these methods with respect to equation parameters R_e, R_m, S_c, mesh sizes h, H and iterative step m are provided. Finally, the efficiency of the proposed methods is confirmed by numerical investigations.
文摘Two algorithms of computing stress increment by using the elasto-plasticity constitutive model are firstly formulated, which are the Euler integration method and the radial return method. Hill＇48 anisotropic yield criterion is used. The Euler integration method can not obtain more accurate computation of the stress increment as the radial return method unless enough subintervals are taken,by which the Euler integration method will take excessive computing time. Without docreasing any accuracy, the radial return method can save much time. Finally, a square cup deep drawing from NUMISHEET＇93 benchmarks is simulated with a self-developed code SheetForrn in order to investigate the accuracy and efficiency of the radial return method.
文摘A mesh-less Refined Integral Algorithm （RIA） of Boundary Element Method （BEM） is proposed to accurately solve the Helmholtz Integral Equation （HIE）.The convergence behavior and the practicability of the method are validated.Computational Fluid Dynamics （CFD）,Finite Element Method （FEM） and RIA are used to predict the propeller excited underwater noise of the submarine hull structure.Firstly the propeller and submarine＇s flows are independently validated,then the self propulsion of the ＂submarine＋propeller＂ system is simulated via CFD and the balanced point of the system is determined as well as the self propulsion factors.Secondly,the transient response of the ＂submarine＋ propeller＂ system is analyzed at the balanced point,and the propeller thrust and torque excitations are calculated.Thirdly the thrust and the torque excitations of the propeller are loaded on the submarine,respectively,to calculate the acoustic response,and the sound power and the main peak frequencies are obtained.Results show that：（1） the thrust mainly excites the submarine axial mode and the high frequency area appears at the two conical-type ends,while the torque mainly excites the circumferential mode and the high frequency area appears at the broadside of the cylindrical section,but with rather smaller sound power and radiation efficiency than the former,（2） the main sound source appears at BPF and 2BPF and comes from the harmonic propeller excitations.So,the main attention should be paid on the thrust excitation control for the sound reduction of the propeller excited submarine structure.
