Finite element analysis of moving contact in mechanically fastened joints

The Finite Element Method (FEM) has made a number of otherwise intractable problems solvable. An important aspect for achieving an economical and accurate solution through FEM is matching the formulation and the computational organisation to the problem. This was realised forcefully in the present case of the solution of a class of moving contact boundary value problems of fastener joints. This paper deals with the problem of changing contact at the pin-hole interface of a fastener joint. Due to moving contact, the stresses and displacements are nonlinear with load. This would, in general, need an interactive-incremental approach for solution. However, by posing the problem in an inverse way, a solution is sought for obtaining loads to suit given contact configuration. Numerical results are given for typical isotropic and composite plates with rigid pins. Two cases of loading are considered: (i) load applied only at the edges of the plate and (ii) load applied at the pin and reacted at a part of the edge of the plate. Load-contact relationships, compliance and stress-patterns are investigated. This paper clearly demonstrates the simplification achieved by a suitable formulation of the problem. The results are of significance to the design and analysis of fastener joints.

Dynamics of dense sheared granular flows. Part 1. Structure and diffusion

Shear flows of inelastic spheres in three dimensions in the Volume fraction range 0.4-0.64 are analysed using event-driven simulations.Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to the line Joining the centres). Here, we have considered both e(t) = +1 and e(t) = e(n) (rough particles) and e(t) =-1 (smooth particles), and the normal coefficient of restitution e(n) was varied in the range 0.6-0.98. Care was taken to avoid inelastic collapse and ensure there are no particle overlaps during the simulation. First, we studied the ordering in the system by examining the icosahedral order parameter Q(6) in three dimensions and the planar order parameter q(6) in the plane perpendicular to the gradient direction. It was found that for shear flows of sufficiently large size, the system Continues to be in the random state, with Q(6) and q(6) close to 0, even for volume fractions between phi = 0.5 and phi = 0.6; in contrast, for a system of elastic particles in the absence of shear, the system orders (crystallizes) at phi = 0.49. This indicates that the shear flow prevents ordering in a system of sufficiently large size. In a shear flow of inelastic particles, the strain rate and the temperature are related through the energy balance equation, and all time scales can be non-dimensionalized by the inverse of the strain rate. Therefore, the dynamics of the system are determined only by the volume fraction and the coefficients of restitution. The variation of the collision frequency with volume fraction and coefficient of estitution was examined. It was found, by plotting the inverse of the collision frequency as a function of volume fraction, that the collision frequency at constant strain rate diverges at a volume fraction phi(ad) (volume fraction for arrested dynamics) which is lower than the random close-packing Volume fraction 0.64 in the absence of shear. The volume fraction phi(ad) decreases as the coefficient of restitution is decreased from e(n) = 1; phi(ad) has a minimum of about 0.585 for coefficient of restitution e(n) in the range 0.6-0.8 for rough particles and is slightly larger for smooth particles. It is found that the dissipation rate and all components of the stress diverge proportional to the collision frequency in the close-packing limit. The qualitative behaviour of the increase in the stress and dissipation rate are well Captured by results derived from kinetic theory, but the quantitative agreement is lacking even if the collision frequency obtained from simulations is used to calculate the pair correlation function used In the theory.

Fatty acids, fibre, carotenoids and tocopherols in relation to glucose metabolism in subjects at high risk for type 2 diabetes : a cross-sectional analysis

