This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-d...This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.展开更多
By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is propose...By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.展开更多
Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and re- sistance to corrosion fatigue, cracking, etc. Co...Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and re- sistance to corrosion fatigue, cracking, etc. Compressive re- sidual stress and dent profile are important factors to eval- uate the effectiveness of shot peening process. In this pa- per, the influence of dimensionless parameters on maximum compressive residual stress and maximum depth of the dent were investigated. Firstly, dimensionless relations of pro- cessing parameters that affect the maximum compressive residual stress and the maximum depth of the dent were de- duced by dimensional analysis method. Secondly, the in- fluence of each dimensionless parameter on dimensionless variables was investigated by the finite element method. Fur- thermore, related empirical formulas were given for each di- mensionless parameter based on the simulation results. Fi- nally, comparison was made and good agreement was found between the simulation results and the empirical formula, which shows that a useful approach is provided in this pa- per for analyzing the influence of each individual parameter.展开更多
Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain met...Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain method and the panicle tracing method, such an algorithm is self-consistent and accurate since the interaction between electromagnetic fields and particles is properly modeled. In the time domain aspect, the generation of multipactor can be easily visualized, which makes it possible to gain a deeper insight into the physical mechanism of this effect. In addition to the classic secondary electron emission model, the measured practical secondary electron yield is used, which increases the accuracy of the algorithm. In order to validate the method, the impedance transformer and ridge waveguide filter are studied. By analyzing the evolution of the secondaries obtained by our method, multipactor thresholds of these components are estimated, which show good agreement with the experimental results. Furthermore, the most sensitive positions where multipactor occurs are determined from the phase focusing phenomenon, which is very meaningful for multipactor analysis and design.展开更多
A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow t...A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM （finite element method） software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.展开更多
The three-dimensional stress distributions in the area surrounding indentation pattern for three different materials, Al2O3, Si3N4 and SiC were analyzed by finite element method(FEM). Those theoretical results were al...The three-dimensional stress distributions in the area surrounding indentation pattern for three different materials, Al2O3, Si3N4 and SiC were analyzed by finite element method(FEM). Those theoretical results were also compared with the experimental ones by Rockwell hardness test. The effect of loading stress on the plastic deformation in specimens, surface was investigated on the assumption of shear strain energy theory by Huber-Mises when the materials were indented. The distributions of nomal stress, shear stress, and Mises stress were analysed with variations of loading conditions. It is clear that the analytical results for the stress distributions, the crack length and its density of probability are in good agreement with the experimental results.展开更多
The estimation of shear strength of rock mass discontinuity is always a focal, but difficult, problem in the field of geotechnical engineering. Considering the disadvantages and limitation of exist- ing estimation met...The estimation of shear strength of rock mass discontinuity is always a focal, but difficult, problem in the field of geotechnical engineering. Considering the disadvantages and limitation of exist- ing estimation methods, a new approach based on the shadow area percentage （SAP） that can be used to quantify surface roughness is proposed in this article. Firstly, by the help of laser scanning technique, the three-dimensional model of the surface of rock discontinuity was established. Secondly, a light source was simulated, and there would be some shadows produced on the model surface. Thirdly, to obtain the value of SAP of each specimen, the shadow detection technique was introduced for use. Fourthly, compared with the result from direct shear testing and based on statistics, an empirical for- mula was found among SAP, normal stress, and shear strength. Data of Yujian （~ River were used as an example, and the following conclusions have been made. （1） In the case of equal normal stress, the peak shear stress is positively proportional to the SAP. （2） The formula for estimating was derived, and the predictions of peak-shear strength made with this equation well agreed with the experimental re- suits obtained in laboratory tests.展开更多
A three-dimensional variational method is proposed to simultaneously retrieve the 3-D atmospheric temperature and moisture profiles from satellite radiance measurements. To include both vertical structure and the hori...A three-dimensional variational method is proposed to simultaneously retrieve the 3-D atmospheric temperature and moisture profiles from satellite radiance measurements. To include both vertical structure and the horizontal patterns of the atmospheric temperature and moisture, an EOF technique is used to decompose the temperature and moisture field in a 3-D space. A number of numerical simulations are conducted and they demonstrate that the 3-D method is less sensitive to the observation errors compared to the 1-D method. When the observation error is more than 2.0 K, to get the best results, the truncation number for the EOF＇s expansion have to be restricted to 2 in the 1-D method, while it can be set as large as 40 in a 3-D method. This results in the truncation error being reduced and the retrieval accuracy being improved in the 3-D method. Compared to the 1-D method, the rms errors of the 3-D method are reduced by 48% and 36% for the temperature and moisture retrievals, respectively. Using the real satellite measured brightness temperatures at 0557 UTC 31 July 2002, the temperature and moisture profiles are retrieved over a region （20°-45°N, 100°- 125°E） and compared with 37 collocated radiosonde observations. The results show that the retrieval accuracy with a 3-D method is significantly higher than those with the 1-D method.展开更多
This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial diffe...This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial differential equations of the above 2D problems is rewritten as an upper triangular differential system. For the associated operator matrix, the existence and the completeness of two normed orthogonal eigenfunction systems in some space are obtained, which belong to the two block operators arising in the operator matrix. Moreover, the general solution to the above 2D problem is given by the eigenfunction expansion method.展开更多
The variation principle is discussed and Rayleigh-Ritz method is proposed for construction of veloci ty field. A kinematically admissible velocity field based on polynomials was appIied to the determina tion of forgin...The variation principle is discussed and Rayleigh-Ritz method is proposed for construction of veloci ty field. A kinematically admissible velocity field based on polynomials was appIied to the determina tion of forging load and deformed buIge profile during upset forging of blocks. Simulation of upsetforging of rectangular blocks under various friction condjtions was performed. Comparison of the computed results with experiments and FEM shows good agreement. It is shown that this techniquecan be used for 3D simulation of metal forming process.展开更多
In this paper, the boundary layer equations (abbreviation BLE) for exterior flow around an obstacle are established using semi-geodesic coordinate system (S-coordinate) based on the curved two dimensional surface of t...In this paper, the boundary layer equations (abbreviation BLE) for exterior flow around an obstacle are established using semi-geodesic coordinate system (S-coordinate) based on the curved two dimensional surface of the obstacle. BLE are nonlinear partial differential equations on unknown normal viscous stress tensor and pressure on the obstacle and the existence of solution of BLE is proved. In addition a dimensional split method for dimensional three Navier-Stokes equations is established by applying several 2D-3C partial differential equations on two dimensional manifolds to approach 3D Navier-Stokes equations. The examples for the exterior flow around spheroid and ellipsoid are presents here.展开更多
This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite differenc...This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.展开更多
To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitr...To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation（BIE） representation solves the two-dimensional convected Helmholtz equation（CHE） and its fundamental solution, which must satisfy a new Sommerfeld radiation condition（SRC） in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green＇s kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.展开更多
The three-dimensional numerical manifold method（NMM） is studied on the basis of two-dimensional numerical manifold method. The three-dimensional cover displacement function is studied. The mechanical analysis and Ha...The three-dimensional numerical manifold method（NMM） is studied on the basis of two-dimensional numerical manifold method. The three-dimensional cover displacement function is studied. The mechanical analysis and Hammer integral method of three-dimensional numerical manifold method are put forward. The stiffness matrix of three-dimensional manifold element is derived and the dissection rules are given. The theoretical system and the numerical realizing method of three-dimensional numerical manifold method are systematically studied. As an example, the cantilever with load on the end is calculated, and the results show that the precision and efficiency are agreeable.展开更多
By introducing a more general auxiliary ordinary differential equation （ODE）, a modified variable separated ordinary differential equation method is presented for solving the （2 ＋ 1）-dimensional sine-Poisson equa...By introducing a more general auxiliary ordinary differential equation （ODE）, a modified variable separated ordinary differential equation method is presented for solving the （2 ＋ 1）-dimensional sine-Poisson equation. As a result, many explicit and exact solutions of the （2 ＋ 1）-dimensional sine-Poisson equation are derived in a simple manner by this technique.展开更多
Principal component analysis and generalized low rank approximation of matrices are two different dimensionality reduction methods. Two different dimensionality reduction algorithms are applied to the L1-CSVM model ba...Principal component analysis and generalized low rank approximation of matrices are two different dimensionality reduction methods. Two different dimensionality reduction algorithms are applied to the L1-CSVM model based on augmented Lagrange method to explore the variation of running time and accuracy of the model in dimensionality reduction space. The results show that the improved algorithm can greatly reduce the running time and improve the accuracy of the algorithm.展开更多
As the traditional forging process has many problems such as low efficiency, high consumption of material and energy, large cylindrical shell rolling is introduced. Large cylindrical shell rolling is a typical rotary ...As the traditional forging process has many problems such as low efficiency, high consumption of material and energy, large cylindrical shell rolling is introduced. Large cylindrical shell rolling is a typical rotary forming technology, and the upper and lower rolls have different radii and speeds. To quickly predict the three-dimensional stresses and eliminate fishtail defect, an improved strip layer method is developed, in which the asymmetry of the upper and lower rolls, non-uniform deformation and stress, as well as the asymmetrical spread on the end surface are considered. The deformation zone is divided into a certain number of layers and strips along the thickness and width, respectively. The transverse displacement model is constructed by polynomial function, in order to increase the computation speed greatly. From the metal plastic mechanics principle, the three-dimensional stress models are established. The genetic algorithm is used for optimization calculation in an industrial experiment example. The results show that the rolling pressure, the normal stresses, the upper and lower friction stress distributions are not similar with those of a general plate rolling. There are two relative maximum values in rolling pressure distribution. The upper and lower longitudinal friction stresses change direction nearby the upper and lower neutral points, respectively. The fishtail profile of spread on the end surface is predicted satisfactorily. The reduction could be helpful to eliminate fishtail defect. The large cylindrical shell rolling example illustrates the calculation results acquired rapidly are good agreements with the finite element simulation and experimental values of previous study. A highly effective and reliable three-dimensional simulation method is proposed for large cylindrical shell rolling and other asymmetrical rolling.展开更多
This paper ix devoted to establishment of the Chebyshev pseudospectral domain de-composition scheme for solving two-dimensional elliptic equation. By the generalized equivalent variatiunal form, we can get the stabili...This paper ix devoted to establishment of the Chebyshev pseudospectral domain de-composition scheme for solving two-dimensional elliptic equation. By the generalized equivalent variatiunal form, we can get the stability and convergence of this new scheme.展开更多
A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream su...A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream surface(curved surface)strip elements along metal flow traces,and the stream surface strip elements were mapped into the corresponding plane strip elements for analysis and computation.The longitudinal distributions of the lateral displacement and the altitudinal displacement of metal were respectively constructed to be a quartic curve and a quadratic curve,of which the lateral distributions were expressed as the third-power spline function,and the altitudinal distributions were fitted in the quadratic curve.From the flow theory of plastic mechanics,the mathematical models of the three-dimensional deformations and stresses of the deformation zone were constructed.Compared with the streamline strip element method proposed by the first author of this paper,the stream surface strip element method takes into account the uneven distributions of stresses and deformations along altitudinal direction,and realizes the precise three-dimensional analysis and computation.The simulation example of continuous hot rolled strip indicates that the method and the model accord with facts and provide a new reliable engineering-computation method for the three-dimensional mechanics simulation of plate and strip rolling process.展开更多
基金supported by the National Natural Science Foundation of China (Grants 11571223, 51404160)Shanxi Province Science Foundation for Youths (Grant 2014021025-1)
文摘This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.
基金supported by the Natural Science Foundation of Zhejiang Province,China(Grant Nos.LY20A010021,LY19A010002,LY20G030025)the Natural Science Founda-tion of Ningbo City,China(Grant Nos.2021J147,2021J235).
文摘By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.
基金supported by the National Natural Science Foun-dation of China (10972228,11002150,and 91016025)the Basic Research Equipment Project of Chinese Academy of Sciences(YZ200930)
文摘Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and re- sistance to corrosion fatigue, cracking, etc. Compressive re- sidual stress and dent profile are important factors to eval- uate the effectiveness of shot peening process. In this pa- per, the influence of dimensionless parameters on maximum compressive residual stress and maximum depth of the dent were investigated. Firstly, dimensionless relations of pro- cessing parameters that affect the maximum compressive residual stress and the maximum depth of the dent were de- duced by dimensional analysis method. Secondly, the in- fluence of each dimensionless parameter on dimensionless variables was investigated by the finite element method. Fur- thermore, related empirical formulas were given for each di- mensionless parameter based on the simulation results. Fi- nally, comparison was made and good agreement was found between the simulation results and the empirical formula, which shows that a useful approach is provided in this pa- per for analyzing the influence of each individual parameter.
基金Project supported by the National Key Laboratory Foundation,China(Grant No.9140C530103110C5301)
文摘Based on the particle-in-cell technology and the secondary electron emission theory, a three-dimensional simulation method for multipactor is presented in this paper. By combining the finite difference time domain method and the panicle tracing method, such an algorithm is self-consistent and accurate since the interaction between electromagnetic fields and particles is properly modeled. In the time domain aspect, the generation of multipactor can be easily visualized, which makes it possible to gain a deeper insight into the physical mechanism of this effect. In addition to the classic secondary electron emission model, the measured practical secondary electron yield is used, which increases the accuracy of the algorithm. In order to validate the method, the impedance transformer and ridge waveguide filter are studied. By analyzing the evolution of the secondaries obtained by our method, multipactor thresholds of these components are estimated, which show good agreement with the experimental results. Furthermore, the most sensitive positions where multipactor occurs are determined from the phase focusing phenomenon, which is very meaningful for multipactor analysis and design.
基金This work was supported by the National Natural Science Foundation of China （No. 50275094）.
文摘A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM （finite element method） software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.
文摘The three-dimensional stress distributions in the area surrounding indentation pattern for three different materials, Al2O3, Si3N4 and SiC were analyzed by finite element method(FEM). Those theoretical results were also compared with the experimental ones by Rockwell hardness test. The effect of loading stress on the plastic deformation in specimens, surface was investigated on the assumption of shear strain energy theory by Huber-Mises when the materials were indented. The distributions of nomal stress, shear stress, and Mises stress were analysed with variations of loading conditions. It is clear that the analytical results for the stress distributions, the crack length and its density of probability are in good agreement with the experimental results.
基金supported by the China Geological Survey (No.1212011014030)the Major State Basic Research Development Program of China (973 Program) (No.2011CB710600)
文摘The estimation of shear strength of rock mass discontinuity is always a focal, but difficult, problem in the field of geotechnical engineering. Considering the disadvantages and limitation of exist- ing estimation methods, a new approach based on the shadow area percentage （SAP） that can be used to quantify surface roughness is proposed in this article. Firstly, by the help of laser scanning technique, the three-dimensional model of the surface of rock discontinuity was established. Secondly, a light source was simulated, and there would be some shadows produced on the model surface. Thirdly, to obtain the value of SAP of each specimen, the shadow detection technique was introduced for use. Fourthly, compared with the result from direct shear testing and based on statistics, an empirical for- mula was found among SAP, normal stress, and shear strength. Data of Yujian （~ River were used as an example, and the following conclusions have been made. （1） In the case of equal normal stress, the peak shear stress is positively proportional to the SAP. （2） The formula for estimating was derived, and the predictions of peak-shear strength made with this equation well agreed with the experimental re- suits obtained in laboratory tests.
基金the 973 Program (Grant No. 2004CB418305)the National Natural Science Foundation of China(Grant No. 40575049).
文摘A three-dimensional variational method is proposed to simultaneously retrieve the 3-D atmospheric temperature and moisture profiles from satellite radiance measurements. To include both vertical structure and the horizontal patterns of the atmospheric temperature and moisture, an EOF technique is used to decompose the temperature and moisture field in a 3-D space. A number of numerical simulations are conducted and they demonstrate that the 3-D method is less sensitive to the observation errors compared to the 1-D method. When the observation error is more than 2.0 K, to get the best results, the truncation number for the EOF＇s expansion have to be restricted to 2 in the 1-D method, while it can be set as large as 40 in a 3-D method. This results in the truncation error being reduced and the retrieval accuracy being improved in the 3-D method. Compared to the 1-D method, the rms errors of the 3-D method are reduced by 48% and 36% for the temperature and moisture retrievals, respectively. Using the real satellite measured brightness temperatures at 0557 UTC 31 July 2002, the temperature and moisture profiles are retrieved over a region （20°-45°N, 100°- 125°E） and compared with 37 collocated radiosonde observations. The results show that the retrieval accuracy with a 3-D method is significantly higher than those with the 1-D method.
基金Project supported by the National Natural Science Foundation of China (No. 10962004)the Special-ized Research Fund for the Doctoral Program of Higher Education of China (No. 20070126002)+1 种基金the Chunhui Program of Ministry of Education of China (No. Z2009-1-01010)the Natural Science Foundation of Inner Mongolia (No. 2009BS0101)
文摘This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial differential equations of the above 2D problems is rewritten as an upper triangular differential system. For the associated operator matrix, the existence and the completeness of two normed orthogonal eigenfunction systems in some space are obtained, which belong to the two block operators arising in the operator matrix. Moreover, the general solution to the above 2D problem is given by the eigenfunction expansion method.
文摘The variation principle is discussed and Rayleigh-Ritz method is proposed for construction of veloci ty field. A kinematically admissible velocity field based on polynomials was appIied to the determina tion of forging load and deformed buIge profile during upset forging of blocks. Simulation of upsetforging of rectangular blocks under various friction condjtions was performed. Comparison of the computed results with experiments and FEM shows good agreement. It is shown that this techniquecan be used for 3D simulation of metal forming process.
文摘In this paper, the boundary layer equations (abbreviation BLE) for exterior flow around an obstacle are established using semi-geodesic coordinate system (S-coordinate) based on the curved two dimensional surface of the obstacle. BLE are nonlinear partial differential equations on unknown normal viscous stress tensor and pressure on the obstacle and the existence of solution of BLE is proved. In addition a dimensional split method for dimensional three Navier-Stokes equations is established by applying several 2D-3C partial differential equations on two dimensional manifolds to approach 3D Navier-Stokes equations. The examples for the exterior flow around spheroid and ellipsoid are presents here.
文摘This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.
基金supported by National Engineering School of Tunis (No.13039.1)
文摘To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation（BIE） representation solves the two-dimensional convected Helmholtz equation（CHE） and its fundamental solution, which must satisfy a new Sommerfeld radiation condition（SRC） in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green＇s kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.
文摘The three-dimensional numerical manifold method（NMM） is studied on the basis of two-dimensional numerical manifold method. The three-dimensional cover displacement function is studied. The mechanical analysis and Hammer integral method of three-dimensional numerical manifold method are put forward. The stiffness matrix of three-dimensional manifold element is derived and the dissection rules are given. The theoretical system and the numerical realizing method of three-dimensional numerical manifold method are systematically studied. As an example, the cantilever with load on the end is calculated, and the results show that the precision and efficiency are agreeable.
基金Project supported by China Postdoctoral Science Foundation （20100481488）, Key Fund Project of Advanced Research of the Weapon Equipment （9140A33040512JB3401）.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10672053)
文摘By introducing a more general auxiliary ordinary differential equation （ODE）, a modified variable separated ordinary differential equation method is presented for solving the （2 ＋ 1）-dimensional sine-Poisson equation. As a result, many explicit and exact solutions of the （2 ＋ 1）-dimensional sine-Poisson equation are derived in a simple manner by this technique.
文摘Principal component analysis and generalized low rank approximation of matrices are two different dimensionality reduction methods. Two different dimensionality reduction algorithms are applied to the L1-CSVM model based on augmented Lagrange method to explore the variation of running time and accuracy of the model in dimensionality reduction space. The results show that the improved algorithm can greatly reduce the running time and improve the accuracy of the algorithm.
基金Supported by National Science and Technology Major Project of China(Grant No.2011ZX04002-101)National Science and Technology Support Plan of China(Grant No.2011BAF15B02)National Natural Science Foundation of China(Grant No.51305388)
文摘As the traditional forging process has many problems such as low efficiency, high consumption of material and energy, large cylindrical shell rolling is introduced. Large cylindrical shell rolling is a typical rotary forming technology, and the upper and lower rolls have different radii and speeds. To quickly predict the three-dimensional stresses and eliminate fishtail defect, an improved strip layer method is developed, in which the asymmetry of the upper and lower rolls, non-uniform deformation and stress, as well as the asymmetrical spread on the end surface are considered. The deformation zone is divided into a certain number of layers and strips along the thickness and width, respectively. The transverse displacement model is constructed by polynomial function, in order to increase the computation speed greatly. From the metal plastic mechanics principle, the three-dimensional stress models are established. The genetic algorithm is used for optimization calculation in an industrial experiment example. The results show that the rolling pressure, the normal stresses, the upper and lower friction stress distributions are not similar with those of a general plate rolling. There are two relative maximum values in rolling pressure distribution. The upper and lower longitudinal friction stresses change direction nearby the upper and lower neutral points, respectively. The fishtail profile of spread on the end surface is predicted satisfactorily. The reduction could be helpful to eliminate fishtail defect. The large cylindrical shell rolling example illustrates the calculation results acquired rapidly are good agreements with the finite element simulation and experimental values of previous study. A highly effective and reliable three-dimensional simulation method is proposed for large cylindrical shell rolling and other asymmetrical rolling.
文摘This paper ix devoted to establishment of the Chebyshev pseudospectral domain de-composition scheme for solving two-dimensional elliptic equation. By the generalized equivalent variatiunal form, we can get the stability and convergence of this new scheme.
基金Sponsored by National Natural Science Foundation of China(50175095)Provincial Natural Science Foundation of Hebei of China(502173)
文摘A new method,the stream surface strip element method,for simulating the three-dimensional deformation of plate and strip rolling process was proposed.The rolling deformation zone was divided into a number of stream surface(curved surface)strip elements along metal flow traces,and the stream surface strip elements were mapped into the corresponding plane strip elements for analysis and computation.The longitudinal distributions of the lateral displacement and the altitudinal displacement of metal were respectively constructed to be a quartic curve and a quadratic curve,of which the lateral distributions were expressed as the third-power spline function,and the altitudinal distributions were fitted in the quadratic curve.From the flow theory of plastic mechanics,the mathematical models of the three-dimensional deformations and stresses of the deformation zone were constructed.Compared with the streamline strip element method proposed by the first author of this paper,the stream surface strip element method takes into account the uneven distributions of stresses and deformations along altitudinal direction,and realizes the precise three-dimensional analysis and computation.The simulation example of continuous hot rolled strip indicates that the method and the model accord with facts and provide a new reliable engineering-computation method for the three-dimensional mechanics simulation of plate and strip rolling process.
