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.

This paper presents the experimental and numerical results of the bending properties of low-height prestressed T-beams. The purpose is to study the bearing capacity, failure state and strain distribution of low-height prestressed T-beams.

First, two 13 m-long full-size test beams were fabricated with different positions of prestressed steel bundles in the span. The load–deflection curves and failure patterns of each test beam were obtained through static load tests. Secondly, the test data were used to validate the finite element model developed to simulate the flexural behavior of low-height prestressed T-beams. Finally, the influence of different parameters (the number of prestressed steel bundles, initial prestress and concrete strength grade) on the flexural performance of the test beams is studied by using a finite element model.

The test results show that when the distance of the prestressed steel beam from the bottom height of the test beam increases from 40 to 120 mm, the cracking load of the test beam decreases from 550.00 to 450.00 kN, reducing by 18.18%, and the ultimate load decreases from 1338.15 to 1227.66 kN, reducing by 8.26%, therefore, the increase of the height of the prestressed steel beam reduces the bearing capacity of the test beam. The numerical simulation results show that when the number of steel bundles increases from 2 to 9, the cracking load increases by 183.60%, the yield load increases by 117.71% and the ultimate load increases by 132.95%. Therefore, the increase in the number of prestressed steel bundles can increase the cracking load, yield load and ultimate load of the test beam. When the initial prestress is from 695 to 1,395 MPa, the cracking load increases by 69.20%, the yield load of the bottom reinforcement increases by 31.61% and the ultimate load increases by 3.97%. Therefore, increasing the initial prestress can increase the cracking load and yield load of the test beam, but it has little effect on the ultimate load. The strength grade of concrete increases from C30 to C80, the cracking load is about 455.00 kN, the yield load is about 850.00 kN and the ultimate load is increased by 4.90%. Therefore, the improvement in concrete strength grade has little influence on the bearing capacity of the test beam.

Based on the experimental study, the bearing capacity of low-height prestressed T-beams with different prestressed steel beam heights is calculated by finite element simulation, and the influence of different parameters on the bearing capacity is discussed. This method not only ensures the accuracy of bearing capacity assessment, but also does not require a large number of samples and has a certain economy. The study of prestressed low-height T-beams is of great significance for understanding the principle and application of prestressed technology. Research on the mechanical behavior and performance of low-height prestressed T beams can provide a scientific basis and technical support for the design and construction of prestressed concrete structures. In addition, the study of prestressed low-height T-beams can also provide a reference for the optimization design and construction of other structural types.

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.

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.

Beginners and even experienced ones have difficulties in completing the structural fire analysis due to numerical difficulties such as convergence errors and singularity and have to spend a lot of time making many repetitive changes on the model. The aim of this article is to highlight the advantages of explicit solver which can eliminate the mentioned difficulties in finite element analysis containing highly nonlinear contacts, clearance between modeled parts at the beginning and large deflections because of high temperature. This article provides important information, especially for researchers and engineers who are new to structural fire analysis.

The finite element method is utilized to achieve mentioned purposes. First, a comparative study is conducted between implicit and explicit solvers by using Abaqus. Then, a validation process is carried out to illustrate the explicit process by using sequentially coupled heat transfer and structural analysis.

Explicit analysis offers an easier solution than implicit analysis for modeling multi-bolted connections under high temperatures. An optimum mesh density for bolted connections is presented to reflect the realistic structural behavior. Presented explicit process with the offered mesh density is used in the validation of an experimental study on multi-bolted splice connection under ISO 834 standard fire curve. A good agreement is achieved.

What makes the study valuable is that the points to be considered in the structural fire analysis are examined and it is a guide that future researchers can benefit from. This is especially true for modeling and analysis of multi-bolted connections in finite element software under high temperatures. The article can help to shorten and even eliminate the iterative debugging phases, which is a problematic and very time-consuming process for many researchers.

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.

When the reinforced concrete beams are reinforced by bonding steel plates to the bottom, excessive use of steel plates will make the reinforced concrete beams become super-reinforced beams, and there are security risks in the actual use of super-reinforced beams. In order to avoid the occurrence of this situation, the purpose of this paper is to study the calculation method of the maximum number of bonded steel plates to reinforce reinforced concrete beams.

First of all, when establishing the limit failure state of the reinforced member, this paper comprehensively considers the role of the tensile steel bar and steel plate and takes the load effect before reinforcement as the negative contribution of the maximum number of bonded steel plates that can be used for reinforcement. Through the definition of the equivalent tensile strength, equivalent elastic modulus and equivalent yield strain of the tensile steel bar and steel plate, a method to determine the relative limit compression zone height of the reinforced member is obtained. Second, based on the maximum ratio of (reinforcement + steel plate), the relative limit compression zone height and the equivalent tensile strength of the tensile steel bar and steel plate of the reinforced member, the calculation method of the maximum number of bonded steel plates is derived. Then, the static load test of the test beam is carried out and the corresponding numerical model is established, and the reliability of the numerical model is verified by comparison. Finally, the accuracy of the calculation method of the maximum number of bonded steel plates is proved by the numerical model.

The numerical simulation results show that when the steel plate width is 800 mm and the thickness is 1–4 mm, the reinforced concrete beam has a delayed yield platform when it reaches the limit state, and the failure mode conforms to the basic stress characteristics of the balanced-reinforced beam. When the steel plate thickness is 5–8 mm, the sudden failure occurs without obvious warning when the reinforced concrete beam reaches the limit state. The failure mode conforms to the basic mechanical characteristics of the super-reinforced beam failure, and the bending moment of the beam failure depends only on the compressive strength of the concrete. The results of the calculation and analysis show that the maximum number of bonded steel plates for reinforced concrete beams in this experiment is 3,487 mm². When the width of the steel plate is 800 mm, the maximum thickness of the steel plate can be 4.36 mm. That is, when the thickness of the steel plate, the reinforced concrete beam is still the balanced-reinforced beam. When the thickness of the steel plate, the reinforced concrete beam will become a super-reinforced beam after reinforcement. The calculation results are in good agreement with the numerical simulation results, which proves the accuracy of the calculation method.

This paper presents a method for calculating the maximum number of steel plates attached to the bottom of reinforced concrete beams. First, based on the experimental research, the failure mode of reinforced concrete beams with different number of steel plates is simulated by the numerical model, and then the result of the calculation method is compared with the result of the numerical simulation to ensure the accuracy of the calculation method of the maximum number of bonded steel plates. And the study does not require a large number of experimental samples, which has a certain economy. The research result can be used to control the number of steel plates in similar reinforcement designs.

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.