Search results

This study implements the fourth-order phase field method (PFM) for modeling fracture in brittle materials. The weak form of the fourth-order PFM requires C¹ basis functions for the crack evolution scalar field in a finite element framework. To address this, non-Sibsonian type shape functions that are nonpolynomial types based on distance measures, are used in the context of natural neighbor shape functions. The capability and efficiency of this method are studied for modeling cracks.

The weak form of the fourth-order PFM is derived from two governing equations for finite element modeling. C⁰ non-Sibsonian shape functions are derived using distance measures on a generalized quad element. Then these shape functions are degree elevated with Bernstein-Bezier (BB) patch to get higher-order continuity (C¹) in the shape function. The quad element is divided into several background triangular elements to apply the Gauss-quadrature rule for numerical integration. Both fourth-order and second-order PFMs are implemented in a finite element framework. The efficiency of the interpolation function is studied in terms of convergence and accuracy for capturing crack topology in the fourth-order PFM.

It is observed that fourth-order PFM has higher accuracy and convergence than second-order PFM using non-Sibsonian type interpolants. The former predicts higher failure loads and failure displacements compared to the second-order model due to the addition of higher-order terms in the energy equation. The fracture pattern is realistic when only the tensile part of the strain energy is taken for fracture evolution. The fracture pattern is also observed in the compressive region when both tensile and compressive energy for crack evolution are taken into account, which is unrealistic. Length scale has a certain specific effect on the failure load of the specimen.

Fourth-order PFM is implemented using C¹ non-Sibsonian type of shape functions. The derivation and implementation are carried out for both the second-order and fourth-order PFM. The length scale effect on both models is shown. The better accuracy and convergence rate of the fourth-order PFM over second-order PFM are studied using the current approach. The critical difference between the isotropic phase field and the hybrid phase field approach is also presented to showcase the importance of strain energy decomposition in PFM.

The purpose of this study is to investigate the post-buckling analysis of functionally graded columns by using three analytical, approximate and numerical methods. A pre-defined function as an initial assumption for the post-buckling path is introduced to solve the differential equation. The finite difference method is used to approximate the lateral deflection of the column based on the differential equation. Moreover, the finite element method is used to derive the tangent stiffness matrix of the column.

The non-linear buckling analysis of functionally graded materials is carried out by using three analytical, finite difference and finite element methods. The elastic deformation and Euler-Bernoulli beam theory are considered to establish the constitutive and kinematics relations, respectively. The governing differential equation of the post-buckling problem is derived through the energy method and the calculus variation.

An incremental iterative solution and the perturbation of the displacement vector at the critical buckling point are performed to determine the post-buckling path. The convergence of the finite element results and the effects of geometric and material characteristics on the post-buckling path are investigated.

The key point of the research is to compare three methods and to detect error sources by considering the derivation process of relations. This comparison shows that a non-incremental solution in the analytical and finite difference methods and an initial assumption in the analytical method lead to an error in results. However, the post-buckling path in the finite element method is traced by the updated tangent stiffness matrix in each load step without any initial limitation.

The purpose of this study is to numerically model the rigid pavement resting on two-parameter soil and to examine its modal parameters.

This study is carried out using a one-dimensional beam element with three rotational and three translational degrees of freedom based on the finite element method. MATLAB programming is used to perform the free vibration analysis of the rigid pavement.

Cyclic frequency and their corresponding mode shapes were determined. It has been investigated how cyclic frequency changes as a result of variations in the thickness, span length of pavement, shear modulus, modulus of subgrade, different boundary conditions and element discretization. Thickness of the pavement and span length has greater effect on the cyclic frequency. Maximum increase of 29.7% is found on increasing the thickness, whereas the cyclic frequency decreases by 63.49% on increasing span length of pavement.

The pavement's free vibration is the sole subject of the current investigation. This study limits for the preliminary design phase of rigid pavements, where a complete three-dimensional finite element analysis is unnecessary. The current approach can be extended to future research using a different method, such as finite element grilling technique, mesh-free technique on reinforced concrete pavements or jointed concrete pavements.

The finite element approach adopted in this paper involves six degrees of freedom for each node. Furthermore, to the best of the authors’ knowledge, no prior study has done seven separate parametric investigations on the modal analysis of rigid pavement resting on two-parameter soil.

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.

The response of the Pine Flat dam–water–foundation rock system is studied by a new described approach (i.e. FE-(FE-TE)-FE). The initial part of study is focused on the time harmonic analysis. In this part, it is possible to compare the transfer functions against corresponding responses obtained by the FE-(FE-HE)-FE approach (referred to as exact method which employs a rigorous fluid hyper-element). Subsequently, the transient analysis is carried out. In that part, it is only possible to compare the results for low and high normalized reservoir length cases. Therefore, the sensitivity of results is controlled due to normalized reservoir length values.

In the present study, dynamic analysis of a typical concrete gravity dam–water–foundation rock system is formulated by the FE-(FE-TE)-FE approach. In this technique, dam and foundation rock are discretized by plane solid finite elements while, water domain near-field region is discretized by plane fluid finite elements. Moreover, the H-W (i.e. Hagstrom–Warburton) high-order condition is imposed at the reservoir truncation boundary. This task is formulated by employing a truncation element at that boundary. It is emphasized that reservoir far-field is excluded from the discretized model.

