The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element ...The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.展开更多
This study presents an experiment of improving the performance of spectral stochastic finite element method using high-order elements. This experiment is implemented through a two-dimensional spectral stochastic finit...This study presents an experiment of improving the performance of spectral stochastic finite element method using high-order elements. This experiment is implemented through a two-dimensional spectral stochastic finite element formulation of an elliptic partial differential equation having stochastic coefficients. Deriving this spectral stochastic finite element formulation couples a two-dimensional deterministic finite element formulation of an elliptic partial differential equation with generalized polynomial chaos expansions of stochastic coefficients. Further inspection of the performance of resulting spectral stochastic finite element formulation with adopting linear and quadratic (9-node or 8-node) quadrilateral elements finds that more accurate standard deviations of unknowns are surprisingly predicted using quadratic quadrilateral elements, especially under high autocorrelation function values of stochastic coefficients. In addition, creating spectral stochastic finite element results using quadratic quadrilateral elements is not unacceptably time-consuming. Therefore, this study concludes that adopting high-order elements can be a lower-cost method to improve the performance of spectral stochastic finite element method.展开更多
In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy o...In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.展开更多
In this paper, we established a finite element (FEM) model to analyze the dynamic characteristics of arch bridges. In this model, the effects of adjustment to the length of a suspender on its geometry stiffness matrix...In this paper, we established a finite element (FEM) model to analyze the dynamic characteristics of arch bridges. In this model, the effects of adjustment to the length of a suspender on its geometry stiffness matrix are stressed. The FEM equations of mechanics characteristics, natural frequency and main mode are set up based on the first order matrix perturbation theory. Applicantion of the proposed model to analyze a real arch bridge proved the improvement in the simulation precision of dynamical characteristics of the arch bridge by considering the effects of suspender length variation.展开更多
The elastoplastic field near crack tips is investigated through finite element simulation.A refined mesh model near the crack tip is proposed. In the mesh refining area, element size continuously varies from the nanom...The elastoplastic field near crack tips is investigated through finite element simulation.A refined mesh model near the crack tip is proposed. In the mesh refining area, element size continuously varies from the nanometer scale to themicrometer scale and the millimeter scale. Graphics of the plastic zone, the crack tip blunting, and the deformed crack tip elements are given in the paper.Based on the curves of stress and plastic strain, closely near the crack tip, the stresssingularity index and the stress intensity factor,as well as the plastic strain singularity index and the plastic strain intensity factor are determined.Thestress and plastic strainsingular index vary with the load, while the dimensions of the stress and the plastic strain intensity factorsdependon the stress and the plastic strain singularity index, respectively. The singular field near the elastoplastic crack tip is characterized by the stress singularity index and the stress intensity factor, or alternativelythe plastic strain singularity index and the plastic strain intensityfactor.At the end of the paper, following Irwin’s concept of fracture mechanics,σδKσδKcriterion andεδQεδQcriterion are proposed.Besides, crack tip angle criterion is also presented.展开更多
Effects of welding current on temperature and velocity fields during gas metal arc welding(GMAW) of commercially pure aluminum were simulated. Equations of conservation of mass, energy and momentum were solved in a th...Effects of welding current on temperature and velocity fields during gas metal arc welding(GMAW) of commercially pure aluminum were simulated. Equations of conservation of mass, energy and momentum were solved in a three-dimensional transient model using FLOW-3 D software. The mathematical model considered buoyancy and surface tension driving forces. Further, effects of droplet heat content and impact force on weld pool surface deformation were added to the model. The results of simulation showed that an increase in the welding current could increase peak temperature and the maximum velocity in the weld pool. The weld pool dimensions and width of the heat-affected zone(HAZ) were enlarged by increasing the welding current. In addition, dimensionless Peclet, Grashof and surface tension Reynolds numbers were calculated to understand the importance of heat transfer by convection and the roles of various driving forces in the weld pool. In order to validate the model, welding experiments were conducted under several welding currents. The predicted weld pool dimensions were compared with the corresponding experimental results, and good agreement between simulation and preliminary test results was achieved.展开更多
Laser welding (LW) becomes one of the most economical high quality joining processes. LW offers the advantage of very controlled heat input resulting in low distortion and the ability to weld heat sensitive components...Laser welding (LW) becomes one of the most economical high quality joining processes. LW offers the advantage of very controlled heat input resulting in low distortion and the ability to weld heat sensitive components. To exploit efficiently the benefits presented by LW, it is necessary to develop an integrated approach to identify and control the welding process variables in order to produce the desired weld characteristics without being forced to use the traditional and fastidious trial and error procedures. The paper presents a study of weld bead geometry characteristics prediction for laser overlap welding of low carbon galvanized steel using 3D numerical modelling and experimental validation. The temperature dependent material properties, metallurgical transformations and enthalpy method constitute the foundation of the proposed modelling approach. An adaptive 3D heat source is adopted to simulate both keyhole and conduction mode of the LW process. The simulations are performed using 3D finite element model on commercial software. The model is used to estimate the weld bead geometry characteristics for various LW parameters, such as laser power, welding speed and laser beam diameter. The calibration and validation of the 3D numerical model are based on experimental data achieved using a 3 kW Nd:Yag laser system, a structured experimental design and confirmed statistical analysis tools. The results reveal that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects on the weld quality. The results show great concordance between predicted and measured values for weld bead geometry characteristics, such as depth of penetration, bead width at the top surface and bead width at the interface between sheets, with an average accuracy greater than 95%.展开更多
A theory invoking concepts from differential geometry of generalized Finsler space in conjunction with diffuse interface modeling is described and implemented in finite element(FE)simulations of dual-phase polycrystal...A theory invoking concepts from differential geometry of generalized Finsler space in conjunction with diffuse interface modeling is described and implemented in finite element(FE)simulations of dual-phase polycrystalline ceramic microstructures.Order parameters accounting for fracture and other structural transformations,notably partial dislocation slip,twinning,or phase changes,are dimensionless entries of an internal state vector of generalized pseudo-Finsler space.Ceramics investigated in computations are a boron carbide-titanium diboride(B4C-TiB2)composite and a diamond-silicon carbide(C-SiC)composite.Deformation mechanisms-in addition to elasticity and cleavage fracture in grains of any phase-include restricted dislocation glide(TiB2 phase),deformation twinning(B4C and-SiC phases),and stress-induced amorphization(B4C phase).The metric tensor of generalized Finsler space is scaled conformally according to dilatation induced by cavitation or other fracture modes and densification induced by phase changes.Simulations of pure shear consider various morphologies and lattice orientations.Effects of microstructure on overall strength of each composite are reported.In B4C-TiB2,minor improvements in shear strength and ductility are observed with an increase in the second phase from 10 to 18%by volume,suggesting that residual stresses or larger-scale crack inhibition may be responsible for toughness gains reported experimentally.In diamond-SiC,a composite consisting of diamond crystals encapsulated in a nano-crystalline SiC matrix shows improved strength and ductility relative to a two-phase composite with isolated bulk SiC grains.展开更多
We suggest an original approach to Lobachevski’s geometry and Hilbert’s Fourth Problem, based on the use of the “mathematics of harmony” and special class of hyperbolic functions, the so-called hyperbolic Fibonacc...We suggest an original approach to Lobachevski’s geometry and Hilbert’s Fourth Problem, based on the use of the “mathematics of harmony” and special class of hyperbolic functions, the so-called hyperbolic Fibonacci l-functions, which are based on the ancient “golden proportion” and its generalization, Spinadel’s “metallic proportions.” The uniqueness of these functions consists in the fact that they are inseparably connected with the Fibonacci numbers and their generalization― Fibonacci l-numbers (l > 0 is a given real number) and have recursive properties. Each of these new classes of hyperbolic functions, the number of which is theoretically infinite, generates Lobachevski’s new geometries, which are close to Lobachevski’s classical geometry and have new geometric and recursive properties. The “golden” hyperbolic geometry with the base (“Bodnar’s geometry) underlies the botanic phenomenon of phyllotaxis. The “silver” hyperbolic geometry with the base ?has the least distance to Lobachevski’s classical geometry. Lobachevski’s new geometries, which are an original solution of Hilbert’s Fourth Problem, are new hyperbolic geometries for physical world.展开更多
This paper presents the solution of coupled radiative transfer equation with heat conduction equation in complex three-dimensional geometries. Due to very different time scales for both physics, the radiative problem ...This paper presents the solution of coupled radiative transfer equation with heat conduction equation in complex three-dimensional geometries. Due to very different time scales for both physics, the radiative problem is considered steady-state but solved at each time iteration of the transient conduction problem. The discrete ordinate method along with the decentered streamline-upwind Petrov-Galerkin method is developed. Since specular reflection is considered on borders, a very accurate algorithm has been developed for calculation of partition ratio coefficients of incident solid angles to the several reflected solid angles. The developed algorithms are tested on a paraboloid-shaped geometry used for example on concentrated solar power technologies.展开更多
基金Project supported by the National Natural Science Foundation of China (Nos. 12102043, 12072375U2241240)the Natural Science Foundation of Hunan Province of China (Nos. 2023JJ40698 and 2021JJ40710)。
文摘The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.
文摘This study presents an experiment of improving the performance of spectral stochastic finite element method using high-order elements. This experiment is implemented through a two-dimensional spectral stochastic finite element formulation of an elliptic partial differential equation having stochastic coefficients. Deriving this spectral stochastic finite element formulation couples a two-dimensional deterministic finite element formulation of an elliptic partial differential equation with generalized polynomial chaos expansions of stochastic coefficients. Further inspection of the performance of resulting spectral stochastic finite element formulation with adopting linear and quadratic (9-node or 8-node) quadrilateral elements finds that more accurate standard deviations of unknowns are surprisingly predicted using quadratic quadrilateral elements, especially under high autocorrelation function values of stochastic coefficients. In addition, creating spectral stochastic finite element results using quadratic quadrilateral elements is not unacceptably time-consuming. Therefore, this study concludes that adopting high-order elements can be a lower-cost method to improve the performance of spectral stochastic finite element method.
文摘In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.
基金Supported by the Key Teacher Foundation of Chongqing University (No. 717411067)
文摘In this paper, we established a finite element (FEM) model to analyze the dynamic characteristics of arch bridges. In this model, the effects of adjustment to the length of a suspender on its geometry stiffness matrix are stressed. The FEM equations of mechanics characteristics, natural frequency and main mode are set up based on the first order matrix perturbation theory. Applicantion of the proposed model to analyze a real arch bridge proved the improvement in the simulation precision of dynamical characteristics of the arch bridge by considering the effects of suspender length variation.
基金The work was supported by the National Natural Science Foundation of China (Grant 11572226).
文摘The elastoplastic field near crack tips is investigated through finite element simulation.A refined mesh model near the crack tip is proposed. In the mesh refining area, element size continuously varies from the nanometer scale to themicrometer scale and the millimeter scale. Graphics of the plastic zone, the crack tip blunting, and the deformed crack tip elements are given in the paper.Based on the curves of stress and plastic strain, closely near the crack tip, the stresssingularity index and the stress intensity factor,as well as the plastic strain singularity index and the plastic strain intensity factor are determined.Thestress and plastic strainsingular index vary with the load, while the dimensions of the stress and the plastic strain intensity factorsdependon the stress and the plastic strain singularity index, respectively. The singular field near the elastoplastic crack tip is characterized by the stress singularity index and the stress intensity factor, or alternativelythe plastic strain singularity index and the plastic strain intensityfactor.At the end of the paper, following Irwin’s concept of fracture mechanics,σδKσδKcriterion andεδQεδQcriterion are proposed.Besides, crack tip angle criterion is also presented.
文摘Effects of welding current on temperature and velocity fields during gas metal arc welding(GMAW) of commercially pure aluminum were simulated. Equations of conservation of mass, energy and momentum were solved in a three-dimensional transient model using FLOW-3 D software. The mathematical model considered buoyancy and surface tension driving forces. Further, effects of droplet heat content and impact force on weld pool surface deformation were added to the model. The results of simulation showed that an increase in the welding current could increase peak temperature and the maximum velocity in the weld pool. The weld pool dimensions and width of the heat-affected zone(HAZ) were enlarged by increasing the welding current. In addition, dimensionless Peclet, Grashof and surface tension Reynolds numbers were calculated to understand the importance of heat transfer by convection and the roles of various driving forces in the weld pool. In order to validate the model, welding experiments were conducted under several welding currents. The predicted weld pool dimensions were compared with the corresponding experimental results, and good agreement between simulation and preliminary test results was achieved.
文摘Laser welding (LW) becomes one of the most economical high quality joining processes. LW offers the advantage of very controlled heat input resulting in low distortion and the ability to weld heat sensitive components. To exploit efficiently the benefits presented by LW, it is necessary to develop an integrated approach to identify and control the welding process variables in order to produce the desired weld characteristics without being forced to use the traditional and fastidious trial and error procedures. The paper presents a study of weld bead geometry characteristics prediction for laser overlap welding of low carbon galvanized steel using 3D numerical modelling and experimental validation. The temperature dependent material properties, metallurgical transformations and enthalpy method constitute the foundation of the proposed modelling approach. An adaptive 3D heat source is adopted to simulate both keyhole and conduction mode of the LW process. The simulations are performed using 3D finite element model on commercial software. The model is used to estimate the weld bead geometry characteristics for various LW parameters, such as laser power, welding speed and laser beam diameter. The calibration and validation of the 3D numerical model are based on experimental data achieved using a 3 kW Nd:Yag laser system, a structured experimental design and confirmed statistical analysis tools. The results reveal that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects on the weld quality. The results show great concordance between predicted and measured values for weld bead geometry characteristics, such as depth of penetration, bead width at the top surface and bead width at the interface between sheets, with an average accuracy greater than 95%.
文摘A theory invoking concepts from differential geometry of generalized Finsler space in conjunction with diffuse interface modeling is described and implemented in finite element(FE)simulations of dual-phase polycrystalline ceramic microstructures.Order parameters accounting for fracture and other structural transformations,notably partial dislocation slip,twinning,or phase changes,are dimensionless entries of an internal state vector of generalized pseudo-Finsler space.Ceramics investigated in computations are a boron carbide-titanium diboride(B4C-TiB2)composite and a diamond-silicon carbide(C-SiC)composite.Deformation mechanisms-in addition to elasticity and cleavage fracture in grains of any phase-include restricted dislocation glide(TiB2 phase),deformation twinning(B4C and-SiC phases),and stress-induced amorphization(B4C phase).The metric tensor of generalized Finsler space is scaled conformally according to dilatation induced by cavitation or other fracture modes and densification induced by phase changes.Simulations of pure shear consider various morphologies and lattice orientations.Effects of microstructure on overall strength of each composite are reported.In B4C-TiB2,minor improvements in shear strength and ductility are observed with an increase in the second phase from 10 to 18%by volume,suggesting that residual stresses or larger-scale crack inhibition may be responsible for toughness gains reported experimentally.In diamond-SiC,a composite consisting of diamond crystals encapsulated in a nano-crystalline SiC matrix shows improved strength and ductility relative to a two-phase composite with isolated bulk SiC grains.
文摘We suggest an original approach to Lobachevski’s geometry and Hilbert’s Fourth Problem, based on the use of the “mathematics of harmony” and special class of hyperbolic functions, the so-called hyperbolic Fibonacci l-functions, which are based on the ancient “golden proportion” and its generalization, Spinadel’s “metallic proportions.” The uniqueness of these functions consists in the fact that they are inseparably connected with the Fibonacci numbers and their generalization― Fibonacci l-numbers (l > 0 is a given real number) and have recursive properties. Each of these new classes of hyperbolic functions, the number of which is theoretically infinite, generates Lobachevski’s new geometries, which are close to Lobachevski’s classical geometry and have new geometric and recursive properties. The “golden” hyperbolic geometry with the base (“Bodnar’s geometry) underlies the botanic phenomenon of phyllotaxis. The “silver” hyperbolic geometry with the base ?has the least distance to Lobachevski’s classical geometry. Lobachevski’s new geometries, which are an original solution of Hilbert’s Fourth Problem, are new hyperbolic geometries for physical world.
文摘This paper presents the solution of coupled radiative transfer equation with heat conduction equation in complex three-dimensional geometries. Due to very different time scales for both physics, the radiative problem is considered steady-state but solved at each time iteration of the transient conduction problem. The discrete ordinate method along with the decentered streamline-upwind Petrov-Galerkin method is developed. Since specular reflection is considered on borders, a very accurate algorithm has been developed for calculation of partition ratio coefficients of incident solid angles to the several reflected solid angles. The developed algorithms are tested on a paraboloid-shaped geometry used for example on concentrated solar power technologies.