Search results

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.

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.

Simulations exist for the prediction of the behaviour of building structural systems under fire, including two-way coupled fire-structure interaction. However, these simulations do not include detailed models of the connections, whereas these connections may impact the overall behaviour of the structure. Therefore, this paper proposes a two-scale method to include screw connections.

The two-scale method consists of (a) a global-scale model that models the overall structural system and (b) a small-scale model to describe a screw connection. Components in the global-scale model are connected by a spring element instead of a modelled screw, and the stiffness of this spring element is predicted by the small-scale model, updated at each load step. For computational efficiency, the small-scale model uses a proprietary technique to model the behaviour of the threads, verified by simulations that model the complete thread geometry, and validated by existing pull-out experiments. For four screw failure modes, load-deformation behaviour and failure predictions of the two-scale method are verified by a detailed system model. Additionally, the two-scale method is validated for a combined load case by existing experiments, and demonstrated for different temperatures. Finally, the two-scale method is illustrated as part of a two-way coupled fire-structure simulation.

It was shown that proprietary ”threaded connection interaction” can predict thread relevant failure modes, i.e. thread failure, shank tension failure, and pull-out. For bearing, shear, tension, and pull-out failure, load-deformation behaviour and failure predictions of the two-scale method correspond with the detailed system model and Eurocode predictions. Related to combined load cases, for a variety of experiments a good correlation has been found between experimental and simulation results, however, pull-out simulations were shown to be inconsistent.

More research is needed before the two-scale method can be used under all conditions. This relates to the failure criteria for pull-out, combined load cases, and temperature loads.

The two-scale method bridges the existing very detailed small-scale screw models with present global-scale structural models, that in the best case only use springs. It shows to be insightful, for it contains a functional separation of scales, revealing their relationships, and it is computationally efficient as it allows for distributed computing. Furthermore, local small-scale non-convergence (e.g. a screw failing) can be handled without convergence problems in the global-scale structural model.

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.

This study aims to derive a novel spatial numerical method based on multidimensional local Taylor series representations for solving high-order advection-diffusion (AD) equations.

The parabolic AD equations are reduced to the nonhomogeneous elliptic system of partial differential equations by utilizing the Chebyshev spectral collocation method (ChSCM) in the temporal variable. The implicit-explicit local differential transform method (IELDTM) is constructed over two- and three-dimensional meshes using continuity equations of the neighbor representations with either explicit or implicit forms in related directions. The IELDTM yields an overdetermined or underdetermined system of algebraic equations solved in the least square sense.

The IELDTM has proven to have excellent convergence properties by experimentally illustrating both h-refinement and p-refinement outcomes. A distinctive feature of the IELDTM over the existing numerical techniques is optimizing the local spatial degrees of freedom. It has been proven that the IELDTM provides more accurate results with far fewer degrees of freedom than the finite difference, finite element and spectral methods.

This study shows the derivation, applicability and performance of the IELDTM for solving 2D and 3D advection-diffusion equations. It has been demonstrated that the IELDTM can be a competitive numerical method for addressing high-space dimensional-parabolic partial differential equations (PDEs) arising in various fields of science and engineering. The novel ChSCM-IELDTM hybridization has been proven to have distinct advantages, such as continuous utilization of time integration and optimized formulation of spatial approximations. Furthermore, the novel ChSCM-IELDTM hybridization can be adapted to address various other types of PDEs by modifying the theoretical derivation accordingly.

The long-span continuous rigid-frame bridges are commonly constructed by the section-by-section symmetrical balance suspension casting method. The deflection of these bridges is increasing over time. Wet joints are a typical construction feature of continuous rigid-frame bridges and will affect their integrity. To investigate the sensitivity of shear surface quality on the mechanical properties of long-span prestressed continuous rigid-frame bridges, a large serviced bridge is selected for analysis.

Its shear surface is examined and classified using the damage measuring method, and four levels are determined statistically based on the core sample integrity, cracking length and cracking depth. Based on the shear-friction theory of the shear surface, a 3D solid element-based finite element model of the selected bridge is established, taking into account factors such as damage location, damage number and damage of the shear surface. The simulated results on the stress distribution of the local segment, the shear surface opening and the beam deflection are extracted and analyzed.

The findings indicate that the main factors affecting the ultimate shear stress and shear strength of the shear surface are size, shear reinforcements, normal stress and friction performance of the shear surface. The connection strength of a single or a few shear surfaces decreases but with little effect on the local stress. Cracking and opening mainly occur at the 1/4 span. Compared with the rigid “Tie” connection, the mid-span deflection of the main span increases by 25.03% and the relative deflection of the section near the shear surface increases by 99.89%. However, when there are penetrating cracks and openings in the shear surface at the 1/2 span, compared with the 1/4 span position, the mid-span deflection of the main span and the relative deflection of the cross-section increase by 4.50%. The deflection of the main span increases with the failure of the shear surface.

These conclusions can guide the analysis of deflection development in long-span prestressed continuous rigid-frame bridges.

The purpose of this paper is to develop an organized structure for damage detection of a cracked cantilever beam using finite element method and experimental method technique.

