Improved hybrid elements for structural analysis

Hybrid elements, which are based on a two-field variational formulation with the displacements and stresses interpolated separately, are known to deliver very high accuracy, and to alleviate to a large extent problems of locking that plague standard displacement-based formulations. The choice of the stress interpolation functions is of course critical in ensuring the high accuracy and robustness of the method. Generally, an attempt is made to keep the stress interpolation to the minimum number of terms that will ensure that the stiffness matrix has no spurious zero-energy modes, since it is known that the stiffness increases with the increase in the number of terms. Although using such a strategy of keeping the number of interpolation terms to a minimum works very well in static problems, it results either in instabilities or fails to converge in transient problems. This is because choosing the stress interpolation functions merely on the basis of removing spurious energy modes can violate some basic principles that interpolation functions should obey. In this work, we address the issue of choosing the interpolation functions based on such basic principles of interpolation theory and mechanics. Although this procedure results in the use of more number of terms than the minimum (and hence in slightly increased stiffness) in many elements, we show that the performance continues to be far superior to displacement-based formulations, and, more importantly, that it also results in considerably increased robustness.

Curved composite beams—optimum lay-up for buckling by ranking

Instability of laminated curved composite beams made of repeated sublaminate construction is studied using finite element method. In repeated sublaminate construction, a full laminate is obtained by repeating a basic sublaminate which has a smaller number of plies. This paper deals with the determination of optimum lay-up for buckling by ranking of such composite curved beams (which may be solid or sandwich). For this purpose, use is made of a two-noded, 16 degress of freedom curved composite beam finite element. The displacements u, v, w of the element reference axis are expressed in terms of one-dimensional first-order Hermite interpolation polynomials, and line member assumptions are invoked in formulation of the elastic stiffness matrix and geometric stiffness matrix. The nonlinear expressions for the strains, occurring in beams subjected to axial, flexural and torsional loads, are incorporated in a general instability analysis. The computer program developed has been used, after extensive checking for correctness, to obtain optimum orientation scheme of the plies in the sublaminate so as to achieve maximum buckling load for typical curved solid/sandwich composite beams.

Frame Soil Interaction And Winkler Model

The Winkler spring model is the most convenient representation of soil support in the domain of linear elasticity for framed structure-soil interaction analyses. The closeness of the analytical results obtained using this model with those corresponding to the elastic half-space continuum has been investigated in the past for foundation beams. The findings, however, are not applicable to framed structures founded on beam or strip footings. Moreover, the past investigations employ the concept of characteristic length which does not adequately account for the stiffness contribution of the superstructure. A framed structure on beam foundation can be described parametrically by the ratios of stiffnesses of superstructure and foundation beams to that of soil. For a practical range of soil allowable pressures, the ranges of these relative stiffness ratios have been established. The present study examines the variation between interactive analyses based on Winkler springs with those using the half-space continuum over these ranges of relative stiffness ratios. The findings enable the analyst to undertake a Winkler spring-based-interaction analysis with knowledge of the likely variation of values with those derived for the more computation-intensive half-space continuum.

Characterization of a large-area PVDF thin film for electro-mechanical and ultrasonic sensing applications

Large-area PVDF thin films have been prepared and characterized for quasi-static and high frequency dynamic strain sensing applications. These films are prepared using hot press method and the piezoelectric phase (beta-phase) has been achieved by thermo-mechanical treatment and poling under DC field. The fabricated films have been characterized for quasi-static strain sensing and the linear strain-voltage relationship obtained is promising. In order to evaluate the ultrasonic sensing properties, a PZT wafer has been used to launch Lamb waves in a metal beam on which the PVDF film sensor is bonded at a distance. The voltage signals obtained from the PVDF films have been compared with another PZT wafer sensor placed on the opposite surface of the beam as a reference signal. Due to higher stiffness and higher thickness of the PZT wafer sensors, certain resonance patterns significantly degrade the sensor sensitivity curves. Whereas, the present results show that the large-area PVDF sensors can be superior with the signal amplitude comparable to that of PZT sensors and with no resonance-induced effect, which is due to low mechanical impedance, smaller thickness and larger area of the PVDF film. Moreover, the developed PVDF sensors are able to capture both A(0) and S-0 modes of Lamb wave, whereas the PZT sensors captures only A(0) mode in the same scale of voltage output. This shows promises in using large-area PVDF films with various surface patterns on structures for distributed sensing and structural health monitoring under quasi-static, vibration and ultrasonic situations. (C) 2010 Elsevier B.V. All rights reserved.

