Experimental validation of a nonlinear μ FE model based on cohesive-frictional plasticity for trabecular bone

Trabecular bone is a porous mineralized tissue playing a major load bearing role in the human body. Prediction of age-related and disease-related fractures and the behavior of bone implant systems needs a thorough understanding of its structure-mechanical property relationships, which can be obtained using microcomputed tomography-based finite element modeling. In this study, a nonlinear model for trabecular bone as a cohesive-frictional material was implemented in a large-scale computational framework and validated by comparison of μFE simulations with experimental tests in uniaxial tension and compression. A good correspondence of stiffness and yield points between simulations and experiments was found for a wide range of bone volume fraction and degree of anisotropy in both tension and compression using a non-calibrated, average set of material parameters. These results demonstrate the ability of the model to capture the effects leading to failure of bone for three anatomical sites and several donors, which may be used to determine the apparent behavior of trabecular bone and its evolution with age, disease, and treatment in the future.

Nonlinear quasi-static finite element simulations predict in vitro strength of human proximal femora assessed in a dynamic sideways fall setup

Osteoporotic proximal femur fractures are caused by low energy trauma, typically when falling on the hip from standing height. Finite element simulations, widely used to predict the fracture load of femora in fall, usually include neither mass-related inertial effects, nor the viscous part of bone's material behavior. The aim of this study was to elucidate if quasi-static non-linear homogenized finite element analyses can predict in vitro mechanical properties of proximal femora assessed in dynamic drop tower experiments. The case-specific numerical models of thirteen femora predicted the strength (R2=0.84, SEE=540 N, 16.2%), stiffness (R2=0.82, SEE=233 N/mm, 18.0%) and fracture energy (R2=0.72, SEE=3.85 J, 39.6%); and provided fair qualitative matches with the fracture patterns. The influence of material anisotropy was negligible for all predictions. These results suggest that quasi-static homogenized finite element analysis may be used to predict mechanical properties of proximal femora in the dynamic sideways fall situation.

Computational Comparison of Continuous and Discontinuous Galerkin Time-Stepping Methods for Nonlinear Initial Value Problems

This article centers on the computational performance of the continuous and discontinuous Galerkin time stepping schemes for general first-order initial value problems in R n , with continuous nonlinearities. We briefly review a recent existence result for discrete solutions from [6], and provide a numerical comparison of the two time discretization methods.

SOMS: SurrOgate MultiStart algorithm for use with nonlinear programming for global optimization

SOMS is a general surrogate-based multistart algorithm, which is used in combination with any local optimizer to find global optima for computationally expensive functions with multiple local minima. SOMS differs from previous multistart methods in that a surrogate approximation is used by the multistart algorithm to help reduce the number of function evaluations necessary to identify the most promising points from which to start each nonlinear programming local search. SOMS’s numerical results are compared with four well-known methods, namely, Multi-Level Single Linkage (MLSL), MATLAB’s MultiStart, MATLAB’s GlobalSearch, and GLOBAL. In addition, we propose a class of wavy test functions that mimic the wavy nature of objective functions arising in many black-box simulations. Extensive comparisons of algorithms on the wavy testfunctions and on earlier standard global-optimization test functions are done for a total of 19 different test problems. The numerical results indicate that SOMS performs favorably in comparison to alternative methods and does especially well on wavy functions when the number of function evaluations allowed is limited.

Nonlinear tests for genetic studies of complex disease

