Search results

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.

This paper aims to estimate the bearing capacity of a surface strip and circular footings lying on layered sand using numerical limit analysis.

Lower and upper bound limit analysis, as well as finite elements and second-order conic programming (SOCP), are used in this analysis. The yield criterion of Mohr-Coulomb is used to model soil behavior. Using this technique, stringent lower and upper bounds on ultimate bearing capacity can be achieved by assuming an associated flow law.

The obtained results indicate that the exact collapse load is typically being bracketed to within 6% about a mean of both the bounds. The obtained results are compared with the existing literature wherever applicable.

To the best of the authors’ knowledge, no study has used lower and upper bound limit analysis, as well as finite elements and SOCP, to estimate the bearing capacity of a surface strip and circular footings lying on layered sand.

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.

This paper aims to investigate the responses of single piles and pile groups due to tunneling-induced ground movements in a two-layered soil system. The analyses mainly focus on the additional single pile responses in terms of bending moment, lateral deflection, axial force, shaft resistance and pile settlement. Subsequently, a series of parametric studies were carried out to better understand the responses of single piles induced by tunneling. To give further understanding regarding the pile groups, a 2 × 2 pile group with two different pile head conditions, namely, free and capped, was considered.

Using the PLAXIS three-dimensional (3D) software, a full 3D numerical modeling is performed to investigate the effects of ground movements caused by tunneling on adjacent pile foundations. The numerical model was validated using centrifuge test data found in the literature. The relevance of the 3D model is also judged by comparison with the 2D plane strain model using the PLAXIS 2D code.

The numerical test results reveal that tunneling induces significant displacements and internal forces in nearby piles. The magnitude and distribution of internal forces depend mainly on the position of the pile toe relative to the tunnel depth and the distance between the pile and the vertical axis of the tunnel. As the volume loss increases from 1% to 3%, the apparent loss of pile capacity increases from 11% to 20%. By increasing the pile length from 0.5 to 1.5 times, the tunnel depth, the maximum pile settlement and lateral deflection decrease by about 63% and 18%, respectively. On the other hand, the maximum bending moment and axial load increase by about 7 and 13 times, respectively. When the pile is located at a distance of 2.5 times the tunnel diameter (Dt), the additional pile responses become insignificant. It was found that an increase in tunnel depth from 1.5Dt to 2.5Dt (with a pile length of 3Dt) increases the maximum lateral deflection by about 420%. Regarding the interaction between tunneling and group of piles, a positive group effect was observed with a significant reduction of the internal forces in rear piles. The maximum bending moment of the front piles was found to be higher than that of the rear piles by about 47%.

Soil is a complex material that shows differently in primary loading, unloading and reloading with stress-dependent stiffness. This general behavior was not possibly being accounted for in simple elastic perfectly plastic Mohr–Coulomb model which is often used to predict the behavior of soils. Thus, in the present study, the more advanced hardening soil model with small-strain stiffness (HSsmall) is used to model the non-linear stress–strain soil behavior. Moreover, unlike previous studies THAT are usually based on the assumption that the soil is homogeneous and using numerical methods by decoupled loadings under plane strain conditions; in this study, the pile responses have been exhaustively investigated in a two-layered soil system using a fully coupled 3D numerical analysis that takes into account the real interactions between tunneling and pile foundations. The paper presents a distinctive set of findings and insights that provide valuable guidance for the design and construction of shield tunnels passing through pile foundations.

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.

Fiber reinforced additive manufacturing (FRAM) is an emerging technology that combines additive manufacturing and composite materials. As a result, design freedom offered by the manufacturing process can be leveraged in design optimization. The purpose of the study is to propose a novel method that improves structural performance by optimizing 3D print orientation of FRAM components.

This work proposes a two-part design optimization method that optimizes 3D global print orientation and topology of a component to improve a structural objective function. The method considers two classes of design variables: (1) print orientation design variables and (2) density-based topology design variables. Print orientation design variables determine a unique 3D print orientation to influence anisotropic material properties. Topology optimization determines an optimal distribution of material within the optimized print orientation.

Two academic examples are used to demonstrate basic behavior of the method in tension and shear. Print orientation and sequential topology optimization improve structural compliance by 90% and 58%, respectively. An industry-level example, an aerospace component, is optimized. The proposed method is used to achieve an 11% and 15% reduction of structural compliance compared to alternative FRAM designs. In addition, compliance is reduced by 43% compared to an equal-mass aluminum design.

Current research surrounding FRAM focuses on the manufacturing process and neglects opportunities to leverage design freedom provided by FRAM. Previous FRAM optimization methods only optimize fiber orientation within a 2D plane and do not establish an optimized 3D print orientation, neglecting exploration of the entire orientation design space.

This study examines the effect of temperature-dependent material models for normal-strength (NSC) and high-strength concrete (HSC) on the thermal analysis of reinforced concrete (RC) walls.