Metal-Insulator Transition in the Hubbard Model on a Triangular Lattice

We discuss the results of an extensive mean-field investigation of the half-filled Hubbard model on a triangular lattice at zero temperature. At intermediate U we find a first-order metal-insulator transition from an incommensurate spiral magnetic metal to a semiconducting state with a commensurate linear spin density wave ordering stabilized by the competition between the kinetic energy and the frustrated nature of the magnetic interaction. At large U the ground state is that of a classical triangular antiferromagnet within our approximation. In the incommensurate spiral metallic phase the Fermi surface has parts in which the wave function renormalization Z is extremely small. The evolution of the Fermi surface and the broadening of the quasi-particle band along with the variation of the plasma frequency and a charge stiffness constant with U/t are discussed.

Spiral states in the square-lattice Hubbard model

We present a variety of physical implications of a mean-field theory for spiral spin-density-wave states in the square-lattice Hubbard model for small deviations from half filling. The phase diagram with the paramagnetic metal, two spiral (semimetallic) states, and ferromagnet is calculated. The momentum distribution function and the (quasiparticle) density of states are discussed. There is a significant broadening of the quasiparticle bands when the antiferromagnetic insulator is doped. The evolution of the Fermi surface and the variation of the plasma frequency and a charge-stiffness constant with U/t and δ are calculated. The connection to results based on the Schwinger-boson-slave-fermion formalism is made.

Combination antiretroviral therapy -associated lipodystrophy : insights into pathogenesis and treatment