Fatty acids, fibre, carotenoids and tocopherols in relation to glucose metabolism in subjects at high risk for type 2 diabetes a cross-sectional analysis Type 2 diabetes (T2D) is a heterogeneous disorder of carbohydrate, lipid and protein metabolism, resulting from genetics, environmental influences and interactions between these. The disease is characterized by insulin resistance, β-cell dysfunction, hepatic glucose overproduction and disordered fat mobilization and storage. The literature on associations between dietary factors and glucose metabolism is inconsistent. One factor behind the discrepant results may be genetic heterogeneity of study populations. Data on nutrient-gene interactions in relation to glucose metabolism are scarce. Thus, investigating high-risk populations and exploring nutrient-gene interactions are essential for improving the understanding of T2D aetiology. Ideally, this information could help to develop prevention programmes that take into account the genetic predisposition to the disease. In this study, associations between measures of glucose metabolism predicting T2D and fatty acids, antioxidative nutrients and fibre were examined in a high-risk population, i.e., in non-diabetic relatives of affected patients. Interactions between the PPARG Pro12Ala polymorphism and fatty acids on glucose metabolism were taken into consideration. This common polymorphism plays an important role in the regulation of glucose metabolism. The inverse associations observed between dietary fibre and insulin resistance are consistent with the prevailing recommendations urging increased intake of fibre to prevent T2D. Beneficial associations observed between the intake of carotenoids and glucose levels stress that a high consumption of vegetables, fruits and berries rich in carotenoids might also play a role in the prevention of T2D. Whether tocopherols have an independent association with glucose metabolism remains questionable. Observed interactions between fatty acids and glucose metabolism suggest that a high intake of palmitic acid is associated with high fasting glucose levels mainly in female Ala allele carriers. Furthermore, the PPARG Pro12Ala polymorphism may modify the metabolic response to dietary marine fat. The beneficial associations of high intake of marine n 3 fatty acids with insulin resistance and glucose levels may be restricted to carriers of the Ala allele. The findings pertain to subjects with a family history of T2D, and the cross-sectional nature of the study precludes inferences about causality. Results nevertheless show that associations of dietary factors with glucose metabolism may be modulated by the genetic makeup of an individual. Additional research is warranted to elucidate the role of probably numerous nutrient-gene interactions, some of which may be sex-specific, in the aetiology of T2D.

Finite element analysis of moving contact in mechanically fastened joints

The Finite Element Method (FEM) has made a number of otherwise intractable problems solvable. An important aspect for achieving an economical and accurate solution through FEM is matching the formulation and the computational organisation to the problem. This was realised forcefully in the present case of the solution of a class of moving contact boundary value problems of fastener joints. This paper deals with the problem of changing contact at the pin-hole interface of a fastener joint. Due to moving contact, the stresses and displacements are nonlinear with load. This would, in general, need an interactive-incremental approach for solution. However, by posing the problem in an inverse way, a solution is sought for obtaining loads to suit given contact configuration. Numerical results are given for typical isotropic and composite plates with rigid pins. Two cases of loading are considered: (i) load applied only at the edges of the plate and (ii) load applied at the pin and reacted at a part of the edge of the plate. Load-contact relationships, compliance and stress-patterns are investigated. This paper clearly demonstrates the simplification achieved by a suitable formulation of the problem. The results are of significance to the design and analysis of fastener joints.

Relaminarization in highly accelerated turbulent boundary layers

The mean flow development in an initially turbulent boundary layer subjected to a large favourable pressure gradient beginning at a point x0 is examined through analyses expected a priori to be valid on either side of relaminarization. The ‘quasi-laminar’ flow in the later stages of reversion, where the Reynolds stresses have by definition no significant effect on the mean flow, is described by an asymptotic theory constructed for large values of a pressure-gradient parameter Λ, scaled on a characteristic Reynolds stress gradient. The limiting flow consists of an inner laminar boundary layer and a matching inviscid (but rotational) outer layer. There is consequently no entrainment to lowest order in Λ−1, and the boundary layer thins down to conserve outer vorticity. In fact, the predictions of the theory for the common measures of boundary-layer thickness are in excellent agreement with experimental results, almost all the way from x0. On the other hand the development of wall parameters like the skin friction suggests the presence of a short bubble-shaped reverse-transitional region on the wall, where neither turbulent nor quasi-laminar calculations are valid. The random velocity fluctuations inherited from the original turbulence decay with distance, in the inner layer, according to inverse-power laws characteristic of quasi-steady perturbations on a laminar flow. In the outer layer, there is evidence that the dominant physical mechanism is a rapid distortion of the turbulence, with viscous and inertia forces playing a secondary role. All the observations available suggest that final retransition to turbulence quickly follows the onset of instability in the inner layer.It is concluded that reversion in highly accelerated flows is essentially due to domination of pressure forces over the slowly responding Reynolds stresses in an originally turbulent flow, accompanied by the generation of a new laminar boundary layer stabilized by the favourable pressure gradient.

