A local pseudo arc-length method(LPALM)for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are ...A local pseudo arc-length method(LPALM)for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.展开更多
A hyperbolic conservation equation can easily generate strong discontinuous solutions such as shock waves and contact discontinuity.By introducing the arc-length parameter,the pseudo arc-length method(PALM)smoothens t...A hyperbolic conservation equation can easily generate strong discontinuous solutions such as shock waves and contact discontinuity.By introducing the arc-length parameter,the pseudo arc-length method(PALM)smoothens the discontinuous solution in the arc-length space.This in turn weakens the singularity of the equation.To avoid constructing a high-order scheme directly in the deformed physical space,the entire calculation process is conducted in a uniform orthogonal arc-length space.Furthermore,to ensure the stability of the equation,the time step is reduced by limiting the moving speed of the mesh.Given that the calculation does not involve the interpolation process of physical quantities after the mesh moves,it maintains a high computational efficiency.The numerical examples show that the algorithm can effectively reduce numerical oscillations while maintaining excellent characteristics such as high precision and high resolution.展开更多
An inverse problem of elastica of a variable-arclength beam subjected to a concentrated load is investigated. The beam is fixed at one end, and can slide freely over a hinge support at the other end. The inverse probl...An inverse problem of elastica of a variable-arclength beam subjected to a concentrated load is investigated. The beam is fixed at one end, and can slide freely over a hinge support at the other end. The inverse problem is to determine the value of the load when the deflection of the action point of the load is given. Based on the elasitca equations and the elliptic integrals, a set of nonlinear equations for the inverse problem are derived, and an analytical solution by means of iterations and Quasi-Newton method is presented. From the results, the relationship between the loads and deflections of the loading point is obtained.展开更多
Arc-length-type and energy-type methods are two main strategies used in structural nonlinear tracing analysis, but the former is widely used due to the explicitness and clarity in conception, as well as the convenienc...Arc-length-type and energy-type methods are two main strategies used in structural nonlinear tracing analysis, but the former is widely used due to the explicitness and clarity in conception, as well as the convenience and reliability in calculation. It is very important to trace the complete load-deflection path in order to know comprehensively the characteristics of structures subjected to loads. Unfortunately, the nonlinear analysis techniques are only workable for tracing the limit-point-type equilibrium path. For the bifurcation-point-type path, most of them cannot secure a satisfactory result. In this paper, main arc-length-type methods are reviewed and compared, and the possible reasons of failures in tracing analysis are briefly discussed. Some improvements are proposed, a displacement perturbation method and a force perturbation method are presented for tracing the bifurcation-point-type paths. Finally, two examples are analyzed to verify the ideas, and some conclusions are drawn with respect to the arc-length-type methods.展开更多
Based on the conventional arc-length method, an improved arc-length method with high-efficiency is proposed. The weighted modifications with respect to the variation of structural stiffness and extra-interpolation mod...Based on the conventional arc-length method, an improved arc-length method with high-efficiency is proposed. The weighted modifications with respect to the variation of structural stiffness and extra-interpolation modification by using the information of known equilibrium points are introduced to improve the incremental arc-length. An approximate expansion method for the accumulated and expected arc-length is used to ensure the convergence at given load levels in large range of applications. Numerical results show that the improved arc-length method has well adaptability and higher efficiency in the post-buckling analysis of plates and shells structures for tracing whole load-deflection path and obtaining the convergence values at any specified load levels.展开更多
We consider the motion of a massive point-like projectile thrown with initial velocity with respect to horizontal in a two-dimensional vertical plane under the influence of gravity in a viscose media. Two different ve...We consider the motion of a massive point-like projectile thrown with initial velocity with respect to horizontal in a two-dimensional vertical plane under the influence of gravity in a viscose media. Two different velocity-dependent resistive media models are considered—linear and quadratic. With an objective to utilizing a Computer Algebra System (CAS), specifically <em>Mathematica</em> [1] numerically we solve the corresponding equations of motions. For a set of compatible parameters characterizing viscose forces graphically we display comparing the trajectories explicitly showing the impact of the models. Utilizing the model-dependent trajectory equations numerically we evaluate their associated arc-lengths. What distinguishes our approach vs. the existing body of work is the notion of the “reverse engineering”. Meaning, utilizing our numeric data we establish their corresponding analytic counter parts. Ultimately, utilizing both outputs numerically and analytically we determine the matching initial projectile angles maximizing their respective arc-lengths.展开更多
Soils with strain-softening behavior — manifesting as a reduction of strength with increasing plastic strain — are commonly found in the natural environment. For slopes in these soils,a progressive failure mechanism...Soils with strain-softening behavior — manifesting as a reduction of strength with increasing plastic strain — are commonly found in the natural environment. For slopes in these soils,a progressive failure mechanism can occur due to a reduction of strength with increasing strain. Finite element method based numerical approaches have been widely performed for simulating such failure mechanism,owning to their ability for tracing the formation and development of the localized shear strain. However,the reliability of the currently used approaches are often affected by poor convergence or significant mesh-dependency,and their applicability is limited by the use of complicated soil models. This paper aims to overcome these limitations by developing a finite element approach using a local arc-length controlled iterative algorithm as the solution strategy. In the proposed finite element approach,the soils are simulated with an elastoplastic constitutive model in conjunction with the Mohr-Coulomb yield function. The strain-softening behavior is represented by a piece-wise linearrelationship between the Mohr-Coulomb strength parameters and the deviatoric plastic strain. To assess the reliability of the proposed finite element approach,comparisons of the numerical solutions obtained by different finite element methods and meshes with various qualities are presented. Moreover,a landslide triggered by excavation in a real expressway construction project is analyzed by the presented finite element approach to demonstrate its applicability for practical engineering problems.展开更多
Nonlinear solution of reinforced concrete structures, particularly complete load-deflection response, requires tracing of the equilibrium path and proper treatment of the limit and bifurcation points. In this regard, ...Nonlinear solution of reinforced concrete structures, particularly complete load-deflection response, requires tracing of the equilibrium path and proper treatment of the limit and bifurcation points. In this regard, ordinary solution techniques lead to instability near the limit points and also have problems in case of snap-through and snap-back. Thus they fail to predict the complete load-displacement response. The arc-length method serves the purpose well in principle, received wide acceptance in finite element analysis, and has been used extensively. However modifications to the basic idea are vital to meet the particular needs of the analysis. This paper reviews some of the recent developments of the method in the last two decades, with particular emphasis on nonlinear finite element analysis of reinforced concrete structures.展开更多
An efficient method for C2 nearly arc-length parameterized curve is presented.An idea of approximation for the arc-length function of parametric curve which interpolates CAD data points is discussed.The parameterizati...An efficient method for C2 nearly arc-length parameterized curve is presented.An idea of approximation for the arc-length function of parametric curve which interpolates CAD data points is discussed.The parameterization is implemented by using parameter transformation.Finally,two numerical examples are given.展开更多
It is regretted that the author corrections requested at the proof stage were not made accurately. There are some incorrect typings in two equations which will lead to inaccurate results if readers perform calculation...It is regretted that the author corrections requested at the proof stage were not made accurately. There are some incorrect typings in two equations which will lead to inaccurate results if readers perform calculations directly展开更多
基金supported by the National Natural Science Foundation of China(11390363 and 11172041)Beijing Higher Education Young Elite Teacher Project(YETP1190)
文摘A local pseudo arc-length method(LPALM)for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.
基金Project supported by the National Natural Science Foundation of China(Nos.11822203 and 12032006)
文摘A hyperbolic conservation equation can easily generate strong discontinuous solutions such as shock waves and contact discontinuity.By introducing the arc-length parameter,the pseudo arc-length method(PALM)smoothens the discontinuous solution in the arc-length space.This in turn weakens the singularity of the equation.To avoid constructing a high-order scheme directly in the deformed physical space,the entire calculation process is conducted in a uniform orthogonal arc-length space.Furthermore,to ensure the stability of the equation,the time step is reduced by limiting the moving speed of the mesh.Given that the calculation does not involve the interpolation process of physical quantities after the mesh moves,it maintains a high computational efficiency.The numerical examples show that the algorithm can effectively reduce numerical oscillations while maintaining excellent characteristics such as high precision and high resolution.
基金The project supported by the National Natural Science Foundation of China(10272011)
文摘An inverse problem of elastica of a variable-arclength beam subjected to a concentrated load is investigated. The beam is fixed at one end, and can slide freely over a hinge support at the other end. The inverse problem is to determine the value of the load when the deflection of the action point of the load is given. Based on the elasitca equations and the elliptic integrals, a set of nonlinear equations for the inverse problem are derived, and an analytical solution by means of iterations and Quasi-Newton method is presented. From the results, the relationship between the loads and deflections of the loading point is obtained.
基金The project supported by the Special Research Fund for Doctor Program of Universities (9424702)
文摘Arc-length-type and energy-type methods are two main strategies used in structural nonlinear tracing analysis, but the former is widely used due to the explicitness and clarity in conception, as well as the convenience and reliability in calculation. It is very important to trace the complete load-deflection path in order to know comprehensively the characteristics of structures subjected to loads. Unfortunately, the nonlinear analysis techniques are only workable for tracing the limit-point-type equilibrium path. For the bifurcation-point-type path, most of them cannot secure a satisfactory result. In this paper, main arc-length-type methods are reviewed and compared, and the possible reasons of failures in tracing analysis are briefly discussed. Some improvements are proposed, a displacement perturbation method and a force perturbation method are presented for tracing the bifurcation-point-type paths. Finally, two examples are analyzed to verify the ideas, and some conclusions are drawn with respect to the arc-length-type methods.
基金Foundation item:the Aeronauties Science Foundation of China(00B01001)
文摘Based on the conventional arc-length method, an improved arc-length method with high-efficiency is proposed. The weighted modifications with respect to the variation of structural stiffness and extra-interpolation modification by using the information of known equilibrium points are introduced to improve the incremental arc-length. An approximate expansion method for the accumulated and expected arc-length is used to ensure the convergence at given load levels in large range of applications. Numerical results show that the improved arc-length method has well adaptability and higher efficiency in the post-buckling analysis of plates and shells structures for tracing whole load-deflection path and obtaining the convergence values at any specified load levels.
文摘We consider the motion of a massive point-like projectile thrown with initial velocity with respect to horizontal in a two-dimensional vertical plane under the influence of gravity in a viscose media. Two different velocity-dependent resistive media models are considered—linear and quadratic. With an objective to utilizing a Computer Algebra System (CAS), specifically <em>Mathematica</em> [1] numerically we solve the corresponding equations of motions. For a set of compatible parameters characterizing viscose forces graphically we display comparing the trajectories explicitly showing the impact of the models. Utilizing the model-dependent trajectory equations numerically we evaluate their associated arc-lengths. What distinguishes our approach vs. the existing body of work is the notion of the “reverse engineering”. Meaning, utilizing our numeric data we establish their corresponding analytic counter parts. Ultimately, utilizing both outputs numerically and analytically we determine the matching initial projectile angles maximizing their respective arc-lengths.
基金funded by the Chinese National Basic Research Program (2010CB731503)
文摘Soils with strain-softening behavior — manifesting as a reduction of strength with increasing plastic strain — are commonly found in the natural environment. For slopes in these soils,a progressive failure mechanism can occur due to a reduction of strength with increasing strain. Finite element method based numerical approaches have been widely performed for simulating such failure mechanism,owning to their ability for tracing the formation and development of the localized shear strain. However,the reliability of the currently used approaches are often affected by poor convergence or significant mesh-dependency,and their applicability is limited by the use of complicated soil models. This paper aims to overcome these limitations by developing a finite element approach using a local arc-length controlled iterative algorithm as the solution strategy. In the proposed finite element approach,the soils are simulated with an elastoplastic constitutive model in conjunction with the Mohr-Coulomb yield function. The strain-softening behavior is represented by a piece-wise linearrelationship between the Mohr-Coulomb strength parameters and the deviatoric plastic strain. To assess the reliability of the proposed finite element approach,comparisons of the numerical solutions obtained by different finite element methods and meshes with various qualities are presented. Moreover,a landslide triggered by excavation in a real expressway construction project is analyzed by the presented finite element approach to demonstrate its applicability for practical engineering problems.
文摘Nonlinear solution of reinforced concrete structures, particularly complete load-deflection response, requires tracing of the equilibrium path and proper treatment of the limit and bifurcation points. In this regard, ordinary solution techniques lead to instability near the limit points and also have problems in case of snap-through and snap-back. Thus they fail to predict the complete load-displacement response. The arc-length method serves the purpose well in principle, received wide acceptance in finite element analysis, and has been used extensively. However modifications to the basic idea are vital to meet the particular needs of the analysis. This paper reviews some of the recent developments of the method in the last two decades, with particular emphasis on nonlinear finite element analysis of reinforced concrete structures.
基金This paper is supported by the science Foundation (01B030)Educational Department of Hunan province
文摘An efficient method for C2 nearly arc-length parameterized curve is presented.An idea of approximation for the arc-length function of parametric curve which interpolates CAD data points is discussed.The parameterization is implemented by using parameter transformation.Finally,two numerical examples are given.
文摘It is regretted that the author corrections requested at the proof stage were not made accurately. There are some incorrect typings in two equations which will lead to inaccurate results if readers perform calculations directly