The Behavior of compressively loaded stiffener runout specimens - Part II: Finite element analysis

The termination of stiffeners in composite aircraft structures give rise to regions of high interlaminar shear and peel stresses as the load in the stiffener is diffused into the skin. This is of particular concern in co-cured composite stiffened structures where there is a relatively low resistance to through-thickness stress components at the skin-stiffener interface. In Part I, experimental results of tested specimens highlighted the influence of local design parameters on their structural response. Indeed some of the observed behavior was unexpected. There is a need to be able to analyse a range of changes in geometry rapidly to allow the analysis to form an integral part of the structural design process.

This work presents the development of a finite element methodology for modelling the failure process of these critical regions. An efficient thick shell element formulation is presented and this element is used in conjuction with the Virtual Crack Closure Technique (VCCT) to predict the crack growth characteristics of the modelled specimens. Three specimens were modelled and the qualitative aspects of crack growth were captured successfully. The shortcomings in the quantitative correlation between the predicted and observed failure loads are discussed. There was evidence to suggest that high through-thickness compressive stresses enhanced the fracture toughness in these critical regions.

Assessment of rolling resistance models in discrete element simulations

Particulate systems are of interest in many disciplines. They are often investigated using the discrete element method because of its capability to investigate particulate systems at the individual particle scale. To model the interaction between two particles and between a particle and a boundary, conventional discrete element models use springs and dampers in both the normal and tangential directions. The significance of particle rotation has been highlighted in both numerical studies and physical experiments. Several researchers have attempted to incorporate a rotational torque to account for the rolling resistance or rolling friction by developing different models. This paper presents a review of the commonly used models for rolling resistance and proposes a more general model. These models are classified into four categories according to their key characteristics. The robustness of these models in reproducing rolling resistance effects arising from different physical situations was assessed by using several benchmarking test cases. The proposed model can be seen to be more general and suitable for modelling problems involving both dynamic and pseudo-static regimes. An example simulation of the formation of a 2D sandpile is also shown. For simplicity, all formulations and examples are presented in 2D form, though the general conclusions are also applicable to 3D systems.

Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials

A numerical method is developed to simulate complex two-dimensional crack propagation in quasi-brittle materials considering random heterogeneous fracture properties. Potential cracks are represented by pre-inserted cohesive elements with tension and shear softening constitutive laws modelled by spatially varying Weibull random fields. Monte Carlo simulations of a concrete specimen under uni-axial tension were carried out with extensive investigation of the effects of important numerical algorithms and material properties on numerical efficiency and stability, crack propagation processes and load-carrying capacities. It was found that the homogeneous model led to incorrect crack patterns and load–displacement curves with strong mesh-dependence, whereas the heterogeneous model predicted realistic, complicated fracture processes and load-carrying capacity of little mesh-dependence. Increasing the variance of the tensile strength random fields with increased heterogeneity led to reduction in the mean peak load and increase in the standard deviation. The developed method provides a simple but effective tool for assessment of structural reliability and calculation of characteristic material strength for structural design.

The development of a surface defect machining method for hard turning processes

Hard turning (HT) is a material removal process employing a combination of a single point cutting tool and high speeds to machine hard ferrous alloys which exhibit hardness values over 45 HRC. In this paper, a surface defect machining (SDM) method for HT is proposed which harnesses the combined advantages of porosity machining and pulsed laser pre-treatment processing. From previous experimental work, this was shown to provide better controllability of the process and improved quality of the machined surface. While the experiments showed promising results, a comprehensive understanding of this new technique could only be achieved through a rigorous, in depth theoretical analysis. Therefore, an assessment of the SDM technique was carried out using both finite element method (FEM) and molecular dynamics (MD) simulations.
FEM modelling was used to compare the conventional HT of AISI 4340 steel (52 HRC) using an Al2O3 insert with the proposed SDM method. The simulations showed very good agreement with the previously published experimental results. Compared to conventional HT, SDM provided favourable machining outcomes, such as reduced shear plane angle, reduced average cutting forces, improved surface roughness, lower residual stresses on the machined surface, reduced tool–chip interface contact length and increased chip flow velocity. Furthermore, a scientific explanation of the improved surface finish was revealed using a state-of-the-art MD simulation model which suggested that during SDM, a combination of both the cutting action and rough polishing action help improve the machined surface finish.

Determination of Power System Response during Small Load Fluctuations

This paper proposes a calculation method to determine power system response during small load perturbations or minor disturbances. The method establishes the initial value of active power transient using traditional reduction technique on admittance matrix, which incorporates voltage variations in the determination. The method examines active power distribution among generators when several loads simultaneously change, and verifies that the superposition principle is applicable for this scenario. The theoretical derivation provided in the paper is validated by numerical simulations using a 3-generator 9-bus benchmark model. The results indicate that the inclusion of voltage variation renders an independent and precise measure of active power response during transient conditions.

Ultimate lateral load resistance of laterally loaded pile