High orders of H-W condition, such as O5-5 considered herein, generate highly accurate responses for both possible excitations under both types of full reflective and absorptive reservoir bottom conditions. It is such that transfer functions are hardly distinguishable from corresponding exact responses obtained through the FE-(FE-HE)-FE approach in time harmonic analyses. This is controlled for both low and high normalized reservoir length cases (L/H = 1 and 3). Moreover, it can be claimed that transient analysis leads practically to exact results (in numerical sense) when one is employing high order H-W truncation element. In other words, the results are not sensitive to reservoir normalized length under these circumstances.

Dynamic analysis of concrete gravity dam–water–foundation rock systems is formulated by a new method. The salient aspect of the technique is that it utilizes H-W high-order condition at the truncation boundary. The method is discussed for all types of excitation and reservoir bottom conditions.

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.

This article aims to present a direct solution to handle linear constraints in finite element (FE) analysis without penalties or the Lagrange multipliers introduced.

First, the system of linear equations corresponding to the linear constraints is solved for the leading variables in terms of the free variables and the constants. Then, the reduced system of equilibrium equations with respect to the free variables is derived from the finite-dimensional virtual work equation. Finally, the algorithm is designed.

The proposed procedure is promising in three typical cases: (1) to enforce displacement constraints in any direction; (2) to implement local refinements by allowing hanging nodes from element subdivision and (3) to treat non-matching grids of distinct parts of the problem domain. The procedure is general and suitable for 3D non-linear analyses.

The algorithm is fitted only to the Galerkin-based numerical methods.

The proposed procedure does not need Lagrange multipliers or penalties. The tangential stiffness matrix of the reduced system of equilibrium equations reserves positive definiteness and symmetry. Besides, many contemporary Galerkin-based numerical methods need to tackle the enforcement of the essential conditions, whose weak forms reduce to linear constraints. As a result, the proposed procedure is quite promising.

The purpose of this study is to develop an efficient numerical procedure for simulating the effect of printing orientation, as one of the primary sources of anisotropy in 3D-printed components, on their fracture properties.

The extended finite element method and the cohesive zone model (XFEM-CZM) are used to develop subroutines for fracture simulation. The ability of two prevalent models, i.e. the continuous-varying fracture properties (CVF) model and the weak plane model (WPM), and a combination of both models (WPM-CVF) are evaluated to capture fracture behavior of the additively manufactured samples. These models are based on the non-local and local forms of the anisotropic maximum tangential stress criterion. The numerical models are assessed by comparing their results with experimental outcomes of 16 different configurations of polycarbonate samples printed using the material extrusion technique.

The results demonstrate that the CVF exaggerates the level of anisotropy, and the WPM cannot detect the mild anisotropy of 3D-printed parts, while the WPM-CVF produces the best results. Additionally, the non-local scheme outperforms the local approach in terms of finite element analysis performance, such as mesh dependency, robustness, etc.

This paper provides a method for modeling anisotropic fracture in 3D-printed objects. A new damage model based on a combination of two prevalent models is offered. Moreover, the developed subroutines for implementing the non-local anisotropic fracture criterion enable a reliable crack propagation simulation in media with varying degrees of complication, such as anisotropy.

The calculation of electromagnetic fields can involve many degrees of freedom (DOFs) to achieve accurate results. The DOFs are directly related to the computational effort of the simulation. The effort is decreased by using the proper generalized decomposition (PGD) and proper orthogonalized decomposition (POD). The purpose of this study is to combine the advantages of both methods. Therefore, a hybrid enrichment strategy is proposed and applied to different electromagnetic formulations.

The POD is an a-priori method, which exploits the solution space by decomposing reference solutions of the field problem. The disadvantage of this method is given by the unknown number of solutions necessary to reconstruct an accurate field representation. The PGD is an a-priori approach, which does not rely on reference solutions, but require much more computational effort than the POD. A hybrid enrichment strategy is proposed, based on building a small POD model and using it as a starting point of the PGD enrichment process.

The hybrid enrichment process is able to accurately approximate the reference system with a smaller computational effort compared to POD and PGD models. The hybrid enrichment process can be combined with the magneto-dynamic T-Ω formulation and the magnetic vector potential formulation to solve eddy current or non-linear problems.

The PGD enrichment process is improved by exploiting a POD. A linear eddy current problem and a non-linear electrical machine simulation are analyzed in terms of accuracy and computational effort. Further the PGD-AV formulation is derived and compared to the PGD-T-Ω reduced order model.

This paper aims to introduce a new numerical approach based on the local weak form and the Petrov–Galerkin idea to numerically simulation of a predator–prey system with two-species, two chemicals and an additional chemotactic influence.

In the first proceeding, the space derivatives are discretized by using the direct meshless local Petrov–Galerkin method. This generates a nonlinear algebraic system of equations. The mentioned system is solved by using the Broyden’s method which this technique is not related to compute the Jacobian matrix.

This current work tries to bring forward a trustworthy and flexible numerical algorithm to simulate the system of predator–prey on the nonrectangular geometries.

The proposed numerical results confirm that the numerical procedure has acceptable results for the system of partial differential equations.