Role of elasto-plastic analysis under cyclic loading in fatigue crack growth studies

Linear Elastic Fracture Mechanics (LEFM) has been widely used in the past for fatigue crack growth studies, but this is acceptable only in situations which are within small scale yielding (SSY). In many practical structural components, conditions of SSY could be violated and one has to look for fracture criteria based on elasto-plastic analysis. Crack closure phenomenon, one of the most striking discoveries based on inelastic deformations during crack growth, has significant effect on fatigue crack growth rate. Numerical simulation of this phenomenon is computationally intensive and involved but has been successfully implemented. Stress intensity factors and strain energy release rates lose their meaning, J-integral (or its incremental) values are applicable only in specific situations, whereas alternate path independent integrals have been proposed in the literature for use with elasto-plastic fracture mechanics (EPFM) based criteria. This paper presents certain salient features of two independent finite element (numerical) studies of relevance to fatigue crack growth, where elasto-plastic analysis becomes significant. These problems can only be handled in the current day computational environment, and would have been only a dream just a few years ago.

The hardnesses of poly(methylmethacrylate)

The paper describes an experimental study of the normal and scratch hardnesses of a poly(methylmethacrylate). The deformations have been introduced using hard steel cones of a range of included cone angles. The influence of the state of interfacial lubrication is examined and rationalized. The observed time dependence of the two types of computed hardness data is compared and the nature of the correlations between these data is evaluated. It is observed that when the imposed strains are modest, say less than 0.2, the scratch hardness and normal hardness deformations produce self consistent data using first order and rather indiscriminate analyses for both types of deformations. At higher levels of imposed strain, a more critical appraisal of the nature of the deformation produced in the two cases is necessary in order to provide mutually consistent hardness values and hence unequivocal rheological characteristics for this polymer.

Durable electrocatalytic-activity of Pt-Au/C cathode in PEMFCs

Longevity remains as one of the central issues in the successful commercialization of polymer electrolyte membrane fuel cells (PEMFCs) and primarily hinges on the durability of the cathode. Incorporation of gold (Au) to platinum (Pt) is known to ameliorate both the electrocatalytic activity and stability of cathode in relation to pristine Pt-cathodes that are currently being used in PEMFCs. In this study, an accelerated stress test (AST) is conducted to simulate prolonged fuel-cell operating conditions by potential cycling the carbon-supported Pt-Au (Pt-Au/C) cathode. The loss in performance of PEMFC with Pt-Au/C cathode is found to be similar to 10% after 7000 accelerated potential-cycles as against similar to 60% for Pt/C cathode under similar conditions. These data are in conformity with the electrochemical surface-area values. PEMFC with Pt-Au/C cathode can withstand > 10 000 potential cycles with very little effect on its performance. X-ray diffraction and transmission electron microscopy studies on the catalyst before and after AST suggest that incorporating Au with Pt helps mitigate aggregation of Pt particles during prolonged fuel-cell operations while X-ray photoelectron spectroscopy reflects that the metallic nature of Pt is retained in the Pt-Au catalyst during AST in comparison to Pt/C that shows a major portion of Pt to be present as oxidic platinum. Field-emission scanning electron microscopy conducted on the membrane electrode assembly before and after AST suggests that incorporating Au with Pt helps mitigating deformations in the catalyst layer.

Optimum design of a helicopter rotor for low vibration using aeroelastic analysis and response surface methods

An aeroelastic analysis based on finite elements in space and time is used to model the helicopter rotor in forward flight. The rotor blade is represented as an elastic cantilever beam undergoing flap and lag bending, elastic torsion and axial deformations. The objective of the improved design is to reduce vibratory loads at the rotor hub that are the main source of helicopter vibration. Constraints are imposed on aeroelastic stability, and move limits are imposed on the blade elastic stiffness design variables. Using the aeroelastic analysis, response surface approximations are constructed for the objective function (vibratory hub loads). It is found that second order polynomial response surfaces constructed using the central composite design of the theory of design of experiments adequately represents the aeroelastic model in the vicinity of the baseline design. Optimization results show a reduction in the objective function of about 30 per cent. A key accomplishment of this paper is the decoupling of the analysis problem and the optimization problems using response surface methods, which should encourage the use of optimization methods by the helicopter industry. (C) 2002 Elsevier Science Ltd. All rights reserved.

Water content - void ratio swell-shrink paths of compacted expansive soils

