Elastic Finite Element Analysis Research Articles

Prediction of laser-welded deformation using artificial neural networks

When predicting welding deformation of the laser-manufactured vehicles and aerospaces, analytical solutions or empirical formulas are not usually accessible in complex problems. Based on the inherent deformation method, a machine learning (ML) approach for predicting welding deformation of welded structures is proposed based on an artificial neural network (ANN). This method is a promising substitute for analytical, empirical, and finite element (FE) solutions due to its accuracy, easy-to-use, efficiency, and universality. First, the outputs of the ANN are determined via dimensionless analysis and comparison of numerical results, which are dimensionally independent. Then, based on the inherent deformation method, the training and validation sets of the ANN are generated through an elastic finite element analysis. At last, the structure of the ANN is determined by analyzing the ANN prediction accuracy with different hidden layers, numbers of neurons, and activation functions. The results show that the ML solutions are in good agreement with the FE results, verifying the effectiveness and generalization ability of the proposed method.

Safety assessment of existing prestressed reinforced concrete bridge decks through different approaches

AbstractThe assessment of existing reinforced concrete (RC) structures, especially, bridges and viaducts, nowadays is one of the most relevant aspects for infrastructural engineers. Due to the old age of a large part of the railways and roadways assets, their safety, especially, with respect to the Ultimate Limit States (ULS) may result compromised. Safety assessment can be performed following different approaches as well as the most sophisticated methodologies based on the application of safety formats for non‐linear finite element (FE) analyses. This study presents the safety assessment at the ULS of existing prestressed RC bridge decks regarding bending moment and axial force failure mechanism adopting different methodologies. Specifically, the safety assessment will be performed through five different approaches (i.e., Courbon theory, elastic FE analysis, and NLFEAs according to the three global safety formats). The decks are modeled using, first, Courbon theory and then, girder FE model composed only of beam elements. The safety assessment is performed by comparing the design resisting moment and design acting moment in the section of interest, therefore, through a local approach. Within the local approach, design values are calculated using the partial safety factors. Afterwards, the decks are modeled with a 3D non‐linear FE model. Three different global safety formats are implemented within a global approach. The application of the different approaches leads to significantly different conclusions about the safety assessment. The two local approaches are the most conservative, therefore, lead to a lower safety margin. The application of global approaches based on global safety formats leads to a safety margin gain compared to local approaches, with significant variability in terms of safety margin gain depending on the specific global safety format.

Structural Performance of Large-Scale 3D-Printed Walls Subjected to Axial Compression Load

Despite the growing incorporation of additive manufacturing in the construction industry, a significant gap persists in the availability of design guidelines and comprehensive structural assessments for 3D-printed components. This research employs experimental methods to evaluate the performance of large-scale hollow unreinforced 3D-printed masonry walls under axial compressive load. Furthermore, a design methodology tailored for 3D-printed masonry walls under compressive loads is proposed in compliance with both American and Canadian masonry standards. This methodology is supported by the application of linear elastic 3D finite element analysis. Moreover, supplementary insights into material behavior are garnered through the examination of small specimens extracted from an additional 3D-printed wall. The results reveal that 3D printing reduces stiffness and compressive strength. The 3D-printed infills played an important role in redistributing stress after cracking. All full-scale 3D-printed masonry walls examined in this research exceeded the performance criteria set by the American masonry standard for load-bearing walls.

Study on Nanomaterials Coated Natural Coir Fibers as Crack Arrestor in Cement Composite