Linkage and association studies are major analytical tools to search for susceptibility genes for complex diseases. With the availability of large collection of single nucleotide polymorphisms (SNPs) and the rapid progresses for high throughput genotyping technologies, together with the ambitious goals of the International HapMap Project, genetic markers covering the whole genome will be available for genome-wide linkage and association studies. In order not to inflate the type I error rate in performing genome-wide linkage and association studies, multiple adjustment for the significant level for each independent linkage and/or association test is required, and this has led to the suggestion of genome-wide significant cut-off as low as 5 × 10 −7. Almost no linkage and/or association study can meet such a stringent threshold by the standard statistical methods. Developing new statistics with high power is urgently needed to tackle this problem. This dissertation proposes and explores a class of novel test statistics that can be used in both population-based and family-based genetic data by employing a completely new strategy, which uses nonlinear transformation of the sample means to construct test statistics for linkage and association studies. Extensive simulation studies are used to illustrate the properties of the nonlinear test statistics. Power calculations are performed using both analytical and empirical methods. Finally, real data sets are analyzed with the nonlinear test statistics. Results show that the nonlinear test statistics have correct type I error rates, and most of the studied nonlinear test statistics have higher power than the standard chi-square test. This dissertation introduces a new idea to design novel test statistics with high power and might open new ways to mapping susceptibility genes for complex diseases. ^

NONLINEAR TIME SERIES/SYSTEM ANALYSIS (STATE-SPACE, DYNAMICS, INTERVENTION, RESILIENCE, KALMAN)

A discussion of nonlinear dynamics, demonstrated by the familiar automobile, is followed by the development of a systematic method of analysis of a possibly nonlinear time series using difference equations in the general state-space format. This format allows recursive state-dependent parameter estimation after each observation thereby revealing the dynamics inherent in the system in combination with random external perturbations.^ The one-step ahead prediction errors at each time period, transformed to have constant variance, and the estimated parametric sequences provide the information to (1) formally test whether time series observations y(,t) are some linear function of random errors (ELEM)(,s), for some t and s, or whether the series would more appropriately be described by a nonlinear model such as bilinear, exponential, threshold, etc., (2) formally test whether a statistically significant change has occurred in structure/level either historically or as it occurs, (3) forecast nonlinear system with a new and innovative (but very old numerical) technique utilizing rational functions to extrapolate individual parameters as smooth functions of time which are then combined to obtain the forecast of y and (4) suggest a measure of resilience, i.e. how much perturbation a structure/level can tolerate, whether internal or external to the system, and remain statistically unchanged. Although similar to one-step control, this provides a less rigid way to think about changes affecting social systems.^ Applications consisting of the analysis of some familiar and some simulated series demonstrate the procedure. Empirical results suggest that this state-space or modified augmented Kalman filter may provide interesting ways to identify particular kinds of nonlinearities as they occur in structural change via the state trajectory.^ A computational flow-chart detailing computations and software input and output is provided in the body of the text. IBM Advanced BASIC program listings to accomplish most of the analysis are provided in the appendix. ^

Analysis of life expectancy using linear and nonlinear models

Life expectancy has consistently increased over the last 150 years due to improvements in nutrition, medicine, and public health. Several studies found that in many developed countries, life expectancy continued to rise following a nearly linear trend, which was contrary to a common belief that the rate of improvement in life expectancy would decelerate and was fit with an S-shaped curve. Using samples of countries that exhibited a wide range of economic development levels, we explored the change in life expectancy over time by employing both nonlinear and linear models. We then observed if there were any significant differences in estimates between linear models, assuming an auto-correlated error structure. When data did not have a sigmoidal shape, nonlinear growth models sometimes failed to provide meaningful parameter estimates. The existence of an inflection point and asymptotes in the growth models made them inflexible with life expectancy data. In linear models, there was no significant difference in the life expectancy growth rate and future estimates between ordinary least squares (OLS) and generalized least squares (GLS). However, the generalized least squares model was more robust because the data involved time-series variables and residuals were positively correlated. ^

La militancia del PCA durante la última dictadura en Argentina : Un análisis sobre la producción, circulación y recepción de la línea partidaria

A partir del golpe de Estado de marzo de 1976 el Partido Comunista argentino (PCA) impulsó la política de "convergencia cívico-militar", una coalición amplia para frenar el avance del "fascismo" o "pinochetismo" dentro de la Junta Militar. Esta posición no era nueva sino que reactualizaba las caracterizaciones elaboradas en los meses previos porque, para el PCA, el golpe no inauguraba una nueva dictadura. En el presente artículo analizamos cómo esta línea, elaborada en el seno del Comité Central, fue recibida, interpretada, elaborada y materializada en prácticas por parte de la militancia comunista. El propósito de esta indagación es contestar por qué los militantes reprodujeron las posturas partidarias, a pesar de que para muchos, más temprano o más tarde, comenzó a resultar evidente que el plan represivo era sostenido por el conjunto de las Fuerzas Armadas y que las caracterizaciones del partido no se ajustaban a la realidad