This paper addresses the behaviour of compacted expansive soils under swell-shrink cycles. Laboratory cyclic swell-shrink tests were conducted on compacted specimens of two expansive soils at surcharge pressures of 6.25, 50.00, and 100.00 kPa. The void ratio and water content of the specimens at several intermediate stages during swelling until the end of swelling and during shrinkage until the end of shrinkage were determined to trace the water content versus void ratio paths with an increasing number of swell-shrink cycles. The test results showed that the swell-shrink path was reversible once the soil reached an equilibrium stage where the vertical deformations during swelling and shrinkage were the same. This usually occurred after about four swell-shrink cycles. The swelling and shrinkage path of each specimen subjected to full swelling - full shrinkage cycles showed an S-shaped curve (two curvilinear portions and a linear portion). However, the swelling and shrinkage path occurred as a part of the S-shaped curve, when the specimen was subjected to full swelling - partial shrinkage cycles. More than 80% of the total volumetric change and more than 50% of the total vertical deformation occurred in the central linear portion of the S-shaped curve. The volumetric change was essentially parallel to the saturation line within a degree of saturation range of 50-80% for the equilibrium cycle. The primary value of the swell-shrink path is to provide information regarding the void ratio change that would occur for a given change in water content for any possible swell-shrink pattern. It is suggested that these swell-shrink paths can be established with a limited number of tests in the laboratory.

Flow of a stratified fluid in a wavy channel

The flow of a stratified fluid in a channel with small and large deformations is investigated. The analogy of this flow with swirling flow in tubes with non-uniform cross-sections is studied. The flow near the wall is blocked when the Froude number takes certain critical values. The possibility of preventing the stagnation zones in the flow field is also discussed

Size-dependent and strain rate effects in quantum dot nanostructures with molecular dynamic simulations

Size and strain rate effects are among several factors which play an important role in determining the response of nanostructures, such as their deformations, to the mechanical loadings. The mechanical deformations in nanostructure systems at finite temperatures are intrinsically dynamic processes. Most of the recent works in this context have been focused on nanowires [1, 2], but very little attention has been paid to such low dimensional nanostructures as quantum dots (QDs). In this contribution, molecular dynamics (MD) simulations with an embedded atom potential method(EAM) are carried out to analyse the size and strain rate effects in the silicon (Si) QDs, as an example. We consider various geometries of QDs such as spherical, cylindrical and cubic. We choose Si QDs as an example due to their major applications in solar cells and biosensing. The analysis has also been focused on the variation in the deformation mechanisms with the size and strain rate for Si QD embedded in a matrix of SiO2 [3] (other cases include SiN and SiC matrices).It is observed that the mechanical properties are the functions of the QD size, shape and strain rate as it is in the case for nanowires [2]. We also present the comparative study resulted from the application of different EAM potentials in particular, the Stillinger-Weber (SW) potential, the Tersoff potentials and the environment-dependent interatomic potential (EDIP) [1]. Finally, based on the stabilized structural properties we compute electronic bandstructures of our nanostructures using an envelope function approach and its finite element implementation.

Design of Geocell Reinforcement for Supporting Embankments on Soft Ground

The methods of design available for geocell-supported embankments are very few. Two of the earlier methods are considered in this paper and a third method is proposed and compared with them. The first method is the slip line method proposed by earlier researchers. The second method is based on slope stability analysis proposed by this author earlier and the new method proposed is based on the finite element analyses. In the first method, plastic bearing failure of the soil was assumed and the additional resistance due to geocell layer is calculated using a non-symmetric slip line field in the soft foundation soil. In the second method, generalpurpose slope stability program was used to design the geocell mattress of required strength for embankment using a composite model to represent the shear strength of geocell layer. In the third method proposed in this paper, geocell reinforcement is designed based on the plane strain finite element analysis of embankments. The geocell layer is modelled as an equivalent composite layer with modified strength and stiffness values. The strength and dimensions of geocell layer is estimated for the required bearing capacity or permissible deformations. These three design methods are compared through a design example.

Geometrically Nonlinear Dynamics of Composite Four-Bar Mechanisms

