Review of continuum, finite element and hybrid techniques in the analysis of stress concentrations in structures

When a uniform flow of any nature is interrupted, the readjustment of the flow results in concentrations and rare-factions, so that the peak value of the flow parameter will be higher than that which an elementary computation would suggest. When stress flow in a structure is interrupted, there are stress concentrations. These are generally localized and often large, in relation to the values indicated by simple equilibrium calculations. With the advent of the industrial revolution, dynamic and repeated loading of materials had become commonplace in engine parts and fast moving vehicles of locomotion. This led to serious fatigue failures arising from stress concentrations. Also, many metal forming processes, fabrication techniques and weak-link type safety systems benefit substantially from the intelligent use or avoidance, as appropriate, of stress concentrations. As a result, in the last 80 years, the study and and evaluation of stress concentrations has been a primary objective in the study of solid mechanics. Exact mathematical analysis of stress concentrations in finite bodies presents considerable difficulty for all but a few problems of infinite fields, concentric annuli and the like, treated under the presumption of small deformation, linear elasticity. A whole series of techniques have been developed to deal with different classes of shapes and domains, causes and sources of concentration, material behaviour, phenomenological formulation, etc. These include real and complex functions, conformal mapping, transform techniques, integral equations, finite differences and relaxation, and, more recently, the finite element methods. With the advent of large high speed computers, development of finite element concepts and a good understanding of functional analysis, it is now, in principle, possible to obtain with economy satisfactory solutions to a whole range of concentration problems by intelligently combining theory and computer application. An example is the hybridization of continuum concepts with computer based finite element formulations. This new situation also makes possible a more direct approach to the problem of design which is the primary purpose of most engineering analyses. The trend would appear to be clear: the computer will shape the theory, analysis and design.

A numerical study of projectile impact on thin aluminium plates

This article deals with a simulation-based Study of the impact of projectiles on thin aluminium plates using LS-DYNA by modelling plates with shell elements and projectiles with solid elements. In order to establish the required modelling criterion in terms of element size for aluminium plates, a convergence Study of residual velocity has been carried Out by varying mesh density in the impact zone. Using the preferred material and meshing criteria arrived at here, extremely good prediction of test residual velocities and ballistic limits given by Gupta et al. (2001) for thin aluminium plates has been obtained. The simulation-based pattern of failure with localized bulging and jagged edge of perforation is similar to the perforation with petalling seen in tests. A number Of simulation-based parametric studies have been carried out and results consistent with published test data have been obtained. Despite the robust correlation achieved against published experimental results, it would be prudent to conduct one's own experiments, for a final correlation via the present modelling procedure and analysis with the explicit LS-DYNTA 970 solver. Hence, a sophisticated ballistic impact testing facility and a high-speed camera have been used to conduct additional tests on grade 1100 aluminium plates of 1 mm thickness with projectiles Of four different nose shapes. Finally, using the developed numerical simulation procedure, an excellent correlation of residual velocity and failure modes with the corresponding test results has been obtained.

CAE-based prediction of projectile residual velocity for impact on single and multi-layered metallic armour plates

The present, paper deals with the CAE-based study Of impact of jacketed projectiles on single- and multi-layered metal armour plates using LS-DYNA. The validation of finite element modelling procedure is mainly based on the mesh convergence study using both shell and solid elements for representing single-layered mild steel target plates. It, is shown that the proper choice of mesh density and the strain rate-dependent material properties are essential for all accurate prediction of projectile residual velocity. The modelling requirements are initially arrived at by correlating against test residual velocities for single-layered mild steel plates of different depths at impact velocities in the ran.-c of approximately 800-870 m/s. The efficacy of correlation is adjudged, in terms of a 'correlation index', defined in the paper: for which values close to unity are desirable. The experience gained for single-layered plates is next; used in simulating projectile impacts on multi-layered mild steel target plates and once again a high degree of correlation with experimental residual velocities is observed. The study is repeated for single- and multi-layered aluminium target plates with a similar level of success in test residual velocity prediction. TO the authors' best knowledge, the present comprehensive study shows in particular for the first time that, with a. proper modelling approach, LS-DYNA can be used with a great degree of confidence in designing perforation-resistant single and multi-layered metallic armour plates.

