Nonlinear fatigue damage model based on the residual strength degradation law

In this paper, a logarithmic expression to describe the residual strength degradation process is developed in order to fatigue test results for normalized carbon steel. The definition and expression of fatigue damage due to symmetrical stress with a constant amplitude are also given. The expression of fatigue damage can also explain the nonlinear properties of fatigue damage. Furthermore, the fatigue damage of structures under random stress is analyzed, and an iterative formula to describe the fatigue damage process is deduced. Finally, an approximate method for evaluating the fatigue life of structures under repeated random stress blocking is presented through various calculation examples.

Failure process analysis of cement under explosive loading with a generalized driven nonlinear threshold model

The microstructural heterogeneity and stress fluctuation play important roles in the failure process of brittle materials. In this paper, a generalized driven nonlinear threshold model with stress fluctuation is presented to study the effects of microstructural heterogeneity on continuum damage evolution. As an illustration, the failure process of cement material under explosive loading is analyzed using the model. The result agrees well with the experimental one, which proves the efficiency of the model.

An observer-based vaccination control law for an SEIR epidemic model based on feedback linearization techniques for nonlinear systems

This paper presents a vaccination strategy for fighting against the propagation of epidemic diseases. The disease propagation is described by an SEIR (susceptible plus infected plus infectious plus removed populations) epidemic model. The model takes into account the total population amounts as a refrain for the illness transmission since its increase makes the contacts among susceptible and infected more difficult. The vaccination strategy is based on a continuous-time nonlinear control law synthesised via an exact feedback input-output linearization approach. An observer is incorporated into the control scheme to provide online estimates for the susceptible and infected populations in the case when their values are not available from online measurement but they are necessary to implement the control law. The vaccination control is generated based on the information provided by the observer. The control objective is to asymptotically eradicate the infection from the population so that the removed-by-immunity population asymptotically tracks the whole one without precise knowledge of the partial populations. The model positivity, the eradication of the infection under feedback vaccination laws and the stability properties as well as the asymptotic convergence of the estimation errors to zero as time tends to infinity are investigated.

Asymptotic scaling in the two-dimensional O(3) Nonlinear sigma model: a Monte Carlo study on parallel computers

We investigate the 2d O(3) model with the standard action by Monte Carlo simulation at couplings β up to 2.05. We measure the energy density, mass gap and susceptibility of the model, and gather high statistics on lattices of size L ≤ 1024 using the Floating Point Systems T-series vector hypercube and the Thinking Machines Corp.'s Connection Machine 2. Asymptotic scaling does not appear to set in for this action, even at β = 2.10, where the correlation length is 420. We observe a 20% difference between our estimate m/Λ^─_(Ms) = 3.52(6) at this β and the recent exact analytical result . We use the overrelaxation algorithm interleaved with Metropolis updates and show that decorrelation time scales with the correlation length and the number of overrelaxation steps per sweep. We determine its effective dynamical critical exponent to be z' = 1.079(10); thus critical slowing down is reduced significantly for this local algorithm that is vectorizable and parallelizable.

We also use the cluster Monte Carlo algorithms, which are non-local Monte Carlo update schemes which can greatly increase the efficiency of computer simulations of spin models. The major computational task in these algorithms is connected component labeling, to identify clusters of connected sites on a lattice. We have devised some new SIMD component labeling algorithms, and implemented them on the Connection Machine. We investigate their performance when applied to the cluster update of the two dimensional Ising spin model.

Finally we use a Monte Carlo Renormalization Group method to directly measure the couplings of block Hamiltonians at different blocking levels. For the usual averaging block transformation we confirm the renormalized trajectory (RT) observed by Okawa. For another improved probabilistic block transformation we find the RT, showing that it is much closer to the Standard Action. We then use this block transformation to obtain the discrete β-function of the model which we compare to the perturbative result. We do not see convergence, except when using a rescaled coupling β_E to effectively resum the series. For the latter case we see agreement for m/ Λ^─_(Ms) at , β = 2.14, 2.26, 2.38 and 2.50. To three loops m/Λ^─_(Ms) = 3.047(35) at β = 2.50, which is very close to the exact value m/ Λ^─_(Ms) = 2.943. Our last point at β = 2.62 disagrees with this estimate however.

SteelConverter and Caltech VirtualShaker: Rapid Nonlinear Cloud-Based Structural Model Conversion and Analysis

STEEL, the Caltech created nonlinear large displacement analysis software, is currently used by a large number of researchers at Caltech. However, due to its complexity, lack of visualization tools (such as pre- and post-processing capabilities) rapid creation and analysis of models using this software was difficult. SteelConverter was created as a means to facilitate model creation through the use of the industry standard finite element solver ETABS. This software allows users to create models in ETABS and intelligently convert model information such as geometry, loading, releases, fixity, etc., into a format that STEEL understands. Models that would take several days to create and verify now take several hours or less. The productivity of the researcher as well as the level of confidence in the model being analyzed is greatly increased.