Realization of first order optical systems using thin lensest

A first order optical system is investigated in full generality within the context of wave optics. The problem is reduced to a study of the ray transfer matrices. The simplest such systems correspond to axially symmetric propagation. Realization of such systems by centrally located lenses separated by finite distances is studied. It is shown that, contrary to the commonly held view, the set of first order systems that can be realized using axially symmetric thin lenses exhausts the entire SL(2, R) group; at most three lenses are needed to realize any element of this group. In particular, the inverse of free propagation can be so realized. Among anisotropic systems it is again shown that every element of the lens group Sp(4, R) can be realized using a finite number of thin lenses.

Forced Oscillations of a Contained Rotating Stratified Fluid

The motion generated by forced oscillations in an incompressible inviscid rotating and/or stratified fluid is examined under linear theory taking the density variation on the inertia terms. The solution consists of numerous internal modes in addition to the mode which oscillates with forcing frequency. Resonance occurs when the forcing frequency is equal to one of the frequencies of the internal modes. Some of these modes grow linearly or exponentially with time rendering the motion unstable and eventually may lead to turbulence. Most of the results discussed here will be missed under Boussinesq approximation.

Dendritic structures produced on solidification of multicomponent aqueous solutions

Dendrite structures of ice produced on undirectional solidification of ternary and quaternary aqueous solutions have been studied. Upon freezing, solutions containing more than one solute produce plate-shaped dendrites of ice. The spacing between dendrites increase linearly with the distance from the chill surface and the square root of local solidification time (or square root of inverse freezing rate) for any fixed composition. For fixed freezing conditions, the dendrite spacings from multicomponent aqueous solutions were a function of the concentrations and diffusion coefficients of the individual solutes. The dendrite spacing produced by freezing of a solution was changed by the addition of a solute different from those already present. If the main diffusion coefficient of the added solute is higher than that of solutes already present, the dendrite spacing is increased and vice versa. The dendrite spacing in multi-component systems increases with the total solute concentration if the constituent solutes are present in equal amounts. The dendrite spacing obtained on freezing of these dilute multicomponent solutions can be expressed by regression equations of the type Image Full-size image (2K) where L is the dendrite spacing in microns, C1, C2 and C3 are concentrations of individual solutes, Θf is the total freezing time and A1 −A8 are constants. A Yates analysis of the dendrite spacings in a factorial design of quaternary solutions indicates that there are strong interactions between individual solutes in regard to their effect on the dendrite spacings. A mass transport analysis has been used to calculate the interdendritic supersaturation ΔC of the individual solutes, the supercooling in the interdendritic liquid ΔT, and the transverse growth velocity of the dendrites, VT. In ternary solutions if two solutes are present in equal amount the supersaturation of the solute with higher main diffusion coefficient is lower, and vice versa. If a solute with higher main diffusion coefficient is added to a binary solution, the interface growth velocity, the interdendritic supersaturation of the base solute and the interdendritic supercooling increase with the quantity of solute added.

EPR studies on donor doped BaTiO3 grain boundary layer ceramic dielectrics

Donor doped BaTiO3 ceramics become insulating5 under controlled conditions with effective dielectric constants >10. The changes in EPR signals indicate that a certain fraction of the donor doped BaTiO3 is cubic even at room temperature and that the cubic fraction increases with the donor content. X-ray powder diffraction data support the EPR results. The coexistence of both the phases over a range of temperature is characteristic of diffused phase transition. The effect of grain size variation on EPR signal intensities indicate that the boundary layers surrounding the grains may constitute the cubic phase as a result of higher Ba-vacancies and donor contents at the grain boundary layer than in the bulk. Since the acceptor states arising from the Ba-vacancies and the impurities are activated in the cubic phase, they capture electrons from the conduction band, rendering the cubic phase electrically more insulating than the semiconductive tetragonal grain interiors. Thus, the cubic grain boundary layers act as effective dielectric media where the field tends to concentrate.

