Estimation of low-velocity impact damage in laminated composite circular plates using nonlinear finite element analysis

Nonlinear finite element analysis is used for the estimation of damage due to low-velocity impact loading of laminated composite circular plates. The impact loading is treated as an equivalent static loading by assuming the impactor to be spherical and the contact to obey Hertzian law. The stresses in the laminate are calculated using a 48 d.o.f. laminated composite sector element. Subsequently, the Tsai-Wu criterion is used to detect the zones of failure and the maximum stress criterion is used to identify the mode of failure. Then the material properties of the laminate are degraded in the failed regions. The stress analysis is performed again using the degraded properties of the plies. The iterative process is repeated until no more failure is detected in the laminate. The problem of a typical T300/N5208 composite [45 degrees/0 degrees/-45 degrees/90 degrees](s) circular plate being impacted by a spherical impactor is solved and the results are compared with experimental and analytical results available in the literature. The method proposed and the computer code developed can handle symmetric, as well as unsymmetric, laminates. It can be easily extended to cover the impact of composite rectangular plates, shell panels and shells.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

Study of Layering Procedures in Finite-Element Analysis of RC Flexural and Torsional Elements

A new postcracking formulation for concrete, along with both implicit and explicit layering procedures, is used in the analysis of reinforced-concrete (RC) flexural and torsional elements. The postcracking formulation accounts for tension stiffening in concrete along the rebar directions, compression softening in cracked concrete based on either stresses or strains, and aggregate interlock based on crack-confining normal stresses. Transverse shear stresses computed using the layering procedures are included in material model considerations that permit the development of inclined cracks through the RC cross section. Examples of a beam analyzed by both the layering techniques, a torsional element, and a column-slab connection region analyzed by the implicit layering procedure are presented here. The study highlights the primary advantages and disadvantages of each layering approach, identifying the class of problems where the application of either procedure is more suitable.

Vortex sheet simulation of a plane �canonical� mixing layer

Discrete vortex simulations of the mixing layer carried out in the past have usually involved large induced velocity fluctuations, and thus demanded rather long time-averaging to obtain satisfactory values of Reynolds stresses and third-order moments. This difficulty has been traced here, in part, to the use of discrete vortices to model what in actuality are continuous vortex sheets. We propose here a novel two-dimensional vortex sheet technique for computing mixing layer flow in the limit of infinite Reynolds number. The method divides the vortex sheet into constant-strength linear elements, whose motions are computed using the Biot-Savart law. The downstream far-field is modelled by a steady vorticity distribution derived by application of conical similarity from the solution obtained in a finite computational domain. The boundary condition on the splitter plate is satisfied rigorously using a doublet sheet. The computed large-scale roll-up of the vortex sheet is qualitatively similar to experimentally obtained shadow-graphs of the plane turbulent mixing layer. The mean streamwise velocity profile and the growth rate agree well with experimental data. The presently computed Reynolds stresses and third-order moments are comparable with experimental and previous vortex-dynamical results, without using any external parameter (such as the vortex core-size) of the kind often used in the latter. The computed autocorrelations are qualitatively similar to experimental results along the top and bottom edges of the mixing layer, and show a well-defined periodicity along the centreline. The accuracy of the present computation is independently established by demonstrating negligibly small changes in the five invariants (including the Hamiltonian) in vortex dynamics.

Lamination-dependent shear deformation models for bending of angle-ply laminated plates

Lamination-dependent shear corrective terms in the analysis of bending of laminated plates are derived from a priori assumed linear thicknesswise distributions for gradients of transverse shear stresses by using CLPT inplane stresses in the two in-plane equilibrium equations of elasticity in each ply. In the development of a general model for angle-ply laminated plates, special cases like cylindrical bending of laminates in either direction, symmetric laminates, cross-ply laminates, antisymmetric angle-ply laminates, homogeneous plates are taken into consideration. Adding these corrective terms to the assumed displacements in (i) Classical Laminate Plate Theory (CLPT) and (ii) Classical Laminate Shear Deformation Theory (CLSDT), two new refined lamination-dependent shear deformation models are developed. Closed form solutions from these models are obtained for antisymmetric angle-ply laminates under sinusoidal load for a type of simply supported boundary conditions. Results obtained from the present models and also from Ren's model (1987) are compared with each other.