Several methods are available for predicting ultimate lateral load resistance of laterally loaded pile. These methods often produce significantly different ultimate lateral resistance. This makes it difficult to select an appropriate method in designing/predicting ultimate lateral resistance of pile. This paper presents a review of two different methods; Meyerh of and Patra & Pise for predicting lateral resistance of pile. Then, the predicted ultimate lateral resistances by these two methods are compared with the experimental results. It is found that Meyerhof's method gives better prediction for single pile with smaller L/d ratio whereas Patra & Pise method gives better predictions for pile groups with higher L/d. Thus, none of these methods can be applicable universally for all possible conditions. Also the parametric study on ultimate lateral resistance revealed that length to diameter ratio, pile spacing, pile configuration in a pile group are important parameters for prediction of lateral load resistance. © 2009 Taylor & Francis Group.

Finite element analysis of interfacial friction and slip in composites with fully unbonded filler

A finite element model is developed to predict the stress-strain behaviour of particulate composites with fully unbonded filler particles. This condition can occur because of the lack of adhesion property of the filler surface. Whilst part of the filler particle is separated from the matrix, another section of filler keeps in contact with the matrix because of the lateral compressive displacement of the matrix. The slip boundary condition is imposed on the section of the interface that remains closed. The states of stress and displacement fields are obtained. The location of any further deformation through crazing or shear band formation is identified. A completely unbonded inclusion with partial slip at a section of the interface reduces the concentration of the stress at the interface significantly. Whereas this might lead to slightly higher strength, it decreases the load transfer efficiency and stiffness of this type of composite.

Evaluation of sources of uncertainties in microtensile bond strength of dental adhesive system for different specimen geometries

Objective: The aim of this research is to use finite element analysis (FEA) to quantify the effect of the sample shape and the imperfections induced during the manufacturing process of samples on the bond strength and modes of failure of dental adhesive systems through microtensile test. Using the FEA prediction for individual parameters effect, estimation of expected variation and spread of the microtensile bond strength results for different sample geometries is made. Methods: The estimated stress distributions for three different sample shapes, hourglass, stick and dumbbell predicted by FEA are used to predict the strength for different fracture modes. Parameters such as the adhesive thickness, uneven interface of the adhesive and composite and dentin, misalignment of axis of loading, the existence of flaws such as induced cracks during shaping the samples or bubbles created during application of the adhesive are considered. Microtensile experiments are performed simultaneously to measure bond strength and modes of failure. These are compared with the FEA results. Results: The relative bonding strength and its standard deviation for the specimens with different geometries measured through the microtensile tests confirm the findings of the FEA. The hourglass shape samples show lower tensile bond strength and standard deviation compared to the stick and dumbbell shape samples. ANOVA analysis confirms no significant difference between dumbbell and stick geometry results, and major differences of these two geometries compared to hourglass shape measured values. Induced flaws in the adhesive and misalignment of the angle of application of load have significant effect on the microtensile bond strength. Using adhesive with higher modulus the differences between the bond strength of the three sample geometries increase. Significance: The result of the research clarifies the importance of the sample geometry chosen in measuring the bond strength. It quantifies the effect of the imperfections on the bond strength for each of the sample geometries through a systematic and all embracing study. The results explain the reasons of the large spread of the microtensile test results reported by various researchers working in different labs and the need for standardization of the test method and sample shape used in evaluation of the dentin-adhesive bonding system. © 2007 Academy of Dental Materials.

Experimental and Computational Approach Investigating Burst Fracture Augmentation Using PMMA and Calcium Phosphate Cements

The aim of the study was to use a computational and experimental approach to evaluate, compare and predict the ability of calcium phosphate (CaP) and poly (methyl methacrylate) (PMMA) augmentation cements to restore mechanical stability to traumatically fractured vertebrae, following a vertebroplasty procedure. Traumatic fractures (n = 17) were generated in a series of porcine vertebrae using a drop-weight method. The fractured vertebrae were imaged using μCT and tested under axial compression. Twelve of the fractured vertebrae were randomly selected to undergo a vertebroplasty procedure using either a PMMA (n = 6) or a CaP cement variation (n = 6). The specimens were imaged using μCT and re-tested. Finite element models of the fractured and augmented vertebrae were generated from the μCT data and used to compare the effect of fracture void fill with augmented specimen stiffness. Significant increases (p <0.05) in failure load were found for both of the augmented specimen groups compared to the fractured group. The experimental and computational results indicated that neither the CaP cement nor PMMA cement could completely restore the vertebral mechanical behavior to the intact level. The effectiveness of the procedure appeared to be more influenced by the volume of fracture filled rather than by the mechanical properties of the cement itself.

Trace Element Abundance and Renal Disease