The study performs an one-at-a-time (OAT) sensitivity analysis to assess the impact of variables defining the constitutive and parametric fire models on the wall's thermal response. Moreover, it extends the sensitivity analysis to a variance-based analysis to assess the effect of constitutive model type, fire model type and constitutive model uncertainty on the RC wall's thermal response variance. The study determines the wall’s thermal behaviour reliability considering the different constitutive models and their uncertainty.

It is found that the impact of the variability in concrete’s conductivity is determined by its temperature-dependent model, which differs for NSC and HSC. Therefore, more testing and improving material modelling are needed. Furthermore, the heating rate of the fire scenario is the dominant factor in deciding fire-resistance performance because it is a causal factor for spalling in HSC walls. And finally the reliability of wall's performance decreased sharply for HSC walls due to the expected spalling of the concrete and loss of cross-section integrity.

Limited studies in the current open literature quantified the impact of constitutive models on the behaviour of RC walls. No studies have examined the effect of material models' uncertainty on wall’s response reliability under fire. Furthermore, the study's results contribute to the ongoing attempts to shape performance-based structural fire engineering.

The need for precise synthesis of customized designs has resulted in the development of advanced coating processes for modern nanomaterials. Achieving accuracy in these processes requires a deep understanding of thermophysical behavior, rheology and complex chemical reactions. The manufacturing flow processes for these coatings are intricate and involve heat and mass transfer phenomena. Magnetic nanoparticles are being used to create intelligent coatings that can be externally manipulated, making them highly desirable. In this study, a Keller box calculation is used to investigate the flow of a coating nanofluid containing a viscoelastic polymer over a circular cylinder.

The rheology of the coating polymer nanofluid is described using the viscoelastic model, while the effects of nanoscale are accounted for by using Buongiorno’s two-component model. The nonlinear PDEs are transformed into dimensionless PDEs via a nonsimilar transformation. The dimensionless PDEs are then solved using the Keller box method.

The transport phenomena are analyzed through a comprehensive parametric study that investigates the effects of various emerging parameters, including thermal radiation, Biot number, Eckert number, Brownian motion, magnetic field and thermophoresis. The results of the numerical analysis, such as the physical variables and flow field, are presented graphically. The momentum boundary layer thickness of the viscoelastic polymer nanofluid decreases as fluid parameter increases. An increase in mixed convection parameter leads to a rise in the Nusselt number. The enhancement of the Brinkman number and Biot number results in an increase in the total entropy generation of the viscoelastic polymer nanofluid.

Intelligent materials rely heavily on the critical characteristic of viscoelasticity, which displays both viscous and elastic effects. Viscoelastic models provide a comprehensive framework for capturing a range of polymeric characteristics, such as stress relaxation, retardation, stretching and molecular reorientation. Consequently, they are a valuable tool in smart coating technologies, as well as in various applications like supercapacitor electrodes, solar collector receivers and power generation. This study has practical applications in the field of coating engineering components that use smart magnetic nanofluids. The results of this research can be used to analyze the dimensions of velocity profiles, heat and mass transfer, which are important factors in coating engineering. The study is a valuable contribution to the literature because it takes into account Joule heating, nonlinear convection and viscous dissipation effects, which have a significant impact on the thermofluid transport characteristics of the coating.

The momentum boundary layer thickness of the viscoelastic polymer nanofluid decreases as the fluid parameter increases. An increase in the mixed convection parameter leads to a rise in the Nusselt number. The enhancement of the Brinkman number and Biot number results in an increase in the total entropy generation of the viscoelastic polymer nanofluid. Increasing the strength of the magnetic field promotes an increase in the density of the streamlines. An increase in the mixed convection parameter results in a decrease in the isotherms and isoconcentration.

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.

In this paper, based on the density functional theory (DFT) and finite element method (FEM), the elastic, buckling and vibrational behaviors of the monolayer bismuthene are studied.

The computed elastic properties based on DFT are used to develop a finite element (FE) model for the monolayer bismuthene in which the Bi-Bi bonds are simulated by beam elements. Furthermore, mass elements are used to model the Bi atoms. The developed FE model is used to compute Young's modulus of monolayer bismuthene. The model is then used to evaluate the buckling force and fundamental natural frequency of the monolayer bismuthene with different geometrical parameters.

Comparing the results of the FEM and DFT, it is shown that the proposed model can predict Young's modulus of the monolayer bismuthene with an acceptable accuracy. It is also shown that the influence of the vertical side length on the fundamental natural frequency of the monolayer bismuthene is not significant. However, vibrational characteristics of the bismuthene are significantly affected by the horizontal side length.

DFT and FEM are used to study the elastic, vibrational and buckling properties of the monolayer bismuthene. The developed model can be used to predict Young's modulus of the monolayer bismuthene accurately. Effect of the vertical side length on the fundamental natural frequency is negligible. However, vibrational characteristics are significantly affected by the horizontal side length.