Lamination-dependent shear deformation models for cylindrical bending of angle-ply laminates

Lamination-dependent shear corrective terms in the analysis of flexure of laminates are derived from a priori assumed linear thicknesswise distributions for gradients of transverse shear stresses and using them in the two in-plane equilibrium equations of elasticity in each ply. Adding these corrective terms to (i) Classical Laminate Plate Theory (CLPT) displacements and (ii) Classical Laminate Shear Deformation Theory (CLSDT) displacements, four new refined lamination-dependent shear deformation models for angle-ply laminates are developed. Performance of these models is evaluated by comparing the results from these models with those from exact elasticity solutions for antisymmetric 2-ply laminates and for 4-ply [15/-15](s) laminates. In general, the model with shear corrective terms based on CLPT and added to CLSDT displacements is sufficient and predicts good estimates, both qualitatively and quantitatively, for all displacements and stresses.

Modelling plane mixing layers using vortex points and sheets

We present here a critical assessment of two vortex approaches (both two-dimensional) to the modelling of turbulent mixing layers. In the first approach the flow is represented by point vortices, and in the second it is simulated as the evolution of a continuous vortex sheet composed of short linear elements or ''panels''. The comparison is based on fresh simulations using approximately the same number of elements in either model, paying due attention in both to the boundary conditions far downstream as well as those on the splitter plate from which the mixing layer issues. The comparisons show that, while both models satisfy the well-known invariants of vortex dynamics approximately to the same accuracy, the vortex panel model, although ultimately not convergent, leads to smoother roll-up and values of stresses and moments that are in closer agreement with the experiment, and has a higher computational efficiency for a given degree of convergence on moments. The point vortex model, while faster for a given number of elements, produces an unsatisfactory roll-up which (for the number of elements used) is rendered worse by the incorporation of the Van der Vooren correction for sheet curvature.

Universal Stability and Temperature Dependent Phase Transformation in Group VIIIB-IB Transition Metal FCC Nanowires

Atomistic simulation of Ag, Al, Au, Cu, Ni, Pd, and Pt FCC metallic nanowires show a universal FCC -> HCP phase transformation below a critical cross-sectional size, which is reported for the first time in this paper. The newly observed HCP structure is also confirmed from previous experimental results. Above the critical cross-sectional size, initial < 100 >/{100} FCC metallic nanowires are found to be metastable. External thermal heating shows the transformation of metastable < 100 >/{100} FCC nanowires into < 110 >/{111} stable configuration. Size dependent metastability/instability is also correlated with initial residual stresses of the nanowire by use of molecular static simulation using the conjugant gradient method at a temperature of 0 K. It is found that a smaller cross-sectional dimension of an initial FCC nanowire shows instability due to higher initial residual stresses, and the nanowire is transformed into the novel HCP structure. The initial residual stress shows reduction with an increase in the cross-sectional size of the nanowires. A size dependent critical temperature is also reported for metastable FCC nanowires using molecular dynamic, to capture the < 110 >/{111} to < 100 >/{100} shape memory and pseudoelasticity.

Analysis, Design, Modeling, and Implementation of an Active Clamp HF Link Converter