This work intends to demonstrate the importance of geometrically nonlinear crosssectional analysis of certain composite beam-based four-bar mechanisms in predicting system dynamic characteristics. All component bars of the mechanism are made of fiber reinforced laminates and have thin rectangular cross-sections. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. They are linked to each other by means of revolute joints. We restrict ourselves to linear materials with small strains within each elastic body (beam). Each component of the mechanism is modeled as a beam based on geometrically nonlinear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and nonlinear 1-D analyses along the four beam reference curves. For thin rectangular cross-sections considered here, the 2-D cross-sectional nonlinearity is overwhelming. This can be perceived from the fact that such sections constitute a limiting case between thin-walled open and closed sections, thus inviting the nonlinear phenomena observed in both. The strong elastic couplings of anisotropic composite laminates complicate the model further. However, a powerful mathematical tool called the Variational Asymptotic Method (VAM) not only enables such a dimensional reduction, but also provides asymptotically correct analytical solutions to the nonlinear cross-sectional analysis. Such closed-form solutions are used here in conjunction with numerical techniques for the rest of the problem to predict multi-body dynamic responses, more quickly and accurately than would otherwise be possible. The analysis methodology can be viewed as a three-step procedure: First, the cross-sectional properties of each bar of the mechanism is determined analytically based on an asymptotic procedure, starting from Classical Laminated Shell Theory (CLST) and taking advantage of its thin strip geometry. Second, the dynamic response of the nonlinear, flexible fourbar mechanism is simulated by treating each bar as a 1-D beam, discretized using finite elements, and employing energy-preserving and -decaying time integration schemes for unconditional stability. Finally, local 3-D deformations and stresses in the entire system are recovered, based on the 1-D responses predicted in the previous step. With the model, tools and procedure in place, we shall attempt to identify and investigate a few problems where the cross-sectional nonlinearities are significant. This will be carried out by varying stacking sequences and material properties, and speculating on the dominating diagonal and coupling terms in the closed-form nonlinear beam stiffness matrix. Numerical examples will be presented and results from this analysis will be compared with those available in the literature, for linear cross-sectional analysis and isotropic materials as special cases.

A Comparative Study of the Formulations for Topology Optimization of Compliant Mechanisms

The topology optimization problem for the synthesis of compliant mechanisms has been formulated in many different ways in the last 15 years, but there is not yet a definitive formulation that is universally accepted. Furthermore, there are two unresolved issues in this problem. In this paper, we present a comparative study of five distinctly different formulations that are reported in the literature. Three benchmark examples are solved with these formulations using the same input and output specifications and the same numerical optimization algorithm. A total of 35 different synthesis examples are implemented. The examples are limited to desired instantaneous output direction for prescribed input force direction. Hence, this study is limited to linear elastic modeling with small deformations. Two design parameterizations, namely, the frame element based ground structure and the density approach using continuum elements, are used. The obtained designs are evaluated with all other objective functions and are compared with each other. The checkerboard patterns, point flexures, the ability to converge from an unbiased uniform initial guess, and the computation time are analyzed. Some observations are noted based on the extensive implementation done in this study. Complete details of the benchmark problems and the results are included. The computer codes related to this study are made available on the internet for ready access.

Instabilities of soft elastic microtubes filled with viscous fluids: Pearls, wrinkles, and sausage strings

A linear stability analysis is presented to study the self-organized instabilities of a highly compliant elastic cylindrical shell filled with a viscous liquid and submerged in another viscous medium. The prototype closely mimics many components of micro-or nanofluidic devices and biological processes such as the budding of a string of pearls inside cells and sausage-string formation of blood vessels. The cylindrical shell is considered to be a soft linear elastic solid with small storage modulus. When the destabilizing capillary force derived from the cross-sectional curvature overcomes the stabilizing elastic and in-plane capillary forces, the microtube can spontaneously self-organize into one of several possible configurations; namely, pearling, in which the viscous fluid in the core of the elastic shell breaks up into droplets; sausage strings, in which the outer interface of the mircrotube deforms more than the inner interface; and wrinkles, in which both interfaces of the thin-walled mircrotube deform in phase with small amplitudes. This study identifies the conditions for the existence of these modes and demonstrates that the ratios of the interfacial tensions at the interfaces, the viscosities, and the thickness of the microtube play crucial roles in the mode selection and the relative amplitudes of deformations at the two interfaces. The analysis also shows asymptotically that an elastic fiber submerged in a viscous liquid is unstable for Y = gamma/(G(e)R) > 6 and an elastic microchannel filled with a viscous liquid should rupture to form spherical cavities (pearling) for Y > 2, where gamma, G(e), and R are the surface tension, elastic shear modulus, and radius, respectively, of the fiber or microchannel.

Synthesis and Characterization of Carbon Nanotube-Polymer Multilayer Structures