It has always been a major goal of Caltech to spread the knowledge created here to other universities. However, due to the complexity of STEEL it was difficult for researchers or engineers from other universities to conduct analyses. While SteelConverter did help researchers at Caltech improve their research, sending SteelConverter and its documentation to other universities was less than ideal. Issues of version control, individual computer requirements, and the difficulty of releasing updates made a more centralized solution preferred. This is where the idea for Caltech VirtualShaker was born. Through the creation of a centralized website where users could log in, submit, analyze, and process models in the cloud, all of the major concerns associated with the utilization of SteelConverter were eliminated. Caltech VirtualShaker allows users to create profiles where defaults associated with their most commonly run models are saved, and allows them to submit multiple jobs to an online virtual server to be analyzed and post-processed. The creation of this website not only allowed for more rapid distribution of this tool, but also created a means for engineers and researchers with no access to powerful computer clusters to run computationally intensive analyses without the excessive cost of building and maintaining a computer cluster.

In order to increase confidence in the use of STEEL as an analysis system, as well as verify the conversion tools, a series of comparisons were done between STEEL and ETABS. Six models of increasing complexity, ranging from a cantilever column to a twenty-story moment frame, were analyzed to determine the ability of STEEL to accurately calculate basic model properties such as elastic stiffness and damping through a free vibration analysis as well as more complex structural properties such as overall structural capacity through a pushover analysis. These analyses showed a very strong agreement between the two softwares on every aspect of each analysis. However, these analyses also showed the ability of the STEEL analysis algorithm to converge at significantly larger drifts than ETABS when using the more computationally expensive and structurally realistic fiber hinges. Following the ETABS analysis, it was decided to repeat the comparisons in a software more capable of conducting highly nonlinear analysis, called Perform. These analyses again showed a very strong agreement between the two softwares in every aspect of each analysis through instability. However, due to some limitations in Perform, free vibration analyses for the three story one bay chevron brace frame, two bay chevron brace frame, and twenty story moment frame could not be conducted. With the current trend towards ultimate capacity analysis, the ability to use fiber based models allows engineers to gain a better understanding of a building’s behavior under these extreme load scenarios.

Following this, a final study was done on Hall’s U20 structure [1] where the structure was analyzed in all three softwares and their results compared. The pushover curves from each software were compared and the differences caused by variations in software implementation explained. From this, conclusions can be drawn on the effectiveness of each analysis tool when attempting to analyze structures through the point of geometric instability. The analyses show that while ETABS was capable of accurately determining the elastic stiffness of the model, following the onset of inelastic behavior the analysis tool failed to converge. However, for the small number of time steps the ETABS analysis was converging, its results exactly matched those of STEEL, leading to the conclusion that ETABS is not an appropriate analysis package for analyzing a structure through the point of collapse when using fiber elements throughout the model. The analyses also showed that while Perform was capable of calculating the response of the structure accurately, restrictions in the material model resulted in a pushover curve that did not match that of STEEL exactly, particularly post collapse. However, such problems could be alleviated by choosing a more simplistic material model.

A nonlinear systems model for the control mechanism of free fatty acid-glucose metabolosm in normal humans

A mathematical model is proposed in this thesis for the control mechanism of free fatty acid-glucose metabolism in healthy individuals under resting conditions. The objective is to explain in a consistent manner some clinical laboratory observations such as glucose, insulin and free fatty acid responses to intravenous injection of glucose, insulin, etc. Responses up to only about two hours from the beginning of infusion are considered. The model is an extension of the one for glucose homeostasis proposed by Charette, Kadish and Sridhar (Modeling and Control Aspects of Glucose Homeostasis. Mathematical Biosciences, 1969). It is based upon a systems approach and agrees with the current theories of glucose and free fatty acid metabolism. The description is in terms of ordinary differential equations. Validation of the model is based on clinical laboratory data available at the present time. Finally procedures are suggested for systematically identifying the parameters associated with the free fatty acid portion of the model.

A Z-scan model for optical nonlinear nanometric films

The Z-scan technique is useful for measuring the nonlinear refractive index of thin films. In conventional Z-scan theories, two effects are often ignored, namely the losses due to the internal multi-interference and the nonlinear absorption inside the sample. Therefore, the theories are restricted to relatively thick films. For films thinner than about 100 nm, the two effects become significant, and thus cannot be ignored. In the present work, we present a Z-scan theory that takes both effects into account. The proposed model calculation is suitable for optical nonlinear films of nanometric thickness. With numerical simulations, we demonstrate dramatic deviations from the conventional Z-scan calculations.

A model for nonlinear viscoelastic mechanical responses of collagenous soft tissues

Experimental observations of the time-dependent mechanical responses of collagenous tissues have demonstrated behavior that deviates from standard treatments of linear or quasi-linear viscoelasticity. In particular, time-dependent deformation can be strongly coupled to strain level, and strain-rate independence can be observed under monotonic loading, even for a tissue with dramatic stress relaxation. It was postulated that this nonlinearity is fundamentally associated with gradual recruitment of individual collagen fibrils during applied mechanical loading. Based on previously observed experimental results for the time-dependent response of collagenous soft tissues, a model is developed to describe the mechanical behavior of these tissues under uniaxial loading. Tissue stresses, under applied strain-controlled loading, are assumed to be a sum of elastic and viscoelastic stress contributions. The relative contributions of elastic and viscoelastic stresses is assumed to vary with strain level, leading to strain- and time-dependent mechanical behavior. The model formulation is examined under conditions of monotonic loading at varying constant strain rates and stress-relaxation at different applied strain levels. The model is compared with experimental data for a membranous biological soft tissue, the amniotic sac, and is found to agree well with experimental results. The limiting behavior of the novel model, at large strains relative to the collagen recruitment, is consistent with the quasi-linear viscoelastic approach. © 2006 Materials Research Society.