Search results

The paper's goal is to examine and illustrate the useful uses of submodeling in finite element modeling for topology optimization and stress analysis. The goal of the study is to demonstrate how submodeling – more especially, a 1D approach – can reliably and effectively produce ideal solutions for challenging structural issues. The paper aims to demonstrate the usefulness of submodeling in obtaining converged solutions for stress analysis and optimized geometry for improved fatigue life by studying a cantilever beam case and using beam formulations. In order to guarantee the precision and dependability of the optimization process, the developed approach will also be validated through experimental testing, such as 3-point bending tests and 3D printing. Using 3D finite element models, the 1D submodeling approach is further validated in the final step, showing a strong correlation with experimental data for deflection calculations.

The authors conducted a literature review to understand the existing research on submodeling and its practical applications in finite element modeling. They selected a cantilever beam case as a test subject to demonstrate stress analysis and topology optimization through submodeling. They developed a 1D submodeling approach to streamline the optimization process and ensure result validity. The authors utilized beam formulations to optimize and validate the outcomes of the submodeling approach. They 3D-printed the optimized models and subjected them to a 3-point bending test to confirm the accuracy of the developed approach. They employed 3D finite element models for submodeling to validate the 1D approach, focusing on specific finite elements for deflection calculations and analyzed the results to demonstrate a strong correlation between the theoretical models and experimental data, showcasing the effectiveness of the submodeling methodology in achieving optimal solutions efficiently and accurately.

The findings of the paper are as follows: 1. The use of submodeling, specifically a 1D submodeling approach, proved to be effective in achieving optimal solutions more efficiently and accurately in finite element modeling. 2. The study conducted on a cantilever beam case demonstrated successful stress analysis and topology optimization through submodeling, resulting in optimized geometry for enhanced fatigue life. 3. Beam formulations were utilized to optimize and validate the outcomes of the submodeling approach, leading to the successful 3D printing and testing of the optimized models through a 3-point bending test. 4. Experimental results confirmed the accuracy and validity of the developed submodeling approach in streamlining the optimization process. 5. The use of 3D finite element models for submodeling further validated the 1D approach, with specific finite elements showing a strong correlation with experimental data in deflection calculations. Overall, the findings highlight the effectiveness of submodeling techniques in achieving optimal solutions and validating results in finite element modeling, stress analysis and optimization processes.

The originality and value of the paper lie in its innovative approach to utilizing submodeling techniques in finite element modeling for structural analysis and optimization. By focusing on the reduction of finite element models and the creation of smaller, more manageable models through submodeling, the paper offers designers a more efficient and accurate way to achieve optimal solutions for complex problems. The study's use of a cantilever beam case to demonstrate stress analysis and topology optimization showcases the practical applications of submodeling in real-world scenarios. The development of a 1D submodeling approach, along with the utilization of beam formulations and 3D printing for experimental validation, adds a novel dimension to the research. Furthermore, the paper's integration of 1D and 3D submodeling techniques for deflection calculations and validation highlights the thoroughness and rigor of the study. The strong correlation between the finite element models and experimental data underscores the reliability and accuracy of the developed approach. Overall, the originality and value of this paper lie in its comprehensive exploration of submodeling techniques, its practical applications in structural analysis and optimization and its successful validation through experimental testing.

The main aim of the current paper is to find a numerical plan for hydraulic fracturing problem with application in extracting natural gases and oil.

First, time discretization is accomplished via Crank-Nicolson and semi-implicit techniques. At the second step, a high-order finite element method using quadratic triangular elements is proposed to derive the spatial discretization. The efficiency and time consuming of both obtained schemes will be investigated. In addition to the popular uniform mesh refinement strategy, an adaptive mesh refinement strategy will be employed to reduce computational costs.

