Search results

The purpose of this paper is to enhance the feasibility of the edge-based smoothed triangular (EST) element, some modifications are made in this study.

First, an efficient strategy based on graph theory is proposed to construct the edge system. Second, the stress is smoothed in global coordinate system based on edge instead of strain, which makes the theory of EST more rigorous and can be easily extended to the situation of multi elements sharing the same edge. Third, the singular degree of freedoms (DOFs) of the nodes linked by edges are restrained in edge local coordinate system, which makes the global stiffness matrix non-singular and can be decomposed successfully.

First, an efficient edge constructing strategy can make EST element more standout. Second, some modifications should be made to EST element to extend it to the situation with multi elements sharing the same edge, so that EST element can deal with the geometrical models with this kind of features. Third, the way to restrain the singular DOFs of EST element must be different from normal isoparametric triangle element, because the stiffness matrix of the smoothing domain is not computed in local coordinate system.

The modified EST element performs stably in engineering analysis including large scale problems and the situation with multi elements sharing the same edge, and the efficiency of edge system constructing is no longer the bottleneck.

The purpose of this paper is to establish a peridynamic method in predicting viscoelastic creep behaviour with recovery stage and to find the suitable numerical parameters of peridynamic method.

A rheological viscoelastic creep constitutive equation including recovery and an elastic peridynamic equation (with integral basis) are examined and used. The elasticity equation within the peridynamic equation is replaced by the viscoelastic equation. A new peridynamic method with two time parameters, i.e. numerical time and viscoelastic real time is designed. The two parameters of peridynamic method, horizon radius and number of nodes per unit volume are studied to get their optimal values. In validating this peridynamic method, comparisons are made between numerical and analytical result and between numerical and experimental data.

The new peridynamic method for viscoelastic creep behaviour is approved by the good matching in numerical-analytical data comparison with difference of < 0.1 per cent and in numerical-experimental data comparison with difference of 4-6 per cent. It can be used for further creep test which may include non-linear viscoelastic behaviour and creep rupture. From this paper, the variation of constants in Burger’s viscoelastic model is also studied and groups of constants values that can simulate solid, fluid and solid-fluid viscoelastic behaviours were obtained. In addition, the numerical peridynamic parameters were also manipulated and examined to achieve the optimal values of the parameters.

The peridynamic model of viscoelastic creep behaviour preferably should have only one time parameter. This can only be done by solving the unstable fluctuation of dynamic results, which is not discussed in this paper. Another limitation is the tertiary region and creep rupture are not included in this paper.

The viscoelastic peridynamic model in this paper can serve as an alternative for conventional numerical simulations in viscoelastic area. This model also is the initial step of developing peridynamic model of viscoelastic creep rupture properties (crack initiation, crack propagation, crack branching, etc.), where this future model has high potential in predicting failure behaviours of any components, tools or structures, and hence increase safety and reduce loss.

The application of viscoelastic creep constitutive model on peridynamic formulation, effect of peridynamic parameters manipulation on numerical result, and optimization of constants of viscoelastic model in simulating three types of viscoelastic creep behaviours.

Pores and shrink holes are unavoidable defects in the die-casting mass production process which may significantly influence the strength, fatigue and fracture behaviour as well as the life span of structures, especially if they are subjected to high static and dynamic loads. Such defects should be considered during the design process or after production, where the defects could be detected with the help of computed tomography (CT) measurements. However, this is usually not done in today's mass production environments. This paper deals with the stress analysis of die-cast structural parts with pores found from CT measurements or that are artificially placed within a structure.

In this paper the authors illustrate two general methodologies to take into account the porosity of die-cast components in the stress analysis. The detailed geometry of a die-cast part including all discontinuities such as pores and shrink holes can be included via STL data provided by CT measurements. The first approach is a combination of the finite element method (FEM) and the finite cell method (FCM), which extends the FEM if the real geometry cuts finite elements. The FCM is only applied in regions with pores. This procedure has the advantage that all simulations with different pore distributions, real or artificial, can be calculated without changing the base finite element mesh. The second approach includes the pore information as STL data into the original CAD model and creates a new adapted finite element mesh for the simulation. Both methods are compared and evaluated for an industrial problem.

The STL data of defects which the authors received from CT measurements could not be directly applied without repairing them. Therefore, for FEM applications an appropriate repair procedure is proposed. The first approach, which combines the FEM with the FCM, the authors have realized within the commercial software tool Abaqus. This combination performs well, which is demonstrated for test examples, and is also applied for a complex industrial project. The developed in-house code still has some limitations which restrict broader application in industry. The second pure FEM-based approach works well without limitations but requires increasing computational effort if many different pore distributions are to be investigated.