The process of inclusion of carbon nanotubes as fibers in cement paste has been proved to have optimistic effect as it enhances the tensile property of cement paste composite. Coir fibers have exceptional mechanical qualities and are thus employed as reinforcement in cement composites. Epoxy resin, which has a high Young’s modulus, is an ideal component for bonding carbon nanotubes (CNTs) to coir fiber. This paper describes a novel kind of nanocomposite made of L-12 epoxy resin and CNTs at the nanolevel, along with coir fibers at the microlevel which operate as crack arrestors. To remove surface contaminants, coir fibers are first treated with sodium hydroxide (NaOH). Epoxy/CNTs polymer coatings were developed at varying CNTs fractions (0.05, 0.1, 0.15, and 0.2 wt.% of cement). Multiwalled CNTs were combined in distilled water, followed by epoxy resin and hardener (9 : 1 v/v) polymer in an ultrasonic sonicator for 90 min to ensure full dispersion of CNTs within the epoxy polymer. This blend is now coated on the treated clustered coir fiber (length 10 cm, 10 strands) and reinforced along the length of a cement composite beam 20 mm × 20 mm × 80 mm in size. Tensile and three-point tests were performed to evaluate the mechanical characteristics of the hybrid composite. The linear elastic finite element analysis is employed to distinguish their failure phenomena via fatigue or fracture behavior. The microstructure behavior and the effect of coating material on the coir fibers were investigated using scanning electron microscope and EDX analysis. The reinforcing impact of nanopolymer coated coir fiber in cement composite diminished the tensile and flexural strength after 0.1 wt.% of CNT fraction but increased the composite’s durability and Young’s modulus. Fourier transform infrared spectroscopy analysis was carried out to assess the chemical interaction between the epoxy/CNTs and the coir fibers. The simulation was performed using ANSYS workbench, and the modeling results were within an acceptable 10% range of the experimental data. Nevertheless, it can be concluded that the hybrid composite is capable of enhancing the composite’s stress and strain capacity by regulating the fracture propagation process at the crack’s end.

Using the Point Method to estimate failure loads in 3D printed graphene-reinforced PLA notched plates

This work estimates failure loads in Fused Filament Fabrication (FFF) printed graphene-reinforced PLA (polylactic acid) plates containing different types of stress risers. With this aim, firstly, several notched plates are tested and conducted to fracture. Then, linear elastic Finite Element (FE) analyses are completed to define the corresponding stress profiles and, finally, the Point Method (PM) is applied to establish the failure criterion. This approach asserts that fracture conditions are achieved when the stress level equates the inherent strength (σ0) at a distance from the notch tip equal to L/2, so both parameters (related to each other through the material fracture toughness, Kmat) have been defined beforehand. The estimations of fracture loads obtained following this approach agree with the experimental results. Thus, the present work demonstrates the accuracy of the PM to estimate failure loads in this 3D printed material.

Feasibility of outrigger structural system for tall timber buildings: A numerical study

Tall timber buildings are prone to serviceability problems (e.g. lateral displacements and wind-induced accelerations). This article investigates the feasibility of timber outrigger system to build up to 20 storeys using linear elastic finite element analysis and relevant European and international standards such as EN 1990, EN 1991-1-1, EN 1991-1-4, EN 1995-1-1 and ISO 10137. The location of up to two outriggers was optimised based on different serviceability criteria. Various outrigger layouts were investigated, and their stiffness demands were discussed. The efficiency of outrigger structures in reducing lateral displacements and wind-induced accelerations is highlighted. A feasibility study considering serviceability requirements was performed using different basic wind velocities, number of storeys and gravity loads. Although the focus of the article is on serviceability requirements, some ultimate limit state considerations (forces and stresses) were briefly discussed. The results demonstrate the feasibility of building up to 16 storeys under medium to high wind velocities.

Deformation and Fracturing Behaviors of Sandstone Disc Containing a Circular Inclusion under Diametrical Compression: Insights from Integrated SG and DIC Experiments

Weak inclusions play a significant role in regulating the deformation and fracturing behavior of rock masses. In this study, sandstone specimens filled with circular inclusions were subjected to diametrical compression until failure. The local strain concentration intensity and cracking process at the rock-inclusion interface were specifically investigated using the integrated strain gauge (SG) and digital image correlation (DIC) techniques. The experimental results reveal that the presence of inclusions introduces more complex deformation and fracturing behaviors in the infilled sandstone specimens. The stiffness of the infilled sandstones, with either low-strength or high-strength inclusions, initially increases and then decreases as the inclusion diameter increases. The effective load-carrying capacity of the infilled sandstone is influenced by the competing effects of inclusion enhancement and hole weakening. Moreover, a more pronounced strain concentration occurs at the inclusion boundary when the inclusion diameter is larger and the inclusion strength is lower. The crack mode in the infilled red sandstone is independent of the strength of the inclusion. For larger hole diameters, fracture planes form along the compression line, whereas interfacial debonding occurs in red sandstone with smaller hole diameters. Finally, the stress distribution characteristics along the loaded diameter and transverse diameter were analyzed using elastic finite element analyses, which can well explain the cracking mechanism of the infilled disc specimens.

Stress and Deformation Analysis of REBCO Pancake Coils With Individual Turn Movement

