Random vibration of a second order non-linear elastic system

A method is presented for obtaining, approximately, the response covariance and probability distribution of a non-linear oscillator under a Gaussian excitation. The method has similarities with the hierarchy closure and the equivalent linearization approaches, but is different. A Gaussianization technique is used to arrive at the output autocorrelation and the input-output cross-correlation. This along with an energy equivalence criterion is used to estimate the response distribution function. The method is applicable in both the transient and steady state response analysis under either stationary or non-stationary excitations. Good comparison has been observed between the predicted and the exact steady state probability distribution of a Duffing oscillator under a white noise input.

A peridynamic theory for linear elastic shells

A state-based peridynamic formulation for linear elastic shells is presented. The emphasis is on introducing, possibly for the first time, a general surface based peridynamic model to represent the deformation characteristics of structures that have one geometric dimension much smaller than the other two. A new notion of curved bonds is exploited to cater for force transfer between the peridynamic particles describing the shell. Starting with the three dimensional force and deformation states, appropriate surface based force, moment and several deformation states are arrived at. Upon application on the curved bonds, such states yield the necessary force and deformation vectors governing the motion of the shell. By incorporating a shear correction factor, the formulation also accommodates analysis of shells that have higher thickness. In order to attain this, a consistent second order approximation to the complementary energy density is considered and incorporated in peridynamics via constitutive correspondence. Unlike the uncoupled constitution for thin shells, a consequence of a first order approximation, constitutive relations for thick shells are fully coupled in that surface wryness influences the in-plane stress resultants and surface strain the moments. Our proposal on the peridynamic shell theory is numerically assessed against simulations on static deformation of spherical and cylindrical shells, that of flat plates and quasi-static fracture propagation in a cylindrical shell. (C) 2016 Elsevier Ltd. All rights reserved.

Bending Of Circular Plate On Elastic-Foundation

The governing differential equation of linear, elastic, thin, circular plate of uniform thickness, subjected to uniformly distributed load and resting on Winkler-Pasternak type foundation is solved using ``Chebyshev Polynomials''. Analysis is carried out using Lenczos' technique, both for simply supported and clamped plates. Numerical results thus obtained by perturbing the differential equation for plates without foundation are compared and are found to be in good agreement with the available results. The effect of foundation on central deflection of the plate is shown in the form of graphs.

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.

A micropolar peridynamic theory in linear elasticity

A state-based micropolar peridynamic theory for linear elastic solids is proposed. The main motivation is to introduce additional micro-rotational degrees of freedom to each material point and thus naturally bring in the physically relevant material length scale parameters into peridynamics. Non-ordinary type modeling via constitutive correspondence is adopted here to define the micropolar peridynamic material. Along with a general three dimensional model, homogenized one dimensional Timoshenko type beam models for both the proposed micropolar and the standard non-polar peridynamic variants are derived. The efficacy of the proposed models in analyzing continua with length scale effects is established via numerical simulations of a few beam and plane-stress problems. (C) 2015 Elsevier Ltd. All rights reserved.

Computation of Thermal Stress Intensity Factors for Bimaterial Interface Cracks Using Domain Integral Method

The concept of domain integral used extensively for J integral has been applied in this work for the formulation of J(2) integral for linear elastic bimaterial body containing a crack at the interface and subjected to thermal loading. It is shown that, in the presence of thermal stresses, the J(k) domain integral over a closed path, which does not enclose singularities, is a function of temperature and body force. A method is proposed to compute the stress intensity factors for bimaterial interface crack subjected to thermal loading by combining this domain integral with the J(k) integral. The proposed method is validated by solving standard problems with known solutions.

Prediction of fatigue strength in plain and reinforced concrete beams

A fatigue crack propagation model for concrete is proposed based on the concepts of fracture mechanics. This model takes into account the loading history, frequency of applied load, and size, effect parameters. Using this model, a method is described based on linear elastic fracture mechanics to assess the residual strength of cracked plain and reinforced concrete (RC) beams. This could be used to predict the residual strength (load carrying capacity) of cracked or damaged plain and reinforced concrete beams at a given level of damage. It has been seen that the fatigue crack propagation rate increases as. the size of plain concrete, beam increases indicating an increase in brittleness. In reinforced concrete (RC) beams, the fracture process becomes stable only when the beam is sufficiently reinforced.

Fracture analysis of stiffened panels under combined tensile, bending, and shear loads

The objective is to present the formulation of numerically integrated modified virtual crack closure integral technique for concentrically and eccentrically stiffened panels for computation of strain-energy release rate and stress intensity factor based on linear elastic fracture mechanics principles. Fracture analysis of cracked stiffened panels under combined tensile, bending, and shear loads has been conducted by employing the stiffened plate/shell finite element model, MQL9S2. This model can be used to analyze plates with arbitrarily located concentric/eccentric stiffeners, without increasing the total number of degrees of freedom, of the plate element. Parametric studies on fracture analysis of stiffened plates under combined tensile and moment loads have been conducted. Based on the results of parametric,studies, polynomial curve fitting has been carried out to get best-fit equations corresponding to each of the stiffener positions. These equations can be used for computation of stress intensity factor for cracked stiffened plates subjected to tensile and moment loads for a given plate size, stiffener configuration, and stiffener position without conducting finite element analysis.

