A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination...In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.展开更多
The present study explored the influence of axial gradation of viscoelastic materials on the dynamic response of the sandwich beam for structural applications.The finite element(FE)formulations are used to model and i...The present study explored the influence of axial gradation of viscoelastic materials on the dynamic response of the sandwich beam for structural applications.The finite element(FE)formulations are used to model and investigate dynamic response of the sandwich beam.The classical beam theory is used to develop the FE formulations and Lagrange's approach is considered to obtain the equations of motion(EOM).FE code is developed and validated with the existing literature and also conducted the convergence study for the developed FE method.Further,the influence of different viscoelastic materials and boundary conditions on the dynamic response of the sandwich beam is investigated.Four different axial gradation configurations of viscoelastic materials are considered for the present work to explore the influence on natural frequency,loss factor and frequency response of the sandwich beam.The modeled axial gradation of viscoelastic material has displayed a considerable impact on the peak vibrational amplitude response of the sandwich beam for all the boundary conditions and these configurations improved the damping capabilities at different configurations for the structural applications.展开更多
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金Project(2017YFC0404802)supported by the National Key R&D Program of ChinaProjects(U1965206,51979143)supported by the National Natural Science Foundation of China。
文摘In the finite element method,the numerical simulation of three-dimensional crack propagation is relatively rare,and it is often realized by commercial programs.In addition to the geometric complexity,the determination of the cracking direction constitutes a great challenge.In most cases,the local stress state provides the fundamental criterion to judge the presence of cracks and the direction of crack propagation.However,in the case of three-dimensional analysis,the coordination relationship between grid elements due to occurrence of cracks becomes a difficult problem for this method.In this paper,based on the extended finite element method,the stress-related function field is introduced into the calculation domain,and then the boundary value problem of the function is solved.Subsequently,the envelope surface of all propagation directions can be obtained at one time.At last,the possible surface can be selected as the direction of crack development.Based on the aforementioned procedure,such method greatly reduces the programming complexity of tracking the crack propagation.As a suitable method for simulating tension-induced failure,it can simulate multiple cracks simultaneously.
基金support from the Department of Science and Technology (DST)file no.ECR/2016/001448 titled“Experimental Investigation of Passive,Semi-active and Active vibration control of Composite Sandwich Structure”funded by Science and Engineering Research Board,Government of India。
文摘The present study explored the influence of axial gradation of viscoelastic materials on the dynamic response of the sandwich beam for structural applications.The finite element(FE)formulations are used to model and investigate dynamic response of the sandwich beam.The classical beam theory is used to develop the FE formulations and Lagrange's approach is considered to obtain the equations of motion(EOM).FE code is developed and validated with the existing literature and also conducted the convergence study for the developed FE method.Further,the influence of different viscoelastic materials and boundary conditions on the dynamic response of the sandwich beam is investigated.Four different axial gradation configurations of viscoelastic materials are considered for the present work to explore the influence on natural frequency,loss factor and frequency response of the sandwich beam.The modeled axial gradation of viscoelastic material has displayed a considerable impact on the peak vibrational amplitude response of the sandwich beam for all the boundary conditions and these configurations improved the damping capabilities at different configurations for the structural applications.