Prediction of solid block masonry prism compressive strength using FE model

Masonry strength is dependent upon characteristics of the masonry unit,the mortar and the bond between them. Empirical formulae as well as analytical and finite element (FE) models have been developed to predict structural behaviour of masonry. This paper is focused on developing a three dimensional non-linear FE model based on micro-modelling approach to predict masonry prism compressive strength and crack pattern. The proposed FE model uses multi-linear stress-strain relationships to model the non-linear behaviour of solid masonry unit and the mortar. Willam-Warnke's five parameter failure theory developed for modelling the tri-axial behaviour of concrete has been adopted to model the failure of masonry materials. The post failure regime has been modelled by applying orthotropic constitutive equations based on the smeared crack approach. Compressive strength of the masonry prism predicted by the proposed FE model has been compared with experimental values as well as the values predicted by other failure theories and Eurocode formula. The crack pattern predicted by the FE model shows vertical splitting cracks in the prism. The FE model predicts the ultimate failure compressive stress close to 85 of the mean experimental compressive strength value.

Robust prediction of residual velocities and ballistic limits of projectiles for impact on thin aluminium plates

This paper deals with the simulation-driven study of the impact of hardened steel projectiles on thin aluminium target plates using explicit finite element analysis as implemented in LS-DYNA. The evaluation of finite element modelling includes a comprehensive mesh convergence study using shell elements for representing target plates and the solid element-based representation of ogivalnosed projectiles. A user-friendly automatic contact detection algorithm is used for capturing interaction between the projectile and the target plate. It is shown that the proper choice of mesh density and strain rate-dependent material properties is crucial as these parameters significantly affect the computed residual velocity. The efficacy of correlation with experimental data is adjudged in terms of a 'correlation index' defined in the present study for which values close to unity are desirable.By simulating laboratory impact tests on thin aluminium plates carried out by earlier investigators, extremely good prediction of experimental ballistic limits has been observed with correlation indices approaching unity. Additional simulation-based parametric studies have been carried out and results consistent with test data have been obtained. The simulation procedures followed in the present study can be applied with confidence in designing thin aluminium armour plates for protection against low calibre projectiles.

Room temperature ferromagnetism in ZnO (core)/graphite (shell) nanowires fabricated by a one-step method

ZnO (core)/graphitic (shell) nanowires were successfully fabricated by a one-step method. Morphology of the as-grown nanowires was studied in detail by scanning electron microscopy, transmission electron microscopy (TEM), and energy dispersive X-ray analysis (EDS). High resolution TEM micrographs and selected area electron diffraction patterns reveal the core/shell morphology of the nanowires that grew along the c-axis of ZnO. EDS study of the nanowires confirms that there are no impurities within the detectable limit. Superconducting quantum interference device magnetometer measurements show room temperature ferromagnetic ordering in these core/shell nanowires. (C) 2010 Elsevier Ltd. All rights reserved.

Multiferroic and magnetoelectric properties of core-shell CoFe2O4@BaTiO3 nanocomposites

Core-shell CoFe2O4@BaTiO3 nanoparticles and nanotubes have been prepared using a combination of solution processing and high temperature calcination. Both the core-shell nanostructures exhibit magnetic and dielectric hysteresis at room temperature and magnetoelectric effect. The dielectric constant of both the nanocomposites decreases upon application of magnetic field. The core-shell nanoparticles exhibit 1.7% change in magnetocapacitance around 134 K at 1 T, while the core-shell nanotubes show a remarkable 4.5% change in magnetocapacitance around 310 K at 2 T.(C) 2010 American Institute of Physics. [doi:10.1063/1.3478231].

Preparation and characterization of silane-stabilized, highly uniform, nanobimetallic Pt-Pd particles in solid and liquid matrixes