A new simulation approach which combines the FEM with the FCM has been developed and implemented into the commercial Abaqus FEM software. This approach the authors have applied to simulate a real engineering die-cast structure with pores. This approach could become a preferred way to consider pores in practical applications, where the porosity can be derived either from CT measurements or are artificially adopted for design purposes. The authors have also shown how pores can be considered in the standard FEM analysis as well.

Failure of the materials occurs once the stress intensity factor (SIF) overtakes the material fracture toughness. At this level, the crack will grow rapidly resulting in unstable crack growth until a complete fracture happens. The SIF calculation of the materials can be conducted by experimental, theoretical and numerical techniques. Prediction of SIF is crucial to ensure safety life from the material failure. The aim of the simulation study is to evaluate the accuracy of SIF prediction using finite element analysis.

The bootstrap resampling method is employed in S-version finite element model (S-FEM) to generate the random variables in this simulation analysis. The SIF analysis studies are promoted by bootstrap S-version Finite Element Model (BootstrapS-FEM). Virtual crack closure-integral method (VCCM) is an important concept to compute the energy release rate and SIF. The semielliptical crack shape is applied with different crack shape aspect ratio in this simulation analysis. The BootstrapS-FEM produces the prediction of SIFs for tension model.

The mean of BootstrapS-FEM is calculated from 100 samples by the resampling method. The bounds are computed based on the lower and upper bounds of the hundred samples of BootstrapS-FEM. The prediction of SIFs is validated with Newman–Raju solution and deterministic S-FEM within 95 percent confidence bounds. All possible values of SIF estimation by BootstrapS-FEM are plotted in a graph. The mean of the BootstrapS-FEM is referred to as point estimation. The Newman–Raju solution and deterministic S-FEM values are within the 95 percent confidence bounds. Thus, the BootstrapS-FEM is considered valid for the prediction with less than 6 percent of percentage error.

The bootstrap resampling method is employed in S-FEM to generate the random variables in this simulation analysis.

It has been usual to prefer an enrichment pattern independent of the mesh when applying singular functions in the Generalized/eXtended finite element method (G/XFEM). This choice, when modeling crack tip singularities through extrinsic enrichment, has been understood as the only way to surpass the typical poor convergence rate obtained with the finite element method (FEM), on uniform or quasi-uniform meshes conforming to the crack. Then, the purpose of this study is to revisit the topological enrichment strategy in the light of a higher-order continuity obtained with a smooth partition of unity (PoU). Aiming to verify the smoothness' impacts on the blending phenomenon, a series of numerical experiments is conceived to compare the two GFEM versions: the conventional one, based on piecewise continuous PoU's, and another which considers PoU's with high-regularity.

The stress approximations right at the crack tip vicinity are qualified by focusing on crack severity parameters. For this purpose, the material forces method originated from the configurational mechanics is employed. Some attempts to improve solution using different polynomial enrichment schemes, besides the singular one, are discussed aiming to verify the transition/blending effects. A classical two-dimensional problem of the linear elastic fracture mechanics (LEFM) is solved, considering the pure mode I and the mixed-mode loadings.

The results reveal that, in the presence of smooth PoU's, the topological enrichment can still be considered as a suitable strategy for extrinsic enrichment. First, because such an enrichment pattern still can treat the crack independently of the mesh and deliver some advantage in terms of convergence rates, under certain conditions, when compared to the conventional FEM. Second, because the topological pattern demands fewer degrees of freedom and impacts conditioning less than the geometrical strategy.

Several outputs are presented, considering estimations for the J–integral and the angle of probable crack advance, this last computed from two different strategies to monitoring blending/transition effects, besides some comments about conditioning. Both h- and p-behaviors are displayed to allow a discussion from different points of view concerning the topological enrichment in smooth GFEM.

The purpose of this paper is to give a review on the newest developments of high-performance finite element methods (FEMs), and exhibit the recent contributions achieved by the authors’ group, especially showing some breakthroughs against inherent difficulties existing in the traditional FEM for a long time.

Three kinds of new FEMs are emphasized and introduced, including the hybrid stress-function element method, the hybrid displacement-function element method for Mindlin–Reissner plate and the improved unsymmetric FEM. The distinguished feature of these three methods is that they all apply the fundamental analytical solutions of elasticity expressed in different coordinates as their trial functions.

