Modern Domain-based Discretization Methods for Damage and Fracture

Standard domain-based discretization methods that have been developed for continuous media are not well suited for treating propagating (or evolving) discontinuities. Indeed, they are approximation methods for the solution of partial differential equations, which are valid on a domain. Discontinuities divide this domain into two or more parts. Conventionally, special interface elements methods are placed a priori between the continuum finite elements to capture discontinuities at locations where they are expected to emerge. More recently, discretization methods have been proposed, which are more flexible than standard finite element methods, while having the potential to capture propagating discontinuities in a robust, efficient and accurate manner. Examples are meshfree methods, finite element methods that exploit the partition-of-unity property of finite element shape functions, and discontinuous Galerkin methods. In this contribution, we shall present an overview of these novel discretization techniques for capturing propagating discontinuities, including a comparison of their similarities and differences.     

Correlated Shadow-Fading in Wireless Networks and its Effect on Call Dropping

We discuss a statistical model to generate correlated shadow-fading patterns for wireless systems in the absence of detailed propagation and landscape information. The currently available autocorrelation models result in anomalous effects that depend on traffic density and mobility, as they propose independent random processes for each mobile. Our approach involves generating a pre-computed fading map with the right marginal distributions and spatial correlations, which avoids inconsistencies such as providing widely differing values for mobiles close to each other. The correlations are introduced via a Gaussian random field, which has a covariance structure that depends on a set of parameters which can be computed from local measurements. The model is efficiently implemented using standard linear-algebra methods. We conclude by describing a simulation experiment to study the effect of correlations on call dropping. The experiment reveals a strong relationship between call dropping and the correlation length of the fading pattern, and indicates circumstances under which dropping may be relatively high.     

Generalized processor sharing with light-tailed and heavy-tailed input

We consider a queue fed by a mixture of light-tailed and heavy-tailed traffic. The two traffic flows are served in accordance with the generalized processor sharing (GPS) discipline. GPS-based scheduling algorithms, such as weighted fair queueing, have emerged as an important mechanism for achieving service differentiation in integrated networks. We derive the asymptotic workload behavior of the light-tailed traffic flow under the assumption that its GPS weight is larger than its traffic intensity. The GPS mechanism ensures that the workload is bounded above by that in an isolated system with the light-tailed flow served in isolation at a constant rate equal to its GPS weight. We show that the workload distribution is in fact asymptotically equivalent to that in the isolated system, multiplied with a certain pre-factor, which accounts for the interaction with the heavy-tailed flow. Specifically, the pre-factor represents the probability that the heavy-tailed flow is backlogged long enough for the light-tailed flow to reach overflow. The results provide crucial qualitative insight in the typical overflow scenario.     

Arterial remodeling after balloon angioplasty or stenting in an atherosclerotic experimental model

The actions of the endogenous ORL1-receptor ligand nociceptin on the membrane properties and synaptic currents in rat periaqueductal gray (PAG) neurons were examined by the use of whole-cell patch-clamp recording in brain slices. Nociceptin produced an outward current in all neurons tested, with an EC50 of 39 +/- 7 nM. The outward current was unaffected by naloxone. Outward currents reversed polarity at -110 +/- 3 mV in 2.5 mM extracellular potassium, and the reversal potential increased when the extracellular potassium concentration was raised (slope = 66.3 mV/log[K+]o mM). Thus, the nociceptin-induced outward current was attributable to an increased K+ conductance. Nociceptin inhibited evoked fast GABAergic (IP-SCs) and glutamatergic (EPSCs) postsynaptic currents and increased paired-pulse facilitation in a subpopulation of PAG neurons. Nociceptin inhibited evoked IPSCs and EPSCs in approximately 50% of neurons throughout the PAG, except in the ventrolateral PAG, where nociceptin inhibited evoked IPSCs in most neurons. Nociceptin decreased the frequency of spontaneous miniature postsynaptic currents (mIPSCs and mEPSCs) in a subpopulation of PAG neurons but had no effect on their amplitude distributions. Thus, nociceptin had a presynaptic inhibitory effect on transmitter release. These findings suggest that nociceptin, via its pre- and postsynaptic actions, has the potential to modulate the analgesic, behavioral, and autonomic functions of the PAG.     

Multidrug resistance protein 1 protects the oropharyngeal mucosal layer and the testicular tubules against drug-induced damage