Leichhardt's vanishing thrips: A remarkable new genus and species of Insecta, Thysanoptera, Phlaeothripinae

Leichhardtithrips evanidus gen. et sp. nov., a wingless, fungus-feeding species of Phlaeothripidae, is described from a single female taken in D'Aguilar National Park, Brisbane. This species has exceptionally reduced sclerotisation rendering systematic comparisons of the new genus difficult. The cylindrical form of the eighth antennal segment is unique in the family, and the compound eyes are not developed on the ventral surface of the head. © The State of Queensland, Queensland Museum 2013.

Methods to determine resistance to anthelmintics when continuing larval development occurs

The in vivo faecal egg count reduction test (FECRT) is the most commonly used test to detect anthelmintic resistance (AR) in gastrointestinal nematodes (GIN) of ruminants in pasture based systems. However, there are several variations on the method, some more appropriate than others in specific circumstances. While in some cases labour and time can be saved by just collecting post-drench faecal worm egg counts (FEC) of treatment groups with controls, or pre- and post-drench FEC of a treatment group with no controls, there are circumstances when pre- and post-drench FEC of an untreated control group as well as from the treatment groups are necessary. Computer simulation techniques were used to determine the most appropriate of several methods for calculating AR when there is continuing larval development during the testing period, as often occurs when anthelmintic treatments against genera of GIN with high biotic potential or high re-infection rates, such as Haemonchus contortus of sheep and Cooperia punctata of cattle, are less than 100% efficacious. Three field FECRT experimental designs were investigated: (I) post-drench FEC of treatment and controls groups, (II) pre- and post-drench FEC of a treatment group only and (III) pre- and post-drench FEC of treatment and control groups. To investigate the performance of methods of indicating AR for each of these designs, simulated animal FEC were generated from negative binominal distributions with subsequent sampling from the binomial distributions to account for drench effect, with varying parameters for worm burden, larval development and drench resistance. Calculations of percent reductions and confidence limits were based on those of the Standing Committee for Agriculture (SCA) guidelines. For the two field methods with pre-drench FEC, confidence limits were also determined from cumulative inverse Beta distributions of FEC, for eggs per gram (epg) and the number of eggs counted at detection levels of 50 and 25. Two rules for determining AR: (1) %reduction (%R) < 95% and lower confidence limit <90%; and (2) upper confidence limit <95%, were also assessed. For each combination of worm burden, larval development and drench resistance parameters, 1000 simulations were run to determine the number of times the theoretical percent reduction fell within the estimated confidence limits and the number of times resistance would have been declared. When continuing larval development occurs during the testing period of the FECRT, the simulations showed AR should be calculated from pre- and post-drench worm egg counts of an untreated control group as well as from the treatment group. If the widely used resistance rule 1 is used to assess resistance, rule 2 should also be applied, especially when %R is in the range 90 to 95% and resistance is suspected.

Location and scale mixtures of Gaussians with flexible tail behaviour: Properties, inference and application to multivariate clustering

The family of location and scale mixtures of Gaussians has the ability to generate a number of flexible distributional forms. The family nests as particular cases several important asymmetric distributions like the Generalized Hyperbolic distribution. The Generalized Hyperbolic distribution in turn nests many other well known distributions such as the Normal Inverse Gaussian. In a multivariate setting, an extension of the standard location and scale mixture concept is proposed into a so called multiple scaled framework which has the advantage of allowing different tail and skewness behaviours in each dimension with arbitrary correlation between dimensions. Estimation of the parameters is provided via an EM algorithm and extended to cover the case of mixtures of such multiple scaled distributions for application to clustering. Assessments on simulated and real data confirm the gain in degrees of freedom and flexibility in modelling data of varying tail behaviour and directional shape.

Focal brachytherapy treatment planning using multi-parametric MRI and biological dose optimisation

