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 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.展开更多
传统有限元法对大坝-不规则地基-库水系统进行建模时存在一定的局限性。基于ABAQUS二次开发接口,实现了比例边界有限元方法(scaled boundary finite element method,SBFEM)与八叉树网格的结合,建立了一种考虑真实地形的高拱坝-不规则地...传统有限元法对大坝-不规则地基-库水系统进行建模时存在一定的局限性。基于ABAQUS二次开发接口,实现了比例边界有限元方法(scaled boundary finite element method,SBFEM)与八叉树网格的结合,建立了一种考虑真实地形的高拱坝-不规则地基-库水系统自动建模方法。利用构建的八叉树比例边界有限元法对某重力坝地震响应进行了数值验证。随后对NG5拱坝系统分别基于平整地基和不规则地基进行线弹性和非线性动力响应分析。结果表明:在地震作用下,相较于简化的平整地基拱坝系统,不规则地基拱坝系统坝顶与坝底横河向相对位移以及第一主应力峰值变化较大,分别增加了73.5%和103.6%;考虑拱坝横缝以及材料非线性的情况下,坝顶与坝底横河向相对位移以及顺河向相对位移和相对速度分别增加了43.9%、32.0%和56.6%,同时边缝的法向开度增加尤为显著,增加了388.9%和381.8%,开度峰值增加了105%,在应力和损伤方面,第一主应力峰值增加了81.6%,损伤较大的区域也沿着坝体底部进行了扩展。展开更多
Variable stiffness composite laminates(VSCLs)are promising in aerospace engineering due to their designable material properties through changing fiber angles and stacking sequences.Aiming to control the thermal postbu...Variable stiffness composite laminates(VSCLs)are promising in aerospace engineering due to their designable material properties through changing fiber angles and stacking sequences.Aiming to control the thermal postbuckling and nonlinear panel flutter motions of VSCLs,a full-order numerical model is developed based on the linear quadratic regulator(LQR)algorithm in control theory,the classical laminate plate theory(CLPT)considering von Kármán geometrical nonlinearity,and the first-order Piston theory.The critical buckling temperature and the critical aerodynamic pressure of VSCLs are parametrically investigated.The location and shape of piezoelectric actuators for optimal control of the dynamic responses of VSCLs are determined through comparing the norms of feedback control gain(NFCG).Numerical simulations show that the temperature field has a great effect on aeroelastic tailoring of VSCLs;the curvilinear fiber path of VSCLs can significantly affect the optimal location and shape of piezoelectric actuator for flutter suppression;the unstable panel flutter and the thermal postbuckling deflection can be suppressed effectively through optimal design of piezoelectric patches.展开更多
Nonlinear static analysis of piezoelectric plates has been carried out using nonlinear finite element method considering electro-mechanical coupling,The geometrical nonlinearity has been taken into account and electri...Nonlinear static analysis of piezoelectric plates has been carried out using nonlinear finite element method considering electro-mechanical coupling,The geometrical nonlinearity has been taken into account and electric potential is assumed to be quadratic across the plate thickness,The governing equations are obtained using potential energy and Hamilton's principle that includes elastic and piezoelectric effects.The finite element model is derived based on constitutive equation of piezoelectric material accounting for coupling between elasticity and electric effect using higher order plate elements,Results are presented for piezoelectric plate under different mechanical boundary conditions,Numerical results for the plate are given in dimensionless graphical forms.Effects of boundary conditions on linear and nonlinear response of the plate are also studied.The numerical results obtained by the present model are in good agreement with the available solutions reported in the literature.展开更多
基金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.
基金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.
文摘传统有限元法对大坝-不规则地基-库水系统进行建模时存在一定的局限性。基于ABAQUS二次开发接口,实现了比例边界有限元方法(scaled boundary finite element method,SBFEM)与八叉树网格的结合,建立了一种考虑真实地形的高拱坝-不规则地基-库水系统自动建模方法。利用构建的八叉树比例边界有限元法对某重力坝地震响应进行了数值验证。随后对NG5拱坝系统分别基于平整地基和不规则地基进行线弹性和非线性动力响应分析。结果表明:在地震作用下,相较于简化的平整地基拱坝系统,不规则地基拱坝系统坝顶与坝底横河向相对位移以及第一主应力峰值变化较大,分别增加了73.5%和103.6%;考虑拱坝横缝以及材料非线性的情况下,坝顶与坝底横河向相对位移以及顺河向相对位移和相对速度分别增加了43.9%、32.0%和56.6%,同时边缝的法向开度增加尤为显著,增加了388.9%和381.8%,开度峰值增加了105%,在应力和损伤方面,第一主应力峰值增加了81.6%,损伤较大的区域也沿着坝体底部进行了扩展。
基金Project(JCYJ20190808175801656)supported by the Science and Technology Innovation Commission of Shenzhen,ChinaProject(2021M691427)supported by Postdoctoral Science Foundation of ChinaProject(9680086)supported by the City University of Hong Kong,China。
文摘Variable stiffness composite laminates(VSCLs)are promising in aerospace engineering due to their designable material properties through changing fiber angles and stacking sequences.Aiming to control the thermal postbuckling and nonlinear panel flutter motions of VSCLs,a full-order numerical model is developed based on the linear quadratic regulator(LQR)algorithm in control theory,the classical laminate plate theory(CLPT)considering von Kármán geometrical nonlinearity,and the first-order Piston theory.The critical buckling temperature and the critical aerodynamic pressure of VSCLs are parametrically investigated.The location and shape of piezoelectric actuators for optimal control of the dynamic responses of VSCLs are determined through comparing the norms of feedback control gain(NFCG).Numerical simulations show that the temperature field has a great effect on aeroelastic tailoring of VSCLs;the curvilinear fiber path of VSCLs can significantly affect the optimal location and shape of piezoelectric actuator for flutter suppression;the unstable panel flutter and the thermal postbuckling deflection can be suppressed effectively through optimal design of piezoelectric patches.
文摘Nonlinear static analysis of piezoelectric plates has been carried out using nonlinear finite element method considering electro-mechanical coupling,The geometrical nonlinearity has been taken into account and electric potential is assumed to be quadratic across the plate thickness,The governing equations are obtained using potential energy and Hamilton's principle that includes elastic and piezoelectric effects.The finite element model is derived based on constitutive equation of piezoelectric material accounting for coupling between elasticity and electric effect using higher order plate elements,Results are presented for piezoelectric plate under different mechanical boundary conditions,Numerical results for the plate are given in dimensionless graphical forms.Effects of boundary conditions on linear and nonlinear response of the plate are also studied.The numerical results obtained by the present model are in good agreement with the available solutions reported in the literature.
