An application of bi-directional evolutionary structural optimisation for optimising energy absorbing structures using a material damage model

The Bi-directional Evolutionary Structural Optimisation (BESO) method is a numerical topology optimisation method developed for use in finite element analysis. This paper presents a particular application of the BESO method to optimise the energy absorbing capability of metallic structures. The optimisation objective is to evolve a structural geometry of minimum mass while ensuring that the kinetic energy of an impacting projectile is reduced to a level which prevents perforation. Individual elements in a finite element mesh are deleted when a prescribed damage criterion is exceeded. An energy absorbing structure subjected to projectile impact will fail once the level of damage results in a critical perforation size. It is therefore necessary to constrain an optimisation algorithm from producing such candidate solutions. An algorithm to detect perforation was implemented within a BESO framework which incorporated a ductile material damage model.

Obtaining an intermediate band photovoltaic material through the Bi insertion in CdTe

Defect interaction can take place in CdTe under Te and Bi rich conditions. We demonstrate in this work through first principles calculations, that this phenomenon allows a Jahn Teller distortion to form an isolated half-filled intermediate band in the host semiconductor band-gap. This delocalized energy band supports the experimental deep level reported in the host band-gap of CdTe at a low bismuth concentration. Furthermore, the calculated optical absorption of CdTe:Bi in this work shows a significant subband-gap absorption that also supports the enhancement of the optical absorption found in the previous experimental results.

Numerical modelling and analysis of the blast performance of laminated glass panels and the influence of material parameters

This paper presents a comprehensive numerical procedure to treat the blast response of laminated glass (LG) panels and studies the influence of important material parameters. Post-crack behaviour of the LG panel and the contribution of the interlayer towards blast resistance are treated. Modelling techniques are validated by comparing with existing experimental results. Findings indicate that the tensile strength of glass considerably influences the blast response of LG panels while the interlayer material properties have a major impact on the response under higher blast loads. Initially, glass panes absorb most of the blast energy, but after the glass breaks, interlayer deforms further and absorbs most of the blast energy. LG panels should be designed to fail by tearing of the interlayer rather than failure at the supports to achieve a desired level of protection. From this aspect, material properties of glass, interlayer and sealant joints play important roles, but unfortunately they are not accounted for in the current design standards. The new information generated in this paper will enhance the capabilities of engineers to better design LG panels under blast loads and use better materials to improve the blast response of LG panels.

Study of effects of dopants on structure of vitreous semiconductors (GeSe3.5)100-xMx (M = Bi, Sb) using high pressure techniques

An investigation of the problem of controlled doping of amorphous chalcogenide semiconductors utilizing a Bridgman anvil high pressure technique, has been undertaken. Bulk amorphous semiconducting materials (GeSe3.5)100-x doped with M = Bi (x = 2, 4, 10) and M = Sb (x = 10) respectively are studied up to a pressure of 100 kbar down to liquid nitrogen temperature, with a view to observe the impurity induced modifications. Measurement of the electrical conductivity of the doped samples under quasi-hydrostatic pressure reveals that the pressure induced effects in lightly doped (2 at % Bi) and heavily doped (x = 4, 10) semiconductors are markedly different. The pressure effects in Sb-doped semiconductors are quite different from those in Bi-doped material.

Interatomic potentials for fusion reactor material simulations