Numerical results show a good agreement between the two schemes as well as the efficiency of the employed techniques to capture acceptable patterns of the model. In central single-crack mode, the experimental results demonstrate that maximal values of displacements in x- and y- directions are 0.1 and 0.08, respectively. They occur around both ends of the line and sides directly next to the line where pressure takes impact. Moreover, the pressure of injected fluid almost gained its initial value, i.e. 3,000 inside and close to the notch. Further, the results for non-central single-crack mode and bifurcated crack mode are depicted. In central single-crack mode and square computational area with a uniform mesh, computational times corresponding to the numerical schemes based on the high order finite element method for spatial discretization and Crank-Nicolson as well as semi-implicit techniques for temporal discretizations are 207.19s and 97.47s, respectively, with 2,048 elements, final time T = 0.2 and time step size τ = 0.01. Also, the simulations effectively illustrate a further decrease in computational time when the method is equipped with an adaptive mesh refinement strategy. The computational cost is reduced to 4.23s when the governed model is solved with the numerical scheme based on the adaptive high order finite element method and semi-implicit technique for spatial and temporal discretizations, respectively. Similarly, in other samples, the reduction of computational cost has been shown.

This is the first time that the high-order finite element method is employed to solve the model investigated in the current paper.

A finite element computational methodology on a curved boundary using an efficient subparametric point transformation is presented. The proposed collocation method uses one-side curved and two-side straight triangular elements to derive exact subparametric shape functions.

Our proposed method builds upon the domain discretization into linear, quadratic and cubic-order elements using subparametric spaces and such a discretization greatly reduces the computational complexity. A unique subparametric transformation for each triangle is derived from the unique parabolic arcs via a one-of-a-kind relationship between the nodal points.

The novel transformation derived in this paper is shown to increase the accuracy of the finite element approximation of the boundary value problem (BVP). Our overall strategy is shown to perform well for the BVP considered in this work. The accuracy of the finite element approximate solution increases with higher-order parabolic arcs.

The proposed collocation method uses one-side curved and two-side straight triangular elements to derive exact subparametric shape functions.

In order to improve the efficiency and reliability of simulation analysis for composite riveting structures in engineering products, a comparative study was conducted on different forms of riveting simulation methods.

Five different rivent simulation models were established using the finite element method, including rigid element CE, flexible element Rbe3 and beam element, and their results were future compared and analyzed.

Under the given technical parameters, the simulation method of Rbe3 (with holes) + beam can meet the analysis requirements of complex engineering products in terms of the rationality of rivet load distribution, calculation error and relatively efficient modeling.

This study proposes a simulation method for the riveting structure of carbon fiber composite materials for engineering applications. This method can satisfy the simulation analysis requirements of transportation vehicles in terms of modeling time, computational efficiency and accuracy. The research can provide technical support for the riveting process and mechanical analysis between carbon fiber composite components in transportation products.

This study aims to propose and numerically assess different ways of discretising a very weak formulation of the Poisson problem.

We use integration by parts twice to shift smoothness requirements to the test functions, thereby allowing low-regularity data and solutions.

Various conforming discretisations are presented and tested, with numerical results indicating good accuracy and stability in different types of problems.

This is one of the first articles to propose and test concrete discretisations for very weak variational formulations in primal form. The numerical results, which include a problem based on real MRI data, indicate the potential of very weak finite element methods for tackling problems with low regularity.

The purpose of this study is to numerically model the rigid pavement resting on Pasternak soil and to examine its various response parameters and stress resultants like deflection, rotation, bending moment and shear force when subjected to aircraft loading.

The study is carried out using a one-dimensional (1D) beam element based on the finite element method (FEM). Each node in this element has three rotational and three translational degrees of freedom (DOF). MATLAB programming is used to perform the static analysis of rigid pavement.

Response parameters and stress resultants of the rigid pavement were determined. The FEM used in this work is validated by two closed-form numerical examples, which are in great accord with previous research findings with a maximum divergence of 4.64%, therefore verifying the finite element approach used in the current study. Additionally, various parametric studies have been carried out to study the variations in response parameters and stress resultants.

The investigation at hand focuses exclusively on the static analysis of the pavement. The study constraints pertaining to the preliminary design phase of rigid pavements are such that a comprehensive three-dimensional finite element analysis is deemed unnecessary.

As limited previous research had performed the static analysis of rigid pavement on Pasternak foundation with 6 DOF. Furthermore, no prior study has done seven separate parametric investigations on the static analysis of rigid pavement.

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.

