Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Optimal patch-leaving behaviour: a case study using the parasitoid Cotesia rubecular

1. Parasitoids are predicted to spend longer in patches with more hosts, but previous work on Cotesia rubecula (Marshall) has not upheld this prediction, Tests of theoretical predictions may be affected by the definition of patch leaving behaviour, which is often ambiguous. 2. In this study whole plants were considered as patches and assumed that wasps move within patches by means of walking or flying. Within-patch and between-patch flights were distinguished based on flight distance. The quality of this classification was tested statistically by examination of log-survivor curves of flight times. 3. Wasps remained longer in patches with higher host densities, which is consistent with predictions of the marginal value theorem (Charnov 1976). tinder the assumption that each flight indicates a patch departure, there is no relationship between host density and leaving tendency. 4. Oviposition influences the patch leaving behaviour of wasps in a count down fashion (Driessen et al. 1995), as predicted by an optimal foraging model (Tenhumberg, Keller & Possingham 2001). 5. Wasps spend significantly longer in the first patch encountered following release, resulting in an increased rate of superparasitism.

Development of a short-term model using petrophysical logging at Century Mine, North Queensland

Pasminco Century Mine has developed a geophysical logging system to provide new data for ore mining/grade control and the generation of Short Term Models for mine planning. Previous work indicated the applicability of petrophysical logging for lithology prediction, however, the automation of the method was not considered reliable enough for the development of a mining model. A test survey was undertaken using two diamond drilled control holes and eight percussion holes. All holes were logged with natural gamma, magnetic susceptibility and density. Calibration of the LogTrans auto-interpretation software using only natural gamma and magnetic susceptibility indicated that both lithology and stratigraphy could be predicted. Development of a capability to enforce stratigraphic order within LogTrans increased the reliability and accuracy of interpretations. After the completion of a feasibility program, Century Mine has invested in a dedicated logging vehicle to log blast holes as well as for use in in-fill drilling programs. Future refinement of the system may lead to the development of GPS controlled excavators for mining ore.

The investigation of factors underlying deficits in self-awareness and self-regulation

Primary objective: To examine a theoretical model which suggests that a contribution of both psychological and neuropsychological factors underlie deficits in self-awareness and self-regulation. Research design: Multivariate design including correlations and analysis of variance (ANOVA). Methods: Sixty-one subjects with acquired brain injury (ABI) were administered standardized measures of self-awareness and self-regulation. Psychological factors included measures of coping-related denial, personality-related denial and personality change. Neuropsychological factors included an estimate of IQ and two measures of executive functioning that assess capacity for volition and purposive behaviour. Main outcomes and results: The findings indicated that the relative contribution of neuropsychological factors to an outcome of deficits in self-awareness and self-regulation had a more direct effect than psychological factors. In general, measures of executive functioning had a direct relationship, while measures of coping-related and personality-related denial had an indirect relationship with measures of self-awareness and self-regulation. Conclusion: The findings highlighted the importance of measuring both neuropsychological and psychological factors and demonstrated that the relative contribution of these variables varies according to different levels of self-awareness and self-regulation.

Genetic sources of covariation among P3(00) and online performance variables in a delayed-response working memory task

Genetic and environmental sources of covariation among the P3(00) and online performance elicited in a delayed-response working memory task, and psychometric IQ assessed by the multidimensional aptitude battery, were examined in an adolescent twin sample. An association between frontal P3 latency and task performance (phenotypic r = -0.33; genotypic r = -0.49) was indicated, with genes (i.e. twin status) accounting for a large part of the covariation ( > 70%). In contrast, genes influencing P3 amplitude mediated only a small part (2%) of the total genetic variation in task performance. While task performance mediated 15% of the total genetic variation in IQ (phenotypic r = 0.22; genotypic r = 0.39) there was no association between P3 latency and IQ or P3 amplitude with IQ. The findings provide some insight into the inter-relationships among psychophysiological, performance and psychometric measures of cognitive ability, and provide support for a levels-of-processing genetic model of cognition where genes act on specific sub-components of cognitive processes.