La militancia del PCA durante la última dictadura en Argentina : Un análisis sobre la producción, circulación y recepción de la línea partidaria

A partir del golpe de Estado de marzo de 1976 el Partido Comunista argentino (PCA) impulsó la política de "convergencia cívico-militar", una coalición amplia para frenar el avance del "fascismo" o "pinochetismo" dentro de la Junta Militar. Esta posición no era nueva sino que reactualizaba las caracterizaciones elaboradas en los meses previos porque, para el PCA, el golpe no inauguraba una nueva dictadura. En el presente artículo analizamos cómo esta línea, elaborada en el seno del Comité Central, fue recibida, interpretada, elaborada y materializada en prácticas por parte de la militancia comunista. El propósito de esta indagación es contestar por qué los militantes reprodujeron las posturas partidarias, a pesar de que para muchos, más temprano o más tarde, comenzó a resultar evidente que el plan represivo era sostenido por el conjunto de las Fuerzas Armadas y que las caracterizaciones del partido no se ajustaban a la realidad

La militancia del PCA durante la última dictadura en Argentina : Un análisis sobre la producción, circulación y recepción de la línea partidaria

A partir del golpe de Estado de marzo de 1976 el Partido Comunista argentino (PCA) impulsó la política de "convergencia cívico-militar", una coalición amplia para frenar el avance del "fascismo" o "pinochetismo" dentro de la Junta Militar. Esta posición no era nueva sino que reactualizaba las caracterizaciones elaboradas en los meses previos porque, para el PCA, el golpe no inauguraba una nueva dictadura. En el presente artículo analizamos cómo esta línea, elaborada en el seno del Comité Central, fue recibida, interpretada, elaborada y materializada en prácticas por parte de la militancia comunista. El propósito de esta indagación es contestar por qué los militantes reprodujeron las posturas partidarias, a pesar de que para muchos, más temprano o más tarde, comenzó a resultar evidente que el plan represivo era sostenido por el conjunto de las Fuerzas Armadas y que las caracterizaciones del partido no se ajustaban a la realidad

Coccolithophore taxa and PCA factor analysis of ODP Site 1233

We provide high-resolution sea surface temperature (SST) and paleoproductivity data focusing on Termination 1. We describe a new method for estimating SSTs based on multivariate statistical analyses performed on modern coccolithophore census data, and we present the first downcore reconstructions derived from coccolithophore assemblages at Ocean Drilling Project (ODP) Site 1233 located offshore Chile. We compare our coccolithophore SST record to alkenone-based SSTs as well as SST reconstructions based on dinoflagellates and radiolaria. All reconstructions generally show a remarkable concordance. As in the alkenone SST record, the Last Glacial Maximum (LGM, 19-23 kyr B.P.) is not clearly defined in our SST reconstruction. After the onset of deglaciation, three major warming steps are recorded: from 18.6 to 18 kyr B.P. (~2.6°C), from 15.7 to 15.3 kyr B.P. (~2.5°C), and from 13 to 11.4 kyr B.P. (~3.4°C). Consistent with the other records from Site 1233 and Antarctic ice core records, we observed a clear Holocene Climatic Optimum (HCO) from ~8-12 kyr B.P. Combining the SST reconstruction with coccolith absolute abundances and accumulation rates, we show that colder temperatures during the LGM are linked to higher coccolithophore productivity offshore Chile and warmer SSTs during the HCO to lower coccolithophore productivity, with indications of weak coastal upwelling. We interpret our data in terms of latitudinal displacements of the Southern Westerlies and the northern margin of the Antarctic Circumpolar Current system over the deglaciation and the Holocene.