The purpose of this study is to simulate the tensile and shear types of fractures using the mixed fracture criteria considering the energy evolution based on the dual bilinear cohesive zone model and investigate the dynamic propagation of tensile and shear fractures induced by an impact load in rock. The propagation of tension and shear at different scales induced by the impact load is also an important aspect of this study.

In this study, based on the well-developed dual bilinear cohesive zone model and combined finite element-discrete element method, the dynamic propagation of tensile and shear fractures induced by the impact load in rock is investigated. Some key technologies, such as the governing partial differential equations, fracture criteria, numerical discretisation and detection and separation, are introduced to form the global algorithm and procedure. By comparing with the tensile and shear fractures induced by the impact load in rock disc in typical experiments, the effectiveness and reliability of the proposed method are well verified.

The dynamic propagation of tensile and shear fractures in the laboratory- and engineering-scale rock disc and rock strata are derived. The influence of mesh sensitivity, impact load velocities and load positions are investigated. The larger load velocities may induce larger fracture width and entire failure. When the impact load is applied near the left support constraint boundary, concentrated shear fractures appear around the loading region, as well as induced shear fracture band, which may induce local instability. The proposed method shows good applicability in studying the propagation of tensile and shear fractures under impact loads.

The proposed method can identify fracture propagation via the stress and energy evolution of rock masses under the impact load, which has potential to be extended into the investigation of the mixed fractures and disturbance of in-situ stresses during dynamic strata mining in deep energy development.

Heterogeneous and micro-structured materials have been the object of multiscale and homogenization techniques aimed at recognizing the elastic properties of the equivalent continuum. The proposed investigation deals with the mechanical characterization of the heterogeneous material structured metamaterials through analyzing the ultimate strength using the limit analysis of the Representative Volume Element (RVE). To get the desired material strength, a novel finite element formulation based on the derivation of self-equilibrated solutions through the finite elements devoted to calculating the lower bound theorem has been implemented together with the limit analysis in Melàn’s formulation.

The finite element formulation is based on discrete mapping of Volterra dislocations in the structure using isoparametric representation. Using standard finite element techniques, the linear operator V, which relates the self-equilibrated internal solicitation to displacement-like nodal parameters, has been built through finite element discretization of displacement and strain.

The proposed work presented an elastic homogenization of the mechanical properties of an elementary cell with a geometry known in the literature, the isotropic truss. The matrix of elastic constants was calculated by subjecting the RVE to numerical load tests, simulated with a commercial FEM calculation code. This step showed the dependence of the isotropy properties, verified with Zener theory, on the density of the RVE. The isotropy condition of the material is only achieved for certain section ratios between body-centered cubic (BCC) and face-centered cubic (FCC), neglecting flexural effects at the nodes. The density that satisfies Zener’s conditions represents the isotropic geomatics of the isotropic truss.

For the isotropic case, the VFEM procedure was used to evaluate the isotropy of the limit domain and was compared with the Mises–Schleicher limit domain. The evaluation of residual ductility and dissipation energy allowed a measurement parameter for the limit anisotropy to be defined. The novelty of the proposal consisted in the formulation of both the linearized and the nonlinear limit locus of the material; hence, it furnished the starting point for further limit analysis of the structures whose elementary volume has been described through the proposed approach.

The purpose of this study is to present a novel extended hybrid analytical method (HAM) that leverages a two-dimensional (2-D) coupling between the semi-analytical Maxwell–Fourier analysis and the finite element method (FEM) in Cartesian coordinates.

The proposed model is applied to flat permanent-magnet linear electrical machines with rotor-dual. The magnetic field solution across the entire machine is established by coupling an exact analytical model (AM), designed for regions with relative magnetic permeability equal to unity, with a FEM in ferromagnetic regions. The coupling between AM and FEM occurs bidirectionally (x, y) along the edges separating teeth regions and their adjacent regions through applied boundary conditions.

The developed HAM yields accurate results concerning the magnetic flux density distribution, cogging force and induced voltage under various operating conditions, including magnetic or geometric parameters. A comparison with hybrid finite-difference and hybrid reluctance network methods demonstrates very satisfactory agreement with 2-D FEM.

The original contribution of this paper lies in establishing a direct coupling between the semi-analytical Maxwell–Fourier analysis and the FEM, particularly at the interface between adjacent regions with differing magnetic parameters.