Mapping Gray's model of personality onto the Eysenck Personality Profiler (EPP)

The aim of this study was to determine how well Gray's model of personality [Gray, J.A. (1982). The neuropsychology of anxiety: an enquiry into the functions of the septo-hippocampal system. Oxfords: Oxford University Press, Gray, J.A. (1987). The psychology of fear and stress. Cambridge: Cambridge University Press], as measured by the Gray Wilson Personality Questionnaire (GWPQ), can provide a full description of personality as measured by the primary scales of the Eysenck Personality Profiler (EPP) and the type scales of the short version or the EPQ-R. Factor analysis of the GWPQ the Anxiety and linpulsivity scales of the EPP and the Learning Styles Questionnaire (LSQ) showed that the GWPQ seemed to measure general activation and inhibition factors, but not the finer features of Gray's theory. When the GWPQ scales were regressed against each scale of the EPP., it was round that they generally provide only a reasonable explanation of the EPP primary scales. It is concluded that the GWPQ measures general propel-ties of Gray's model, that the linpulsivity and Anxiety scales of the EPP also scent related to the GWPQ scales, and that Gray's model of personality provides only a partial explanation of personality in general. (C) 2002 Published by Elsevier Science Ltd. All rights reserved.

Review and future directions in the modelling and control of continuous drum granulation

Many granulation plants operate well below design capacity, suffering from high recycle rates and even periodic instabilities. This behaviour cannot be fully predicted using the present models. The main objective of the paper is to provide an overview of the current status of model development for granulation processes and suggest future directions for research and development. The end-use of the models is focused on the optimal design and control of granulation plants using the improved predictions of process dynamics. The development of novel models involving mechanistically based structural switching methods is proposed in the paper. A number of guidelines are proposed for the selection of control relevant model structures. (C) 2002 Published by Elsevier Science B.V.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.

The role of melanocortin-1 receptor polymorphism in skin cancer risk phenotypes

We have examined melanocortin-1 receptor (MC1R) variant allele frequencies in the general population and in a collection of adolescent dizygotic and monozygotic twins to determine statistical associations of pigmentation phenotypes with increased skin cancer risk. This included hair and skin color, freckling, mole count and sun exposed skin reflectance. Nine variants were studied and designated as either strong R (OR = 63; 95% CI 32-140) or weak r (OR = 5; 95% CI 3-11) red hair alleles. Penetrance of each MC1R variant allele was consistent with an allelic model where effects were multiplicative for red hair but additive for skin reflectance. To assess the interaction of the brown eye color gene BEY2/OCA2 on the phenotypic effects of variant MC1R alleles we imputed OCA2 genotype in the twin collection. A modifying effect of OCA2 on MC1R variant alleles was seen on constitutive skin color, freckling and mole count. In order to study the individual effects of these variants on pigmentation phenotype we have established a series of human primary melanocyte strains genotyped for the MC1R receptor. These include strains which are MC1R wild-type consensus, variant heterozygotes, and homozygotes for strong R alleles Arg151Cys and Arg160Trp. Ultrastructural analysis demonstrated that only consensus strains contained stage III and IV melanosomes in their terminal dendrites whereas Arg151Cys and Arg160Trp homozygous strains contained only immature stage I and II melanosomes. Such genetic association studies combined with the functional analysis of MC1R variant alleles in melanocytic cells should provide a link in understanding the association between pigmentary phototypes and skin cancer risk.

Linking physiological and architectural models of cotton