We develop lightweight, multilayer materials composed of alternating layers of poly dimethyl siloxane (PDMS) polymer and vertically aligned carbon nanotube (CNT) arrays, and characterize their mechanical response in compression. The CNT arrays used In the assembly are synthesized with graded mechanical properties along their thickness, and their use enables the creation of multilayer structures with low density (0.12-0.28 g/cm(3)). We test the mechanical response of structures composed of different numbers of CNT layers partially embedded in PDMS polymer, under quasi-static and dynamic loading. The resulting materials exhibit a hierarchical, fibrous structure with unique mechanical properties: They can sustain large compressive deformations (up to similar to 0.8 strain) with a nearly complete recovery and present strain localization in selected sections of the materials. Energy absorption, as determined by the hysteresis observed In stress-strain curves, is found to be at least 3 orders of magnitude larger than that of natural and synthetic cellular materials of comparable density. Conductive bucky paper Is Included within the polymer interlayers. This allows the measurement of resistance variation as a function of applied stress, showing strong correlation with the observed strain localization In compression.

Effect of shear deformation on interfacial stresses in plated beams subjected to arbitrary loading

Bonding a fibre reinforced polymer (FRP) composite or metallic plate to the soffit of a reinforced concrete (RC), timber or metallic beam can significantly increase its strength and other aspects of structural performance. These hybrid beams are often found to fail due to premature debonding of the plate from the original beam in a brittle manner. This has led to the development of many analytical solutions over the last two decades to quantify the interfacial shear and normal stresses between the adherends. The adherends are subjected to axial, bending and shear deformations. However, most analytical solutions have neglected the influence of shear deformation of the adherends. For the few solutions which consider this effect in an approximate manner, their applicability is limited to one or two specific load cases. This paper presents a general analytical solution for the interfacial stresses in plated beams under an arbitrary loading with the shear deformation of the adherends duly considered. The shear stress distribution is assumed to be parabolic through the depth of the adherends in predicting the interfacial shear stress and Timoshenko's beam theory is adopted in predicting interfacial normal stress to account for the shear deformation. The solution is applicable to a beam of arbitrary prismatic cross-section bonded symmetrically or asymmetrically with a thin or thick plate, both having linear elastic material properties. The effect of shear deformation is illustrated through an example beam. The influence of material and geometric parameters of the adherends and adhesive on the interfacial stress concentrations at the plate end is discussed. (C) 2011 Elsevier Ltd. All rights reserved.

Macromolecular Transport Through Porous Arterial Walls

Arteries are heterogeneous, composite structures that undergo large cyclic deformations during blood transport. Presence, build-up and consequent rupture of blockages in blood vessels, called atherosclerotic plaques, lead to disruption in the blood flow that can eventually be fatal. Abnormal lipid profile and hypertension are the main risk factors for plaque progression. Treatments span from pharmacological methods, to minimally invasive balloon angioplasty and stent procedures, and finally to surgical alternatives. There is a need to understand arterial disease progression and devise methods to detect, control, treat and manage arterial disease through early intervention. Local delivery through drug eluting stents also provide an attractive option for maintaining vessel integrity and restoring blood flow while releasing controlled amount of drug to reduce and alleviate symptoms. Development of drug eluting stents is hence interesting albeit challenging because it requires an integration of knowledge of mechanical properties with material transport of drug through the arterial wall to produce a desired biochemical effect. Although experimental models are useful in studying such complex multivariate phenomena, numerical models of mass transport in the vessel have proved immensely useful to understand and delineate complex interactions between chemical species, physical parameters and biological variables. The goals of this review are to summarize literature based on studies of mass transport involving low density lipoproteins in the arterial wall. We also discuss numerical models of drug elution from stents in layered and porous arterial walls that provide a unique platform that can be exploited for the design of novel drug eluting stents.

Nonlinear viscoelasticity of freestanding and polymer-anchored vertically aligned carbon nanotube foams

Vertical arrays of carbon nanotubes (VACNTs) show unique mechanical behavior in compression, with a highly nonlinear response similar to that of open cell foams and the ability to recover large deformations. Here, we study the viscoelastic response of both freestanding VACNT arrays and sandwich structures composed of a VACNT array partially embedded between two layers of poly(dimethylsiloxane) (PDMS) and bucky paper. The VACNTs tested are similar to 2 mm thick foams grown via an injection chemical vapor deposition method. Both freestanding and sandwich structures exhibit a time-dependent behavior under compression. A power-law function of time is used to describe the main features observed in creep and stress-relaxation tests. The power-law exponents show nonlinear viscoelastic behavior in which the rate of creep is dependent upon the stress level and the rate of stress relaxation is dependent upon the strain level. The results show a marginal effect of the thin PDMS/bucky paper layers on the viscoelastic responses. At high strain levels (epsilon - 0.8), the peak stress for the anchored CNTs reaches similar to 45 MPa, whereas it is only similar to 15MPa for freestanding CNTs, suggesting a large effect of PDMS on the structural response of the sandwich structures. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.3699184]