The element energy projection (EEP) method for computation of super- convergent resulting in a one-dimensional finite element method (FEM) is successfully used to self-adaptive FEM analysis of various linear probl...The element energy projection (EEP) method for computation of super- convergent resulting in a one-dimensional finite element method (FEM) is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton's method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation (ODE) of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.展开更多
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 new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization...A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.展开更多
The structural damage identification through modal data often leads to solving a set of linear equations. Special numerical treatment is sometimes required for an accurate and stable solution owing to the ill conditio...The structural damage identification through modal data often leads to solving a set of linear equations. Special numerical treatment is sometimes required for an accurate and stable solution owing to the ill conditioning of the equations. Based on the singular value decomposition (SVD) of the coefficient matrix, an error based truncation algorithm is proposed in this paper. By rejection of selected small singular values, the influence of noise can be reduced. A simply-supported beam is used as a simulation example to compare the results to other methods. Illustrative numerical examples demonstrate the good efficiency and stability of the algorithm in the nondestructive identification of structural damage through modal data.展开更多
The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathe...The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathematical descriptions,namely,limit state functions of failure modes.Several problems are to be solved in the use of traditional methods for gravity dams.One is how to define the limit state function really reflecting the mechanical mechanism of the failure mode;another is how to understand the relationship among failure modes and enable the probability of the whole structure to be determined.Performing SFP analysis for a gravity dam system is a challenging task.This work proposes a novel nonlinear finite-element-based SFP analysis method for gravity dams.Firstly,reasonable nonlinear constitutive modes for dam concrete,concrete/rock interface and rock foundation are respectively introduced according to corresponding mechanical mechanisms.Meanwhile the response surface(RS) method is used to model limit state functions of main failure modes through the Monte Carlo(MC) simulation results of the dam-interface-foundation interaction finite element(FE) analysis.Secondly,a numerical SFP method is studied to compute the probabilities of several failure modes efficiently by simple matrix integration operations.Then,the nonlinear FE-based SFP analysis methodology for gravity dams considering correlated failure modes with the additional sensitivity analysis is proposed.Finally,a comprehensive computational platform for interfacing the proposed method with the open source FE code Code Aster is developed via a freely available MATLAB software tool(FERUM).This methodology is demonstrated by a case study of an existing gravity dam analysis,in which the dominant failure modes are identified,and the corresponding performance functions are established.Then,the dam failure probability of the structural system is obtained by the proposed method considering the correlation relationship of main failure modes on the basis of the mechanical mechanism analysis with the MC-FE simulations.展开更多
Anonlinearfinite elementprogram takinginto accountofthe nonlinearshearstress strainrelationisused to simulatethe progressive damage ofcompositelaminates with a circularhole subjectedtotensileloading.Stress, strain,...Anonlinearfinite elementprogram takinginto accountofthe nonlinearshearstress strainrelationisused to simulatethe progressive damage ofcompositelaminates with a circularhole subjectedtotensileloading.Stress, strain,and ultimate strength ofthe laminate are analyzed. Comparison between calculated result and test data showsthatthe nonlinear approach gives more accurateresult,especiallyfororthotropiclaminate .Parametricstudy showsthatthe variation ofnonlinearfactorhas more effectontheresultfor orthotropic platethanforquasi isotropic plate.展开更多
The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling ...The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.展开更多
We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is ...We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.展开更多
In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element...In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.展开更多
This paper presents the application of anisotropic damage theory to the study of forming limit diagram of A12024T3 aluminum alloy sheet. In the prediction of limiting strains of the aluminum sheet structure, a finite ...This paper presents the application of anisotropic damage theory to the study of forming limit diagram of A12024T3 aluminum alloy sheet. In the prediction of limiting strains of the aluminum sheet structure, a finite element cell model has been constructed. The cell model consists of two phases, the aluminum alloy matrix and the intermetallic cluster. The material behavior of the aluminum alloy matrix is described with a fully coupled elasto-plastic damage constitutive equation. The intermetallic cluster is assumed to be elastic and brittle. By varying the stretching ratio, the limiting strains of the sheet under biaxial stretching have been predicted by using the necking criterion proposed. The prediction is in good agreement with the experimental findings. Moreover, the finite element cell model can provide information for understanding the microscopic damage mechanism of the aluminum alloy. Over-estimation of the limit strains may result if the effect of material damage is ignored in the sheet metal forming study.展开更多
Traditional fracture analysis is based on fracture mechanics and damage mechanics. They focus on the propagation of the fracture. However, their propagation criterions are not easily applied in practice and the curren...Traditional fracture analysis is based on fracture mechanics and damage mechanics. They focus on the propagation of the fracture. However, their propagation criterions are not easily applied in practice and the current analysis is limited in planar problem. This paper presents a new theory that the occurrence of the unbalanced force (derived from the Deformation Reinforcement Theory) could be the criterion of the initiation of the fracture, and the distribution area and propagation of the unbalanced force could be the indication of the fracture propagation direction. By aggregate analysis with Stress Intensity Factor (SIF) criterion, the unbalanced force actually is the opposite external load that is the SIF difference incurred between the external loads and permitted by the structure. Numerical simulation and physical experiments on pre-fracture cuboid rock specimens proved that the occurrence of the unbalanced force could be the initiation of the fracture. Mesh size dependence was also considered by analysis of different mesh size finite element gravity dam models. Furthermore, the theory was applied to the feasibility analysis of the Baihetan arch dam together with physical experiments in order to evaluate the fracture propagation of dam heel. The results show that it is an effective way to use unbalanced force to analyze the fracture initiation and propagation when performing 3-dimensional nonlinear FEM calculation.展开更多
The offshore reinforced concrete structures are always subject to cyclic load, such as wave load.In this paper a new finite element analysis model is developed to analyze the stress and strain state of reinforced conc...The offshore reinforced concrete structures are always subject to cyclic load, such as wave load.In this paper a new finite element analysis model is developed to analyze the stress and strain state of reinforced concrete structures including offshore concrete structures, subject to any number of the cyclic load. On the basis of the anal ysis of the experimental data,this model simplifies the number of cycles-total cyclic strain curve of concrete as three straight line segments,and it is assumed that the stress-strain curves of different cycles in each segment are the same, thus the elastoplastic analysis is only needed for the first cycle of each segment, and the stress or strain corresponding to any number of cycles can be obtained by superposition of stress or strain obtained by the above e lastoplastic analysis based on the cyclic numbers in each segment.This model spends less computer time,and can obtain the stress and strain states of the structures after any number of cycles.The endochronic-damage and ideal offshore concrete platform subject to cyclic loading are experimented and analyzed by the finite element method based on the model proposed in this paper. The results between the experiment and the finite element analysis are in good agreement,which demonstrates the validity and accuracy of the proposed model.展开更多
A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as th...A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.展开更多
The nonlinear analysis of reinforced concrete rectangular slabs undermonotonic transverse loads is performed by finite element method.The layered rectangu-lar element with 4 nodes and 20 degrees of freedom is develope...The nonlinear analysis of reinforced concrete rectangular slabs undermonotonic transverse loads is performed by finite element method.The layered rectangu-lar element with 4 nodes and 20 degrees of freedom is developed,in whichbending-stretching coupling effect is taken into account.An orthotropic equivalentuniaxial stress-strain constitutive model of concrete is used.A program is worked out andused to calculate two reinforced concrete slabs.The results of calculation are in goodconformity with the corresponding test results.In addition,the influence of tension stif-fening effect of cracked concrete on the results of calculation is discussed.展开更多
In order to analyze bellows effectively and practically, the finite_element_displacement_perturbation method (FEDPM) is proposed for the geometric nonlinear behaviors of shells of revolution subjected to pure bending ...In order to analyze bellows effectively and practically, the finite_element_displacement_perturbation method (FEDPM) is proposed for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes. The formulations are mainly based upon the idea of perturba_ tion that the nodal displacement vector and the nodal force vector of each finite element are expanded by taking root_mean_square value of circumferential strains of the shells as a perturbation parameter. The load steps and the iteration times are not as arbitrary and unpredictable as in usual nonlinear analysis. Instead, there are certain relations between the load steps and the displacement increments, and no need of iteration for each load step. Besides, in the formulations, the shell is idealized into a series of conical frusta for the convenience of practice, Sander's nonlinear geometric equations of moderate small rotation are used, and the shell made of more than one material ply is also considered.展开更多
Some theoretical methods have been reported to deal with nonlinear problems of composite materials but the accuracy is not so good. In the meantime, a lot of linear problems are difficult to be managed by the theoreti...Some theoretical methods have been reported to deal with nonlinear problems of composite materials but the accuracy is not so good. In the meantime, a lot of linear problems are difficult to be managed by the theoretical methods. The present study aims to use the developed method, the random microstructure finite element method, to deal with these nonlinear problems. In this paper, the random microstructure finite element method is used to deal with all three kinds of nonlinear property problems of composite materials. The analyzed results suggest the influences of the nonlinear phenomena on the effective properties of composite materials are significant and the random microstructure finite element method is an effective tool to investigate the nonlinear problems.展开更多
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.展开更多
In this paper a nonlinear response of a fixed offshore platform under the combined forces of waves, wind and sea currents is presented. Wave force acting on the elements is calculated using Morison equation. Hydrodyna...In this paper a nonlinear response of a fixed offshore platform under the combined forces of waves, wind and sea currents is presented. Wave force acting on the elements is calculated using Morison equation. Hydrodynamic loads on horizontal and vertical tubular members and the dynamic response of offshore fixed platform coupled with distribution of displacement, axial force, and bending moment along the base of the platform for regular and severe cases have been investigated. The structure must be able maintain production in a one-year wave return period condition and also to be able to continue with one hundred-year storm return period. The results of this study show that bending moment values with a one-year wave return period condition for the base platform and junction of platform to deck are 70 percent and 59 percent, respectively more than bending moment with a one-year wave return period. The direction of wave and wind hit has significant effects on the shift platform response, also nonlinear response is important for the safe design and operation of offshore structures.展开更多
The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite ...The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.展开更多
In this paper, the finite element method using vector potential in applications to 2D nonlinear eddy current field is discussed. The authors use the equivalent magnetic energy method to deal with magnetization curve o...In this paper, the finite element method using vector potential in applications to 2D nonlinear eddy current field is discussed. The authors use the equivalent magnetic energy method to deal with magnetization curve of ferromagnetic material,and present the formulation of 2D nonlinear eddy current field.With this method the authors analyze the eddy current field in an induction ladle furnace and the force distribution in the charge (molten metal),and plot the corresponding curves.展开更多
基金supported by the National Natural Science Foundation of China(Nos.51378293,51078199,50678093,and 50278046)the Program for Changjiang Scholars and the Innovative Research Team in University of China(No.IRT00736)
文摘The element energy projection (EEP) method for computation of super- convergent resulting in a one-dimensional finite element method (FEM) is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton's method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation (ODE) of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.
基金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.
基金supported by a Major Research Project in Higher Education Institutions in Henan Province,with Project Number 23A560015.
文摘A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.
文摘The structural damage identification through modal data often leads to solving a set of linear equations. Special numerical treatment is sometimes required for an accurate and stable solution owing to the ill conditioning of the equations. Based on the singular value decomposition (SVD) of the coefficient matrix, an error based truncation algorithm is proposed in this paper. By rejection of selected small singular values, the influence of noise can be reduced. A simply-supported beam is used as a simulation example to compare the results to other methods. Illustrative numerical examples demonstrate the good efficiency and stability of the algorithm in the nondestructive identification of structural damage through modal data.
基金Projects(51409167,51139001,51179066)supported by the National Natural Science Foundation of ChinaProjects(201401022,201501036)supported by the Ministry of Water Resources Public Welfare Industry Research Special Fund,ChinaProjects(GG201532,GG201546)supported by the Scientific and Technological Research for Water Conservancy,Henan Province,China
文摘The structural system failure probability(SFP) is a valuable tool for evaluating the global safety level of concrete gravity dams.Traditional methods for estimating the failure probabilities are based on defined mathematical descriptions,namely,limit state functions of failure modes.Several problems are to be solved in the use of traditional methods for gravity dams.One is how to define the limit state function really reflecting the mechanical mechanism of the failure mode;another is how to understand the relationship among failure modes and enable the probability of the whole structure to be determined.Performing SFP analysis for a gravity dam system is a challenging task.This work proposes a novel nonlinear finite-element-based SFP analysis method for gravity dams.Firstly,reasonable nonlinear constitutive modes for dam concrete,concrete/rock interface and rock foundation are respectively introduced according to corresponding mechanical mechanisms.Meanwhile the response surface(RS) method is used to model limit state functions of main failure modes through the Monte Carlo(MC) simulation results of the dam-interface-foundation interaction finite element(FE) analysis.Secondly,a numerical SFP method is studied to compute the probabilities of several failure modes efficiently by simple matrix integration operations.Then,the nonlinear FE-based SFP analysis methodology for gravity dams considering correlated failure modes with the additional sensitivity analysis is proposed.Finally,a comprehensive computational platform for interfacing the proposed method with the open source FE code Code Aster is developed via a freely available MATLAB software tool(FERUM).This methodology is demonstrated by a case study of an existing gravity dam analysis,in which the dominant failure modes are identified,and the corresponding performance functions are established.Then,the dam failure probability of the structural system is obtained by the proposed method considering the correlation relationship of main failure modes on the basis of the mechanical mechanism analysis with the MC-FE simulations.
文摘Anonlinearfinite elementprogram takinginto accountofthe nonlinearshearstress strainrelationisused to simulatethe progressive damage ofcompositelaminates with a circularhole subjectedtotensileloading.Stress, strain,and ultimate strength ofthe laminate are analyzed. Comparison between calculated result and test data showsthatthe nonlinear approach gives more accurateresult,especiallyfororthotropiclaminate .Parametricstudy showsthatthe variation ofnonlinearfactorhas more effectontheresultfor orthotropic platethanforquasi isotropic plate.
基金The Construction S&T Project of the Department of Transportation of Sichuan Province(Grant No.2023A02)the National Natural Science Foundation of China(No.52109135).
文摘The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.
基金supported by the National Natural Science Foundation of China(Nos.12132001 and 52192632)。
文摘We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.
文摘In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.
基金Project supported by the Research Committee of The Hong Kong Polytechnic University (No.G-YX34).
文摘This paper presents the application of anisotropic damage theory to the study of forming limit diagram of A12024T3 aluminum alloy sheet. In the prediction of limiting strains of the aluminum sheet structure, a finite element cell model has been constructed. The cell model consists of two phases, the aluminum alloy matrix and the intermetallic cluster. The material behavior of the aluminum alloy matrix is described with a fully coupled elasto-plastic damage constitutive equation. The intermetallic cluster is assumed to be elastic and brittle. By varying the stretching ratio, the limiting strains of the sheet under biaxial stretching have been predicted by using the necking criterion proposed. The prediction is in good agreement with the experimental findings. Moreover, the finite element cell model can provide information for understanding the microscopic damage mechanism of the aluminum alloy. Over-estimation of the limit strains may result if the effect of material damage is ignored in the sheet metal forming study.
基金supported by the National Natural Science Foundation of China (Grant No. 50709014)China National Funds for Distinguished Young Scientists (Grant No. 50925931)State Key Laboratory of Hydroscience and Engineering of China (Grant No. 2008-TC-2)
文摘Traditional fracture analysis is based on fracture mechanics and damage mechanics. They focus on the propagation of the fracture. However, their propagation criterions are not easily applied in practice and the current analysis is limited in planar problem. This paper presents a new theory that the occurrence of the unbalanced force (derived from the Deformation Reinforcement Theory) could be the criterion of the initiation of the fracture, and the distribution area and propagation of the unbalanced force could be the indication of the fracture propagation direction. By aggregate analysis with Stress Intensity Factor (SIF) criterion, the unbalanced force actually is the opposite external load that is the SIF difference incurred between the external loads and permitted by the structure. Numerical simulation and physical experiments on pre-fracture cuboid rock specimens proved that the occurrence of the unbalanced force could be the initiation of the fracture. Mesh size dependence was also considered by analysis of different mesh size finite element gravity dam models. Furthermore, the theory was applied to the feasibility analysis of the Baihetan arch dam together with physical experiments in order to evaluate the fracture propagation of dam heel. The results show that it is an effective way to use unbalanced force to analyze the fracture initiation and propagation when performing 3-dimensional nonlinear FEM calculation.
文摘The offshore reinforced concrete structures are always subject to cyclic load, such as wave load.In this paper a new finite element analysis model is developed to analyze the stress and strain state of reinforced concrete structures including offshore concrete structures, subject to any number of the cyclic load. On the basis of the anal ysis of the experimental data,this model simplifies the number of cycles-total cyclic strain curve of concrete as three straight line segments,and it is assumed that the stress-strain curves of different cycles in each segment are the same, thus the elastoplastic analysis is only needed for the first cycle of each segment, and the stress or strain corresponding to any number of cycles can be obtained by superposition of stress or strain obtained by the above e lastoplastic analysis based on the cyclic numbers in each segment.This model spends less computer time,and can obtain the stress and strain states of the structures after any number of cycles.The endochronic-damage and ideal offshore concrete platform subject to cyclic loading are experimented and analyzed by the finite element method based on the model proposed in this paper. The results between the experiment and the finite element analysis are in good agreement,which demonstrates the validity and accuracy of the proposed model.
文摘A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.
文摘The nonlinear analysis of reinforced concrete rectangular slabs undermonotonic transverse loads is performed by finite element method.The layered rectangu-lar element with 4 nodes and 20 degrees of freedom is developed,in whichbending-stretching coupling effect is taken into account.An orthotropic equivalentuniaxial stress-strain constitutive model of concrete is used.A program is worked out andused to calculate two reinforced concrete slabs.The results of calculation are in goodconformity with the corresponding test results.In addition,the influence of tension stif-fening effect of cracked concrete on the results of calculation is discussed.
文摘In order to analyze bellows effectively and practically, the finite_element_displacement_perturbation method (FEDPM) is proposed for the geometric nonlinear behaviors of shells of revolution subjected to pure bending moments or lateral forces in one of their meridional planes. The formulations are mainly based upon the idea of perturba_ tion that the nodal displacement vector and the nodal force vector of each finite element are expanded by taking root_mean_square value of circumferential strains of the shells as a perturbation parameter. The load steps and the iteration times are not as arbitrary and unpredictable as in usual nonlinear analysis. Instead, there are certain relations between the load steps and the displacement increments, and no need of iteration for each load step. Besides, in the formulations, the shell is idealized into a series of conical frusta for the convenience of practice, Sander's nonlinear geometric equations of moderate small rotation are used, and the shell made of more than one material ply is also considered.
基金This work is supported by the National Natural Science Foundation of China under the Grant 19772037 and 19902014
文摘Some theoretical methods have been reported to deal with nonlinear problems of composite materials but the accuracy is not so good. In the meantime, a lot of linear problems are difficult to be managed by the theoretical methods. The present study aims to use the developed method, the random microstructure finite element method, to deal with these nonlinear problems. In this paper, the random microstructure finite element method is used to deal with all three kinds of nonlinear property problems of composite materials. The analyzed results suggest the influences of the nonlinear phenomena on the effective properties of composite materials are significant and the random microstructure finite element method is an effective tool to investigate the nonlinear problems.
文摘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.
文摘In this paper a nonlinear response of a fixed offshore platform under the combined forces of waves, wind and sea currents is presented. Wave force acting on the elements is calculated using Morison equation. Hydrodynamic loads on horizontal and vertical tubular members and the dynamic response of offshore fixed platform coupled with distribution of displacement, axial force, and bending moment along the base of the platform for regular and severe cases have been investigated. The structure must be able maintain production in a one-year wave return period condition and also to be able to continue with one hundred-year storm return period. The results of this study show that bending moment values with a one-year wave return period condition for the base platform and junction of platform to deck are 70 percent and 59 percent, respectively more than bending moment with a one-year wave return period. The direction of wave and wind hit has significant effects on the shift platform response, also nonlinear response is important for the safe design and operation of offshore structures.
基金This work was supported in part by the National Science Foundation under grant DMS-1620288。
文摘The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.
文摘In this paper, the finite element method using vector potential in applications to 2D nonlinear eddy current field is discussed. The authors use the equivalent magnetic energy method to deal with magnetization curve of ferromagnetic material,and present the formulation of 2D nonlinear eddy current field.With this method the authors analyze the eddy current field in an induction ladle furnace and the force distribution in the charge (molten metal),and plot the corresponding curves.