The new FEMs show advantages from both analytical and numerical approaches. All the models exhibit outstanding capacity for resisting various severe mesh distortions, and even perform well when other models cannot work. Some difficulties in the history of FEM are also broken through, such as the limitations defined by MacNeal’s theorem and the edge-effect problems of Mindlin–Reissner plate.

These contributions possess high value for solving the difficulties in engineering computations, and promote the progress of FEM.

This paper aims to present a computational framework to generate numeric enrichment functions for two-dimensional problems dealing with single/multiple local phenomenon/phenomena. The two-scale generalized/extended finite element method (G/XFEM) approach used here is based on the solution decomposition, having global- and local-scale components. This strategy allows the use of a coarse mesh even when the problem produces complex local phenomena. For this purpose, local problems can be defined where these local phenomena are observed and are solved separately by using fine meshes. The results of the local problems are used to enrich the global one improving the approximate solution.

The implementation of the two-scale G/XFEM formulation follows the object-oriented approach presented by the authors in a previous work, where it is possible to combine different kinds of elements and analyses models with the partition of unity enrichment scheme. Beside the extension of the G/XFEM implementation to enclose the global–local strategy, the imposition of different boundary conditions is also generalized.

The generalization done for boundary conditions is very important, as the global–local approach relies on the boundary information transferring process between the two scales of the analysis. The flexibility for the numerical analysis of the proposed framework is illustrated by several examples. Different analysis models, element formulations and enrichment functions are used, and the accuracy, robustness and computational efficiency are demonstrated.

This work shows a generalize imposition of different boundary conditions for global–local G/XFEM analysis through an object-oriented implementation. This generalization is very important, as the global–local approach relies on the boundary information transferring process between the two scales of the analysis. Also, solving multiple local problems simultaneously and solving plate problems using global–local G/XFEM are other contributions of this work.

This paper aims to improve the finite element modelling of transformers subjected to DC excitation, by including core joint details.

Geomagnetically induced currents (GICs) or leakage DC can cause part-cycle, half wave saturation of a power transformer’s core. Practical measurements and finite element matrix (FEM) simulation were carried out using three laboratory-scale, untanked single-phase four limb transformers resembling real power transformers in terms of the core steel and parallel winding assemblies. “Equivalent air gaps” at the joints, based on AC measurements, were applied to the FEM models for simultaneous AC and DC excitation.

Measurements confirm that introducing equivalent air gaps at the joints improves the FEM simulation of transformers carrying DC.

The FEM simulations based on the laboratory transformers are exemplary, showing the difference between modelling core joints as solid or including equivalent air gaps. They show that, for more representative results, laboratory transformers used for research should have mitred core joints (like power transformers).

This research shows why joint details are important in FEM models for analysing transformer core saturation in the presence of DC/GICs. Extending this, other core structures of power transformers with mitred joints should improve the understanding of the leakage flux during half-wave saturation.

Of particular interest is the ability of the extended finite element method (XFEM) to capture transient solution and motion of phase boundaries without adaptive remeshing or moving-mesh algorithms for a physically nonlinear phase change problem. The paper aims to discuss this issue.

The XFEM is applied to solve nonlinear transient problems with a phase change. Thermal conductivity and volumetric heat capacity are assumed to be dependent on temperature. The nonlinearities in the governing equations make it necessary to employ an effective iterative approach to solve the problem. The Newton-Raphson method is used and the incremental discrete XFEM equations are derived.

The robustness and utility of the method are demonstrated on several one-dimensional benchmark problems.

The novel procedure based on the XFEM is developed to solve physically nonlinear phase change problems.

The purpose of this paper is to evaluate some numerical integration strategies used in generalized (G)/extended finite element method (XFEM) to solve linear elastic fracture mechanics problems. A range of parameters are here analyzed, evidencing how the numerical integration error and the computational efficiency are improved when particularities from these examples are properly considered.

Numerical integration strategies were implemented in an existing computational environment that provides a finite element method and G/XFEM tools. The main parameters of the analysis are considered and the performance using such strategies is compared with standard integration results.

Known numerical integration strategies suitable for fracture mechanics analysis are studied and implemented. Results from different crack configurations are presented and discussed, highlighting the necessity of alternative integration techniques for problems with singularities and/or discontinuities.

This study presents a variety of fracture mechanics examples solved by G/XFEM in which the use of standard numerical integration with Gauss quadratures results in loss of precision. It is discussed the behaviour of subdivision of elements and mapping of integration points strategies for a range of meshes and cracks geometries, also featuring distorted elements and how they affect strain energy and stress intensity factors evaluation for both strategies.