Vertical rigidity of the space self adaptive 530 high rigidity mill is calculated by applying the boundary element method (BEM) of three dimension elastic contact problem,which can update the existed deforming s...Vertical rigidity of the space self adaptive 530 high rigidity mill is calculated by applying the boundary element method (BEM) of three dimension elastic contact problem,which can update the existed deforming separation calculating theory and corresponding methods of material mechanics,elastic mechanics and finite element method.The method has less hypotheses and stronger synthesis in contact type calculating model.The advantages of the method are high calculating rate,high calculating accuracy,etc..展开更多
The attitude optimal control problem (OCP) of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the ...The attitude optimal control problem (OCP) of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a dis- crete nonlinear programming (NLP) problem using the Chebyshev-Gauss pseudospectral method (CGPM). Solutions of the NLP problem can be obtained using the sequential quadratic programming (SQP) algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss (CG) points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform (FFT). To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and con- trol variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.展开更多
This publication is a revised version of the previous article. Seismic rigidity method despite its widespread use is the object of harsh criticism from scientists who oppose it to the methodology and results of seismo...This publication is a revised version of the previous article. Seismic rigidity method despite its widespread use is the object of harsh criticism from scientists who oppose it to the methodology and results of seismological registration of earthquakes and microseisms. The article substantiates the original approach based on the solution of the direct problem of seismic microzonation for the model of real soil thickness. A new formula of the seismic rigidity method is proposed, taking into account the lithological, hydrogeological and spectral features of the soil mass, as well as the position of the new seismic scale of the SSI. The formula was tested on the example of the correct description of the features of macroseismic effects on the territory of Leninakan at the Spitak earthquake in 1988. Linear estimates according to the formula of seismic rigidity in the seismic microzoning area represent changes in seismic intensity in the most contrast way. It is shown that the real estimates of seismic intensity under strong seismic effects (by I > VII degree) will not exceed those given by the formula of the seismic rigidity method.展开更多
Based on the theoretical model of rigidity correlation method, the study on application was carried out with Chinese Liyuan face rockfill dam as example. The linear relation equations between the rockfill rigidity and...Based on the theoretical model of rigidity correlation method, the study on application was carried out with Chinese Liyuan face rockfill dam as example. The linear relation equations between the rockfill rigidity and density measured by pit method were established, and the regression performance and accuracy of rigidity correlation method were analyzed by calculating the inversion values of density. The results show that the regression equations of rigidity correlation method are high significant so as to work out the rockfill density precisely;rigidity correlation method is used for density inversion of rockfill with minor error and namely high accuracy, which is proper with satisfactory results.展开更多
3-D rigid visco-plastic finite element method (FEM) is used in the analysisof metal forming processes, including strip and plate rolling, shape rolling, slab edging, specialstrip rolling. The shifted incomplete Choles...3-D rigid visco-plastic finite element method (FEM) is used in the analysisof metal forming processes, including strip and plate rolling, shape rolling, slab edging, specialstrip rolling. The shifted incomplete Cholesky decomposition of the stiffness matrix with thesolution of the equations for velocity increment by the conjugate gradient method is combined. Thistechnique, termed the shifted ICCG method, is then employed to solve the slab edging problem. Theperformance of this algorithm in terms of the number of iterations, friction variation, shiftedparameter psi and the results of simulation for processing parameters are analysed. Numerical testsand application of this technique verify the efficiency and stability of the shifted ICCG method inthe analysis of slab edging.展开更多
An optimal motion planning scheme based on the quasi-Newton method is proposed for a rigid spacecraft with two momentum wheels. A cost functional is introduced to incorporate the control energy, the final state errors...An optimal motion planning scheme based on the quasi-Newton method is proposed for a rigid spacecraft with two momentum wheels. A cost functional is introduced to incorporate the control energy, the final state errors and the constraints on states. The motion planning for determining control inputs to minimize the cost functional is formulated as a nonlinear optimal control problem. Using the control parametrization, one can transform the infinite dimensional optimal control problem to a finite dimensional one that is solved via the quasi-Newton methods for a feasible trajectory which satisfies the nonholonomic constraint. The optimal motion planning scheme was applied to a rigid spacecraft with two momentum wheels. The simulation results show the effectiveness of the proposed optimal motion planning scheme.展开更多
A conventional complex variable boundary integral equation (CVBIE) in plane elasticity is provided. After using the Somigliana identity between a particular fundamental stress field and a physical stress field, an a...A conventional complex variable boundary integral equation (CVBIE) in plane elasticity is provided. After using the Somigliana identity between a particular fundamental stress field and a physical stress field, an additional integral equality is obtained. By adding both sides of this integral equality to both sides of the conventional CVBIE, the amended boundary integral equation (BIE) is obtained. The method based on the discretization of the amended BIE is called the amended influence matrix method. With this method, for the Neumann boundary value problem (BVP) of an interior region, a unique solution for the displacement can be obtained. Several numerical examples are provided to prove the efficiency of the suggested method.展开更多
Applications of certain multi-parameter acceleration techniques used with themodified New-ton-Raphson (mN-R) methods to solve the nonlinear equations arising from rigid-plasticfinite element analysis are investigated....Applications of certain multi-parameter acceleration techniques used with themodified New-ton-Raphson (mN-R) methods to solve the nonlinear equations arising from rigid-plasticfinite element analysis are investigated. New modified multi-parameter techniques, developed fromCrisfield's multi-parameter methods, are utilized to solve these nonlinear equations. The numericalperformance of these techniques is compared with the standard Newton-Raphson method (sN-R),Crisfield's single parameter method (C1), Crisfield's two parameter method (C2) and Crisfield'sthree parameter method (C3). The new techniques do not involve additional residual force calculationand require little extra computational effort. In addition, they are more robust and efficient thanother existing acceleration techniques.展开更多
Common compliant joints generally have limited range of motion, reduced fatigue life and high stress concentration. To overcome these shortcomings, periodically corrugated cantilever beam is applied to design complian...Common compliant joints generally have limited range of motion, reduced fatigue life and high stress concentration. To overcome these shortcomings, periodically corrugated cantilever beam is applied to design compliant joints. Basic corrugated beam unit is modeled by using pseudo-rigid-body method. The trajectory and deformation behavior of periodically corrugated cantilever beam are estimated by the transformation of coordinate and superposition of the deformation of corrugated beam units. Finite element analysis(FEA) is carried out on corrugated cantilever beam to estimate the accuracy of the pseudo-rigid-body model. Results show that the kinetostatic behaviors obtained by this method, which has a relative error less than 6%, has good applicability and corrugated cantilever beam has the characteristics of a large range of motion and high mechanical strength. The corrugated cantilever beam is then applied to design a flexible rotational joint to obtain a larger angle output. The paper proposes a pseudo-rigid-body model for corrugated cantilever beam and designed a flexible rotational joint with large angle output.展开更多
According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the...According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the Thermo-Parameter Method(TPM)and Parametric Variational principles(PVP)are used to reduce the computational effort while maintaining the accuracy of solutions.A better solution isalso obtained in this paper.展开更多
The distinct element method(DEM)has been used successfully for the dynamic analysis of rigid block sys- tems.One of many difficulties associated with DEM is modeling of damping.In this paper,new procedures are propose...The distinct element method(DEM)has been used successfully for the dynamic analysis of rigid block sys- tems.One of many difficulties associated with DEM is modeling of damping.In this paper,new procedures are proposed for the damping modeling and its numerical implementation in distinct element analysis of rigid muhi-block systems.The stiff- ness proportional damping is constructed for the prescribed damping ratio,based on the non-zero fundamental frequency ef- fective during the time interval while the boundary conditions remain essentially constant.At this time interval,the funda- mental frequency can be estimated without complete eigenvalue analysis.The damping coefficients will vary while the damp- ing ratio remains the same throughout the entire analysis.A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation phases.These procedures were implemented in the development of the dis- tinet element method for the dynamic analyses of piled multi-block systems.The analysis results |or the single-block and two-block systems were in a good agreement with the analytic predictions.Applications to the seismic analyses of piled four- block systems revealed that the new procedures can make a significant difference and may lead to much-improved results.展开更多
The compressibility properties of systems consisting of generic rotating rigid triangles are analyzed and discussed. It is shown that these systems which are usually associated with auxeticity can exhibit strongly ani...The compressibility properties of systems consisting of generic rotating rigid triangles are analyzed and discussed. It is shown that these systems which are usually associated with auxeticity can exhibit strongly anisotropic properties for certain conformations, which may give rise to the anomalous property of negative linear compressibility (NLC), that is, the system with particular geometry will expand along one direction when loaded hydrostatically. It is also shown that through carefully choosing the geometric features (i.e. the dimensions and the alignment of the rotating triangles as well as the angles between them) and the direction along which the linear compressibility is measured, one may control the magnitude and range of the NLC. All this provides a novel but effective method of manufacturing the systems which can be tailored to achieve particular values of NLC to fit particular practical applications.展开更多
Springback is one of important factors influencing the forming quality of numerical control (NC) bending of thin-walled tube. In this paper, a numerical-analytic method for springback angle prediction of the process...Springback is one of important factors influencing the forming quality of numerical control (NC) bending of thin-walled tube. In this paper, a numerical-analytic method for springback angle prediction of the process was put forward. The method is based on springback angle model derived using analytic method and simulation results from three-dimensional (3D) rigid-plastic finite element method (FEM). The method is validated through comparison with experimental results. The features of the method are as follows: (1) The method is high in efficiency because it combines advantages of rigid-plastic FEM and analytic method. (2) The method is satisfactory in accuracy, since the field variables used in the model is resulting from 3D rigid-plastic FEM solution, and the effects both of axial force and strain neutral axis shift have been included. (3) Research on multi-factor effects can be carried out using the method due to its advantage inheriting from rigid-plastic FEM. The method described here is also of general significance to other bending processes.展开更多
We proposed a mesh-free method, the called node-based smoothed point interpolation method(NS-PIM),for dynamic analysis of rotating beams. A gradient smoothing technique is used, and the requirements on the consistence...We proposed a mesh-free method, the called node-based smoothed point interpolation method(NS-PIM),for dynamic analysis of rotating beams. A gradient smoothing technique is used, and the requirements on the consistence of the displacement functions are further weakened. In static problems, the beams with three types of boundary conditions are analyzed, and the results are compared with the exact solution, which shows the effectiveness of this method and can provide an upper bound solution for the deflection.This means that the NS-PIM makes the system soften. The NS-PIM is then further extended for solving a rigid-flexible coupled system dynamics problem, considering a rotating flexible cantilever beam. In this case, the rotating flexible cantilever beam considers not only the transverse deformations,but also the longitudinal deformations. The rigid-flexible coupled dynamic equations of the system are derived via employing Lagrange’s equations of the second type. Simulation results of the NS-PIM are compared with those obtained using finite element method(FEM) and assumed mode method. It is found that compared with FEM, the NS-PIM has anti-ill solving ability under the same calculation conditions.展开更多
The paper presents two methods for the formulation of free vibration analysis of collecting electrodes of precipitators.The first,called the hybrid finite element method, combines the finit element method used for cal...The paper presents two methods for the formulation of free vibration analysis of collecting electrodes of precipitators.The first,called the hybrid finite element method, combines the finit element method used for calculations of spring deformations with the rigid finite element method used to reflect mass and geometrical features,which is called the hybrid finite element method.As a result,a model with a diagonal mass matrix is obtained.Due to a specific geometry of the electrodes,which are long plates of complicated shapes,the second method proposed is the strip method which is a semi-analytical method.The strip method allows us to formulate the equations of motion with a considerably smaller number of generalized coordinates.Results of numerical calculations obtained by both methods are compared with those obtained using commercial software like ANSYS and ABAQUS.Good compatibility of results is achieved.展开更多
Previous work examined the effect of the attached stiffness matrix terms on stability of an elastic beam undergoing prescribed large overall motion. The aim of the present work is to extend the nonlinear formulations ...Previous work examined the effect of the attached stiffness matrix terms on stability of an elastic beam undergoing prescribed large overall motion. The aim of the present work is to extend the nonlinear formulations to an elastic beam with free large overall motion. Based on initial stress method, the nonlinear coupling equations of elastic beams are obtained with free large overall motion and the attached stiffness matrix is derived by solving sub-static formulation. The angular velocity and the tip deformation of the elastic pendulum are calculated. The analytical results show that the simulation results of the present model are tabled and coincide with the one-order approximate model. It is shown that the simulation results accord with energy conservation principle.展开更多
The problem of the elastic interaction between a screw dislocation and a three-phase circular inclusion with interracial rigid lines (anti-cracks) is investigated. An efficient and concise method for the complex mul...The problem of the elastic interaction between a screw dislocation and a three-phase circular inclusion with interracial rigid lines (anti-cracks) is investigated. An efficient and concise method for the complex multiply connected region is developed, with which explicit series form solutions of the complex potentials in the matrix, and the interphase layer and inclusion regions are derived. Based on the complex potentials, the image force on the screw dislocation is then calculated by using the Peach-Koehler formula. The equilibrium position of the dislocation is discussed in detail for various rigid line geometries, interphase layer thicknesses and material property combinations. The main results show that the interracial rigid lines exert a significant perturbation effect on the motion of the screw dislocation near the circular inclusion surrounded by an interphase layer.展开更多
To avoid mesh distortion and iterative remeshing in mesh-based numerical analysis,a meshless approach based on element free Galerkin (EFG) method is applied to the metal forming analysis of ring compression. Discrete ...To avoid mesh distortion and iterative remeshing in mesh-based numerical analysis,a meshless approach based on element free Galerkin (EFG) method is applied to the metal forming analysis of ring compression. Discrete equations are formulated upon the moving least-squares (MLS) approximation and modified Markov variational principles for rigid-plastic/ rigid-viscoplastic (RP/RVP) material models. The penalty function is used for the incompressible condition without volumetric locking. Based on the axisymmetric mechanical model,ring tests with different friction coefficients are studied. The deformed nodal configurations and shaded contours of equivalent strains are shown by developed meshless post processor. The comparison of meshless and finite element (FE) results validates the feasibility and accuracy for meshless method to simulate metal forming process.展开更多
文摘Vertical rigidity of the space self adaptive 530 high rigidity mill is calculated by applying the boundary element method (BEM) of three dimension elastic contact problem,which can update the existed deforming separation calculating theory and corresponding methods of material mechanics,elastic mechanics and finite element method.The method has less hypotheses and stronger synthesis in contact type calculating model.The advantages of the method are high calculating rate,high calculating accuracy,etc..
基金supported by the National Natural Science Foundation of China(No.11472058)
文摘The attitude optimal control problem (OCP) of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a dis- crete nonlinear programming (NLP) problem using the Chebyshev-Gauss pseudospectral method (CGPM). Solutions of the NLP problem can be obtained using the sequential quadratic programming (SQP) algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss (CG) points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform (FFT). To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and con- trol variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.
文摘This publication is a revised version of the previous article. Seismic rigidity method despite its widespread use is the object of harsh criticism from scientists who oppose it to the methodology and results of seismological registration of earthquakes and microseisms. The article substantiates the original approach based on the solution of the direct problem of seismic microzonation for the model of real soil thickness. A new formula of the seismic rigidity method is proposed, taking into account the lithological, hydrogeological and spectral features of the soil mass, as well as the position of the new seismic scale of the SSI. The formula was tested on the example of the correct description of the features of macroseismic effects on the territory of Leninakan at the Spitak earthquake in 1988. Linear estimates according to the formula of seismic rigidity in the seismic microzoning area represent changes in seismic intensity in the most contrast way. It is shown that the real estimates of seismic intensity under strong seismic effects (by I > VII degree) will not exceed those given by the formula of the seismic rigidity method.
文摘Based on the theoretical model of rigidity correlation method, the study on application was carried out with Chinese Liyuan face rockfill dam as example. The linear relation equations between the rockfill rigidity and density measured by pit method were established, and the regression performance and accuracy of rigidity correlation method were analyzed by calculating the inversion values of density. The results show that the regression equations of rigidity correlation method are high significant so as to work out the rockfill density precisely;rigidity correlation method is used for density inversion of rockfill with minor error and namely high accuracy, which is proper with satisfactory results.
基金supported by Huo Yingdong Young Teachers Foundation,Ministry of State Education of ChinaNational Natural Science Foundation of China(No.59904003).
文摘3-D rigid visco-plastic finite element method (FEM) is used in the analysisof metal forming processes, including strip and plate rolling, shape rolling, slab edging, specialstrip rolling. The shifted incomplete Cholesky decomposition of the stiffness matrix with thesolution of the equations for velocity increment by the conjugate gradient method is combined. Thistechnique, termed the shifted ICCG method, is then employed to solve the slab edging problem. Theperformance of this algorithm in terms of the number of iterations, friction variation, shiftedparameter psi and the results of simulation for processing parameters are analysed. Numerical testsand application of this technique verify the efficiency and stability of the shifted ICCG method inthe analysis of slab edging.
基金Project supported by the National Natural Science Foundation of China (No. 10372014).
文摘An optimal motion planning scheme based on the quasi-Newton method is proposed for a rigid spacecraft with two momentum wheels. A cost functional is introduced to incorporate the control energy, the final state errors and the constraints on states. The motion planning for determining control inputs to minimize the cost functional is formulated as a nonlinear optimal control problem. Using the control parametrization, one can transform the infinite dimensional optimal control problem to a finite dimensional one that is solved via the quasi-Newton methods for a feasible trajectory which satisfies the nonholonomic constraint. The optimal motion planning scheme was applied to a rigid spacecraft with two momentum wheels. The simulation results show the effectiveness of the proposed optimal motion planning scheme.
文摘A conventional complex variable boundary integral equation (CVBIE) in plane elasticity is provided. After using the Somigliana identity between a particular fundamental stress field and a physical stress field, an additional integral equality is obtained. By adding both sides of this integral equality to both sides of the conventional CVBIE, the amended boundary integral equation (BIE) is obtained. The method based on the discretization of the amended BIE is called the amended influence matrix method. With this method, for the Neumann boundary value problem (BVP) of an interior region, a unique solution for the displacement can be obtained. Several numerical examples are provided to prove the efficiency of the suggested method.
文摘Applications of certain multi-parameter acceleration techniques used with themodified New-ton-Raphson (mN-R) methods to solve the nonlinear equations arising from rigid-plasticfinite element analysis are investigated. New modified multi-parameter techniques, developed fromCrisfield's multi-parameter methods, are utilized to solve these nonlinear equations. The numericalperformance of these techniques is compared with the standard Newton-Raphson method (sN-R),Crisfield's single parameter method (C1), Crisfield's two parameter method (C2) and Crisfield'sthree parameter method (C3). The new techniques do not involve additional residual force calculationand require little extra computational effort. In addition, they are more robust and efficient thanother existing acceleration techniques.
基金supported by National Natural Science Foundation of China(Grant Nos.51205134,91223201)Doctoral Fund of Ministry of Education of China(Grant No.20120172120001)+2 种基金Research Project of State Key Laboratory of Mechanical System and Vibration of China(Grant No.MSV201405)Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme(GDUPS,2010)Fundamental Research Funds for the Central Universities(Grant No.2013ZM012)
文摘Common compliant joints generally have limited range of motion, reduced fatigue life and high stress concentration. To overcome these shortcomings, periodically corrugated cantilever beam is applied to design compliant joints. Basic corrugated beam unit is modeled by using pseudo-rigid-body method. The trajectory and deformation behavior of periodically corrugated cantilever beam are estimated by the transformation of coordinate and superposition of the deformation of corrugated beam units. Finite element analysis(FEA) is carried out on corrugated cantilever beam to estimate the accuracy of the pseudo-rigid-body model. Results show that the kinetostatic behaviors obtained by this method, which has a relative error less than 6%, has good applicability and corrugated cantilever beam has the characteristics of a large range of motion and high mechanical strength. The corrugated cantilever beam is then applied to design a flexible rotational joint to obtain a larger angle output. The paper proposes a pseudo-rigid-body model for corrugated cantilever beam and designed a flexible rotational joint with large angle output.
基金The project supported by National Natural Science Foundation of China
文摘According to the lower-bound theorem of limit analysis the Rigid Finite Element Meth-od(RFEM)is applied to structural limit analysis and the linear programmings for limit analysis are deducedin this paper.Moreover,the Thermo-Parameter Method(TPM)and Parametric Variational principles(PVP)are used to reduce the computational effort while maintaining the accuracy of solutions.A better solution isalso obtained in this paper.
文摘The distinct element method(DEM)has been used successfully for the dynamic analysis of rigid block sys- tems.One of many difficulties associated with DEM is modeling of damping.In this paper,new procedures are proposed for the damping modeling and its numerical implementation in distinct element analysis of rigid muhi-block systems.The stiff- ness proportional damping is constructed for the prescribed damping ratio,based on the non-zero fundamental frequency ef- fective during the time interval while the boundary conditions remain essentially constant.At this time interval,the funda- mental frequency can be estimated without complete eigenvalue analysis.The damping coefficients will vary while the damp- ing ratio remains the same throughout the entire analysis.A new numerical procedure is developed to prevent unnecessary energy loss that can occur during the separation phases.These procedures were implemented in the development of the dis- tinet element method for the dynamic analyses of piled multi-block systems.The analysis results |or the single-block and two-block systems were in a good agreement with the analytic predictions.Applications to the seismic analyses of piled four- block systems revealed that the new procedures can make a significant difference and may lead to much-improved results.
基金Project supported by the National Natural Science Foundation of China(Grant No.51475208)
文摘The compressibility properties of systems consisting of generic rotating rigid triangles are analyzed and discussed. It is shown that these systems which are usually associated with auxeticity can exhibit strongly anisotropic properties for certain conformations, which may give rise to the anomalous property of negative linear compressibility (NLC), that is, the system with particular geometry will expand along one direction when loaded hydrostatically. It is also shown that through carefully choosing the geometric features (i.e. the dimensions and the alignment of the rotating triangles as well as the angles between them) and the direction along which the linear compressibility is measured, one may control the magnitude and range of the NLC. All this provides a novel but effective method of manufacturing the systems which can be tailored to achieve particular values of NLC to fit particular practical applications.
基金This work was supported by the National Natural Science Foundation of China for Distinguished Young Scholars (Grant No. 50225518)the Teaching and Research Award Program for 0utstanding Young Teachers in Higher Education Institution of M0E, PRCthe Aeronautical Science Foundation of China (Grant No. 04H53057).
文摘Springback is one of important factors influencing the forming quality of numerical control (NC) bending of thin-walled tube. In this paper, a numerical-analytic method for springback angle prediction of the process was put forward. The method is based on springback angle model derived using analytic method and simulation results from three-dimensional (3D) rigid-plastic finite element method (FEM). The method is validated through comparison with experimental results. The features of the method are as follows: (1) The method is high in efficiency because it combines advantages of rigid-plastic FEM and analytic method. (2) The method is satisfactory in accuracy, since the field variables used in the model is resulting from 3D rigid-plastic FEM solution, and the effects both of axial force and strain neutral axis shift have been included. (3) Research on multi-factor effects can be carried out using the method due to its advantage inheriting from rigid-plastic FEM. The method described here is also of general significance to other bending processes.
基金the support from the National Natural Science Foundation of China (Grants 11272155, 11132007, and 11502113)the Fundamental Research Funds for Central Universities (Grant 30917011103)the China Scholarship Council for one year study at the University of Cincinnati
文摘We proposed a mesh-free method, the called node-based smoothed point interpolation method(NS-PIM),for dynamic analysis of rotating beams. A gradient smoothing technique is used, and the requirements on the consistence of the displacement functions are further weakened. In static problems, the beams with three types of boundary conditions are analyzed, and the results are compared with the exact solution, which shows the effectiveness of this method and can provide an upper bound solution for the deflection.This means that the NS-PIM makes the system soften. The NS-PIM is then further extended for solving a rigid-flexible coupled system dynamics problem, considering a rotating flexible cantilever beam. In this case, the rotating flexible cantilever beam considers not only the transverse deformations,but also the longitudinal deformations. The rigid-flexible coupled dynamic equations of the system are derived via employing Lagrange’s equations of the second type. Simulation results of the NS-PIM are compared with those obtained using finite element method(FEM) and assumed mode method. It is found that compared with FEM, the NS-PIM has anti-ill solving ability under the same calculation conditions.
基金Research is financed from the project NR03-0036-04/2008
文摘The paper presents two methods for the formulation of free vibration analysis of collecting electrodes of precipitators.The first,called the hybrid finite element method, combines the finit element method used for calculations of spring deformations with the rigid finite element method used to reflect mass and geometrical features,which is called the hybrid finite element method.As a result,a model with a diagonal mass matrix is obtained.Due to a specific geometry of the electrodes,which are long plates of complicated shapes,the second method proposed is the strip method which is a semi-analytical method.The strip method allows us to formulate the equations of motion with a considerably smaller number of generalized coordinates.Results of numerical calculations obtained by both methods are compared with those obtained using commercial software like ANSYS and ABAQUS.Good compatibility of results is achieved.
基金supported by the National Natural Science Foundation of China (11132007)
文摘Previous work examined the effect of the attached stiffness matrix terms on stability of an elastic beam undergoing prescribed large overall motion. The aim of the present work is to extend the nonlinear formulations to an elastic beam with free large overall motion. Based on initial stress method, the nonlinear coupling equations of elastic beams are obtained with free large overall motion and the attached stiffness matrix is derived by solving sub-static formulation. The angular velocity and the tip deformation of the elastic pendulum are calculated. The analytical results show that the simulation results of the present model are tabled and coincide with the one-order approximate model. It is shown that the simulation results accord with energy conservation principle.
基金Project supported by the National Natural Science Foundation of China (No.10472030).
文摘The problem of the elastic interaction between a screw dislocation and a three-phase circular inclusion with interracial rigid lines (anti-cracks) is investigated. An efficient and concise method for the complex multiply connected region is developed, with which explicit series form solutions of the complex potentials in the matrix, and the interphase layer and inclusion regions are derived. Based on the complex potentials, the image force on the screw dislocation is then calculated by using the Peach-Koehler formula. The equilibrium position of the dislocation is discussed in detail for various rigid line geometries, interphase layer thicknesses and material property combinations. The main results show that the interracial rigid lines exert a significant perturbation effect on the motion of the screw dislocation near the circular inclusion surrounded by an interphase layer.
文摘To avoid mesh distortion and iterative remeshing in mesh-based numerical analysis,a meshless approach based on element free Galerkin (EFG) method is applied to the metal forming analysis of ring compression. Discrete equations are formulated upon the moving least-squares (MLS) approximation and modified Markov variational principles for rigid-plastic/ rigid-viscoplastic (RP/RVP) material models. The penalty function is used for the incompressible condition without volumetric locking. Based on the axisymmetric mechanical model,ring tests with different friction coefficients are studied. The deformed nodal configurations and shaded contours of equivalent strains are shown by developed meshless post processor. The comparison of meshless and finite element (FE) results validates the feasibility and accuracy for meshless method to simulate metal forming process.