Fusion energy is a clean and safe solution for the intricate question of how to produce non-polluting and sustainable energy for the constantly growing population. The fusion process does not result in any harmful waste or green-house gases, since small amounts of helium is the only bi-product that is produced when using the hydrogen isotopes deuterium and tritium as fuel. Moreover, deuterium is abundant in seawater and tritium can be bred from lithium, a common metal in the Earth's crust, rendering the fuel reservoirs practically bottomless. Due to its enormous mass, the Sun has been able to utilize fusion as its main energy source ever since it was born. But here on Earth, we must find other means to achieve the same. Inertial fusion involving powerful lasers and thermonuclear fusion employing extreme temperatures are examples of successful methods. However, these have yet to produce more energy than they consume. In thermonuclear fusion, the fuel is held inside a tokamak, which is a doughnut-shaped chamber with strong magnets wrapped around it. Once the fuel is heated up, it is controlled with the help of these magnets, since the required temperatures (over 100 million degrees C) will separate the electrons from the nuclei, forming a plasma. Once the fusion reactions occur, excess binding energy is released as energetic neutrons, which are absorbed in water in order to produce steam that runs turbines. Keeping the power losses from the plasma low, thus allowing for a high number of reactions, is a challenge. Another challenge is related to the reactor materials, since the confinement of the plasma particles is not perfect, resulting in particle bombardment of the reactor walls and structures. Material erosion and activation as well as plasma contamination are expected. Adding to this, the high energy neutrons will cause radiation damage in the materials, causing, for instance, swelling and embrittlement. In this thesis, the behaviour of a material situated in a fusion reactor was studied using molecular dynamics simulations. Simulations of processes in the next generation fusion reactor ITER include the reactor materials beryllium, carbon and tungsten as well as the plasma hydrogen isotopes. This means that interaction models, {\it i.e. interatomic potentials}, for this complicated quaternary system are needed. The task of finding such potentials is nonetheless nearly at its end, since models for the beryllium-carbon-hydrogen interactions were constructed in this thesis and as a continuation of that work, a beryllium-tungsten model is under development. These potentials are combinable with the earlier tungsten-carbon-hydrogen ones. The potentials were used to explain the chemical sputtering of beryllium due to deuterium plasma exposure. During experiments, a large fraction of the sputtered beryllium atoms were observed to be released as BeD molecules, and the simulations identified the swift chemical sputtering mechanism, previously not believed to be important in metals, as the underlying mechanism. Radiation damage in the reactor structural materials vanadium, iron and iron chromium, as well as in the wall material tungsten and the mixed alloy tungsten carbide, was also studied in this thesis. Interatomic potentials for vanadium, tungsten and iron were modified to be better suited for simulating collision cascades that are formed during particle irradiation, and the potential features affecting the resulting primary damage were identified. Including the often neglected electronic effects in the simulations was also shown to have an impact on the damage. With proper tuning of the electron-phonon interaction strength, experimentally measured quantities related to ion-beam mixing in iron could be reproduced. The damage in tungsten carbide alloys showed elemental asymmetry, as the major part of the damage consisted of carbon defects. On the other hand, modelling the damage in the iron chromium alloy, essentially representing steel, showed that small additions of chromium do not noticeably affect the primary damage in iron. Since a complete assessment of the response of a material in a future full-scale fusion reactor is not achievable using only experimental techniques, molecular dynamics simulations are of vital help. This thesis has not only provided insight into complicated reactor processes and improved current methods, but also offered tools for further simulations. It is therefore an important step towards making fusion energy more than a future goal.

Behaviour of Bi-adhesive Joints

The use of relatively low modulus adhesive at the ends of overlap in a bi-adhesive bondline of a bonded joint can reduce the stress concentration significantly and, therefore, potentially lead to higher strength of the joint. This study presents the two-dimensional and three-dimensional nonlinear (geometric and material) finite element analyses of adhesively bonded single lap joints having modulus-graded bondline under monotonic loading conditions. The adhesives were modelled as an elasto-plastic multi-linear material, while the substrates were regarded as both linear elastic and bi-linear elasto-plastic material. The computational simulations have been performed to investigate the bondline behaviour by studying the stress and strain distributions both at the mid-plane as well as at the interface of the bondline. It has been observed that the static strength is higher for joints with bi-adhesive bondlines compared to those with single adhesives in bondline. Higher joint strength has also been observed for optimum bi-adhesive bondline ratio through parametric studies. Effects of load level, and bondline thickness on stress distribution in the bi-adhesive bondline have also been studied. 3D analysis results reveal the existence of complex multi-axial stress/strain state at the ends of the overlap in the bondline which cannot be observed in 2D plane strain analysis. About 1/3rd of the width of the joint from the free edge in the width direction has 3D stress state, especially in the compliant adhesive of the bondline. Magnitudes of longitudinal and lateral stress/strain components are comparable to peel stress/strain components. It has also been analytically shown that the in-plane global stiffness of the joint remains unaffected by modulus gradation of the bondline adhesive. (C) Koninklijke Brill NV, Leiden, 2010.

Effect of prestrain on fatigue crack propagation in mild steel containing non-metallic inclusions

The effect of tensile prestrain on fatigue crack propagation behaviour of commercial mild steel with significant amount of stringer inclusions has been studied. In prestrained materials the usual stable stage II crack growth region is preceded by a phase wherein a retardation in crack growth rate occurs. No such behaviour is observed in annealed material. The amount of retardation is found to increase with increase in prestrain. A mechanism for the observed retardation in crack growth rate is also presented.

Effects of temper level on the dependence of fatigue crack growth threshold and crack closure on the prior austenitic grain size

An attempt has been made to systematically investigate the effects of microstructural parameters, such as the prior austenite grain size (PAGS), in influencing the resistance to fatigue crack growth (FCG) in the near-threshold region under three different temper levels in a quenched and tempered high-strength steel. By austenitizing at various temperatures, the PAGS was varied from about 0.7 to 96 μm. The microstructures with these grain sizes were tempered at 200 °C, 400 °C, and 530 °C and tested for fatigue thresholds and crack closure. It has been found that, in general, three different trends in the dependence of both the total threshold stress intensity range, ΔK th , and the intrinsic threshold stress intensity range, ΔK eff, th , on the PAGS are observable. By considering in detail the factors such as cyclic stress-strain behavior, environmental effects on FCG, and embrittlement during tempering, the present observations could be rationalized. The strong dependence of ΔK th and ΔK eff, th on PAGS in microstructures tempered at 530 °C has been primarily attributed to cyclic softening and thereby the strong interaction of the crack tip deformation field with the grain boundary. On the other hand, a less strong dependence of ΔK th and ΔK eff, th on PAGS is suggested to be caused by the cyclic hardening behavior of lightly tempered microstructures occurring in 200 °C temper. In both microstructures, crack closure influenced near-threshold FCG (NTFCG) to a significant extent, and its magnitude was large at large grain sizes. Microstructures tempered at the intermediate temperatures failed to show a systematic variation of ΔKth and ΔKeff, th with PAGS. The mechanisms of intergranular fracture vary between grain sizes in this temper. A transition from “microstructure-sensitive” to “microstructure-insensitive” crack growth has been found to occur when the zone of cyclic deformation at the crack tip becomes more or less equal to PAGS. Detailed observations on fracture morphology and crack paths corroborate the grain size effects on fatigue thresholds and crack closure.