Introduction: Combination antiretroviral therapy (cART) has decreased morbidity and mortality of individuals infected with human immunodeficiency virus type 1 (HIV-1). Its use, however, is associated with adverse effects which increase the patients risk of conditions such as diabetes and coronary heart disease. Perhaps the most stigmatizing side effect is lipodystrophy, i.e., the loss of subcutaneous adipose tissue (SAT) in the face, limbs and trunk while fat accumulates intra-abdominally and dorsocervically. The pathogenesis of cART-associated lipodystrophy is obscure. Nucleoside reverse transcriptase inhibitors (NRTI) have been implicated to cause lipoatrophy via mitochondrial toxicity. There is no known effective treatment for cART-associated lipodystrophy during unchanged antiretroviral regimen in humans, but in vitro data have shown uridine to abrogate NRTI-induced toxicity in adipocytes. Aims: To investigate whether i) cART or lipodystrophy associated with its use affect arterial stiffness; ii) lipoatrophic SAT is inflamed compared to non-lipoatrophic SAT; iii) abdominal SAT from patients with compared to those without cART-associated lipoatrophy differs with respect to mitochondrial DNA (mtDNA) content, adipose tissue inflammation and gene expression, and if NRTIs stavudine and zidovudine are associated with different degree of changes; iv) lipoatrophic abdominal SAT differs from preserved dorsocervical SAT with respect to mtDNA content, adipose tissue inflammation and gene expression in patients with cART-associated lipodystrophy and v) whether uridine can revert lipoatrophy and the associated metabolic disturbances in patients on stavudine or zidovudine based cART. Subjects and methods: 64 cART-treated patients with (n=45) and without lipodystrophy/-atrophy (n=19) were compared cross-sectionally. A marker of arterial stiffness, heart rate corrected augmentation index (AgIHR), was measured by pulse wave analysis. Body composition was measured by magnetic resonance imaging and dual-energy X-ray absorptiometry, and liver fat content by proton magnetic resonance spectroscopy. Gene expression and mtDNA content in SAT were assessed by real-time polymerase chain reaction and microarray. Adipose tissue composition and inflammation were assessed by histology and immunohistochemistry. Dorsocervical and abdominal SAT were studied. The efficacy and safety of uridine for the treatment of cART-associated lipoatrophy were evaluated in a randomized, double-blind, placebo-controlled 3-month trial in 20 lipoatrophic cART-treated patients. Results: Duration of antiretroviral treatment and cumulative exposure to NRTIs and protease inhibitors, but not the presence of cART-associated lipodystrophy, predicted AgIHR independent of age and blood pressure. Gene expression of inflammatory markers was increased in SAT of lipodystrophic as compared to non-lipodystrophic patients. Expression of genes involved in adipogenesis, triglyceride synthesis and glucose disposal was lower and of those involved in mitochondrial biogenesis, apoptosis and oxidative stress higher in SAT of patients with than without cART-associated lipoatrophy. Most changes were more pronounced in stavudine-treated than in zidovudine-treated individuals. Lipoatrophic SAT had lower mtDNA than SAT of non-lipoatrophic patients. Expression of inflammatory genes was lower in dorsocervical than in abdominal SAT. Neither depot had characteristics of brown adipose tissue. Despite being spared from lipoatrophy, dorsocervical SAT of lipodystrophic patients had lower mtDNA than the phenotypically similar corresponding depot of non-lipodystrophic patients. The greatest difference in gene expression between dorsocervical and abdominal SAT, irrespective of lipodystrophy status, was in expression of homeobox genes that regulate transcription and regionalization of organs during embryonal development. Uridine increased limb fat and its proportion of total fat, but had no effect on liver fat content and markers of insulin resistance. Conclusions: Long-term cART is associated with increased arterial stiffness and, thus, with higher cardiovascular risk. Lipoatrophic abdominal SAT is characterized by inflammation, apoptosis and mtDNA depletion. As mtDNA is depleted even in non-lipoatrophic dorsocervical SAT, lipoatrophy is unlikely to be caused directly by mtDNA depletion. Preserved dorsocervical SAT of patients with cART-associated lipodystrophy is less inflamed than their lipoatrophic abdominal SAT, and does not resemble brown adipose tissue. The greatest difference in gene expression between dorsocervical and abdominal SAT is in expression of transcriptional regulators, homeobox genes, which might explain the differential susceptibility of these adipose tissue depots to cART-induced toxicity. Uridine is able to increase peripheral SAT in lipoatrophic patients during unchanged cART.

Temperature-Sensitive Src-Transformed Mdck Cells nn 3d Cultures as a Model of Malignant Transformation of Epithelial Cells

