We have deduced incremental harmonic balance an iteration scheme in the (IHB) method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial va...We have deduced incremental harmonic balance an iteration scheme in the (IHB) method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence.展开更多
Single-point incremental forming (SPIF) is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation a...Single-point incremental forming (SPIF) is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation analysis on forming process becomes an important and useful method for the planning of shell products, the choice of material, the design of the forming process and the planning of the forming tool. Using solid brick elements, the finite element method(FEM) model of truncated pyramid was established. Based on the theory of anisotropy and assumed strain formulation, the SPIF processes with different parameters were simulated. The resulted comparison between the simulations and the experiments shows that the FEM model is feasible and effective. Then, according to the simulated forming process, the deformation pattern of SPIF can be summarized as the combination of plane-stretching deformation and bending deformation. And the study about the process parameters' impact on deformation shows that the process parameter of interlayer spacing is a dominant factor on the deformation. Decreasing interlayer spacing, the strain of one step decreases and the formability of blank will be improved. With bigger interlayer spacing, the plastic deformation zone increases and the forming force will be bigger.展开更多
The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered...The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities.展开更多
Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices cal...Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy.展开更多
The usage of steel-mixed composite beams is quite extensive today.During an event of constructing steel-mixed composite bridges,the incremental launching construction method is generally adopted.This paper mainly anal...The usage of steel-mixed composite beams is quite extensive today.During an event of constructing steel-mixed composite bridges,the incremental launching construction method is generally adopted.This paper mainly analyzes the force of incremental launching construction on a steel-concrete continuous beam bridge,the classification of incremental launching construction,the application of incremental launching construction in steel-mixed composite beams,the temporary facilities existing in incremental launching construction as well as their existing problems.Lastly,the analysis of the stress of composite beams in incremental launching construction is described by using the reference for the construction of mixed composite continuous beam bridges provided.展开更多
By the improvement of Riks’and Crisfield’s arc-length method,the adaptiveparameter incremental method is preasted for predicting the snapping response ofstructures. Its justification is fulfilled. Finally,the effect...By the improvement of Riks’and Crisfield’s arc-length method,the adaptiveparameter incremental method is preasted for predicting the snapping response ofstructures. Its justification is fulfilled. Finally,the effectiveness of this method isdemonstrated by solving the snapping response of spherical caps subjected to centrallydistributed pressures.展开更多
Accelerated proximal gradient methods have recently been developed for solving quasi-static incremental problems of elastoplastic analysis with some different yield criteria.It has been demonstrated through numerical ...Accelerated proximal gradient methods have recently been developed for solving quasi-static incremental problems of elastoplastic analysis with some different yield criteria.It has been demonstrated through numerical experiments that these methods can outperform conventional optimization-based approaches in computational plasticity.However,in literature these algorithms are described individually for specific yield criteria,and hence there exists no guide for application of the algorithms to other yield criteria.This short paper presents a general form of algorithm design,independent of specific forms of yield criteria,that unifies the existing proximal gradient methods.Clear interpretation is also given to each step of the presented general algorithm so that each update rule is linked to the underlying physical laws in terms of mechanical quantities.展开更多
In the incremental sheet forming (ISF) process, springback is a very important factor that affects the quality of parts. Predicting and controlling springback accurately is essential for the design of the toolpath f...In the incremental sheet forming (ISF) process, springback is a very important factor that affects the quality of parts. Predicting and controlling springback accurately is essential for the design of the toolpath for ISF. A three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those from the experiment. The springback angle was found to be in accordance with the experimental result, proving the FEM to be effective. A coupled artificial neural networks (ANN) and finite element method technique was developed to simulate and predict springback responses to changes in the processing parameters. A particle swarm optimization (PSO) algorithm was used to optimize the weights and thresholds of the neural network model. The neural network was trained using available FEM simulation data. The results showed that a more accurate prediction of s!oringback can be acquired using the FEM-PSONN model.展开更多
In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the pr...In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those of experiment. The results of numerical simulations, such as the strain history and distribution, the stress state and distribution, sheet thickness distribution, etc, were discussed in details, and the influences of process parameters on these results were also analyzed. The simulated results of the radial strain and the thickness distribution are in good agreement with experimental results. The simulations reveal that the deformation is localized around the tool and constantly remains close to a plane strain state. With decreasing depth step, increasing tool diameter and wall inclination angle, the axial stress reduces, leading to less thinning and more homogeneous plastic strain and thickness distribution. During ISF, the plastic strain increases stepwise under the action of the tool. Each increase in plastic strain is accompanied by hydrostatic pressure, which explains why obtainable deformation using ISF exceeds the forming limits of conventional sheet forming.展开更多
A novel incremental nonlinear detection algorithm is presented for Multiple-Input Multiple-Output (MIMO) system. In this algorithm, the received data at multiple receiver antennas are nonlinearly mapped and then sum...A novel incremental nonlinear detection algorithm is presented for Multiple-Input Multiple-Output (MIMO) system. In this algorithm, the received data at multiple receiver antennas are nonlinearly mapped and then summed with weights. The weight coefficients are incrementally computed to avoid direct computation of the inverse of a matrix, which greatly reduce the computational complexity. Simulation and comparison show that the proposed algorithm can obtain better performance of Bit Error Rate (BER) than linear Minimum Mean Square Error (MMSE).展开更多
Warp yarns and weft yarns of plain woven fabric are the principal axes of mate- rial of fabric. They are orthogonal in their original con?guration, but are obliquely crisscross in deformed con?guration in general. I...Warp yarns and weft yarns of plain woven fabric are the principal axes of mate- rial of fabric. They are orthogonal in their original con?guration, but are obliquely crisscross in deformed con?guration in general. In this paper the expressions of incremental components of strain tensor are derived, the non-linear model of woven fabric is linearized physically and its geometric non-linearity survives. The convenience of determining the total deformation is shown by the choice of the coordinate system of the principal axes of the material, with the convergence of the incremental methods illustrated by examples. This incremental model furnishes a basis for numerical simulations of fabric draping and wrinkling based on the micro-mechanical model of fabric.展开更多
A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successful...A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successfully developed for the analysis of truss,beam,frame,and 2D continua problems. In these analyses,LIMcan provide more precise stress results and less computational time consumption compared with displacement-based FEM. The plate element was based on the Mindlin-Reissner plate theory which took into account the transverse shear effects.Numerical examples were presented to study its performance including accuracy and convergence behavior,and the results were compared with the results have been obtained from the displacementbased quadrilateral plate elements and the analytical solutions. The4NQP13 element can analyze the moderately thick plates and the thin plates using LIMand is free from spurious zero energy modes and free from shear locking for thin plate analysis.展开更多
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.展开更多
A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed ...A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed by using the moving least square approximation,and the discrete governing equations for elasto-plastic material are constructed with the direct collo- cation method.The boundary conditions are also imposed by collocation.The method established is a truly meshless one,as it does not need any mesh,either for the purpose of interpolation of the solution variables,or for the purpose of construction of the discrete equations.It is simply formu- lated and very efficient,and no post-processing procedure is required to compute the derivatives of the unknown variables,since the solution from this method based on the moving least square approximation is already smooth enough.Numerical examples are given to verify the accuracy of the meshless method proposed for elasto-plasticity analysis.展开更多
Elasto-plastic consolidation is one of the classic coupling questions in geomechanics. To solve this problem, an elasto-plastic constitutive model is derived based on the numerical modeling method. The model is applie...Elasto-plastic consolidation is one of the classic coupling questions in geomechanics. To solve this problem, an elasto-plastic constitutive model is derived based on the numerical modeling method. The model is applied to Blot's consolidation theory. Incremental governing partial differential equations are established using this method. According to the stress path, the decoupling condition of these equations is discussed. Based on these conditions, an incremental diffusion equation and uncoupling governing equations are presented. The method is then applied to numerical analyses of three examples. The results show that (1) the effect of the stress path should be taken into account in the simulation of the soil consolidation question; (2) this decoupling method can predict the evolvement of pore water pressure; (3) the settlement using cam-clay model is less than that using numerical model because of dilatancy.展开更多
The nonlinear quasi-conforming FEM is presented based on the basic concept of the quasi- -conforming finite element. First, the incremental principle of stationary potential energy is discussed, Then, the formulation ...The nonlinear quasi-conforming FEM is presented based on the basic concept of the quasi- -conforming finite element. First, the incremental principle of stationary potential energy is discussed, Then, the formulation process of the nonlinear quasi-conforming FEM is given. Lastly, two computational examples of shells are given.展开更多
The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quad...The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures.展开更多
Current building design codes allow the appearance of structural and nonstructural damage under design basis earthquakes.The research regarding probabilistic seismic loss estimation in domestic building structure is u...Current building design codes allow the appearance of structural and nonstructural damage under design basis earthquakes.The research regarding probabilistic seismic loss estimation in domestic building structure is urgent.The evaluation in this paper is based on a 11-story reinforced concrete office building,incremental dynamic analysis(IDA)is conducted in Perform 3D program using models capable to simulate all possible limit states up to collapse.Next,the probability distribution of post-earthquake casualties,rebuild costs repair costs and business downtime loss are calculated in PACT software for the studied building considering the modified component vulnerability groups and population flow models.The evaluation procedure can also shed light on other types of buildings in China.For non-typical functional building structures,this article proposes to build a finite element model of structural components and to classify the vulnerability groups based on the construction drawings,and to supply and improve the vulnerability library of appendages in FEMA P-58 according to the actual situation.In this way,the application scope of building seismic performance evaluation can be expanded.展开更多
基金supported by the National Natural Science Foundation of China (10772202)Doctoral Program Foundation of Ministry of Education of China (20050558032)Guangdong Province Natural Science Foundation (07003680, 05003295)
文摘We have deduced incremental harmonic balance an iteration scheme in the (IHB) method using the harmonic balance plus the Newton-Raphson method. Since the convergence of the iteration is dependent upon the initial values in the iteration, the convergent region is greatly restricted for some cases. In this contribution, in order to enlarge the convergent region of the IHB method, we constructed the zeroth-order deformation equation using the homotopy analysis method, in which the IHB method is employed to solve the deformation equation with an embedding parameter as the active increment. Taking the Duffing and the van der Pol equations as examples, we obtained the highly accurate solutions. Importantly, the presented approach renders a convenient way to control and adjust the convergence.
基金supported by National Natural Science Foundation of China(No. 50175034).
文摘Single-point incremental forming (SPIF) is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation analysis on forming process becomes an important and useful method for the planning of shell products, the choice of material, the design of the forming process and the planning of the forming tool. Using solid brick elements, the finite element method(FEM) model of truncated pyramid was established. Based on the theory of anisotropy and assumed strain formulation, the SPIF processes with different parameters were simulated. The resulted comparison between the simulations and the experiments shows that the FEM model is feasible and effective. Then, according to the simulated forming process, the deformation pattern of SPIF can be summarized as the combination of plane-stretching deformation and bending deformation. And the study about the process parameters' impact on deformation shows that the process parameter of interlayer spacing is a dominant factor on the deformation. Decreasing interlayer spacing, the strain of one step decreases and the formability of blank will be improved. With bigger interlayer spacing, the plastic deformation zone increases and the forming force will be bigger.
基金Project supported by the Ph. D. Programs Foundation of Ministry of Education of China (No.20050558032) the Natural Science Foundation of Guangdong Province of China (No.05003295) the Foundation of Sun Yat-sen University Advanced Research Center (No.06M8) the Young Teacher Scientific Research Foundation of Sun Sat-sen University (No.1131011)
文摘The incremental harmonic balance method was extended to analyze the flutter of systems with multiple structural strong nonlinearities. The strongly nonlinear cubic plunging and pitching stiffness terms were considered in the flutter equations of two-dimensional airfoil. First, the equations were transferred into matrix form, then the vibration process was divided into the persistent incremental processes of vibration moments. And the expression of their solutions could be obtained by using a certain amplitude as control parameter in the harmonic balance process, and then the bifurcation, limit cycle flutter phenomena and the number of harmonic terms were analyzed. Finally, numerical results calculated by the Runge-Kutta method were given to verify the results obtained by the proposed procedure. It has been shown that the incremental harmonic method is effective and precise in the analysis of strongly nonlinear flutter with multiple structural nonlinearities.
基金the National Natural Science Foundation of China(Nos.11702215 and11972277)the Natural Science Basic Research Plan in Shaanxi Province of China(Nos.2017JQ5062 and 2018JQ1029)。
文摘Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy.
文摘The usage of steel-mixed composite beams is quite extensive today.During an event of constructing steel-mixed composite bridges,the incremental launching construction method is generally adopted.This paper mainly analyzes the force of incremental launching construction on a steel-concrete continuous beam bridge,the classification of incremental launching construction,the application of incremental launching construction in steel-mixed composite beams,the temporary facilities existing in incremental launching construction as well as their existing problems.Lastly,the analysis of the stress of composite beams in incremental launching construction is described by using the reference for the construction of mixed composite continuous beam bridges provided.
文摘By the improvement of Riks’and Crisfield’s arc-length method,the adaptiveparameter incremental method is preasted for predicting the snapping response ofstructures. Its justification is fulfilled. Finally,the effectiveness of this method isdemonstrated by solving the snapping response of spherical caps subjected to centrallydistributed pressures.
文摘Accelerated proximal gradient methods have recently been developed for solving quasi-static incremental problems of elastoplastic analysis with some different yield criteria.It has been demonstrated through numerical experiments that these methods can outperform conventional optimization-based approaches in computational plasticity.However,in literature these algorithms are described individually for specific yield criteria,and hence there exists no guide for application of the algorithms to other yield criteria.This short paper presents a general form of algorithm design,independent of specific forms of yield criteria,that unifies the existing proximal gradient methods.Clear interpretation is also given to each step of the presented general algorithm so that each update rule is linked to the underlying physical laws in terms of mechanical quantities.
基金Project(50175034) supported by the National Natural Science Foundation of China
文摘In the incremental sheet forming (ISF) process, springback is a very important factor that affects the quality of parts. Predicting and controlling springback accurately is essential for the design of the toolpath for ISF. A three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those from the experiment. The springback angle was found to be in accordance with the experimental result, proving the FEM to be effective. A coupled artificial neural networks (ANN) and finite element method technique was developed to simulate and predict springback responses to changes in the processing parameters. A particle swarm optimization (PSO) algorithm was used to optimize the weights and thresholds of the neural network model. The neural network was trained using available FEM simulation data. The results showed that a more accurate prediction of s!oringback can be acquired using the FEM-PSONN model.
基金Project(50175034) supported by the National Natural Science Foundation of China
文摘In order to investigate the process of incremental sheet forming (ISF) through both experimental and numerical approaches, a three-dimensional elasto-plastic finite element model (FEM) was developed to simulate the process and the simulated results were compared with those of experiment. The results of numerical simulations, such as the strain history and distribution, the stress state and distribution, sheet thickness distribution, etc, were discussed in details, and the influences of process parameters on these results were also analyzed. The simulated results of the radial strain and the thickness distribution are in good agreement with experimental results. The simulations reveal that the deformation is localized around the tool and constantly remains close to a plane strain state. With decreasing depth step, increasing tool diameter and wall inclination angle, the axial stress reduces, leading to less thinning and more homogeneous plastic strain and thickness distribution. During ISF, the plastic strain increases stepwise under the action of the tool. Each increase in plastic strain is accompanied by hydrostatic pressure, which explains why obtainable deformation using ISF exceeds the forming limits of conventional sheet forming.
文摘A novel incremental nonlinear detection algorithm is presented for Multiple-Input Multiple-Output (MIMO) system. In this algorithm, the received data at multiple receiver antennas are nonlinearly mapped and then summed with weights. The weight coefficients are incrementally computed to avoid direct computation of the inverse of a matrix, which greatly reduce the computational complexity. Simulation and comparison show that the proposed algorithm can obtain better performance of Bit Error Rate (BER) than linear Minimum Mean Square Error (MMSE).
基金Project supported by the National Natural Science Foundation of China (No. 10272079).
文摘Warp yarns and weft yarns of plain woven fabric are the principal axes of mate- rial of fabric. They are orthogonal in their original con?guration, but are obliquely crisscross in deformed con?guration in general. In this paper the expressions of incremental components of strain tensor are derived, the non-linear model of woven fabric is linearized physically and its geometric non-linearity survives. The convenience of determining the total deformation is shown by the choice of the coordinate system of the principal axes of the material, with the convergence of the incremental methods illustrated by examples. This incremental model furnishes a basis for numerical simulations of fabric draping and wrinkling based on the micro-mechanical model of fabric.
基金National Natural Science Foundation of China(No.10872128)
文摘A force-based quadrilateral plate element( 4NQP13) for the analysis of the plate bending problems using large increment method( LIM) was proposed. The LIM, a force-based finite element method( FEM),has been successfully developed for the analysis of truss,beam,frame,and 2D continua problems. In these analyses,LIMcan provide more precise stress results and less computational time consumption compared with displacement-based FEM. The plate element was based on the Mindlin-Reissner plate theory which took into account the transverse shear effects.Numerical examples were presented to study its performance including accuracy and convergence behavior,and the results were compared with the results have been obtained from the displacementbased quadrilateral plate elements and the analytical solutions. The4NQP13 element can analyze the moderately thick plates and the thin plates using LIMand is free from spurious zero energy modes and free from shear locking for thin plate analysis.
基金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.
基金Project supported by the National Natural Science Foundation of China(No.10172052).
文摘A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed by using the moving least square approximation,and the discrete governing equations for elasto-plastic material are constructed with the direct collo- cation method.The boundary conditions are also imposed by collocation.The method established is a truly meshless one,as it does not need any mesh,either for the purpose of interpolation of the solution variables,or for the purpose of construction of the discrete equations.It is simply formu- lated and very efficient,and no post-processing procedure is required to compute the derivatives of the unknown variables,since the solution from this method based on the moving least square approximation is already smooth enough.Numerical examples are given to verify the accuracy of the meshless method proposed for elasto-plasticity analysis.
文摘Elasto-plastic consolidation is one of the classic coupling questions in geomechanics. To solve this problem, an elasto-plastic constitutive model is derived based on the numerical modeling method. The model is applied to Blot's consolidation theory. Incremental governing partial differential equations are established using this method. According to the stress path, the decoupling condition of these equations is discussed. Based on these conditions, an incremental diffusion equation and uncoupling governing equations are presented. The method is then applied to numerical analyses of three examples. The results show that (1) the effect of the stress path should be taken into account in the simulation of the soil consolidation question; (2) this decoupling method can predict the evolvement of pore water pressure; (3) the settlement using cam-clay model is less than that using numerical model because of dilatancy.
文摘The nonlinear quasi-conforming FEM is presented based on the basic concept of the quasi- -conforming finite element. First, the incremental principle of stationary potential energy is discussed, Then, the formulation process of the nonlinear quasi-conforming FEM is given. Lastly, two computational examples of shells are given.
文摘The suitability of six higher order root solvers is examined for solving the nonlinear equilibrium equations in large deformation analysis of structures.The applied methods have a better convergence rate than the quadratic Newton-Raphson method.These six methods do not require higher order derivatives to achieve a higher convergence rate.Six algorithms are developed to use the higher order methods in place of the Newton-Raphson method to solve the nonlinear equilibrium equations in geometrically nonlinear analysis of structures.The higher order methods are applied to both continuum and discrete problems(spherical shell and dome truss).The computational cost and the sensitivity of the higher order solution methods and the Newton-Raphson method with respect to the load increment size are comparatively investigated.The numerical results reveal that the higher order methods require a lower number of iterations that the Newton-Raphson method to converge.It is also shown that these methods are less sensitive to the variation of the load increment size.As it is indicated in numerical results,the average residual reduces in a lower number of iterations by the application of the higher order methods in the nonlinear analysis of structures.
基金This research has been supported by the National Natural ScienceFoundation of China (Grant No. 51778135 )the Natural Science Foundation of JiangsuProvince (Grant No. BK20160207)+1 种基金Aeronautical Science Foundation of China (GrantNo. 20130969010)the Fundamental Research Funds for the Central Universities andPostgraduate Research & Practice Innovation Program of Jiangsu Province (Grant No.KYCX18_0113 and KYLX16_0253).
文摘Current building design codes allow the appearance of structural and nonstructural damage under design basis earthquakes.The research regarding probabilistic seismic loss estimation in domestic building structure is urgent.The evaluation in this paper is based on a 11-story reinforced concrete office building,incremental dynamic analysis(IDA)is conducted in Perform 3D program using models capable to simulate all possible limit states up to collapse.Next,the probability distribution of post-earthquake casualties,rebuild costs repair costs and business downtime loss are calculated in PACT software for the studied building considering the modified component vulnerability groups and population flow models.The evaluation procedure can also shed light on other types of buildings in China.For non-typical functional building structures,this article proposes to build a finite element model of structural components and to classify the vulnerability groups based on the construction drawings,and to supply and improve the vulnerability library of appendages in FEMA P-58 according to the actual situation.In this way,the application scope of building seismic performance evaluation can be expanded.