Highly uniform, stable nanobimetallic dispersions are prepared in a single si ep in the form of sols, gels, and monoliths, using organically modified silicates as the matrix and the stabilizer. The Pt-Pd bimetallic dispersions are characterized by W-vis, TEM, SEM, and XRD measurements. The evolution of silicate was followed by IR spectroscopy. XPS and CO adsorption studies reveal that the structure of the particles consists of a palladium core and a platinum shell. Electrocatalysis of ascorbic acid oxidation has been demonstrated using thin films of silicate containing the nanobimetal particles on a glassy carbon electrode.

Using Water as a Design Element in Crystal Engineering. Host-Guest Compounds of Hydrated 3,5-Dihydroxybenzoic Acid

Thirteen host guest compounds of 3,5-dihydroxybenzoic acid (DHBA) have been structurally characterized. Water molecules occupy the peripheries of a hexagonal void, created with DHBA molecules, and act as ``hooks'' to connect the guest molecules with the host-framework via hydrogen bonding. The ``water hook'' is an OH group acting as a donor. Consequently, the guest molecules were chosen so that they contain good hydrogen bond acceptor functionalities. A number of multicomponent hydrates were isolated with stoichiometries (DHBA)(x)(H2O). (guest),. Of these, compounds with the following as guests were obtained as crystals that were good enough for single crystal work: ethyl acetate (EtOAc), diethyl oxalate, dimethyl oxalate, di(n-propyl) oxalate, diethyl malonate, diethyl succinate, chloroacetonitrile, N,N-dimethyl formamide (DMF), acetone, dimethyl sulfoxide (DMSO), 1-propanol, and 2-butanol. From 2-butanol, a hemihydrate, (DHBA)(2)(H2O), was also obtained concomitantly. Further to guest stabilization, water acts as a good mediator of effective crystal packing and also determines the topology of the host framework. En the present series of compounds, the role of water is wide ranging, and it is not easy to classify it specifically as a host or as a guest.

Shell-model calculation of defect energies in alkali halides employing crystal-independent interionic potential parameters

Shell model calculation of defect energies in alkali halides have been carried out using the ion-dependent, crystal-independent potential parameters of Sangster and Atwood (1978). Results indicate that appreciable differences exist between barrier heights for migration of cations and anions. While barrier heights for cations are generally lower than for anions in alkali halides of NaCl structure, the opposite is true in alkali halides of CsCl structure.

Synthesis and characterization of Pd(0), PdS, and Pd@PdO core-shell nanoparticles by solventless thermolysis of a Pd-thiolate cluster

Colloids of palladium nanoparticles have been prepared by the solvated metal atom dispersion (SMAD) method. The as-prepared Pd colloid consists of particles with an average diameter of 2.8 +/- 0.1 nm. Digestive ripening of the as-prepared Pd colloid, a process involving refluxing the as-prepared colloid at or near the boiling point of the solvent in the presence of a passivating agent, dodecanethiol resulted in a previously reported Pd-thiolate cluster, Pd(SC12H25)(2)](6) but did not render the expected narrowing down of the particle size distribution. Solventless thermolysis of the Pd-thiolate complex resulted in various Pd systems such as Pd(0), PdS, and Pd@PdO core-shell nanoparticles thus demonstrating its versatility. These I'd nanostructures have been characterized using high-resolution electron microscopy and powder X-ray diffraction methods. (C) 2010 Elsevier Inc. All rights reserved.

A numerical study of projectile impact on mild steel armour plates

The present article deals with the development of a finite element modelling approach for the prediction of residual velocities of hard core ogival-nose projectiles following normal impact on mild steel target plates causing perforation. The impact velocities for the cases analysed are in the range 818–866.3 m/s. Assessment of finite element modelling and analysis includes a comprehensive mesh convergence study using shell elements for representing target plates and solid elements for jacketed projectiles with a copper sheath and a rigid core. Dynamic analyses were carried out with the explicit contact-impact LS-DYNA 970 solver. It has been shown that proper choice of element size and strain rate-based material modelling of target plate are crucial for obtaining test-based residual velocity.The present modelling procedure also leads to realistic representation of target plate failure and projectile sheath erosion during perforation, and confirms earlier observations that thermal effects are not significant for impact problems within the ordnance range. To the best of our knowledge, any aspect of projectile failure or degradation obtained in simulation has not been reported earlier in the literature. The validated simulation approach was applied to compute the ballistic limits and to study the effects of plate thickness and projectile diameter on residual velocity, and trends consistent with experimental data for similar situations were obtained.