In recent years, rare-earth barium copper oxide (REBCO) pancake coils have shown excellent performances in high magnetic field applications. The no-insulation (NI) winding technique enables to enhance the thermal stability and current density of REBCO pancake coils. Still, further developments of REBCO coils are expected; e.g., the mechanical deteriorations of REBCO tapes caused by electromagnetic forces have been reported. The screening currents are considered as one of the causes to damage the REBCO tapes because of a nonuniform current distribution. Meanwhile, the electromagnetic forces on non-impregnant NI REBCO coils, whose windings can be deformed separately with turns and moved easily, are complicated, and the detailed mechanism has not yet been clarified. Commonly, the circumferential movement/deformation of the winding is considered as a rigid body. In this study, we investigate the stresses and displacements due to the electromagnetic forces on NI REBCO pancake coils in high magnetic field using a current simulation and a 2D elastic finite element analysis in the radial and circumferential directions; i.e., each turn can individually move in the radial and circumferential directions considering the spiral winding structure. The simulation results show the circumferential movement/deformation of windings and the nonuniform hoop stress due to the electromagnetic force.

Assessment of residual stresses in welded T-joints using contour method

Welded T-joints are commonly used in manufacturing of stiffened panels or crossbeams. However, welding of web-to-flange joint induces residual stresses and distortions resulting in significant challenges during both manufacturing and service. These issues can ultimately impact the structural performance, including fatigue and stability, and lifetime of the weldment. The purpose of current paper is determining the distribution and magnitude of manufacturing-induced longitudinal residual stresses in welded T-joints. Influence of hot rolling, thermal cutting, prestressing, welding, etc. are all included using the so-called contour method. Deformed surfaces are measured using coordinate measuring machine after wire electrical discharge machining of specimens. An inverse linear elastic finite element analysis is carried out for assessing residual stresses after data smoothing. Two-dimensional longitudinal residual stress patterns are evaluated in representative cross-sections clearly indicating high stress peaks near the welds, even exceeding yield strength of the base material. In addition, longitudinal membrane residual stresses are determined based on through-thickness stresses. In addition, a simplified residual stress model, satisfying the requirement for equilibrium of internal forces and bending moments, is developed for welded T-joints with double-sided fillet welds based on contour method results.

Investigation of anchorage plate opening on bursting force

Three-dimensional (3D) Finite Element analyses were performed to investigate transverse stresses and bursting forces at the anchorage zone in post-tensioned concrete beams. The behavior of the concrete beam was investigated through linear elastic Finite Element analyses. To validate the Numerical modeling results, an experimental study was performed on a beam and results are compared in terms of ultimate strength. The effect of duct hole and eccentricity are investigated by varying diameter and eccentric distances on the concrete beam. Furthermore, the results are validated for bursting stress and bursting force. Altogether, the equation of bursting force in IS 1343:2012 hasn’t considered the effect of duct hole and eccentricity, the results are compared with IS 1343:2012, and an equation is presented considering the effect of duct and eccentricity.

Data-driven artificial neural network for elastic plastic stress and strain computation for notched bodies

An integration of artificial neural network (ANN) and finite element (FE) analysis was developed to predict elastic–plastic stress and strains at notch locations based on a linear elastic FE analysis solution. Two FE models were created for different materials and multiaxial load cases to produce hypothetical elastic and elastic–plastic stress and strain datasets. One model was based on an elastic state, and the other was based on an elastic–plastic state. The ANN was trained with the elastic stress data from the linear elastic FE model as input and the elastic–plastic stress–strain data from the nonlinear elastic–plastic FE model as output. The dataset was divided into three groups: training, verification, and testing data. The ANN was trained using the training data, evaluated using the verification data, and tested for generalizability using the testing data. The results showed that the proposed methodology can predict elastic–plastic stress and strain fields for notched bodies under multiaxial loadings accurately and efficiently using only the elastic FE analysis solution.

A linear elastic finite element approach to fatigue life estimation for defect laden materials

This work assumed that the stress-life behaviors of defect-laden materials converge at a certain low life cycle and attempted to estimate fatigue life by applying relative knock-downs to the endurance limit of a reference stress-life data. Linear elastic finite element analysis on small cracks initiated from defects was performed to evaluate the effect of defect size, aspect ratio, and location on endurance limit. While for internal defects size and shape are both influential, size and depth were found to be more important for surface defects. These findings were validated against the fatigue lives of laser powder bed fused Ti-6Al-4V specimens.

