A thermo-elastic contact problem of functionally graded materials (FGMs) rotating brake disk with different pure brake pad areas under temperature dependent material properties is solved by Finite Element Method (FEM). The properties of brake disk change gradually from metal to ceramic by power-law distribution along the radial direction from the inner to the outer surface. Areas of the pure pad are changing while the vertical force is constant. The ratio of brake pad thickness to FGMs brake disk thickness is assumed 0.66. Two sources of thermal loads are considered: (1) Heat generation between the pad and brake disk due to contact friction, and (2) External thermal load due to a constant temperature at inner and outer surfaces. Mechanical responses of FGMs disk are compared with several pad contact areas. The results for temperature-dependent and temperature-independent material properties are investigated and presented. The results show that the absolute value of the shear stress in temperature-dependent material can be greater than that for temperature-independent material. The radial stress for some specific grading index (n = 1.5) is compressive near the inner surface for double contact while it is tensile for a single contact. It is concluded that the radial strain for some specific value of grading index (n = 1) is lower than other FGMs and pure double side contact brake disks.
Nowadays, development of cable-stayed bridges is increasing around the world. The mitigation of seismic forces to these bridges are obligatory to prevent damages or failure of its structural members. Herein, this paper aimed to determine the near-fault ground motion effect on an existing cablestayed bridge equipped with lead-rubber bearing. In this context, Shipshaw cable-stayed bridge is selected as the case study. The selected bridge has a span of 183.2 m composite deck and 43 m height of steel tower. 2D finite element models of the non-isolated and base isolated bridges are modelled by using SAP2000. Three different near-fault ground motions which are Tabas 1978, Cape Mendocino 1992 and Kobe 1995 were subjected to the 2D FEM models in order to determine the seismic behaviour of the bridge. The near-fault ground motions were applied to the bridge in the longitudinal direction. Nonlinear dynamic analysis was performed to determine the dynamic responses of the bridge. Comparison of dynamic response of nonisolated and base isolated bridge under three different near-fault ground motions were conducted. The results obtained from numerical analyses of the bridge showed that the isolation system lengthened the period of bridge and minimised deck displacement, base shear and base moment of the bridge. It is concluded that the isolation system significantly reduced the destructive effects of near-fault ground motions on the bridge.
Invention of milling combined laser sintering system (MLSS) is able to reduce the mould manufacturing time and improve the mould accuracy. Thus, more study is needed to increase the understanding for the laser sintered material machining characteristic to gain benefit from the invention of MLSS. This paper clarified the analysis of laser sintered material machinability with the application of Finite Element Method (FEM). Mild steel AISI1055 was applied in developing the Finite Element model in this study due to its popularity in machinability test and adequate level of data availability. 2D orthogonal cutting was employed on edge design tools with updated Lagrangian coupled thermo mechanical plane strain model. Adaptive meshing, tool edge radius and various types of friction models were assigned to obtain efficient simulations and precise cutting results. Cutting force and cutting-edge temperature estimated by Finite Element Method are validated against corresponding experimental values by previous researchers. In the study, cutting force increases when radial depth increases and lowest error acquired when the shear friction factor of 0.8 was applied. Machining simulation for laser sintered materials estimated lower cutting force compared with mild steel AISI1055 due to lower Young modulus. Higher cutting temperature estimated for machining simulation laser sintered material compared with machining simulation mild steel AISI1055 due to its low thermal conductivity.
The tensile strength of intact rock materials has been determined by indirect method more frequent than the direct method.
The most commonly used indirect method is Brazilian test. Stress and deformability undergo during the test reflected by
geometry shape of the samples with respect to the different diameter ratio. This study focuses on influence of geometry
shape in solid and ring disk with different diameter ratio on the stress distribution and deformations within sandstone
subjected to indirect tensile loading by Brazilian test. Then, the finite element method in RS2 software was utilised to
simulate and gain in depth understanding the behaviour of Brazilian test. The analysis shown that the maximum tensile
strength in a ring disk with diameter ratio of 0.1 is three times higher than in solid disk. Meanwhile, as the diameter
ratio of ring disk increases, it produces lower tensile strength. The numerical simulation also has successfully illustrated
the shear failure which observed near the loading platen of solid disk during Brazilian test. The finite element analysis
utilised in this research has successfully enables the stress distribution and deformation behaviour of the rock under
tension to be studied closely
A large number of shallow buried tunnels are built in the city nowadays and the special strata such as large upper-soft and lower-hard ground often encountered. Deformation control of strata is the focus issue related to the construction safety. Based on Dalian metro Hing Street station with the classical geological condition of upper-soft and lower-hard ground, this paper fully used a combined control method including six different support measures to control the deformation of surrounding rock. 3D finite element model was setup to analyze the construction effect of combined control measures and the monitoring in-site was carried out to verify the deformation control effect of combined control method. It shows that the maximum surface subsidence value is gradually reduced with the support measures gradually increasing. In the case of various supports the maximum sedimentation value is 2.67 cm, which is 42. 1% lower than that of not using control method and the control effect is obvious. In addition, it can be seen that the two-layer initial support and additional large arch foot have the best effect on controlling the ground surface settlement with reduction of 11.7% and 20.2%, respectively. The research results can provide practical experience for the construction of such tunnels, and guide the design and construction of the tunnel in the future.
The development degree of fissure water in underground rock is a great trouble to the construction of railway tunnel, which will cause a series of environmental geological problems. Take the surrounding rock-section of the typical red clay in Lvliang-Mt. railway tunnel below the underground water level as an example, several aspects about the red clay surrounding rock will be researched, including pore water pressure, volume moisture content, stress of surrounding rock, vault subsidence and horizontal convergence through the field monitoring. Taking into account the importance of railway tunnel engineering, the large shear test of red clay was carried out at the construction site specially and the reliable situ shear strength parameters of surrounding rock will be obtained. These investigations and field tests helped to do a series of work: Three dimensional finite element numerical model of railway tunnel will be established, the deformation law of the red clay surrounding rock will be investigated, respectively, for the water-stress coupling effect and without considering it, the variation of the pore water pressure during excavation, the influence degree about the displacement field and stress field of water-stress coupling on red clay-rock will be discussed and the mechanism of the surrounding rock deformation will be submitted. Finally, the paper puts forward the feasible drainage scheme of the surrounding rock and the tunnel cathode. The geological environment safety of tunnel construction is effectively protected.
This study aimed to develop a three-dimensional finite element model of a functionally graded femoral prosthesis. The model consisted of a femoral prosthesis created from functionally graded materials (FGMs), cement, and femur. The hip prosthesis was composed of FGMs made of titanium alloy, chrome-cobalt, and hydroxyapatite at volume fraction gradient exponents of 0, 1, and 5, respectively. The stress was measured on the femoral prosthesis, cement, and femur. Stress on the neck of the femoral prosthesis was not sensitive to the properties of the constituent material. However, stress on the stem and cement decreased proportionally as the volume fraction gradient exponent of the FGM increased. Meanwhile, stress became uniform on the cement mantle layer. In addition, stress on the femur in the proximal part increased and a high surface area of the femoral part was involved in absorbing the stress. As such, the stress-shielding area decreased. The results obtained in this study are significant in the design and longevity of new prosthetic devices because FGMs offer the potential to achieve stress distribution that more closely resembles that of the natural bone in the femur.
The mechanical properties of tissue scaffolds are essential in providing stability for tissue repair and growth. Thus, the ability of scaffolds to withstand specific loads is crucial for scaffold design. Most research on scaffold pores focuses on grids with pore size and gradient structure, and many research models are based on scaffolding with vertically arranged holes. However, little attention is paid to the influence of the distribution of holes on the mechanical properties of the scaffold. To address this gap, this research investigates the effect of pore distribution on the mechanical properties of tissue scaffolds. The study involves four types of scaffold designs with regular and staggered pore arrangements and porosity ranging from 30% to 80%. Finite element analysis (FEA) was used to compare the mechanical properties of different scaffold designs, with von-Mises stress distribution maps generated for each scaffold. The results show that scaffolds with regular vertical holes exhibit a more uniform stress distribution and better mechanical performance than those with irregular holes. In contrast, the scaffold with a staggered arrangement of holes had a higher probability of stress concentration. The study emphasized the importance of balancing porosity and strength in scaffold design.
This paper investigated the static behaviour of glass fibre reinforced polymer (GFRP) built-up hollow and concrete filled built-up beams tested under four-point bending with a span-to-depth ratio of 1.67, therefore focusing their shear performance. Two parameters considered for hollow sections were longitudinal web stiffener and strengthening at the web-flange junction. The experimental results indicated that the GFRP hollow beams failed by web crushing at supports; therefore, the longitudinal web stiffener has an insignificant effect on improving the maximum load. Strengthening web-flange junctions using rectangular hollow sections increased the maximum load by 47%. Concrete infill could effectively prevent the web crushing, and it demonstrated the highest load increment of 162%. The concrete filled GFRP composite beam failed by diagonal tension in the lightweight concrete core. The finite element models adopting Hashin damage criteria yielded are in good agreement with the experimental results in terms of maximum load and failure mode. Based on the numerical study, the longitudinal web stiffener could prevent the web buckling of the slender GFRP beam and improved the maximum load by 136%. The maximum load may be further improved by increasing the thickness of the GFRP section and the size of rectangular hollow sections used for strengthening. It was found that the bond-slip at the concrete-GFRP interface affected the shear resistance of concrete-GFRP composite beam.
The mechanical behavior of graphene/polymer interfaces in the graphene-reinforced epoxy nanocomposite is one of the factors that dictates the deformation and damage response of the nanocomposites. In this study, hybrid molecular dynamic (MD) and finite element (FE) simulations of a graphene/polymer nanocomposite are developed to characterize the elastic-damage behavior of graphene/polymer interfaces under a tensile separation condition. The MD results show that the graphene/epoxy interface behaves in the form of elastic-softening exponential regressive law. The FE results verify the adequacy of the cohesive zone model in accurate prediction of the interface damage behavior. The graphene/epoxy cohesive interface is characterized by normal stiffness, tensile strength, and fracture energy of 5 × 10-8 (aPa·nm-1), 9.75 × 10-10 (nm), 2.1 × 10-10 (N·nm-1) respectively, that is followed by an exponential regressive law with the exponent, α = 7.74. It is shown that the commonly assumed bilinear softening law of the cohesive interface could lead up to 55% error in the predicted separation of the interface.
The main purpose of this research is to design a high-fatigue performance hoop wrapped compressed natural gas (CNG) composite cylinder. To this end, an optimization algorithm was presented as a combination of finite element simulation (FES) and response surface analysis (RSA). The geometrical model was prepared as a variable wall-thickness following the experimental measurements. Next, transient dynamic analysis was performed subjected to the refueling process, including the minimum and maximum internal pressures of 20 and 200 bar, respectively. The time histories of stress tensor components were extracted in the critical region. Furthermore, RSA was utilized to investigate the interaction effects of various polymer composite shell manufacturing process parameters (thickness and fiber angle) on the fatigue life of polymer composite CNG pressure tank (type-4). In the optimization procedure, four parameters including wall-thickness of the composite shell in three different sections of the CNG tank and fiber angle were considered as input variables. In addition, the maximum principal stress of the component was considered as the objective function. Eventually, the fatigue life of the polymer composite tank was calculated using stress-based failure criterion. The results indicated that the proposed new design (applying optimal parameters) leads to improve the fatigue life of the polymer composite tank with polyethylene liner about 2.4 times in comparison with the initial design.
This paper aims to propose a temperature-dependent cohesive model to predict the delamination of dissimilar metal-composite material hybrid under Mode-I and Mode-II delamination. Commercial nonlinear finite element (FE) code LS-DYNA was used to simulate the material and cohesive model of hybrid aluminium-glass fibre-reinforced polymer (GFRP) laminate. For an accurate representation of the Mode-I and Mode-II delamination between aluminium and GFRP laminates, cohesive zone modelling with bilinear traction separation law was implemented. Cohesive zone properties at different temperatures were obtained by applying trends of experimental results from double cantilever beam and end notched flexural tests. Results from experimental tests were compared with simulation results at 30, 70 and 110 °C to verify the validity of the model. Mode-I and Mode-II FE models compared to experimental tests show a good correlation of 5.73% and 7.26% discrepancy, respectively. Crack front stress distribution at 30 °C is characterised by a smooth gradual decrease in Mode-I stress from the centre to the edge of the specimen. At 70 °C, the entire crack front reaches the maximum Mode-I stress with the exception of much lower stress build-up at the specimen's edge. On the other hand, the Mode-II stress increases progressively from the centre to the edge at 30 °C. At 70 °C, uniform low stress is built up along the crack front with the exception of significantly higher stress concentrated only at the free edge. At 110 °C, the stress distribution for both modes transforms back to the similar profile, as observed in the 30 °C case.
The stiffness response or load-deformation/displacement behavior is the most important mechanical behavior that frequently being utilized for validation of the mathematical-physical models representing the mechanical behavior of solid objects in numerical method, compared to actual experimental data. This numerical study aims to investigate the linear-nonlinear stiffness behavior of carbon fiber-reinforced polymer (CFRP) composites at material and structural levels, and its dependency to the sets of individual/group elastic and damage model parameters. In this regard, a validated constitutive damage model, elastic-damage properties as reference data, and simulation process, that account for elastic, yielding, and damage evolution, are considered in the finite element model development process. The linear-nonlinear stiffness responses of four cases are examined, including a unidirectional CFRP composite laminate (material level) under tensile load, and also three multidirectional composite structures under flexural loads. The result indicated a direct dependency of the stiffness response at the material level to the elastic properties. However, the stiffness behavior of the composite structures depends both on the structural configuration, geometry, lay-ups as well as the mechanical properties of the CFRP composite. The value of maximum reaction force and displacement of the composite structures, as well as the nonlinear response of the structures are highly dependent not only to the mechanical properties, but also to the geometry and the configuration of the structures.
Composite structures are made of multidirectional (MD) fiber-reinforced polymer (FRP) composite laminates, which fail due to multiple damages in matrix, interface, and fiber constituents at different scales. The yield point of a unidirectional FRP composite is assumed as the lamina strength limit representing the damage initiation phenomena, while yielding of MD composites in structural applications are not quantified due to the complexity of the sequence of damage evolutions in different laminas dependent on their angle and specification. This paper proposes a new method to identify the yield point of MD composite structures based on the evolution of the damage dissipation energy (DDE). Such a characteristic evolution curve is computed using a validated finite element model with a mesoscale damage-based constitutive model that accounts for different matrix and fiber failure modes in angle lamina. The yield point of composite structures is identified to correspond to a 5% increase in the initial slope of the DDE evolution curve. The yield points of three antisymmetric MD FRP composite structures under flexural loading conditions are established based on Hashin unidirectional (UD) criteria and the energy-based criterion. It is shown that the new energy concept provides a significantly larger safe limit of yield for MD composite structures compared to UD criteria, in which the accumulation of energy dissipated due to all damage modes is less than 5% of the fracture energy required for the structural rupture.
The determination of the myocardium's tissue properties is important in constructing functional finite element (FE) models of the human heart. To obtain accurate properties especially for functional modeling of a heart, tissue properties have to be determined in vivo. At present, there are only few in vivo methods that can be applied to characterize the internal myocardium tissue mechanics. This work introduced and evaluated an FE inverse method to determine the myocardial tissue compressibility. Specifically, it combined an inverse FE method with the experimentally-measured left ventricular (LV) internal cavity pressure and volume versus time curves. Results indicated that the FE inverse method showed good correlation between LV repolarization and the variations in the myocardium tissue bulk modulus K (K = 1/compressibility), as well as provided an ability to describe in vivo human myocardium material behavior. The myocardium bulk modulus can be effectively used as a diagnostic tool of the heart ejection fraction. The model developed is proved to be robust and efficient. It offers a new perspective and means to the study of living-myocardium tissue properties, as it shows the variation of the bulk modulus throughout the cardiac cycle.
Bio-composites of oil palm empty fruit bunch (OPEFB) fibres and polycaprolactones (PCL) with a thickness of 1 mm were prepared and characterized. The composites produced from these materials are low in density, inexpensive, environmentally friendly, and possess good dielectric characteristics. The magnitudes of the reflection and transmission coefficients of OPEFB fibre-reinforced PCL composites with different percentages of filler were measured using a rectangular waveguide in conjunction with a microwave vector network analyzer (VNA) in the X-band frequency range. In contrast to the effective medium theory, which states that polymer-based composites with a high dielectric constant can be obtained by doping a filler with a high dielectric constant into a host material with a low dielectric constant, this paper demonstrates that the use of a low filler percentage (12.2%OPEFB) and a high matrix percentage (87.8%PCL) provides excellent results for the dielectric constant and loss factor, whereas 63.8% filler material with 36.2% host material results in lower values for both the dielectric constant and loss factor. The open-ended probe technique (OEC), connected with the Agilent vector network analyzer (VNA), is used to determine the dielectric properties of the materials under investigation. The comparative approach indicates that the mean relative error of FEM is smaller than that of NRW in terms of the corresponding S21 magnitude. The present calculation of the matrix/filler percentages endorses the exact amounts of substrate utilized in various physics applications.
Deep beams are commonly used in tall buildings, offshore structures, and foundations. According to many codes and standards, strut-and-tie model (STM) is recommended as a rational approach for deep beam analyses. This research focuses on the STM recommended by ACI 318-11 and AASHTO LRFD and uses experimental results to modify the strut effectiveness factor in STM for reinforced concrete (RC) deep beams. This study aims to refine STM through the strut effectiveness factor and increase result accuracy. Six RC deep beams with different shear span to effective-depth ratios (a/d) of 0.75, 1.00, 1.25, 1.50, 1.75, and 2.00 were experimentally tested under a four-point bending set-up. The ultimate shear strength of deep beams obtained from non-linear finite element modeling and STM recommended by ACI 318-11 as well as AASHTO LRFD (2012) were compared with the experimental results. An empirical equation was proposed to modify the principal tensile strain value in the bottle-shaped strut of deep beams. The equation of the strut effectiveness factor from AASHTTO LRFD was then modified through the aforementioned empirical equation. An investigation on the failure mode and crack propagation in RC deep beams subjected to load was also conducted.
This paper proposes an alternative approach to extract transformer's winding parameters of resistance (R), inductance (L), capacitance (C) and conductance (G) based on Finite Element Method (FEM). The capacitance and conductance were computed based on Fast Multiple Method (FMM) and Method of Moment (MoM) through quasi-electrostatics approach. The AC resistances and inductances were computed based on MoM through quasi-magnetostatics approach. Maxwell's equations were used to compute the DC resistances and inductances. Based on the FEM computed parameters, the frequency response of the winding was obtained through the Bode plot function. The simulated frequency response by FEM model was compared with the simulated frequency response based on the Multi-conductor Transmission Line (MTL) model and the measured frequency response of a 33/11 kV, 30 MVA transformer. The statistical indices such as Root Mean Square Error (RMSE) and Absolute Sum of Logarithmic Error (ASLE) were used to analyze the performance of the proposed FEM model. It is found that the simulated frequency response by FEM model is quite close to measured frequency response at low and mid frequency regions as compared to simulated frequency response by MTL model based on RMSE and ASLE analysis.
Effects of different surface textures on the interface shear strength, interface slip, and failure modes of the concrete-to-concrete bond are examined through finite element numerical model and experimental methods in the presence of the horizontal load with 'push-off' technique under different normal stresses. Three different surface textures are considered; smooth, indented, and transversely roughened to finish the top surfaces of the concrete bases. In the three-dimensional modeling via the ABAQUS solver, the Cohesive Zone Model (CZM) is used to simulate the interface shear failure. It is observed that the interface shear strength increases with the applied normal stress. The transversely roughened surface achieves the highest interface shear strength compared with those finished with the indented and smooth approaches. The smooth and indented surfaces are controlled by the adhesive failure mode while the transversely roughened surface is dominated by the cohesive failure mode. Also, it is observed that the CZM approach can accurately model the interface shear failure with 3-29% differences between the modeled and the experimental test findings.
This paper investigates the behaviour of C-shaped and L-shaped angle shear connectors embedded in solid concrete slabs. An effective finite element model is proposed to simulate the push out tests of these shear connectors that encompass nonlinear material behaviour, large displacement and damage plasticity. The finite element models are validated against test results. Parametric studies using this nonlinear model are performed to investigate the variations in concrete strength and connector dimensions. The finite element analyses also confirm the test results that increasing the length of shear connector increases their shear strength proportionately. It is observed that the maximum stress in L-shaped angle connectors takes place in the weld attachment to the beam, whereas in the C-shaped angle connectors, it is in the attached leg. The location of maximum concrete compressive damage is rendered in each case. Finally, a new equation for prediction of the shear capacity of C-shaped angle connectors is proposed.