A Smooth Finite Element Method Based on Reproducing Kernel DMS-Splines

The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials There is thus a case for combining these advantages in a so-called hybrid scheme or a `smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C-p (p >= 1) continuity One such recent attempt, a NURBS based parametric bridging method (Shaw et al 2008b), uses polynomial reproducing, tensor-product non-uniform rational B-splines (NURBS) over a typical FE mesh and relies upon a (possibly piecewise) bijective geometric map between the physical domain and a rectangular (cuboidal) parametric domain The present work aims at a significant extension and improvement of this concept by replacing NURBS with DMS-splines (say, of degree n > 0) that are defined over triangles and provide Cn-1 continuity across the triangle edges This relieves the need for a geometric map that could precipitate ill-conditioning of the discretized equations Delaunay triangulation is used to discretize the physical domain and shape functions are constructed via the polynomial reproduction condition, which quite remarkably relieves the solution of its sensitive dependence on the selected knotsets Derivatives of shape functions are also constructed based on the principle of reproduction of derivatives of polynomials (Shaw and Roy 2008a) Within the present scheme, the triangles also serve as background integration cells in weak formulations thereby overcoming non-conformability issues Numerical examples involving the evaluation of derivatives of targeted functions up to the fourth order and applications of the method to a few boundary value problems of general interest in solid mechanics over (non-simply connected) bounded domains in 2D are presented towards the end of the paper

Flow-Acoustic Performance of Perforated Element Mufflers With Application to Design

A fairly comprehensive computer program incorporating explicit expressions for the four-pole parameters of concentric-tube resonators, plug mufflers, and three-duct cross-flow perforated elements has been used for parametric studies. The parameters considered are hole diameter, the center-to-center distance between consecutive holes (which decides porosity), the incoming mean flow Mach number, the area expansion ratio, the number of partitions of chambers within a given overall shell length, and the relative lengths of these partitions or chambers, all normalized with respect to the exhaust pipe diameter. Transmission loss has been plotted as a function of a normalized frequency parameter. Additionally, the effect of the tail pipe length on insertion loss for an anechoic source has also been studied. These studies have been supplemented by empirical expressions for the normalized static pressure drop for different types of perforated-element mufflers developed from experimental observations.

The Reproducing Kernel DMS-FEM: 3D Shape Functions and Applications to Linear Solid Mechanics

We propose a family of 3D versions of a smooth finite element method (Sunilkumar and Roy 2010), wherein the globally smooth shape functions are derivable through the condition of polynomial reproduction with the tetrahedral B-splines (DMS-splines) or tensor-product forms of triangular B-splines and ID NURBS bases acting as the kernel functions. While the domain decomposition is accomplished through tetrahedral or triangular prism elements, an additional requirement here is an appropriate generation of knotclouds around the element vertices or corners. The possibility of sensitive dependence of numerical solutions to the placements of knotclouds is largely arrested by enforcing the condition of polynomial reproduction whilst deriving the shape functions. Nevertheless, given the higher complexity in forming the knotclouds for tetrahedral elements especially when higher demand is placed on the order of continuity of the shape functions across inter-element boundaries, we presently emphasize an exploration of the triangular prism based formulation in the context of several benchmark problems of interest in linear solid mechanics. In the absence of a more rigorous study on the convergence analyses, the numerical exercise, reported herein, helps establish the method as one of remarkable accuracy and robust performance against numerical ill-conditioning (such as locking of different kinds) vis-a-vis the conventional FEM.