Dynamic growth of tensile cracks by ductile and brittle fracture mechanisms in a viscoplastic material

In this paper, a finite element analysis of steady-state dynamic crack growth under Mode I, plane strain, small-scale yielding conditions is performed in a rate dependent plastic material characterized by the over-stress model. The main objective of the paper is to obtain theoretically the dependence of dynamic fracture toughness on crack speed. Crack propagation due to a ductile (micro-void) mechanism or a brittle (cleavage) mechanism, as well as transition from one mode to another are considered. The conversion from ductile to brittle has been observed experimentally but has received very little attention using analytical methods. Local fracture criteria based on strains and stresses are used to describe ductile and brittle fracture mechanisms. The results obtained in this paper are in general agreement with micro-structural observations of mode conversion during fracture initiation. Finally, the particular roles played by material rate sensitivity and inertia are examined in some detail.

Composite material and piezoelectric coefficient uncertainty effects on structural health monitoring using feedback control gains as damage indicators

This article addresses uncertainty effect on the health monitoring of a smart structure using control gain shifts as damage indicators. A finite element model of the smart composite plate with surface-bonded piezoelectric sensors and actuators is formulated using first-order shear deformation theory and a matrix crack model is integrated into the finite element model. A constant gain velocity/position feedback control algorithm is used to provide active damping to the structure. Numerical results show that the response of the structure is changed due to matrix cracks and this change can be compensated by actively tuning the feedback controller. This change in control gain can be used as a damage indicator for structural health monitoring. Monte Carlo simulation is conducted to study the effect of material uncertainty on the damage indicator by considering composite material properties and piezoelectric coefficients as independent random variables. It is found that the change in position feedback control gain is a robust damage indicator.

Re-examination of crack opening model of interface fracture

The complex singularity associated with a crack at the interface between two dissimilar, isotropic and homogeneous materials leads to mathematical artefacts, such as stress oscillations and crack face interpenetrations in the vicinity of the crack tip. To avoid these unrealistic features, Sinclair (Sinclair GB. On the stress singularity at an interface crack. International Journal of Fracture 1980;16(2):111-9) assumed a finite crack opening angle (COA) such that the singularity lambda became real equal to 1/2. This paper extends the COA model by considering real singularities not necessarily equal to 1/2. When COA is 0 degrees: the interface crack singularity is complex with a real part equal to 1/2. On increasing COA, the imaginary part of the singularity decreases and becomes zero at a threshold value of COA; at this point, the singularity is a real, repeated value. A further increase in COA results in a pair of real singularities. Different crack opening configurations and material combinations are studied, and results presented for threshold COAs and associated values of singularity. Stress analyses for these three regimes: (a) complex, (b) real pair and (c) real repeated singularities, are reported. It is seen that additional complexities are present in the last case. Typical results for stress fields are also included for comparing with standard fields. (C) 1999 Elsevier Science Ltd. All rights reserved.

Solution combustion synthesis of gamma(L)-Bi(2)MoO(6) and photocatalytic activity under solar radiation

The facile method of solution combustion was used to synthesize gamma(L)-Bi(2)MoO(6). The material was crystallized in a purely crystalline orthorhombic phase with sizes varying from 300 to 500 nm. Because the band gap was 2.51 eV, the degradation of wide variety of cationic and anionic dyes was investigated under solar radiation. Despite the low surface area (< 1 m(2)/g) of the synthesized material, gamma(L)-Bi(2)MoO(6) showed high photocatalytic activity under solar radiation due to its electronic and morphological properties. (C) 2011 Elsevier Ltd. All rights reserved.

The effect of T-stress on plane strain dynamic crack growth in elastic-plastic materials

In this paper, the effects of T -stress on steady, dynamic crack growth in an elastic-plastic material are examined using a modified boundary layer formulation. The analyses are carried out under mode I, plane strain conditions by employing a special finite element procedure based on moving crack tip coordinates. The material is assumed to obey the J (2) flow theory of plasticity with isotropic power law hardening. The results show that the crack opening profile as well as the opening stress at a finite distance from the tip are strongly affected by the magnitude and sign of the T -stress at any given crack speed. Further, it is found that the fracture toughness predicted by the analyses enhances significantly with negative T -stress for both ductile and cleavage mode of crack growth.