Adhesively-bonded Patch Repair with Composites

Adhesively-bonded composite patch repairs over cracked or corrosion-damaged metallic aircraft structures have shown great promise for extending life of ageing structures. This study presents the numerical investigation into the interface behaviour of adhesively-bonded cracked aluminum alloy substrate patched with fibre-reinforced composite material. The adhesive is modelled as an elasto-plastic bilinear material to characterise the debond behaviour, while the defective substrate is regarded as linear elastic continuum. Two typical patch shapes were selected based on information available in the literature. Geometric and material nonlinear analyses for square and octagonal patches were performed to capture peel and shear stresses developed between the substrate and the patch to examine the possibility of interface delamination/debonding. Parametric studies on adhesive thickness and patch thickness were carried out to predict their infuence on damage tolerance of repaired structures.

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.

Adhesively-bonded Patch Repair with Composites

Adhesively-bonded composite patch repairs over cracked or corrosion-damaged metallic aircraft structures have shown great promise for extending life of ageing structures. This study presents the numerical investigation into the interface behaviour of adhesively-bonded cracked aluminum alloy substrate patched with fibre-reinforced composite material. The adhesive is modelled as an elasto-plastic bilinear material to characterise the debond behaviour, while the defective substrate is regarded as linear elastic continuum. Two typical patch shapes were selected based on information available in the literature. Geometric and material nonlinear analyses for square and octagonal patches were performed to capture peel and shear stresses developed between the substrate and the patch to examine the possibility of interface delamination/debonding. Parametric studies on adhesive thickness and patch thickness were carried out to predict their infuence on damage tolerance of repaired structures.

Integrating strong and weak discontinuities without integration subcells and example applications in an XFEM/GFEM framework

Partition of unity methods, such as the extended finite element method, allows discontinuities to be simulated independently of the mesh (Int. J. Numer. Meth. Engng. 1999; 45:601-620). This eliminates the need for the mesh to be aligned with the discontinuity or cumbersome re-meshing, as the discontinuity evolves. However, to compute the stiffness matrix of the elements intersected by the discontinuity, a subdivision of the elements into quadrature subcells aligned with the discontinuity is commonly adopted. In this paper, we use a simple integration technique, proposed for polygonal domains (Int. J. Nuttier Meth. Engng 2009; 80(1):103-134. DOI: 10.1002/nme.2589) to suppress the need for element subdivision. Numerical results presented for a few benchmark problems in the context of linear elastic fracture mechanics and a multi-material problem show that the proposed method yields accurate results. Owing to its simplicity, the proposed integration technique can be easily integrated in any existing code. Copyright (C) 2010 John Wiley & Sons, Ltd.

An ensemble of B-DNA dinucleotide geometries lead to characteristic nucleosomal DNA structure and provide plasticity required for gene expression

Background: A nucleosome is the fundamental repeating unit of the eukaryotic chromosome. It has been shown that the positioning of a majority of nucleosomes is primarily controlled by factors other than the intrinsic preference of the DNA sequence. One of the key questions in this context is the role, if any, that can be played by the variability of nucleosomal DNA structure. Results: In this study, we have addressed this question by analysing the variability at the dinucleotide and trinucleotide as well as longer length scales in a dataset of nucleosome X-ray crystal structures. We observe that the nucleosome structure displays remarkable local level structural versatility within the B-DNA family. The nucleosomal DNA also incorporates a large number of kinks. Conclusions: Based on our results, we propose that the local and global level versatility of B-DNA structure may be a significant factor modulating the formation of nucleosomes in the vicinity of high-plasticity genes, and in varying the probability of binding by regulatory proteins. Hence, these factors should be incorporated in the prediction algorithms and there may not be a unique `template' for predicting putative nucleosome sequences. In addition, the multimodal distribution of dinucleotide parameters for some steps and the presence of a large number of kinks in the nucleosomal DNA structure indicate that the linear elastic model, used by several algorithms to predict the energetic cost of nucleosome formation, may lead to incorrect results.

A general procedure for modified crack closure integral in 3D problems with cracks

In linear elastic fracture mechanics (LEFM), Irwin's crack closure integral (CCI) is one of the signficant concepts for the estimation of strain energy release rates (SERR) G, in individual as well as mixed-mode configurations. For effective utilization of this concept in conjunction with the finite element method (FEM), Rybicki and Kanninen [Engng Fracture Mech. 9, 931 938 (1977)] have proposed simple and direct estimations of the CCI in terms of nodal forces and displacements in the elements forming the crack tip from a single finite element analysis instead of the conventional two configuration analyses. These modified CCI (MCCI) expressions are basically element dependent. A systematic derivation of these expressions using element stress and displacement distributions is required. In the present work, a general procedure is given for the derivation of MCCI expressions in 3D problems with cracks. Further, a concept of sub-area integration is proposed which facilitates evaluation of SERR at a large number of points along the crack front without refining the finite element mesh. Numerical data are presented for two standard problems, a thick centre-cracked tension specimen and a semi-elliptical surface crack in a thick slab. Estimates for the stress intensity factor based on MCCI expressions corresponding to eight-noded brick elements are obtained and compared with available results in the literature.