STRESS ANALYSIS OF TEETH, PERIODONTAL TISSUES, AND MAXILLA BONES DUE TO BRUXISM

This paper investigates the stress distribution on the teeth and jaw during bruxism that may provide additional insights to the consequences of bruxism, splint design considerations, a basis for possible therapy protocols, and a baseline for future bruxism studies. A three-dimensional (3D) solid model of a complete jaw comprised of teeth, periodontal ligaments, and maxilla bone is reconstructed based on a Cone Beam Computed Tomography (CBCT) scan of a patient. A finite element model of the jaw is then constructed using the geometry imported from the solid model that has been reworked to remove imperfections during reconstruction and to allow proper contact between teeth. A linear elastic finite element analyses were then performed to simulate the bruxism phenomenon which includes biting and grinding. Numerical analysis is conducted only to the maxilla, with forces are placed on the teeth’s surface where contact occurs. The analysis is conducted for all stages of bruxism and possible movements patterns. Results suggest that high stresses occur during grinding, indicating a strong correlation between bruxism and dental issues such as tooth wear and alveolar resorption.

Artificial neural networks-based [formula omitted]-integral prediction for cracked bodies under elasto-plastic deformation state –monotonic loading

Artificial neural networks (ANNs) integrated with a finite element (FE)-based equivalent domain integral method are developed to compute J-integral at the vicinity of crack tips through a time-efficient approach. Robust ANN models are trained to establish nonlinear relationships between FE predicted elastic and elasto-plastic stress, strain, and displacement fields of stainless steel (SS304). Subsequently, elastic–plastic J-integral can be determined by using only elastic FE analysis solution rather than computationally expensive elasto-plastic FE analysis solution. The results show that well-trained ANN models can efficiently and accurately determine J-integral around the crack tips on the basis of numerical elastic FE analysis solution.

Study on thermal deformation of hybrid printed circuit heat exchanger for advanced nuclear reactor

The intermediate heat exchanger (IHX) is one of the key components in advanced nuclear reactor power generation system. Printed circuit heat exchanger (PCHE) with mini-channels is a good candidate for the IHX. But the whole PCHE with tens of thousands of mini-channels is hard to be modeled and calculated at high temperature. Thus, the homogenization method, which can calculate the thermal deformation of the whole hybrid mini-channel PCHE, is conducted in this study and both the elastic and elastic-plastic finite element analysis methods are used. Here, the performance of the PCHE is assumed to be isotropic in the simulation. The results show that the macroscopic thermal deformation of the whole PCHE follows a gradually downward trend from hot side to cold side, while obvious thermal stress concentration is observed around the junctions between the PCHE core and cover plate. The stress and strain calculated by elastic and elastic-plastic methods under different constraints show quite different patterns. For the thermal deformation analysis of PCHE under weak constraint or small temperature difference conditions, the stress and strain results by elastic analysis and elastic-plastic analysis are similar so that the elastic analysis can be used. Under the strong condition, the stress of PCHE calculated by elastic analysis is much higher than that of elastic-plastic analysis, while the strain is lower than that of elastic-plastic analysis. The relative deviations of stress are 100.883% and 120.786% along path 1 and path 2 respectively. And the relative deviations of strain are 6.987% and 10.139% along path 1 and path 2 respectively. The elastic-plastic analysis should be used instead of the elastic analysis.

Characteristics of compressive stress around dowel joint in concrete pavement system

ABSTRACT Although Friberg’s theoretical analysis is frequently used to analyse the dowel bar-concrete interaction in jointed plain concrete pavement (JPCP) system, the compressive stress concentration at top and bottom of the dowel slot cannot be explained with this theory. In this paper, two modifications including the elastic deformation of concrete support and the stress distribution factor considering the sinusoidal contact stress distribution are proposed to improve the Friberg’s theoretical analysis. Modifications are then validated against elastic finite element analysis (FEA) results in terms of the contact stress distribution, the maximum compressive stress as well as the range of compression zone. To further investigate the maximum compressive stress at the joint surface, relevant parameters including the dowel bar diameter, the modulus of elasticity of concrete as well as the joint width are studied through elastic finite element analysis (FEA). Apart from proposing referenced moduli of dowel support, a close relationship between the modified modulus of dowel support and dowel bar diameter is also found and verified. Finally, differences between the Friberg’s theoretical analysis and the modified method are summarised in a flow chart at the end of the paper.