We describe a novel approach to treatment planning for focal brachytherapy utilizing a biologically based inverse optimization algorithm and biological imaging to target an ablative dose at known regions of significant tumour burden and a lower, therapeutic dose to low risk regions.

Suomen maankuoren kivilajikoostumuksen tulkinta seismisistä aineistoista

Aim of this study is to investigate composition of the crust in Finland using seismic wide-angle velocity models and laboratory measurements on P- and S-wave velocities of different rock types. The velocities adopted from wide-angle velocity models were compared with laboratory velocities of different rock types corrected for the crustal PT conditions in the study area. The wide-angle velocity models indicate that the P-wave velocity does not only increase step-wise at boundaries of major crustal layers, but there is also gradual increase of velocity within the layers. On the other hand, the laboratory measurements of velocities indicate that no single rock type is able to provide the gradual downward increasing trends. Thus, there must be gradual vertical changes in rock composition. The downward increase of velocities indicates that the composition of the crust becomes gradually more mafic with increasing depth. Even though single rock types cannot simulate the wide-angle model velocities, it can be done with a mixture of rock types. There are a large number of rock type mixtures giving the correct P-wave velocities. Therefore, the inverse solution of rock types and their proportions from velocities is a non-unique problem if only P-wave velocities is available. Amount of the possible rock type mixtures can be limitted using S-wave velocities, reflection seismic results and other geological and geophysical results of the study area. Crustal model FINMIX-2 is presented in this study and it suggest that the crustal velocity profiles can be simulated with rock type mixtures, where the upper crust consists of felsic gneisses and granitic-granodioritic rocks with a minor contribution of quartzite, amphibolite and diabase. In the middle crust the amphibolite proportion increases. The lower crust consists of tonalitic gneiss, mafic garnet granulite, hornblendite, pyroxenite and minor mafic eclogite. This composition model is in agreement with deep crustal kimberlite-hosted xenolith data in eastern Finland and reflectivity of the FIRE (Finnish Reflection Experiment). According to FINMIX-2 model the Moho is deeper and the crustal composition is a more mafic than an average global continental model would suggest. Composition models of southern Finland are quite similar than FINMIX-2 model. However, there are minor differencies between the models, which indicates areal differences of composition. Models of northern Finland shows that the crustal thickness is smaller than southern Finland and composition of the upper crust is different. Density profiles calculated from the lithological models suggest that there is practically no density contrast at Moho in areas of the high-velocity lower crust. This implies that crustal thickness in the central Fennoscandian Shield may have been controlled by the densities of the lower crustal and upper mantle rocks.

Investigation of the effects of extracellular osmotic pressure on morphology and mechanical 1 properties of individual chondrocyte

It has been demonstrated that most cells of the body respond to osmotic pressure in a systematic manner. The disruption of the collagen network in the early stages of osteoarthritis causes an increase in water content of cartilage which leads to a reduction of pericellular osmolality in chondrocytes distributed within the extracellular environment. It is therefore arguable that an insight into the mechanical properties of chondrocytes under varying osmotic pressure would provide a better understanding of chondrocyte mechanotransduction and potentially contribute to knowledge on cartilage degeneration. In this present study, the chondrocyte cells were exposed to solutions with different osmolality. Changes in their dimensions and mechanical properties were measured over time. Atomic Force Microscopy (AFM) was used to apply load at various strain-rates and the force-time curves were logged. The thin-layer elastic model was used to extract the elastic stiffness of chondrocytes at different strain-rates and at different solution osmolality. In addition, the porohyperelastic (PHE) model was used to investigate the strain-rate dependent responses under the loading and osmotic pressure conditions. The results revealed that the hypo-osmotic external environment increased chondrocyte dimensions and reduced Young’s modulus of the cells at all strain-rates tested. In contrast, the hyper-osmotic external environment reduced dimensions and increased Young’s modulus. Moreover, by using the PHE model coupled with inverse FEA simulation, we established that the hydraulic permeability of chondrocytes increased with decreasing extracellular osmolality which is consistent with previous work in the literature. This could be due to a higher intracellular fluid volume fraction with lower osmolality.