The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess v...The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.展开更多
A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FE...A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FEM is described and compared with ES-FEM.A practical modification of IES-FEM is then introduced that used the technique employed by ES-FEM for the nodal strain calculation.The differences in the strain computation among ES-FEM,IES-FEM,and FEM are then discussed.The modified IES-FEM exhibited superior performance in displacement and a slight advantage in stress compared to FEM using the same mesh according to the results obtained from both the regular and irregular elements.The robustness of the IES-FEM to severely deformed meshes was also verified.展开更多
The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES...The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.展开更多
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 work presents a locking-free smoothed finite element method(S-FEM)for the simulation of soft matter modelled by the equations of quasi-incompressible hyperelasticity.The proposed method overcomes well-known issue...This work presents a locking-free smoothed finite element method(S-FEM)for the simulation of soft matter modelled by the equations of quasi-incompressible hyperelasticity.The proposed method overcomes well-known issues of standard finite element methods(FEM)in the incompressible limit:the over-estimation of stiffness and sensitivity to severely distorted meshes.The concepts of cell-based,edge-based and node-based S-FEMs are extended in this paper to three-dimensions.Additionally,a cubic bubble function is utilized to improve accuracy and stability.For the bubble function,an additional displacement degree of freedom is added at the centroid of the element.Several numerical studies are performed demonstrating the stability and validity of the proposed approach.The obtained results are compared with standard FEM and with analytical solutions to show the effectiveness of the method.展开更多
The aim of this work is to employ a modified cell-based smoothed finite element method(S-FEM)for topology optimization with the domain discretized with arbitrary polygons.In the present work,the linear polynomial basi...The aim of this work is to employ a modified cell-based smoothed finite element method(S-FEM)for topology optimization with the domain discretized with arbitrary polygons.In the present work,the linear polynomial basis function is used as the weight function instead of the constant weight function used in the standard S-FEM.This improves the accuracy and yields an optimal convergence rate.The gradients are smoothed over each smoothing domain,then used to compute the stiffness matrix.Within the proposed scheme,an optimum topology procedure is conducted over the smoothing domains.Structural materials are distributed over each smoothing domain and the filtering scheme relies on the smoothing domain.Numerical tests are carried out to pursue the performance of the proposed optimization by comparing convergence,efficiency and accuracy.展开更多
A four-node quadrilateral shell element with smoothed membrane-bending based on Mindlin-Reissner theory is proposed. The element is a combination of a plate bending and membrane element. It is based on mixed interpola...A four-node quadrilateral shell element with smoothed membrane-bending based on Mindlin-Reissner theory is proposed. The element is a combination of a plate bending and membrane element. It is based on mixed interpolation where the bending and membrane stiffness matrices are calculated on the boundaries of the smoothing cells while the shear terms are approximated by independent interpolation functions in natural coordinates. The proposed element is robust, computationally inexpensive and free of locking. Since the integration is done on the element boundaries for the bending and membrane terms, the element is more accurate than the MITC4 element for distorted meshes. This will be demonstrated for several numerical examples.展开更多
Formulation and numerical evaluation of a novel twice-interpolation finite element method （TFEM） is presented for solid mechanics problems. In this method, the trial function for Galerkin weak form is constructed th...Formulation and numerical evaluation of a novel twice-interpolation finite element method （TFEM） is presented for solid mechanics problems. In this method, the trial function for Galerkin weak form is constructed through two stages of consecutive interpolation. The primary interpolation follows exactly the same procedure of standard FEM and is further reproduced according to both nodal values and averaged nodal gradients obtained from primary interpolation. The trial functions thus constructed have continuous nodal gradients and contain higher order polynomial without increasing total freedoms. Several benchmark examples and a real dam problem are used to examine the TFEM in terms of accuracy and convergence. Compared with standard FEM, TFEM can achieve significantly better accuracy and higher convergence rate, and the continuous nodal stress can be obtained without any smoothing operation. It is also found that TFEM is insensitive to the quality of the elemental mesh. In addition, the present TFEM can treat the incompressible material without any modification.展开更多
Optical full-field measurement methods are now widely applied in various domains.In general,the displacement fields can be directly obtained from the measurement,however in mechanical analysis strain fields are prefer...Optical full-field measurement methods are now widely applied in various domains.In general,the displacement fields can be directly obtained from the measurement,however in mechanical analysis strain fields are preferred.To extract strain fields from noisy displacement fields is always a challenging topic.In this study,a finite element method for smoothing displacement fields and calculating strain fields is proposed.An experimental test case on a holed aluminum specimen under tension is applied to validate this method.The heterogeneous displacement fields are measured by digital image correlation(DIC).By this proposed method,the result shows that the measuring noise on experimental displacement fields can be successfully removed,and strain fields can be reconstructed in the arbitrary area.展开更多
The smoothed finite element method (S-FEM) was originated by G R Liu by combining some meshfree techniques with the well-established standard finite element method (FEM). It has a family of models carefully designed w...The smoothed finite element method (S-FEM) was originated by G R Liu by combining some meshfree techniques with the well-established standard finite element method (FEM). It has a family of models carefully designed with innovative types of smoothing domains. These models are found having a number of important and theoretically profound properties. This article first provides a concise and easy-to-follow presentation of key formulations used in the S-FEM. A number of important properties and unique features of S-FEM models are discussed in detail, including 1) theoretically proven softening effects;2) upper-bound solutions;3) accurate solutions and higher convergence rates;4) insensitivity to mesh distortion;5) Jacobian?free;6) volumetric-locking-free;and most importantly 7) working well with triangular and tetrahedral meshes that can be automatically generated. The S-FEM is thus ideal for automation in computations and adaptive analyses, and hence has profound impact on Al-assisted modeling and simulation. Most importantly, one can now purposely design an S-FEM model to obtain solutions with special properties as wish, meaning that S-FEM offers a framework for design numerical models with desired properties. This novel concept of numerical model demand may drastically change the landscape of modeling and simulation. Future directions of research are also provided.展开更多
The paper examines three selective schemes for the smoothed finite element method （SFEM） which was formulated by incorporating a cell-wise strain smoothing operation into the standard compatible finite element meth...The paper examines three selective schemes for the smoothed finite element method （SFEM） which was formulated by incorporating a cell-wise strain smoothing operation into the standard compatible finite element method （FEM）. These selective SFEM schemes were formulated based on three selective integration FEM schemes with similar properties found between the number of smoothing cells in the SFEM and the number of Gaussian integration points in the FEM. Both scheme 1 and scheme 2 are free of nearly incompressible locking, but scheme 2 is more general and gives better results than scheme 1. In addition, scheme 2 can be applied to anisotropic and nonlinear situations, while scheme 1 can only be applied to isotropic and linear situations. Scheme 3 is free of shear locking. This scheme can be applied to plate and shell problems. Results of the numerical study show that the selective SFEM schemes give more accurate results than the FEM schemes.展开更多
Free-surface flows, especially those associated with fluid-structure interactions(FSIs), pose challenging problems in numerical simulations. The authors of this work recently developed a smoothed particle element meth...Free-surface flows, especially those associated with fluid-structure interactions(FSIs), pose challenging problems in numerical simulations. The authors of this work recently developed a smoothed particle element method(SPEM) to simulate FSIs. In this method, both the fluid and solid regions are initially modeled using a smoothed finite element method(S-FEM) in a Lagrangian frame, whereas the fluid regions undergoing large deformations are adaptively converted into particles and modeled with an improved smoothed particle hydrodynamics(SPH) method. This approach greatly improves computational accuracy and efficiency because of the advantages of the S-FEM in efficiently treating solid/fluid regions showing small deformations and the SPH method in effectively modeling moving interfaces. In this work, we further enhance the efficiency of the SPEM while effectively capturing local fluid information by introducing a multi-resolution technique to the SPEM and developing an effective approach to treat multi-resolution element-particle interfaces. Various numerical examples demonstrate that the multiresolution SPEM can significantly reduce the computational cost relative to the original version with a constant resolution.Moreover, the novel approach is effective in modeling various incompressible flow problems involving FSIs.展开更多
An effective hybrid optimization method is proposed by integrating an adaptive Kriging(A-Kriging)into an improved partial swarm optimization algorithm(IPSO)to give a so-called A-Kriging-IPSO for maximizing the bucklin...An effective hybrid optimization method is proposed by integrating an adaptive Kriging(A-Kriging)into an improved partial swarm optimization algorithm(IPSO)to give a so-called A-Kriging-IPSO for maximizing the buckling load of laminated composite plates(LCPs)under uniaxial and biaxial compressions.In this method,a novel iterative adaptive Kriging model,which is structured using two training sample sets as active and adaptive points,is utilized to directly predict the buckling load of the LCPs and to improve the efficiency of the optimization process.The active points are selected from the initial data set while the adaptive points are generated using the radial random-based convex samples.The cell-based smoothed discrete shear gap method(CS-DSG3)is employed to analyze the buckling behavior of the LCPs to provide the response of adaptive and input data sets.The buckling load of the LCPs is maximized by utilizing the IPSO algorithm.To demonstrate the efficiency and accuracy of the proposed methodology,the LCPs with different layers(2,3,4,and 10 layers),boundary conditions,aspect ratios and load patterns(biaxial and uniaxial loads)are investigated.The results obtained by proposed method are in good agreement with the literature results,but with less computational burden.By applying adaptive radial Kriging model,the accurate optimal resultsebased predictions of the buckling load are obtained for the studied LCPs.展开更多
A strain smoothing formulation for NURBS (non-uniform rational B-splines) based isogeometric finite element analysis is presented. This approach is formulated within the framework of assumed strain methods and strain ...A strain smoothing formulation for NURBS (non-uniform rational B-splines) based isogeometric finite element analysis is presented. This approach is formulated within the framework of assumed strain methods and strain smoothing operations. The strain smoothing is defined through strain averaging in the element sub-domains which are subsequently used for numerical integration of the Galerkin weak form. This formulation satisfies the orthogonality condition of the assumed strain methods. Meanwhile the present formulation totally avoids the gradient computation of the rational NURBS basis functions in the formulation of stiffness matrix. A transformation method is employed to accurately enforce the displacement boundary conditions. Numerical results demonstrate that the present formation gives very satisfactory solution accuracy simultaneously with improved computational efficiency.展开更多
近年来,部分特高压换流变压器在运行中相继发生电弧短路故障并引发爆炸、起火事故,严重威胁直流系统安全稳定运行。充油设备故障冲击下的结构失效机制尚不明确,且缺乏成熟的数值计算方法,制约故障防爆技术的发展。基于此,该文提出一套...近年来,部分特高压换流变压器在运行中相继发生电弧短路故障并引发爆炸、起火事故,严重威胁直流系统安全稳定运行。充油设备故障冲击下的结构失效机制尚不明确,且缺乏成熟的数值计算方法,制约故障防爆技术的发展。基于此,该文提出一套适用于高能电弧故障冲击的结构失效仿真计算方法。首先,建立有限腔体内油中电弧能量持续注入的气泡动力学模型,准确描述故障气泡的脉动膨胀行为;其次,提出自适应有限元-光滑粒子流体动力学(finite element method-smoothed particle hydrodynamics,FEM-SPH)耦合方法,利用SPH粒子继承失效前的物理信息参与FEM计算;进行不同能量、不同位置的电弧故障仿真计算,获得换流变压器结构的薄弱区域及其破裂行为,复现了油箱结构失效行为。研究发现,油箱顶盖两侧及侧壁转角接缝位置容易发生应力集中现象。一旦结构出现裂缝,将在极短时间内沿着应力集中方向快速发展,最终导致壁面整体撕裂。可知,该方法揭示的结构失效行为可为改进变压器设计和提高设备安全性提供依据。展开更多
基金part of the TPS projecta Vied-Newton PhD scholarship+1 种基金a Dixon scholarship from Imperial College London,UKthe Dean’s Fund from Imperial College London for financial support(2017-2020)。
文摘The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.
基金the National Natural Science Foundation of China(No.11672238)the 111 Project(No.BP0719007)the Shaanxi Province Natural Science Foundation(No.2020JZ-06)for the financial support.
文摘A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FEM is described and compared with ES-FEM.A practical modification of IES-FEM is then introduced that used the technique employed by ES-FEM for the nodal strain calculation.The differences in the strain computation among ES-FEM,IES-FEM,and FEM are then discussed.The modified IES-FEM exhibited superior performance in displacement and a slight advantage in stress compared to FEM using the same mesh according to the results obtained from both the regular and irregular elements.The robustness of the IES-FEM to severely deformed meshes was also verified.
基金funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under Grant number 107.02-2019.330。
文摘The main purpose of this paper is to present numerical results of static bending and free vibration of functionally graded porous(FGP) variable-thickness plates by using an edge-based smoothed finite element method(ES-FEM) associate with the mixed interpolation of tensorial components technique for the three-node triangular element(MITC3), so-called ES-MITC3. This ES-MITC3 element is performed to eliminate the shear locking problem and to enhance the accuracy of the existing MITC3 element. In the ES-MITC3 element, the stiffness matrices are obtained by using the strain smoothing technique over the smoothing domains formed by two adjacent MITC3 triangular elements sharing an edge. Materials of the plate are FGP with a power-law index(k) and maximum porosity distributions(U) in the forms of cosine functions. The influences of some geometric parameters, material properties on static bending, and natural frequency of the FGP variable-thickness plates are examined in detail.
文摘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.
基金Changkye Lee and Jurng-Jae Yee would like to thank the support by Basic Science Research Program through the National Research Foundation(NRF)funded by Korea through Ministry of Education(No.2016R1A6A1A03012812).
文摘This work presents a locking-free smoothed finite element method(S-FEM)for the simulation of soft matter modelled by the equations of quasi-incompressible hyperelasticity.The proposed method overcomes well-known issues of standard finite element methods(FEM)in the incompressible limit:the over-estimation of stiffness and sensitivity to severely distorted meshes.The concepts of cell-based,edge-based and node-based S-FEMs are extended in this paper to three-dimensions.Additionally,a cubic bubble function is utilized to improve accuracy and stability.For the bubble function,an additional displacement degree of freedom is added at the centroid of the element.Several numerical studies are performed demonstrating the stability and validity of the proposed approach.The obtained results are compared with standard FEM and with analytical solutions to show the effectiveness of the method.
基金support by Basic Science Research Program through the National Research Foundation(NRF)funded by Korea Ministry of Education(No.2016R1A6A1A0312812).
文摘The aim of this work is to employ a modified cell-based smoothed finite element method(S-FEM)for topology optimization with the domain discretized with arbitrary polygons.In the present work,the linear polynomial basis function is used as the weight function instead of the constant weight function used in the standard S-FEM.This improves the accuracy and yields an optimal convergence rate.The gradients are smoothed over each smoothing domain,then used to compute the stiffness matrix.Within the proposed scheme,an optimum topology procedure is conducted over the smoothing domains.Structural materials are distributed over each smoothing domain and the filtering scheme relies on the smoothing domain.Numerical tests are carried out to pursue the performance of the proposed optimization by comparing convergence,efficiency and accuracy.
文摘A four-node quadrilateral shell element with smoothed membrane-bending based on Mindlin-Reissner theory is proposed. The element is a combination of a plate bending and membrane element. It is based on mixed interpolation where the bending and membrane stiffness matrices are calculated on the boundaries of the smoothing cells while the shear terms are approximated by independent interpolation functions in natural coordinates. The proposed element is robust, computationally inexpensive and free of locking. Since the integration is done on the element boundaries for the bending and membrane terms, the element is more accurate than the MITC4 element for distorted meshes. This will be demonstrated for several numerical examples.
基金supported by the National Natural Science Foundation of China(50474053,50475134 and 50675081)the 863 project (2007AA042142)
文摘Formulation and numerical evaluation of a novel twice-interpolation finite element method （TFEM） is presented for solid mechanics problems. In this method, the trial function for Galerkin weak form is constructed through two stages of consecutive interpolation. The primary interpolation follows exactly the same procedure of standard FEM and is further reproduced according to both nodal values and averaged nodal gradients obtained from primary interpolation. The trial functions thus constructed have continuous nodal gradients and contain higher order polynomial without increasing total freedoms. Several benchmark examples and a real dam problem are used to examine the TFEM in terms of accuracy and convergence. Compared with standard FEM, TFEM can achieve significantly better accuracy and higher convergence rate, and the continuous nodal stress can be obtained without any smoothing operation. It is also found that TFEM is insensitive to the quality of the elemental mesh. In addition, the present TFEM can treat the incompressible material without any modification.
基金supported by the National Basic Research Program of China("973"Project,Grant No.2010CB631005,2011CB606105)the National Natural Science Foundation of China(Grant No.10625209,10732080,90916010)+2 种基金China Postdoctoral Science Foundation(Grant No.20090460335)Specialized Research Fund for the Doctoral Program of Higher Education(Grant No.20090002110048)the opening funds from the State Key Laboratory of Explosion Science and Technology(KFJJ10-18Y)
文摘Optical full-field measurement methods are now widely applied in various domains.In general,the displacement fields can be directly obtained from the measurement,however in mechanical analysis strain fields are preferred.To extract strain fields from noisy displacement fields is always a challenging topic.In this study,a finite element method for smoothing displacement fields and calculating strain fields is proposed.An experimental test case on a holed aluminum specimen under tension is applied to validate this method.The heterogeneous displacement fields are measured by digital image correlation(DIC).By this proposed method,the result shows that the measuring noise on experimental displacement fields can be successfully removed,and strain fields can be reconstructed in the arbitrary area.
文摘The smoothed finite element method (S-FEM) was originated by G R Liu by combining some meshfree techniques with the well-established standard finite element method (FEM). It has a family of models carefully designed with innovative types of smoothing domains. These models are found having a number of important and theoretically profound properties. This article first provides a concise and easy-to-follow presentation of key formulations used in the S-FEM. A number of important properties and unique features of S-FEM models are discussed in detail, including 1) theoretically proven softening effects;2) upper-bound solutions;3) accurate solutions and higher convergence rates;4) insensitivity to mesh distortion;5) Jacobian?free;6) volumetric-locking-free;and most importantly 7) working well with triangular and tetrahedral meshes that can be automatically generated. The S-FEM is thus ideal for automation in computations and adaptive analyses, and hence has profound impact on Al-assisted modeling and simulation. Most importantly, one can now purposely design an S-FEM model to obtain solutions with special properties as wish, meaning that S-FEM offers a framework for design numerical models with desired properties. This novel concept of numerical model demand may drastically change the landscape of modeling and simulation. Future directions of research are also provided.
文摘The paper examines three selective schemes for the smoothed finite element method （SFEM） which was formulated by incorporating a cell-wise strain smoothing operation into the standard compatible finite element method （FEM）. These selective SFEM schemes were formulated based on three selective integration FEM schemes with similar properties found between the number of smoothing cells in the SFEM and the number of Gaussian integration points in the FEM. Both scheme 1 and scheme 2 are free of nearly incompressible locking, but scheme 2 is more general and gives better results than scheme 1. In addition, scheme 2 can be applied to anisotropic and nonlinear situations, while scheme 1 can only be applied to isotropic and linear situations. Scheme 3 is free of shear locking. This scheme can be applied to plate and shell problems. Results of the numerical study show that the selective SFEM schemes give more accurate results than the FEM schemes.
基金supported by the National Numerical Wind Tunnel Project (Grant No. NNW2019ZT2-B02)the National Natural Science Foundation of China (Grant Nos. 12032002,51779003,and 11902005)the SinoGerman Mobility Programme (Grant No. M-0210)。
文摘Free-surface flows, especially those associated with fluid-structure interactions(FSIs), pose challenging problems in numerical simulations. The authors of this work recently developed a smoothed particle element method(SPEM) to simulate FSIs. In this method, both the fluid and solid regions are initially modeled using a smoothed finite element method(S-FEM) in a Lagrangian frame, whereas the fluid regions undergoing large deformations are adaptively converted into particles and modeled with an improved smoothed particle hydrodynamics(SPH) method. This approach greatly improves computational accuracy and efficiency because of the advantages of the S-FEM in efficiently treating solid/fluid regions showing small deformations and the SPH method in effectively modeling moving interfaces. In this work, we further enhance the efficiency of the SPEM while effectively capturing local fluid information by introducing a multi-resolution technique to the SPEM and developing an effective approach to treat multi-resolution element-particle interfaces. Various numerical examples demonstrate that the multiresolution SPEM can significantly reduce the computational cost relative to the original version with a constant resolution.Moreover, the novel approach is effective in modeling various incompressible flow problems involving FSIs.
基金Vietnam National Foundation for Science and Technology Development(NAFOSTED)under Grant number 107.02-2019.330.
文摘An effective hybrid optimization method is proposed by integrating an adaptive Kriging(A-Kriging)into an improved partial swarm optimization algorithm(IPSO)to give a so-called A-Kriging-IPSO for maximizing the buckling load of laminated composite plates(LCPs)under uniaxial and biaxial compressions.In this method,a novel iterative adaptive Kriging model,which is structured using two training sample sets as active and adaptive points,is utilized to directly predict the buckling load of the LCPs and to improve the efficiency of the optimization process.The active points are selected from the initial data set while the adaptive points are generated using the radial random-based convex samples.The cell-based smoothed discrete shear gap method(CS-DSG3)is employed to analyze the buckling behavior of the LCPs to provide the response of adaptive and input data sets.The buckling load of the LCPs is maximized by utilizing the IPSO algorithm.To demonstrate the efficiency and accuracy of the proposed methodology,the LCPs with different layers(2,3,4,and 10 layers),boundary conditions,aspect ratios and load patterns(biaxial and uniaxial loads)are investigated.The results obtained by proposed method are in good agreement with the literature results,but with less computational burden.By applying adaptive radial Kriging model,the accurate optimal resultsebased predictions of the buckling load are obtained for the studied LCPs.
基金supported by the National Natural Science Foundation of China (Grant No. 10972188)the Program for New Century Excellent Talents in University from China Ministry of Education (Grant No. NCET-09-0678)the Fundamental Research Funds for the Central Universities of China (Grant No. 2010121073)
文摘A strain smoothing formulation for NURBS (non-uniform rational B-splines) based isogeometric finite element analysis is presented. This approach is formulated within the framework of assumed strain methods and strain smoothing operations. The strain smoothing is defined through strain averaging in the element sub-domains which are subsequently used for numerical integration of the Galerkin weak form. This formulation satisfies the orthogonality condition of the assumed strain methods. Meanwhile the present formulation totally avoids the gradient computation of the rational NURBS basis functions in the formulation of stiffness matrix. A transformation method is employed to accurately enforce the displacement boundary conditions. Numerical results demonstrate that the present formation gives very satisfactory solution accuracy simultaneously with improved computational efficiency.
文摘近年来,部分特高压换流变压器在运行中相继发生电弧短路故障并引发爆炸、起火事故,严重威胁直流系统安全稳定运行。充油设备故障冲击下的结构失效机制尚不明确,且缺乏成熟的数值计算方法,制约故障防爆技术的发展。基于此,该文提出一套适用于高能电弧故障冲击的结构失效仿真计算方法。首先,建立有限腔体内油中电弧能量持续注入的气泡动力学模型,准确描述故障气泡的脉动膨胀行为;其次,提出自适应有限元-光滑粒子流体动力学(finite element method-smoothed particle hydrodynamics,FEM-SPH)耦合方法,利用SPH粒子继承失效前的物理信息参与FEM计算;进行不同能量、不同位置的电弧故障仿真计算,获得换流变压器结构的薄弱区域及其破裂行为,复现了油箱结构失效行为。研究发现,油箱顶盖两侧及侧壁转角接缝位置容易发生应力集中现象。一旦结构出现裂缝,将在极短时间内沿着应力集中方向快速发展,最终导致壁面整体撕裂。可知,该方法揭示的结构失效行为可为改进变压器设计和提高设备安全性提供依据。
