This paper presents the elastic and plastic deformation of the steel helmet with coldextrusion moulding. The plastic streamline of the plastic mould-making process for ellipse thinplate is described. The distribution ...This paper presents the elastic and plastic deformation of the steel helmet with coldextrusion moulding. The plastic streamline of the plastic mould-making process for ellipse thinplate is described. The distribution of slip-line is established based on the plastic streamline. Theextrusion force of plastic moulding of the steel helmet is calculated by using of slip-line method.Furthermore, an applied example is given.展开更多
The axi-symmetric active earth pressure for layered backfills was investigated using the slip line method. Due to different soil properties, the interfaces of the soil layers were considered as discontinuity surface. ...The axi-symmetric active earth pressure for layered backfills was investigated using the slip line method. Due to different soil properties, the interfaces of the soil layers were considered as discontinuity surface. Accordingly, the change of the major principal direction was obtained and a new computation scheme was proposed to deal with the discontinuity in the calculation, finally the slip line method was extended to layered backfills. Results indicate that, the major principal direction, as well as the earth pressure, has a finite jump on passing the soil interfaces. The magnitude of the jump depends on the soil properties mainly. Generally, the active earth pressure when a strong layer is overlying a weak layer is much larger than that when the weak layer is lying on the strong layer. The present solution can be reduced to plane strain case and it has been compared with Rankine's and Coulomb's results, and a good agreement is observed.展开更多
By using speetrally stable targets, the empirical line （EL） method was tested to correct the multispectral IKONOS imagery acquired over Putuo Mountain, Zhejiang, China. A series of calibration targets, which were sp...By using speetrally stable targets, the empirical line （EL） method was tested to correct the multispectral IKONOS imagery acquired over Putuo Mountain, Zhejiang, China. A series of calibration targets, which were spectrally stable over time, were selected to establish the linear predicted equation. Subsequently, a series of spectrally stable validation targets were selected to assess the accuracy of the equations. And, validation targets, which were speetrally unstable over time, were used to test the feasibility of using the EL method to calibrate the archival remotely sensed data. Ground reflectance measurements for each target were made using an ASD FieldSpec spectroradiometer. A Trimble GeoXTTM GPS unit with sub-meter accuracy was used to estimate the target position accurately. Linear regression equations for four tKONOS bands were derived. The coefficients of determination for the blue, green, and red bands were all greater than 0.9800 and it was 0.9697 for the near infrared band. It was concluded that reasonable results could be obtained by using speetrally stable targets.展开更多
This paper retrospectively analyzes the complicated rectovaginal fistula in 2 cases after the successful treatment. Through literatures collected, with the difficulties and contradictions of the treatment of complicat...This paper retrospectively analyzes the complicated rectovaginal fistula in 2 cases after the successful treatment. Through literatures collected, with the difficulties and contradictions of the treatment of complicated rectovaginal fistula as the starting point, the currently used hanging line method, the advantage and deficiency of surgical treatment in recent years, and the change of treatment method were reviewed. We wish to explore which current treatments could be a better choice.展开更多
Slip line method for sintered powder materials under condition of axial symmetry is proposed based on the simplified yield condition of sintered powder materials and Haar-von Karman perfect plastic criterion. The equa...Slip line method for sintered powder materials under condition of axial symmetry is proposed based on the simplified yield condition of sintered powder materials and Haar-von Karman perfect plastic criterion. The equations of slip line and stress along slip line are derived, and numerical solutions are given. Deformation load in closed die upsetting of sintered copper cylinder is calculated by slip line method, and theoretical solutions are compared with experimental results.展开更多
To extract quantitative biophysical parameters such as leaf biomass and leaf chlorophyll concentration from the remotely sensed imagery, the effect of atmospheric attenuation must be removed. The refined empirical li...To extract quantitative biophysical parameters such as leaf biomass and leaf chlorophyll concentration from the remotely sensed imagery, the effect of atmospheric attenuation must be removed. The refined empirical line （REL） method was used to calibrate the IKONOS multispectral imagery. The IKONOS digital numbers （DN） were converted to the at-satellite reflectance, then the linear relation between at-satellite reflectance and surface spectral reflectance （ρA） was derived from six bright targets of known reflectance in the image, and modelled estimates of the image reflectance at ρA=0. Validation targets were used to test the feasibility of REL method. The mean relative errors between ρA retrieved from IKONOS image using REL method and ground-measured ρA were 11%, 13%, 3% and 5% in the IKONOS blue, green, red and near-infrared （NIR） respectively. When dark targets are unavailable or measurement of dark target is inconvenient, the REL method was most crucial for retrieving surface spectral reflectance. The REL offers a simple approach for quantitative retrieval of biophysical parameters from IKONOS imagery.展开更多
This paper puts forward a two-parameter family of nonlinear conjugate gradient(CG)method without line search for solving unconstrained optimization problem.The main feature of this method is that it does not rely on a...This paper puts forward a two-parameter family of nonlinear conjugate gradient(CG)method without line search for solving unconstrained optimization problem.The main feature of this method is that it does not rely on any line search and only requires a simple step size formula to always generate a sufficient descent direction.Under certain assumptions,the proposed method is proved to possess global convergence.Finally,our method is compared with other potential methods.A large number of numerical experiments show that our method is more competitive and effective.展开更多
In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake d...In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence,this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES)method.Based on the dynamic k-equation large-eddy simulation(LES),this method uses a precursor method to generate atmospheric inflow turbulence,models the tower and nacelle wakes,and improves the body force projection method based on an anisotropic Gaussian distribution function.For these three improvements,three wind tunnel experiments are used to validate the numerical accuracy of this method.The results show that the numerical results calculated in the far-wake region can reflect the characteristics of typical onshore and offshore wind conditions compared with the experimental results.After modeling the tower and nacelle wakes,the wake velocity distribution is consistent with the experimental result.The radial migration velocity of the tip vortex calculated by the improved blade body force distribution model is 0.32 m/s,which is about 6%different from the experimental value and improves the prediction accuracy of the tip vortex radial movement.The method proposed in this paper is very helpful for wind turbine wake dynamic analysis and wind farm power prediction.展开更多
Rolling force and rolling moment are prime process parameter of external spline cold rolling. However, the precise theoretical formulae of rolling force and rolling moment are still very fewer, and the determination o...Rolling force and rolling moment are prime process parameter of external spline cold rolling. However, the precise theoretical formulae of rolling force and rolling moment are still very fewer, and the determination of them depends on experience. In the present study, the mathematical models of rolling force and rolling moment are established based on stress field theory of slip-line. And the isotropic hardening is used to improve the yield criterion. Based on MATLAB program language environment, calculation program is developed according to mathematical models established. The rolling force and rolling moment could be predicted quickly via the calculation program, and then the reliability of the models is validated by FEM. Within the range of module of spline m=0.5-1.5 mm, pressure angle of reference circle α=30.0°-45.0°, and number of spline teeth Z=19-54, the rolling force and rolling moment in rolling process （finishing rolling is excluded） are researched by means of virtualizing orthogonal experiment design. The results of the present study indicate that： the influences of module and number of spline teeth on the maximum rolling force and rolling moment in the process are remarkable; in the case of pressure angle of reference circle is little, module of spline is great, and number of spline teeth is little, the peak value of rolling force in rolling process may appear in the midst of the process; the peak value of rolling moment in rolling process appears in the midst of the process, and then oscillator weaken to a stable value. The results of the present study may provide guidelines for the determination of power of the motor and the design of hydraulic system of special machine, and provide basis for the farther researches on the precise forming process of external spline cold rolling.展开更多
The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good ac...The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.展开更多
Many methods are proposed to deal with the type synthesis of parallel kinematic mechanisms（PKMs）, but most of them are less intuitive to some extent. Thus, to propose a concise and intuitive type synthesis method fo...Many methods are proposed to deal with the type synthesis of parallel kinematic mechanisms（PKMs）, but most of them are less intuitive to some extent. Thus, to propose a concise and intuitive type synthesis method for engineering application is a very challenging issue, which should be further studied in the field. Grassmann line geometry, which can investigate the dimensions of spatial line-clusters in a concise way, is taken as the mathematic foundation. Atlas method is introduced to visually describe the degrees of freedom（DOFs） and constraints of a mechanism, and the dual rule is brought in to realize the mutual conversion of the freedom-space and constraint-space. Consequently, a systematic method based on Grassmann line geometry and Atlas method is generated and the entire type synthesis process is presented. Three type 4-DOF PKMs, i.e., 1T3R, 2T2R and 3T1R（T： translational DOF; R： rotational DOF）, are classified according to the different combinations of the translational DOFs and rotational DOFs. The type synthesis of 4-DOF PKMs is carried out and the possible configurations are thoroughly investigated. Some new PKMs with useful functions are generated during this procedure. The type synthesis method based on Grassmann line geometry and Atlas method is intuitive and concise, and can reduce the complexity of the PKMs＇ type synthesis. Moreover, this method can provide theoretical guidance for other PKMs＇ type synthesis and engineering application. A novel type synthesis method is proposed, which solves the existing methods＇ problems in terms of complicated, not intuitive and unsuitable for practical application.展开更多
This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two node...This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two nodes is briefly presented on the basis of the total Lagrangian formulation. The static and dynamic equilibrium equations of mooring lines are established. An incremental-iterative method is used to determine the initial static equilibrium state of cable systems under the action of self weights, buoyancy and current. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method, and examine the effect of various parameters.展开更多
The Finite Element Method of Lines （FEMOL） is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to...The Finite Element Method of Lines （FEMOL） is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to a set of examples such as edge-crack plate, the central-crack plate, the plate with cracks emanating from a hole under tensile or under combination loads of tensile and bending. Their dimensionless stress distribution, the stress intensify factor （SIF） and crack opening displacement （COD） are obtained, and comparison with known solutions by other methods are reported. It is found that a good accuracy is achieved by FEMOL. The method is successfully modified to remarkably increase the accuracy and reduce convergence difficulties. So it is a very useful and new tool in studying fracture mechanics problems.展开更多
In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs...In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.展开更多
In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-differenc...In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-difference and finite element methods to approximate spatial derivatives, this new technique does not require a mesh in the problem domain, and a set of scattered nodes provided by initial data is required for the solution of the problem using some radial basis functions. Accuracy of the method is assessed in terms of the error norms L2, L∞ and the three invariants C1, C2, C3. Numerical experiments are performed to demonstrate the accuracy and easy implementation of this method for the three classes of time-dependent nonlinear coupled partial differential equations.展开更多
The finite element method (FEM) and the boundary element method (BEM) are often adopted. However, they are not convenient to spatially vary thermal properties of functionally graded material (FGM). Therefore, the meth...The finite element method (FEM) and the boundary element method (BEM) are often adopted. However, they are not convenient to spatially vary thermal properties of functionally graded material (FGM). Therefore, the method of lines (MOL) is introduced to solve the temperature field of FGM. The basic idea of the method is to semi-discretize the governing equation into a system of ordinary differential equations (ODEs) defined on discrete lines by means of the finite difference method. The temperature field of FGM can be obtained by solving the ODEs. The functions of thermal properties are directly embodied in these equations and these properties are not discretized in the domain. Thus, difficulty of FEM and BEM is overcome by the method. As a numerical example, the temperature field of a plane problem is analyzed for FGMs through varying thermal conductivity coefficient by the MOL.展开更多
The method of lines based on Hu Hai-chang 's theory for the vibration and stability of moderate thick plates is developed. The standard nonlinear ordinary differential equation (ODE) system for natural frequencies...The method of lines based on Hu Hai-chang 's theory for the vibration and stability of moderate thick plates is developed. The standard nonlinear ordinary differential equation (ODE) system for natural frequencies and critical load is given by use of ODE techniques, and then any indicated eigenvalue could be obtained directly from ODE solver by employing the so-called initial eigenfunction technique instead of the mode orthogonality condition. Numerical examples show that the present method is very effective and reliable.展开更多
The method of lines is applied to the boundary-value problem for third order partial differential equation. Explicit expression and order of convergence for the approximate solution are obtained.
In this paper, we solve chiral nonlinear Schrodinger equation (CNSE) numerically. Two numerical methods are derived using the explicit Runge-Kutta method of order four and the linear multistep method (Predictor-Correc...In this paper, we solve chiral nonlinear Schrodinger equation (CNSE) numerically. Two numerical methods are derived using the explicit Runge-Kutta method of order four and the linear multistep method (Predictor-Corrector method of fourth order). The resulting schemes of fourth order accuracy in spatial and temporal directions. The CNSE is non-integrable and has two kinds of soliton solutions: bright and dark soliton. The exact solutions and the conserved quantities of CNSE are used to display the efficiency and robustness of the numerical methods we derived. Interaction of two bright solitons for different parameters is also displayed.展开更多
Based on the sub-region generalized variationM principle, a sub-region mixed version of the newly-developed semi-analytical ＇finite element method of lines＇ （FEMOL） is proposed in this paper for accurate and effic...Based on the sub-region generalized variationM principle, a sub-region mixed version of the newly-developed semi-analytical ＇finite element method of lines＇ （FEMOL） is proposed in this paper for accurate and efficient computation of stress intensity factors （SIFs） of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations （ODEs） and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.展开更多
文摘This paper presents the elastic and plastic deformation of the steel helmet with coldextrusion moulding. The plastic streamline of the plastic mould-making process for ellipse thinplate is described. The distribution of slip-line is established based on the plastic streamline. Theextrusion force of plastic moulding of the steel helmet is calculated by using of slip-line method.Furthermore, an applied example is given.
基金the National Natural Science Foundation of China (No. 50679041)the Shanghai Leading Academic Discipline Project (No. B208)
文摘The axi-symmetric active earth pressure for layered backfills was investigated using the slip line method. Due to different soil properties, the interfaces of the soil layers were considered as discontinuity surface. Accordingly, the change of the major principal direction was obtained and a new computation scheme was proposed to deal with the discontinuity in the calculation, finally the slip line method was extended to layered backfills. Results indicate that, the major principal direction, as well as the earth pressure, has a finite jump on passing the soil interfaces. The magnitude of the jump depends on the soil properties mainly. Generally, the active earth pressure when a strong layer is overlying a weak layer is much larger than that when the weak layer is lying on the strong layer. The present solution can be reduced to plane strain case and it has been compared with Rankine's and Coulomb's results, and a good agreement is observed.
基金Project supported by the National Natural Science Foundation of China (No.40171065)the National High Technology Research and Development Program of China (Nos.2002AA130010-2-7 and 2003AA131020-04-06).
文摘By using speetrally stable targets, the empirical line （EL） method was tested to correct the multispectral IKONOS imagery acquired over Putuo Mountain, Zhejiang, China. A series of calibration targets, which were spectrally stable over time, were selected to establish the linear predicted equation. Subsequently, a series of spectrally stable validation targets were selected to assess the accuracy of the equations. And, validation targets, which were speetrally unstable over time, were used to test the feasibility of using the EL method to calibrate the archival remotely sensed data. Ground reflectance measurements for each target were made using an ASD FieldSpec spectroradiometer. A Trimble GeoXTTM GPS unit with sub-meter accuracy was used to estimate the target position accurately. Linear regression equations for four tKONOS bands were derived. The coefficients of determination for the blue, green, and red bands were all greater than 0.9800 and it was 0.9697 for the near infrared band. It was concluded that reasonable results could be obtained by using speetrally stable targets.
文摘This paper retrospectively analyzes the complicated rectovaginal fistula in 2 cases after the successful treatment. Through literatures collected, with the difficulties and contradictions of the treatment of complicated rectovaginal fistula as the starting point, the currently used hanging line method, the advantage and deficiency of surgical treatment in recent years, and the change of treatment method were reviewed. We wish to explore which current treatments could be a better choice.
文摘Slip line method for sintered powder materials under condition of axial symmetry is proposed based on the simplified yield condition of sintered powder materials and Haar-von Karman perfect plastic criterion. The equations of slip line and stress along slip line are derived, and numerical solutions are given. Deformation load in closed die upsetting of sintered copper cylinder is calculated by slip line method, and theoretical solutions are compared with experimental results.
基金Project supported by the Hi-Tech Research and Development Pro-gram (863) of China (Nos. 2002AA130010-2-7 and 2003AA131020- 04-06) and the National Natural Science Foundation of China (No. 40171065)
文摘To extract quantitative biophysical parameters such as leaf biomass and leaf chlorophyll concentration from the remotely sensed imagery, the effect of atmospheric attenuation must be removed. The refined empirical line （REL） method was used to calibrate the IKONOS multispectral imagery. The IKONOS digital numbers （DN） were converted to the at-satellite reflectance, then the linear relation between at-satellite reflectance and surface spectral reflectance （ρA） was derived from six bright targets of known reflectance in the image, and modelled estimates of the image reflectance at ρA=0. Validation targets were used to test the feasibility of REL method. The mean relative errors between ρA retrieved from IKONOS image using REL method and ground-measured ρA were 11%, 13%, 3% and 5% in the IKONOS blue, green, red and near-infrared （NIR） respectively. When dark targets are unavailable or measurement of dark target is inconvenient, the REL method was most crucial for retrieving surface spectral reflectance. The REL offers a simple approach for quantitative retrieval of biophysical parameters from IKONOS imagery.
基金Supported by 2023 Inner Mongolia University of Finance and Economics,General Scientific Research for Universities directly under Inner Mon‐golia,China (NCYWT23026)2024 High-quality Research Achievements Cultivation Fund Project of Inner Mongolia University of Finance and Economics,China (GZCG2479)。
文摘This paper puts forward a two-parameter family of nonlinear conjugate gradient(CG)method without line search for solving unconstrained optimization problem.The main feature of this method is that it does not rely on any line search and only requires a simple step size formula to always generate a sufficient descent direction.Under certain assumptions,the proposed method is proved to possess global convergence.Finally,our method is compared with other potential methods.A large number of numerical experiments show that our method is more competitive and effective.
基金Project supported by the National Key Research and Development Program of China(Nos.2019YFE0192600,2017YFE0132000,and 2019YFB1503700)the National Natural Science Foundation of China(Nos.51761135012 and 11872248)。
文摘In a large wind farm,the wakes of upstream and downstream wind turbines can interfere with each other,affecting the overall power output of the wind farm.To further improve the numerical accuracy of the turbine wake dynamics under atmosphere turbulence,this work proposes some improvements to the actuator line-large-eddy simulation(AL-LES)method.Based on the dynamic k-equation large-eddy simulation(LES),this method uses a precursor method to generate atmospheric inflow turbulence,models the tower and nacelle wakes,and improves the body force projection method based on an anisotropic Gaussian distribution function.For these three improvements,three wind tunnel experiments are used to validate the numerical accuracy of this method.The results show that the numerical results calculated in the far-wake region can reflect the characteristics of typical onshore and offshore wind conditions compared with the experimental results.After modeling the tower and nacelle wakes,the wake velocity distribution is consistent with the experimental result.The radial migration velocity of the tip vortex calculated by the improved blade body force distribution model is 0.32 m/s,which is about 6%different from the experimental value and improves the prediction accuracy of the tip vortex radial movement.The method proposed in this paper is very helpful for wind turbine wake dynamic analysis and wind farm power prediction.
基金supported by National Natural Science Foundation of China (Grant No. 50675145)Shanxi Provincial Key Project of Science and Technology of China (Grant No. 2006031147)+1 种基金Shanxi Provincial Innovation Project for Graduate Students of China (Grant No. 20061027)Shanxi Provincial Key Project for Studied-abroad Returnee of China
文摘Rolling force and rolling moment are prime process parameter of external spline cold rolling. However, the precise theoretical formulae of rolling force and rolling moment are still very fewer, and the determination of them depends on experience. In the present study, the mathematical models of rolling force and rolling moment are established based on stress field theory of slip-line. And the isotropic hardening is used to improve the yield criterion. Based on MATLAB program language environment, calculation program is developed according to mathematical models established. The rolling force and rolling moment could be predicted quickly via the calculation program, and then the reliability of the models is validated by FEM. Within the range of module of spline m=0.5-1.5 mm, pressure angle of reference circle α=30.0°-45.0°, and number of spline teeth Z=19-54, the rolling force and rolling moment in rolling process （finishing rolling is excluded） are researched by means of virtualizing orthogonal experiment design. The results of the present study indicate that： the influences of module and number of spline teeth on the maximum rolling force and rolling moment in the process are remarkable; in the case of pressure angle of reference circle is little, module of spline is great, and number of spline teeth is little, the peak value of rolling force in rolling process may appear in the midst of the process; the peak value of rolling moment in rolling process appears in the midst of the process, and then oscillator weaken to a stable value. The results of the present study may provide guidelines for the determination of power of the motor and the design of hydraulic system of special machine, and provide basis for the farther researches on the precise forming process of external spline cold rolling.
文摘The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.
基金supported by National Natural Science Foundation of China(Grant No.51135008)National Basic Research Program of China(973 Program,Grant No.2013CB035400)China Postdoctoral Science Foundation(Grant Nos.2012M520256,2013T60107)
文摘Many methods are proposed to deal with the type synthesis of parallel kinematic mechanisms（PKMs）, but most of them are less intuitive to some extent. Thus, to propose a concise and intuitive type synthesis method for engineering application is a very challenging issue, which should be further studied in the field. Grassmann line geometry, which can investigate the dimensions of spatial line-clusters in a concise way, is taken as the mathematic foundation. Atlas method is introduced to visually describe the degrees of freedom（DOFs） and constraints of a mechanism, and the dual rule is brought in to realize the mutual conversion of the freedom-space and constraint-space. Consequently, a systematic method based on Grassmann line geometry and Atlas method is generated and the entire type synthesis process is presented. Three type 4-DOF PKMs, i.e., 1T3R, 2T2R and 3T1R（T： translational DOF; R： rotational DOF）, are classified according to the different combinations of the translational DOFs and rotational DOFs. The type synthesis of 4-DOF PKMs is carried out and the possible configurations are thoroughly investigated. Some new PKMs with useful functions are generated during this procedure. The type synthesis method based on Grassmann line geometry and Atlas method is intuitive and concise, and can reduce the complexity of the PKMs＇ type synthesis. Moreover, this method can provide theoretical guidance for other PKMs＇ type synthesis and engineering application. A novel type synthesis method is proposed, which solves the existing methods＇ problems in terms of complicated, not intuitive and unsuitable for practical application.
基金supported by the National Natural Science Foundation of China (Grant No.11072052)the National High Technology Research and Development Program of China (863 Program,Grant No.2006AA09A109-3)
文摘This study has focused on developing numerical procedures for the static and dynamic nonlinear analysis of mooring lines. A geometrically nonlinear finite element method using isoparametric cable element with two nodes is briefly presented on the basis of the total Lagrangian formulation. The static and dynamic equilibrium equations of mooring lines are established. An incremental-iterative method is used to determine the initial static equilibrium state of cable systems under the action of self weights, buoyancy and current. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method, and examine the effect of various parameters.
文摘The Finite Element Method of Lines （FEMOL） is a semi-analytic approach and takes a position between FEM and analytic methods. First, FEMOL in Fracture Mechanics is presented in detail. Then, the method is applied to a set of examples such as edge-crack plate, the central-crack plate, the plate with cracks emanating from a hole under tensile or under combination loads of tensile and bending. Their dimensionless stress distribution, the stress intensify factor （SIF） and crack opening displacement （COD） are obtained, and comparison with known solutions by other methods are reported. It is found that a good accuracy is achieved by FEMOL. The method is successfully modified to remarkably increase the accuracy and reduce convergence difficulties. So it is a very useful and new tool in studying fracture mechanics problems.
文摘In this study, the method of lines (MOLs) with higher order central difference approximation method coupled with the classical fourth order Runge-Kutta (RK(4,4)) method is used in solving shallow water equations (SWEs) in Cartesian coordinates to foresee water levels associated with a storm accurately along the coast of Bangladesh. In doing so, the partial derivatives of the SWEs with respect to the space variables were discretized with 5-point central difference, as a test case, to obtain a system of ordinary differential equations with time as an independent variable for every spatial grid point, which with initial conditions were solved by the RK(4,4) method. The complex land-sea interface and bottom topographic details were incorporated closely using nested schemes. The coastal and island boundaries were rectangularized through proper stair step representation, and the storing positions of the scalar and momentum variables were specified according to the rules of structured C-grid. A stable tidal regime was made over the model domain considering the effect of the major tidal constituent, M2 along the southern open boundary of the outermost parent scheme. The Meghna River fresh water discharge was taken into account for the inner most child scheme. To take into account the dynamic interaction of tide and surge, the generated tidal regime was introduced as the initial state of the sea, and the surge was then made to come over it through computer simulation. Numerical experiments were performed with the cyclone April 1991 to simulate water levels due to tide, surge, and their interaction at different stations along the coast of Bangladesh. Our computed results were found to compare reasonable well with the limited observed data obtained from Bangladesh Inland Water Transport Authority (BIWTA) and were found to be better in comparison with the results obtained through the regular finite difference method and the 3-point central difference MOLs coupled with the RK(4,4) method with regard to the root mean square error values.
文摘In this paper a meshless method of lines is proposed for the numerical solution of time-dependent nonlinear coupled partial differential equations. Contrary to mesh oriented methods of lines using the finite-difference and finite element methods to approximate spatial derivatives, this new technique does not require a mesh in the problem domain, and a set of scattered nodes provided by initial data is required for the solution of the problem using some radial basis functions. Accuracy of the method is assessed in terms of the error norms L2, L∞ and the three invariants C1, C2, C3. Numerical experiments are performed to demonstrate the accuracy and easy implementation of this method for the three classes of time-dependent nonlinear coupled partial differential equations.
文摘The finite element method (FEM) and the boundary element method (BEM) are often adopted. However, they are not convenient to spatially vary thermal properties of functionally graded material (FGM). Therefore, the method of lines (MOL) is introduced to solve the temperature field of FGM. The basic idea of the method is to semi-discretize the governing equation into a system of ordinary differential equations (ODEs) defined on discrete lines by means of the finite difference method. The temperature field of FGM can be obtained by solving the ODEs. The functions of thermal properties are directly embodied in these equations and these properties are not discretized in the domain. Thus, difficulty of FEM and BEM is overcome by the method. As a numerical example, the temperature field of a plane problem is analyzed for FGMs through varying thermal conductivity coefficient by the MOL.
基金The project supported by the Pioneer Fundation of Tongji University
文摘The method of lines based on Hu Hai-chang 's theory for the vibration and stability of moderate thick plates is developed. The standard nonlinear ordinary differential equation (ODE) system for natural frequencies and critical load is given by use of ODE techniques, and then any indicated eigenvalue could be obtained directly from ODE solver by employing the so-called initial eigenfunction technique instead of the mode orthogonality condition. Numerical examples show that the present method is very effective and reliable.
文摘The method of lines is applied to the boundary-value problem for third order partial differential equation. Explicit expression and order of convergence for the approximate solution are obtained.
文摘In this paper, we solve chiral nonlinear Schrodinger equation (CNSE) numerically. Two numerical methods are derived using the explicit Runge-Kutta method of order four and the linear multistep method (Predictor-Corrector method of fourth order). The resulting schemes of fourth order accuracy in spatial and temporal directions. The CNSE is non-integrable and has two kinds of soliton solutions: bright and dark soliton. The exact solutions and the conserved quantities of CNSE are used to display the efficiency and robustness of the numerical methods we derived. Interaction of two bright solitons for different parameters is also displayed.
基金Project supported by the National Natural Sciences Foundation of China(Nos.59525813 and 19872066)the Cardiff Advanced Chinese Engineering Centre of Cardiff University.
文摘Based on the sub-region generalized variationM principle, a sub-region mixed version of the newly-developed semi-analytical ＇finite element method of lines＇ （FEMOL） is proposed in this paper for accurate and efficient computation of stress intensity factors （SIFs） of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations （ODEs） and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.