Despite the strong influence of plant architecture on crop yield, most crop models either ignore it or deal with it in a very rudimentary way. This paper demonstrates the feasibility of linking a model that simulates the morphogenesis and resultant architecture of individual cotton plants with a crop model that simulates the effects of environmental factors on critical physiological processes and resulting yield in cotton. First the varietal parameters of the models were made concordant. Then routines were developed to allocate the flower buds produced each day by the crop model amongst the potential positions generated by the architectural model. This allocation is done according to a set of heuristic rules. The final weight of individual bolls and the shedding of buds and fruit caused by water, N, and C stresses are processed in a similar manner. Observations of the positions of harvestable fruits, both within and between plants, made under a variety of agronomic conditions that had resulted in a broad range of plant architectures were compared to those predicted by the model with the same environmental inputs. As illustrated by comparisons of plant maps, the linked models performed reasonably well, though performance of the fruiting point allocation and shedding algorithms could probably be improved by further analysis of the spatial relationships of retained fruit. (C) 2002 Elsevier Science Ltd. All rights reserved.

On the economic analysis of response to preventive measures

There are many circumstances in which the effectiveness of preventive measures depends to a large extent on the compliance of the patient in changing his or her behavior or lifestyle. It is shown how economic techniques can be used (i) to describe the rationale of individuals and predict their behavior (Section 2); and (ii) to assess preventive measures that, by requiring a change of conduct, imply "costs" to the individual due to a decline in the quality of life (Appendix). Cigarette smoking and coronary heart disease are used as an illustration. While the analysis of Section 2 uses graphical techniques, a simple textbook-type of lifetime utility model with a mathematical emphasis is used in the Appendix. It is also shown that techniques often used to assess health care programs such as the QALYs (Quality-Adjusted Life Years) are inappropriate to the evaluation of preventive programs aiming at behavioral changes. Finally, topics that call for further research are indicated.

Surgical correction of scoliosis: Numerical analysis and optimization of the procedure

A previously developed model is used to numerically simulate real clinical cases of the surgical correction of scoliosis. This model consists of one-dimensional finite elements with spatial deformation in which (i) the column is represented by its axis; (ii) the vertebrae are assumed to be rigid; and (iii) the deformability of the column is concentrated in springs that connect the successive rigid elements. The metallic rods used for the surgical correction are modeled by beam elements with linear elastic behavior. To obtain the forces at the connections between the metallic rods and the vertebrae geometrically, non-linear finite element analyses are performed. The tightening sequence determines the magnitude of the forces applied to the patient column, and it is desirable to keep those forces as small as possible. In this study, a Genetic Algorithm optimization is applied to this model in order to determine the sequence that minimizes the corrective forces applied during the surgery. This amounts to find the optimal permutation of integers 1, ... , n, n being the number of vertebrae involved. As such, we are faced with a combinatorial optimization problem isomorph to the Traveling Salesman Problem. The fitness evaluation requires one computing intensive Finite Element Analysis per candidate solution and, thus, a parallel implementation of the Genetic Algorithm is developed.

The social adjustment of patients with schizophrenia: implications to the mental health policy in Brazil

A sample (n=124) of schizophrenic patients from a defined catchment area of the city os S.Paulo, Brazil, who had been consecutively admitted to hospital, was assessed for psychopathological status and social adjustment levels. Sociodemographic, socio-economic and occupational characteristics were recorded: almost 30% of the subjects had no occupation and received no social benefit, more than two-thirds had a monthly per capita income of US$ 100.00 or less. Sixty-five percent presented with Schneiderian firstrank symptoms. Nearly half the sample showed poor or very poor social adjustment in the month prior to admission. The most affected areas of social functioning were participation in the household activities, work and social withdrawal. The current mental health policy of promoting extra-mural care as an alternative to the previous hospital-based model will then mean the investment in a network of new community-based services, that give effective treatment and support to patients and their families. The need of further research into the current picture of mental disorders in the country is stressed.

Optimal insurance with counterparty default risk

We study the design of optimal insurance contracts when the insurer can default on its obligations. In our model default arises endogenously from the interaction of the insurance premium, the indemnity schedule and the insurer’s assets. This allows us to understand the joint effect of insolvency risk and background risk on efficient contracts. The results may shed light on the aggregate risk retention sched- ules observed in catastrophe reinsurance markets, and can assist in the design of (re)insurance programs and guarantee funds.