Due to presence of cracks the dynamic characteristics of structure change. The change in dynamic behavior has been used as one of the criteria of fault diagnosis for structures. Major characteristics of the structure which undergo change due to presence of crack are: natural frequencies, the amplitude responses due to vibration and the mode shapes. Therefore, an attempt has been made to formulate a smart technique for minimizing the amplitude of vibration for crack cantilever beam structures. In the analysis both single and double cracks are taken into account.

The results of the active vibration control experiments proved that piezoelectric sensor/actuator pair is an effective sensor and actuator configuration for active vibration control to reduce the amplitude of vibration for closed-loop system.

It is necessary that structures must safely work during its service life, but damages initiate a breakdown period on the structures which directly affect the industrial growth. It is a recognized fact that dynamic behavior of structures changes due to presence of crack. It has been observed that the presence of cracks in structures or in machine members leads to operational problem as well as premature failure.

1) The OE-MLPG penalty meshfree method is developed to solve cracked structure.(2) Smartening the numerical meshfree method by combining the particle swarm optimization (PSO) optimization algorithms and Voronoi computational geometric algorithm. (3). Selection of base functions, finding optimal penalty factor and distribution of appropriate nodal points to the accuracy of calculation in the meshless local Petrov–Galekrin (MLPG) meshless method.

Using appropriate shape functions and distribution of nodal points in local domains and sub-domains and choosing an approximation or interpolation method has an effective role in the application of meshless methods for the analysis of computational fracture mechanics problems, especially problems with geometric discontinuity and cracks. In this research, computational geometry technique, based on the Voronoi diagram (VD) and Delaunay triangulation and PSO algorithm, are used to distribute nodal points in the sub-domain of analysis (crack line and around it on the crack plane).

By doing this process, the problems caused by too closeness of nodal points in computationally sensitive areas that exist in general methods of nodal point distribution are also solved. Comparing the effect of the number of sentences of basic functions and their order in the definition of shape functions, performing the mono-objective PSO algorithm to find the penalty factor, the coefficient, convergence, arrangement of nodal points during the three stages of VD implementation and the accuracy of the answers found indicates, the efficiency of V-E-MLPG method with Ns = 7 and ß = 0.0037–0.0075 to estimation of 3D-stress intensity factors (3D-SIFs) in computational fracture mechanics.

The present manuscript is a continuation of the studies (Ref. [33]) carried out by the authors, about; feasibility assessment, improvement and solution of challenges, introduction of more capacities and capabilities of the numerical MLPG method have been used. In order to validate the modeling and accuracy of calculations, the results have been compared with the findings of reference article [34] and [35].

This study aims to elucidate the role of curved walls in the presence of identical mass of porous bed with identical heating at a wall for two heating objectives: enhancement of heat transfer to fluid saturated porous beds and reduction of entropy production for thermal and flow irreversibilities.

Two heating configurations have been proposed: Case 1: isothermal heating at bottom straight wall with cold side curved walls and Case 2: isothermal heating at left straight wall with cold horizontal curved walls. Galerkin finite element method is used to obtain the streamfunctions and heatfunctions associated with local entropy generation terms.

The flow and thermal maps show significant variation from Case 1 to Case 2 arrangements. Case 1 configuration may be the optimal strategy as it offers larger heat transfer rates at larger values of Darcy number, Da_m. However, Case 2 may be the optimal strategy as it provides moderate heat transfer rates involving savings on entropy production at larger values of Da_m. On the other hand, at lower values of Da_m (Da_m ≤ 10⁻³), Case 1 or 2 exhibits almost similar heat transfer rates, while Case 1 is preferred for savings of entropy production.

The concave wall is found to be effective to enhance heat transfer rates to promote convection, while convex wall exhibits reduction of entropy production rate. Comparison between Case 1 and Case 2 heating strategies enlightens efficient heating strategies involving concave or convex walls for various values of Da_m.

A novel method for modelling permanent magnets is investigated based on numerical approximations with rational functions. This study aims to introduce the AAA algorithm and other recently developed, cutting-edge mathematical tools, which provide outstandingly fast and accurate numerical computation of potentials and vector fields.

First, the AAA algorithm is briefly introduced along with its main variants and other advanced mathematical tools involved in the modelling. Then, the analysis of a circular Halbach array with a one-pole pair is carried out by means of the AAA-least squares method, focusing on vector potential and flux density in the bore and validating results by means of classic finite element software. Finally, the investigation is completed by a finite difference analysis.

AAA methods for field analysis prove to be strikingly fast and accurate. Results are in excellent agreement with those provided by the finite element model, and the very good agreement with those from finite differences suggests future improvements. They are also easy programming; the MATLAB code is less than 200 lines. This indicates they can provide an effective tool for rapid analysis.

AAA methods in magnetostatics are novel, but their extension to analogous physical problems seems straightforward. Being a meshless method, it is unlikely that local non-linearities can be considered. An aspect of particular interest, left for future research, is the capability of handling inhomogeneous domains, i.e. solving general interface problems.

The authors use cutting-edge mathematical tools for the modelling of complex physical objects in magnetostatics.