The equilibrium between cell proliferation, differentiation, and apoptosis is crucial for maintaining homeostasis in epithelial tissues. In order for the epithelium to function properly, individual cells must gain normal structural and functional polarity. The junctional proteins have an important role both in binding the cells together and in taking part in cell signaling. Cadherins form adherens junctions. Cadherins initiate the polarization process by first recognizing and binding the neighboring cells together, and then guiding the formation of tight junctions. Tight junctions form a barrier in dividing the plasma membranes to apical and basolateral membrane domains. In glandular tissues, single layered and polarized epithelium is folded into tubes or spheres, in which the basal side of the epithelial layer faces the outer basal membrane, and the apical side the lumen. In carcinogenesis, the differentiated architecture of an epithelial layer is disrupted. Filling of the luminal space is a hallmark of early epithelial tumors in tubular and glandular structures. In order for the transformed tumor cells to populate the lumen, enhanced proliferation as well as inhibition of apoptosis is required. Most advances in cancer biology have been achieved by using two-dimensional (2D) cell culture models, in which the cells are cultured on flat surfaces as monolayers. However, the 2D cultures are limited in their capacity to recapitulate the structural and functional features of tubular structures and to represent cell growth and differentiation in vivo. The development of three-dimensional (3D) cell culture methods enables the cells to grow and to be studied in a more natural environment. Despite the wide use of 2D cell culture models and the development of novel 3D culture methods, it is not clear how the change of the dimensionality of culture conditions alters the polarization and transformation process and the molecular mechanisms behind them. Src is a well-known oncogene. It is found in focal and adherens junctions of cultured cells. Active src disrupts cell-cell junctions and interferes with cell-matrix binding. It promotes cell motility and survival. Src transformation in 2D disrupts adherens junctions and the fibroblastic phenotype of the cells. In 3D, the adherens junctions are weakened, and in glandular structures, the lumen is filled with nonpolarized vital cells. Madin-Darby canine kidney (MDCK) cells are an epithelial cell type commonly used as a model for cell polarization. Its-src-transformed variants are useful model systems for analyzing the changes in cell morphology, and they play a role in src-induced malignant transformation. This study investigates src-transformed cells in 3D cell cultures as a model for malignant transformation. The following questions were posed. Firstly: What is the role of the composition and stiffness of the extracellular matrix (ECM) on the polarization and transformation of ts v-src MDCK cells in 3D cell cultures? Secondly: How do the culture conditions affect gene expression? What is the effect of v-src transformation in 2D and in 3D cell models? How does the shift from 2D to 3D affect cell polarity and gene expression? Thirdly: What is the role of survivin and its regulator phosphatase and tensin homolog protein (PTEN) in cell polarization and transformation, and in determining cell fate? How does their expression correlate with impaired mitochondrial function in transformed cells? In order to answer the above questions, novel methods of culturing and monitoring cells had to be created: novel 3D methods of culturing epithelial cells were engineered, enabling real time monitoring of a polarization and transformation process, and functional testing of 3D cell cultures. Novel 3D cell culture models and imaging techniques were created for the study. Attention was focused especially on confocal microscopy and live-cell imaging. Src-transformation disturbed the polarization of the epithelium by disrupting cell adhesion, and sensitized the cells to their environment. With active src, the morphology of the cell cluster depended on the composition and stiffness of the matrix. Gene expression studies revealed a broader impact of src transformation than mere continuous activity of src-kinase. In 2D cultures, src transformation altered the expression of immunological, actin cytoskeleton and extracellular matrix (ECM). In 3D, the genes regulating cell division, inhibition of apoptosis, cell metabolism, mitochondrial function, actin cytoskeleton and mechano-sensing proteins were altered. Surprisingly, changing the culture conditions from 2D to 3D affected also gene expression considerably. The microarray hit survivin, an inhibitor of apoptosis, played a crucial role in the survival and proliferation of src-transformed cells.

Analysis of continuum using boundary compatibility conditions of integrated force method

Static and vibration problems of an indeterminate continuum are traditionally analyzed by the stiffness method. The force method is more or less non-existent for such problems. This situation is primarily due to the incomplete state of development of the compatibility conditions which are essential for the analysis of indeterminate structures by the flexibility method. The understanding of the Compatibility Conditions (CC) has been substantially augmented. Based on the understanding of CC, a novel formulation termed the Integrated Force Method (IFM) has been established. In this paper IFM has been extended for the static and vibration analyses of a continuum. The IFM analysis is illustrated taking three examples: 1. (1) rectangular plate in flexure 2. (2) analysis of a cantilevered dam 3. (3) free vibration analysis of a beam. From the examples solved it is observed that the force response of an indeterminate continuum with mixed boundary conditions can be generated by IFM without any reference to displacements in the field or on the boundary. Displacements if required can be calculated by back substitution.

Influence of an additional discrete elastic support on the stability of a Leipholz column

The stability characteristics of a conservatively loaded structure are expected to improve if additional supports are provided to the structure. The same, however, may not be said of a non-conservatively loaded structure; several factors, such as the location and stiffness of supports, type of structure and loading, have a significant influence on the stability characteristics. The influence of an arbitrarily located elastic support on the stability characteristics of a Leipholz column is examined.

Method of Analysis if Uniaxial Deformation Data