This paper deals with the system oriented analysis, design, modeling, and implementation of active clamp HF link three phase converter. The main advantage of the topology is reduced size, weight, and cost of the isolation transformer. However, violation of basic power conversion rules due to presence of the leakage inductance in the HF transformer causes over voltage stresses across the cycloconverter devices. It makes use of the snubber circuit necessary in such topologies. The conventional RCD snubbers are dissipative in nature and hence inefficient. The efficiency of the system is greatly improved by using regenerative snubber or active clamp circuit. It consists of an active switching device with an anti-parallel diode and one capacitor to absorb the energy stored in the leakage inductance of the isolation transformer and to regenerate the same without affecting circuit performance. The turn on instant and duration of the active device are selected such that it requires simple commutation requirements. The time domain expressions for circuit dynamics, design criteria of the snubber capacitor with two conflicting constrains (over voltage stress across the devices and the resonating current duration), the simulation results based on generalized circuit model and the experimental results based on laboratory prototype are presented.

Microstructural stability and superplasticity in an electrodeposited nanocrystalline Ni-P alloy

Stabilization of nanocrystalline grain sizes by second phase particles can facilitate superplasticity at high strain rates and/or low temperatures. A metastable single phase nano-Ni-P alloy prepared by electrodeposition, with a grain size of similar to 6 nm, transforms to a nanoduplex structure at T> 673 K, with similar to 4 vol.% Ni3P particles at triple junctions and within Ni grains. The nanoduplex microstructure is reasonably stable up to 777 K, and the growth of Ni grains occurs in a coupled manner with the growth of Ni3P particles such that the ratio of the two mean sizes (Z) is essentially constant. High temperature tests for a grain size of 290 nm reveal superplastic behavior with an optimum elongation to failure of 810% at a strain rate of 7 x 10(-4) s(-1) and a relatively low temperature of 777 K. Superplastic deformation enhances both grain growth and the ratio Z, implying that grain boundary sliding (GBS) significantly influences the microstructural dynamics. Analysis of the deformation processes suggests that superplasticity is associated with GBS controlled by the overcoming of intragranular particles by dislocations, so that deformation is independent of the grain size. The nano-Ni-P alloy exhibits lower ductility than nano-Ni due to concurrent cavitation caused by higher stresses. (C) 2011 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Unsteady axisymmetric stagnation-point flow of a viscous fluid on a cylinder

The unsteady viscous flow in the vicinity of an axisymmetric stagnation point of an infinite circular cylinder is investigated when both the free stream velocity and the velocity of the cylinder vary arbitrarily with time. The cylinder moves either in the same direction as that of the free stream or in the opposite direction. The flow is initially (t = 0) steady and then at t > 0 it becomes unsteady. The semi-similar solution of the unsteady Navier-Stokes equations has been obtained numerically using an implicit finite-difference scheme. Also the self-similar solution of the Navier-Stokes equations is obtained when the velocity of the cylinder and the free stream velocity vary inversely as a linear function of time. For small Reynolds number, a closed form solution is obtained. When the Reynolds number tends to infinity, the Navier-Stokes equations reduce to those of the two-dimensional stagnation-point flow. The shear stresses corresponding to stationary and the moving cylinder increase with the Reynolds number. The shear stresses increase with time for the accelerating flow but decrease with increasing time for the decelerating flow. For the decelerating case flow reversal occurs in the velocity profiles after a certain instant of time. (C) 1999 Elsevier Science Ltd. All rights reserved.

Tracking and erosion resistance of silicone rubber nanocomposites under positive and negative dc voltages

ASTM D2303 standard provides a method for evaluating the tracking and erosion resistance of polymeric insulators under ac voltages. In this paper, the above method has been extended for evaluating the performance of the insulators under dc stresses. Tests were conducted on polymeric silicone rubber (SR) insulators under positive and negative dc stresses. Micron sized Alumina trihydrate (uATH) and nano sized Alumina (nALU) were used as fillers in SR matrix to improve the resistance to tracking and erosion. Results suggest that SR composites perform better under negative dc than under positive dc voltages. Eroded mass and leakage current data support the above result. Samples with low concentration of nano alumina fillers performed on par with the samples with large loadings of uATH.

Critical stress intensity factors for cracked hollow pipes under transient thermal loads