The multidrug resistance protein 1 (MRP1) gene encodes a transporter protein that helps to protect cells against xenobiotics. Elevated levels of MRP1 in tumor cells can result in active extrusion of a wide range of (anticancer) drugs with different cellular targets, a phenomenon called multidrug resistance (MDR). To explore the protective function of the mouse mrp1 protein during drug treatment, we investigated the toxicity caused by the anticancer drug etoposide-phosphate (ETOPOPHOS) in mice lacking the mrp1 gene (mrp1(-/-) mice). We show here that the lack of mrp1 protein results in increased etoposide-induced damage to the mucosa of the oropharyngeal cavity and to the seminiferous tubules of the testis. The high concentrations of mrp1 that we find in the basal layers of the oropharyngeal mucosa and in the basal membrane of the Sertoli cells in the testis apparently protect wild-type mice against this tissue damage. We also find drug-induced polyuria in mrp1(-/-) mice, which correlates with the presence of mrp1 protein in the urinary collecting tubules, the major site of kidney water reabsorption. Our results indicate that specific inhibitors of MRP1 used to reverse MDR, in combination with carcinostatic drugs transported by MRP1, might lead to drug-induced mucositis, (temporary) infertility, and diabetes insipidus.     

The use of the Hoffman yield criterion in finite element analysis of anisotropic composites

A return mapping algorithm has been developed for the Hoffman yield function of anisotropic plasticity. The accuracy of the algorithm has been assessed by means of iso-error maps for trial stress increments in the deviatoric and volumetric plane. A tangent operator that is consistent with the developed integration algorithm has been formulated. The Hoffman model has been applied to a plate structure and to two shell structures.     

Analysis of fracture and delamination in laminates using 3D numerical modelling

The static failure behaviour of the fibre-metal laminate GLARE is examined using 3D finite element simulations. The configuration analysed is a centre-cracked tensile specimen composed of two aluminium layers sandwiching a cross-plied, fibre-epoxy layer. The crack and delamination growths are simulated by means of interface elements equipped with a mixed-mode damage model. The mode-mixity is derived from an energy criterion typically used in linear elastic fracture mechanics studies. The damage kinetic law is rate-dependent, in order to simulate rate effects during interfacial delamination and to avoid numerical convergence problems due to crack bifurcations. The numerical implementation of the interface damage model is based on a backward Euler approach. In the boundary value problem studied, the failure responses of GLARE specimens containing elastic aluminium layers and elasto-plastic aluminium layers are compared. The development of plastic deformations in the aluminium layers stabilizes the effective failure response, and increases the residual strength of the laminate. For a ‘quasi-brittle’ GLARE specimen with elastic aluminium layers, the residual strength is governed by the toughness for interfacial delamination, and is in close correspondence with the residual strength obtained from a closed-form expression derived from energy considerations. Conversely, for a ‘ductile’ GLARE specimen with elasto-plastic aluminium layers, the residual strength is also determined by the relation between the fracture strength and the yield strength of the aluminium. The amount of constraint by the horizontal displacements at the vertical specimen edges has a moderate to small influence on the residual strength. Furthermore, the ultimate laminate strength is lower for a larger initial crack length, and shows to be in good correspondence with experimental values.     

A unified framework for concrete damage and fracture models including size effects

A unified approach is given for isotropic and anisotropic damage formulations for concrete fracture. The formulation encompasses the classical fixed and rotating smeared crack models, but also more refined approaches based on the microplane concept. Higher-order strain gradients are introduced to avoid the boundary value problem from becoming ill-posed at the onset of softening. Again, this has been accomplished in a unified setting. By analysing a geometrically identical anchor bolt of three different sizes it is shown that higher-order strain gradients also introduce a proper size effect in the model. This revised version was published online in July 2006 with corrections to the Cover Date.     

Excimer laser induced bubble: dimensions, theory, and implications for laser angioplasty

BACKGROUND AND OBJECTIVE: Previous studies have demonstrated that during Xenon-Chloride excimer laser ablation of tissue, rapidly expanding and imploding bubbles (diameter < 3 mm), predominantly containing water vapor, are formed. These short lived bubbles (life time < 300 microseconds) induce mechanical damage in adjacent tissue. In the present study, a theoretical analysis of the volume of vaporized water is correlated with measured bubble volumes formed in hemoglobin solution. STUDY DESIGN/MATERIALS AND METHODS: The dimensions of the rapidly expanding and imploding vapor bubble induced by the XeCl excimer laser pulses (308 nm, 115 ns), delivered via a 300, 550, or 950 microns diameter monofiber in 16% w/v hemoglobin solution (at 37 degrees C), were measured. RESULTS: Theoretical analysis and the experimental data correlated well (correlation coefficient r = 0.97). The diameter of excimer laser induced bubbles increased with increasing pulse energy. For a given radiant exposure, the bubble size was decreased by either decreasing the fiber tip area or by decreasing the absorption coefficient of the hemoglobin solution. CONCLUSION: We conclude that, for a wide range of conditions, theory agrees well with experimental data. Thus, during delivery of excimer laser pulses in blood, bubble dimensions can be reduced by flushing with saline or by reduction of the area radiated with each laser pulse, for example, by pulse multiplexing or using a smaller multifiber catheter.