A comprehensive scheme for analysing uniaxial deformation data, taking into account the finite stiffness of the testing machine is presented. Equations relevant to tension and stress relaxation tests carried out under cross head speed control, and to creep testing under constant load, are described. For the first two cases, the implications of not using gauge length extensometry but relying upon cross head displacement for inferring specimen extension, and the role of uncertainty in machine stiffness are also examined. The final section touches upon the extension of the present scheme to account for specimen anelasticity.

Parametric instability of stochastic columns

Columns which have stochastically distributed Young's modulus and mass density and are subjected to deterministic periodic axial loadings are considered. The general case of a column supported on a Winkler elastic foundation of random stiffness and also on discrete elastic supports which are also random is considered. Material property fluctuations are modeled as independent one-dimensional univariate homogeneous real random fields in space. In addition to autocorrelation functions or their equivalent power spectral density functions, the input random fields are characterized by scale of fluctuations or variance functions for their second order properties. The foundation stiffness coefficient and the stiffnesses of discrete elastic supports are treated to constitute independent random variables. The system equations of boundary frequencies are obtained using Bolotin's method for deterministic systems. Stochastic FEM is used to obtain the discrete system with random as well as periodic coefficients. Statistical properties of boundary frequencies are derived in terms of input parameter statistics. A complete covariance structure is obtained. The equations developed are illustrated using a numerical example employing a practical correlation structure.

A monolithic strategy for fluid-structure interaction problems

In this work, we present a new monolithic strategy for solving fluid-structure interaction problems involving incompressible fluids, within the context of the finite element method. This strategy, similar to the continuum dynamics, conserves certain properties, and thus provides a rational basis for the design of the time-stepping strategy; detailed proofs of the conservation of these properties are provided. The proposed algorithm works with displacement and velocity variables for the structure and fluid, respectively, and introduces no new variables to enforce velocity or traction continuity. Any existing structural dynamics algorithm can be used without change in the proposed method. Use of the exact tangent stiffness matrix ensures that the algorithm converges quadratically within each time step. An analytical solution is presented for one of the benchmark problems used in the literature, namely, the piston problem. A number of benchmark problems including problems involving free surfaces such as sloshing and the breaking dam problem are used to demonstrate the good performance of the proposed method. Copyright (C) 2010 John Wiley & Sons, Ltd.

Dense granular flow down an inclined plane A comparison between the hard particle model and soft particle simulations

The granular flow down an inclined plane is simulated using the discrete element (DE) technique to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model First, we examine the average coordination number for the particles in the flow down an inclined plane using the DE technique using the linear contact model with and without friction, and the Hertzian contact model with friction The simulations show that the average coordination number decreases below 1 for values of the spring stiffness corresponding to real materials, such as sand and glass, even when the angle of inclination is only 10 larger than the angle of repose Additional measures of correlations in the system, such as the fraction of particles with multibody contact, the force ratio (average ratio of the magnitudes of the largest and the second largest force on a particle), and the angle between the two largest forces on the particle, show no evidence of force chains or other correlated motions in the system An analysis of the bond-orientational order parameter indicates that the flow is in the random state, as in event-driven (ED) simulations V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are compared with the theory V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] and are found to be in quantitative agreement In addition, we also conduct a comparison of the collision frequency and the distribution of the precollisional relative velocities of particles in contact The strong correlation effects exhibited by these two quantities in event-driven simulations V Kumaran, J Fluid Mech 632, 145 (2009)] are also found in the DE simulations (C) 2010 American Institute of Physics doi 10 1063/1 3504660]

Exact Solution of a One-Dimensional Multicomponent Lattice Gas with Hyperbolic Interaction

We present the exact solution to a one-dimensional multicomponent quantum lattice model interacting by an exchange operator which falls off as the inverse sinh square of the distance. This interaction contains a variable range as a parameter and can thus interpolate between the known solutions for the nearest-neighbor chain and the inverse-square chain. The energy, susceptibility, charge stiffness, and the dispersion relations for low-lying excitations are explicitly calculated for the absolute ground state, as a function of both the range of the interaction and the number of species of fermions.