Finite element analyses of a long hollow cylinder having an axisymmetric circumferential internal edge crack, subjected to convective cooling on the inner surface are performed. The transient thermal stress intensity factor is estimated using a domain version of the J-integral method. The effect of the thickness of the cylinder, crack length, and heat transfer coefficient on the stress intensity factor history are studied. The variations of critical normalized stress intensity factor with crack length-to-thickness ratio for different parameters are presented. The results show that if a small inner surface crack begins to grow, its stress intensity factor will increase with increase in crack length, reach a maximum, and then begin to drop. Based on the results, a fracture-based design methodology for cracked hollow pipes under transient thermal loads is discussed.

An assessment of creep deformation and fracture behavior of 2.25Cr-1Mo similar and dissimilar weld joints

The evaluation of the creep deformation and fracture behavior of a 2.25Cr-1Mo steel base metal, a 2.25Cr-1Mo/2.25Cr-1Mo similar weld joint, and a 2.25Cr-1Mo/Alloy 800 dissimilar weld joint at 823 K over a stress range of 90 to WO MPa has been carried out. The specimens for creep testing were taken from single-V weld pads fabricated by a shielded metal arc-welding process using 2.25Cr-1Mo steel (for similar-joint) and INCONEL 182 (for dissimilar-joint) electrodes. The weld pads were subsequently given a postweld hear treatment (PWHT) of 973 K for I hour. The microstructure and microhardness of the weld joints were evaluated in the as-welded, postweld heat-treated, and creep-tested conditions. The heat-affected zone (HAZ) of similar weld joint consisted of bainite in the coarse-prior-austenitic-grain (CPAG) region near the fusion line, followed by bainite in the fine-prior-austenitic-grain (FPAG) and intercritical regions merging with the unaffected base metal. In addition to the HAZ structures in the 2.25Cr-1Mo steel, the dissimilar weld joint displayed a definite INCONEL/2.25Cr-1Mo weld interface structure present either as a sharp line or as a diffuse region. A hardness trough was observed in the intercritical region of the HAZ in both weld joints, while a maxima in hardness was seen at the weld interface of the dissimilar weld joint. Both weld joints exhibited significantly lower rupture lives compared to the 2.25Cr-1Mo base metal. The dissimilar weld joint exhibited poor rupture life compared to the similar weld joint, at applied stresses lower than 130 MPa. In both weld joints, the strain distribution across the specimen gage length during creep testing varied significantly. During creep testing, localization of deformation occurred in the intercritical HAZ. In the similar weld joint, at all stress levels investigated, and in the dissimilar weld joint, at stresses greater than or equal to 150 MPa, the creep failure occulted in the intercritical HAZ. The fracture occurred by transgranular mode with a large number of dimples. At stresses below 150 MPa, the failure in the dissimilar weld joint occurred in the CPAG HAZ near to the weld interface. The failure occurred by extensive intergranular creep cavity formation.

Analysis of deformation mechanisms in superplastic yttria-stabilized tetragonal zirconia

There have been extensive experimental observations of changes in the apparent rate controlling creep parameters in studies on superplastic materials. The three most common explanations associated with these changes in the stress exponent, n, the activation energy Q and the inverse grain size exponent, p involve the effect of concurrent grain growth, the operation of a threshold stress or transitions in creep mechanisms. Each of these factors may influence experimental creep data in a similar manner. Therefore, a careful analysis of the consequences of all three factors must involve the development of a consistent set of experimental observations in order to adequately distinguish the effects of each. This paper discusses the role of concurrent grain growth, a threshold stress and transitions in creep mechanisms in superplastic materials. Specific attention is given to the analysis of data on superplastic yttria-stabilized zirconia ceramics for which an increase in n has been observed at low applied stresses. It is demonstrated that neither concurrent grain growth nor a threshold stress can account for all the relevant experimental observations in this material. It is concluded that the changes in rate controlling creep parameters are associated with the operation of two distinct sequential mechanisms as part of a grain boundary sliding process.