Search results

This study aims to introduce the hybrid finite element (FE) – meshfree method and multiscale variational principle into the traditional mixed FE formulation, leading to a stable mixed formulation for incompressible linear elasticity which circumvents the need to satisfy inf-sup condition.

Using the hybrid FE–meshfree method, the displacement and pressure are interpolated conveniently with the same order so that a continuous pressure field can be obtained with low-order elements. The multiscale variational principle is then introduced into the Galerkin form to obtain stable and convergent results.

The present method is capable of overcoming volume locking and does not exhibit unphysical oscillations near the incompressible limit. Moreover, there are no extra unknowns introduced in the present method because the fine-scale unknowns are eliminated using the static condensation technique, and there is no need to evaluate any user-defined stability parameter as the classical stabilization methods do. The shape functions constructed in the present model possess continuous derivatives at nodes, which gives a continuous and more precise stress field with no need of an additional smooth process. The shape functions in the present model also possess the Kronecker delta property, so that it is convenient to impose essential boundary conditions.

The proposed model can be implemented easily. Its convergence rates and accuracy in displacement, energy and pressure are even comparable to those of second-order mixed elements.

The purpose of this paper is to propose a new three-node triangular element in the framework of the numerical manifold method (NMM), which is designated by Trig3-MLScns.

The formulation uses the improved parametric shape functions of classical triangular elements (Trig3-0) to construct the partition of unity (PU) and the moving least square (MLS) interpolation method to construct the local approximation function.

Compared with the classical three-node element (Trig3-0), the Trig3-MLScns element has a higher order of approximations, much better accuracy and continuous nodal stress. Moreover, the linear dependence problem associated with many PU-based methods with high-order approximations is eliminated in the present element. A number of numerical examples indicate the high accuracy and robustness of the Trig3-MLScns element.

The proposed element inherits the individual merits of the NMM and the MLS.

The goal of structural optimization is to find the best possible configuration that minimizes the objective function and satisfies a set of constraints. Here we present a method based on the evolutionary structural optimization method, where the quality of the solution is improved by avoiding the chain‐like sets of elements which are sources of potential kinematic instabilities, and by including local error estimators. Both of these enhancements are employed to activate refining the mesh so as to obtain accurate and stable solutions as the volume removal proceeds. Several related contributions of Professor E. Hinton are cited.

This study presents finite element (FE) and generalized regression neural network (GRNN) approaches for modeling multiple crack growth problems and predicting crack-growth directions under the influence of multiple crack parameters.

To determine the crack-growth direction in aluminum specimens, multiple crack parameters representing some degree of crack propagation complexity, including crack length, inclination angle, offset and distance, were examined. FE method models were developed for multiple crack growth simulations. To capture the complex relationships among multiple crack-growth variables, GRNN models were developed as nonlinear regression models. Six input variables and one output variable comprising 65 training and 20 test datasets were established.

The FE model could conveniently simulate the crack-growth directions. However, several multiple crack parameters could affect the simulation accuracy. The GRNN offers a reliable method for modeling the growth of multiple cracks. Using 76% of the total dataset, the NN model attained an R2 value of 0.985.

The models are presented for static multiple crack growth problems. No material anisotropy is observed.

In practical crack-growth analyses, the NN approach provides significant benefits and savings.

The proposed GRNN model is simple to develop and accurate. Its performance was superior to that of other NN models. This model is also suitable for modeling multiple crack growths with arbitrary geometries. The proposed GRNN model demonstrates its prediction capability with a simpler learning process, thus producing efficient multiple crack growth predictions and assessments.

In finite element analysis of pile driving, the nodes of the finite element mesh are the most important locations for output stresses. Especially at the pile‐soil interface, it is essential to obtain accurate nodal stresses. Several global and local stress smoothing methods available in the literature were reviewed and examined. Global methods are found to be computationally expensive, so results obtained from several local stress smoothing methods are compared. It is shown that accurate nodal stresses can be obtained by approximating the stress distribution inside four‐element patches by a polynomial with order equal to the order of the shape functions. Equally good results can be obtained by approximating the stress distribution inside each element by a bilinear surface. When a method taking into account both equilibrium and boundary conditions was applied, a set of ill‐conditioned matrices was produced for the four‐element patches. Such methods are therefore not recommended.

HAVING discussed in the standard longhand notation the main ideas and methods for the calculation of redundant structures on the basis of forces as unknowns we now turn our attention to the matrix formulation of the analysis. Consider a system consisting of s structural elements with a total number n of redundancies which may be forces (stresses), moments or any generalized forces. We select a basic system by ‘cutting’ a number r of redundancies where r<n. Thus, the simple idea of a statically determinate basic system (r=n) is but a particular case of our investigations.

A finite element model for numerical simulation of non‐steady but continuous chip formation under orthogonal cutting conditions is described. The problem is treated as coupled thermo‐mechanical. A velocity approach has been adopted for the proposed solution. The computational algorithm takes care of dynamic contact conditions and makes use of an automatic remeshing procedure. The results of simulation yield complete history of chip initiation and growth as well as distributions of strain rate, strain, stress and temperature. The paper includes a detailed presentation of computational results for an illustrative case.

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 address error‐controlled adaptive finite element (FE) method for thin and thick plates. A procedure is presented for determining the most suitable plate model (among available hierarchical plate models) for each particular FE of the selected mesh, that is provided as the final output of the mesh adaptivity procedure.

The model adaptivity procedure can be seen as an appropriate extension to model adaptivity for linear elastic plates of so‐called equilibrated boundary traction approach error estimates, previously proposed for 2D/3D linear elasticity. Model error indicator is based on a posteriori element‐wise computation of improved (continuous) equilibrated boundary stress resultants, and on a set of hierarchical plate models. The paper illustrates the details of proposed model adaptivity procedure for choosing between two most frequently used plate models: the one of Kirchhoff and the other of Reissner‐Mindlin. The implementation details are provided for a particular case of the discrete Kirchhoff quadrilateral four‐node plate FE and the corresponding Reissner‐Mindlin quadrilateral with the same number of nodes. The key feature for those elements that they both provide the same quality of the discretization space (and thus the same discretization error) is the one which justifies uncoupling of the proposed model adaptivity from the mesh adaptivity.

Several numerical examples are presented in order to illustrate a very satisfying performance of the proposed methodology in guiding the final choice of the optimal model and mesh in analysis of complex plate structures.

The paper confirms that one can make an automatic selection of the most appropriate plate model for thin and thick plates on the basis of proposed model adaptivity procedure.

This article aims to develop an accurate and efficient meshfree Galerkin method based on the strain smoothing technique for linear elastic continuous and fracture problems.

This paper proposed a generalized linear smoothed meshfree method (LSMM), in which the compatible strain is reconstructed by the linear smoothed strains. Based on the idea of the weighted residual method and employing three linearly independent weight functions, the linear smoothed strains can be created easily in a smoothing domain. Using various types of basic functions, LSMM can solve the linear elastic continuous and fracture problems in a unified way.

On the one hand, the LSMM inherits the properties of high efficiency and stability from the stabilized conforming nodal integration (SCNI). On the other hand, the LSMM is more accurate than the SCNI, because it can produce continuous strains instead of the piece-wise strains obtained by SCNI. Those excellent performances ensure that the LSMM has the capability to precisely track the crack propagation problems. Several numerical examples are investigated to verify the accurate, convergence rate and robustness of the present LSMM.

This study provides an accurate and efficient meshfree method for simulating crack growth.