Directions of high cycle fatigue cracks emanating from circular notches studied by optical profilometry

Current models for predicting the fatigue endurance of notched solids use the stresses along a straight line, beginning at the notch root, as a simplification of the real crack propagation path. In this work, the experimental crack paths for hollow notched samples were analysed through different microscopy techniques, with the objective of establishing high cycle fatigue crack growth directions in a mild steel. Fully reversed tension–compression fatigue tests (R=−1) of thin-walled tube specimens with a passing-through hole were carried out. The crack paths observed in the outer cylindrical surface were studied in each case, with special attention to the crack initiation point and the crack direction along the first grains. Moreover, the analysis of the fracture surfaces allowed the same analysis to be performed to determine the internal crack paths. It was observed that the crack initiation point was close to the maximum principal stress point at the hole contour as obtained from linear elastic finite element analysis, and the crack direction in its initiation was generally close to Mode I direction, contrary to the conventionally accepted 45∘ crack growth direction.

Parametric optimization of corner radius in hexagonal honeycombs under in-plane compression

What is the optimum corner radius for a regular hexagonal honeycomb in the context of in-plane, quasistatic compression loading? Drawing inspiration from social insect hexagonal cell nests, where a non-zero corner radius appears to be a design feature, a 400-point design of experiments study is conducted using 2D plane strain Finite Element Analysis (FEA) to study the influence of cell size, beam thickness and corner radius on effective modulus and maximum corner stress in the honeycomb cells. An experimental study is conducted to examine these relationships beyond the bounds of the small deformation, linear elastic FEA analysis. The study finds that corner radii always increase the effective modulus of the honeycomb, even after accounting for the additional mass associated with the corner fillet. A key finding of this work is that while corner radii also reduces the maximum corner stress, there exists a clear optimum, beyond which stresses rise again as the corner radius increases in magnitude. This optimum corner radius is shown to be a function of the beam thickness (t) to beam length (l) ratio (t/l), with the optimum value increasing with increasing t/l. The experimental study shows that the presence of a corner radius shifts the failure mechanism from nodal fracture to plastic hinging for honeycombs with thick beams, and may have benefits for energy absorption applications. This work makes the case for the treatment of the corner radius as an independent design feature for optimization in the wider context of cellular materials, as well as has implications for the study of the geometry of insect nests.

Minimum thickness of solid concrete floor slab with column line beams to avoid walking vibration discomfort

AbstractIn this study, an attempt has been made to evaluate the minimum slab thickness of reinforced concrete (RC) solid floor slab with column‐line beams to minimize vibration discomfort when subjected to human walking‐induced excitation. An elastic time history finite element analysis methodology has been adopted for the study. An intermediate floor panel of a multistory building has been modeled using appropriate material properties, element dimensioning and stiffness modifiers. The developed model was subjected to a walking excitation using a time history function. The finite element modeling and time history analysis scheme has been validated against experimental study described in the design guide by Fanella and Mota. The study was then focused on studying the walking‐induced vibration characteristics of the floor slab using the span length and slab thickness as the study parameters. The acceleration data obtained for the central point of the slab was subjected to fast fourier transform (FFT) to convert the response to the frequency domain. Comparisons were then made to obtain the minimum slab thickness for each span length using both frequency and acceleration criteria following ACI serviceability limit and the design guide limit by Fanella and Mota. It was observed that a combination of the two limits were needed to determine the minimum slab thickness. The findings will help engineers to design slabs with adequate thickness to minimize human walking‐induced vibration related discomfort.

Bounded Elastic Moduli Multiplier Technique for Limit Loads: Part 1—Theory

Abstract Statically admissible stress distributions are necessary to evaluate lower bound limit loads. Over the last three decades, several methods have been postulated to obtain these distributions using iterative elastic finite element analyses. Some of the pioneering techniques are the reduced modulus, r-node, elastic compensation, and linear matching methods, to mention a few. A new method, called the bounded elastic moduli multiplier technique (BEMMT), is proposed and the theoretical underpinnings thereof are explained in this paper. BEMMT demonstrates greater robustness, more generality, and better stress distributions, consistently leading to lower-bound limit loads that are closer to elastoplastic finite element analysis estimates. It also leads us to question the validity of the prevailing power law-based stationary stress distributions. Part 2 of this paper offers several case studies to validate this claim.