This research investigates the relationship between elevated trace elements in soils, stream sediments and stream water and the prevalence of Chronic Kidney Disease (CKD). The study uses a collaboration of datasets provided from the UK Renal Registry Report (UKRR) on patients with renal diseases requiring treatment including Renal Replacement Therapy (RRT), the soil geochemical dataset for Northern Ireland provided by the Tellus Survey, Geological Survey of Northern Ireland (GSNI) and the bioaccessibility of Potentially Toxic Elements (PTEs) from soil samples which were obtained from the Unified Barge Method (UBM). The relationship between these factors derives from the UKRR report which highlights incidence rates of renal impaired patients showing regional variation with cases of unknown aetiology. Studies suggest a potential cause of the large variation and uncertain aetiology is associated with underlying environmental factors such as the oral bioaccessibility of trace elements in the gastrointestinal tract.
As previous research indicates that long term exposure is related to environmental factors, Northern Ireland is ideally placed for this research as people traditionally live in the same location for long periods of time. Exploratory data analysis and multivariate analyses are used to examine the soil, stream sediments and stream water geochemistry data for a range of key elements including arsenic, lead, cadmium and mercury identified from a review of previous renal disease literature. The spatial prevalence of patients with long term CKD is analysed on an area basis. Further work includes cluster analysis to detect areas of low or high incidences of CKD that are significantly correlated in space, Geographical Weighted Regression (GWR) and Poisson kriging to examine locally varying relationship between elevated concentrations of PTEs and the prevalence of CKD.

Study of localized damage in composite laminates using micro-macro approach

A robust multiscale scheme referred to as micro–macro method has been developed for the prediction of localized damage in fiber reinforced composites and implemented in a finite element framework. The micro–macro method is based on the idea of partial homogenization of a structure. In this method, the microstructural details are included in a small region of interest in the structure and the rest is modeled as a homogeneous continuum. The solution to the microstructural fields is then obtained on solving the two different domains, simultaneously. This method accurately predicts local stress fields in stress concentration regions and is computationally efficient as compared with the solution of a full scale microstructural model. This scheme has been applied to obtain localized damage at high and low stress zones of a V-notched rail shear specimen. The prominent damage mechanisms under shear loading, namely, matrix cracking and interfacial debonding, have been modeled using Mohr–Coulomb plasticity and traction separation law, respectively. The average stress at the notch has been found to be 44% higher than the average stresses away from the notch for a 90 N shear load. This stress rise is a direct outcome of the geometry of the notch.

Transonic Nonlinear Aeroelastic Simulations Using An Harmonic Balance Method

This paper presents an approach to compute transonic Limit Cycle O
scillations using a coupled Harmonic Balance formulation based on the Euler equations for fluid dynamics and finite element models. The paper will investigate the role of aerodynamic (shocks) and structural nonlinearities in driving the limit cycle behaviour. Part icular attention will be given to nonlinear interactions for subcritical LCOs. The Aero elastic Harmonic Balance formulation, allows for solutions of the coupled structural dynamics and CFD system at a reduced cost.

A new method for transmission loss allocation considering the circulating currents between generators

A new method is presented for transmission loss allocation based on the separation of transmission loss caused by load and the loss due to circulating currents between generators. The theoretical basis for and derivation of the loss formulae are presented using simple systems. The concept is then extended to a general power system using the Ybus model. Details of the application of the proposed method to a typical power system are presented along with results from the IEEE 30 bus test system. The results from both the small system and the standard IEEE test system demonstrate the validity of the proposed method.

Finite element modeling of debonding failures in FRP-strengthened RC beams: A dynamic approach

This paper is concerned with the finite element simulation of debonding failures in FRP-strengthened concrete beams. A key challenge for such simulations is that common solution techniques such as the Newton-Raphson method and the arc-length method often fail to converge. This paper examines the effectiveness of using a dynamic analysis approach in such FE simulations, in which debonding failure is treated as a dynamic problem and solved using an appropriate time integration method. Numerical results are presented to show that an appropriate dynamic approach effectively overcomes the convergence problem and provides accurate predictions of test results.

Experimental and discrete element modelling of cohesive iron ore fines

In this study, the behaviour of iron ore fines with varying levels of adhesion was investigated using a confined compression test and a uniaxial test. The uniaxial test was conducted using the semi-automated uniaxial EPT tester in which the cohesive strength of a bulk solid is evaluated from an unconfined compression test following a period of consolidation to a pre-defined vertical stress. The iron ore fines were also tested by measuring both the vertical and circumferential strains on the cylindrical container walls under vertical loading in a separate confined compression tester - the K0 tester, to determine the lateral pressure ratio. Discrete Element Method simulations of both experiments were carried out and the predictions were compared with the experimental observations. A recently developed DEM contact model for cohesive solids, an Elasto-Plastic Adhesive model, was used. This particle contact model uses hysteretic non-linear loading and unloading paths and an adhesion parameter which is a function of the maximum contact overlap. The model parameters for the simulations are phenomenologically based to reproduce the key bulk characteristics exhibited by the solid. The simulation results show a good agreement in capturing the stress history dependent behaviour depicted by the flow function of the cohesive iron ore fines while also providing a reasonably good match for the lateral pressure ratio observed during the confined compression K0 tests. This demonstrates the potential for the DEM model to be used in the simulation of bulk handling applications.