Predicting breeding values for milk yield of Guzera (Bos indicus) cows using random regression models

Given the importance of Guzera breeding programs for milk production in the tropics, the objective of this study was to compare alternative random regression models for estimation of genetic parameters and prediction of breeding values. Test-day milk yields records (TDR) were collected monthly, in a maximum of 10 measurements. The database included 20,524 records of first lactation from 2816 Guzera cows. TDR data were analyzed by random regression models (RRM) considering additive genetic, permanent environmental and residual effects as random and the effects of contemporary group (CG), calving age as a covariate (linear and quadratic effects) and mean lactation curve as fixed. The genetic additive and permanent environmental effects were modeled by RRM using Wilmink, All and Schaeffer and cubic B-spline functions as well as Legendre polynomials. Residual variances were considered as heterogeneous classes, grouped differently according to the model used. Multi-trait analysis using finite-dimensional models (FDM) for testday milk records (TDR) and a single-trait model for 305-days milk yields (default) using the restricted maximum likelihood method were also carried out as further comparisons. Through the statistical criteria adopted, the best RRM was the one that used the cubic B-spline function with five random regression coefficients for the genetic additive and permanent environmental effects. However, the models using the Ali and Schaeffer function or Legendre polynomials with second and fifth order for, respectively, the additive genetic and permanent environmental effects can be adopted, as little variation was observed in the genetic parameter estimates compared to those estimated by models using the B-spline function. Therefore, due to the lower complexity in the (co)variance estimations, the model using Legendre polynomials represented the best option for the genetic evaluation of the Guzera lactation records. An increase of 3.6% in the accuracy of the estimated breeding values was verified when using RRM. The ranks of animals were very close whatever the RRM for the data set used to predict breeding values. Considering P305, results indicated only small to medium difference in the animals' ranking based on breeding values predicted by the conventional model or by RRM. Therefore, the sum of all the RRM-predicted breeding values along the lactation period (RRM305) can be used as a selection criterion for 305-day milk production. (c) 2014 Elsevier B.V. All rights reserved.

The Fokker-Planck equation for a bistable potential

The Fokker-Planck equation is studied through its relation to a Schrodinger-type equation. The advantage of this combination is that we can construct the probability distribution of the Fokker-Planck equation by using well-known solutions of the Schrodinger equation. By making use of such a combination, we present the solution of the Fokker-Planck equation for a bistable potential related to a double oscillator. Thus, we can observe the temporal evolution of the system describing its dynamic properties such as the time tau to overcome the barrier. By calculating the rates k = 1/tau as a function of the inverse scaled temperature 1/D, where D is the diffusion coefficient, we compare the aspect of the curve k x 1/D, with the ones obtained from other studies related to four different kinds of activated process. We notice that there are similarities in some ranges of the scaled temperatures, where the different processes follow the Arrhenius behavior. We propose that the type of bistable potential used in this study may be used, qualitatively, as a simple model, whose rates share common features with the rates of some single rate-limited thermally activated processes. (C) 2014 Elsevier B.V. All rights reserved.

Information-entropic measure of energy-degenerate kinks in two-field models

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Lumped element modeling of operational structures by inverting the mobility models

This paper is concerned with what a source precisely sees when it drives a receiver such as a continuous structural object. An equivalent lumped element system consisting of masses, springs and dampers is developed to visually represent the operational structural dynamics of a single-input structure at the driving point. The development is solely based on the mobility model of the driving point response. The mobility model is mathematically inverted to give the impedance model that is suitable for lumped element modeling. The two types of structures studied are unconstrained inertial objects and constrained resilient objects. The lumped element systems presented suggest a new view to dynamics that a single-input flexible structure in operation can be decomposed into the two subsystems: a base system of single degree of freedom (or of a mass for an inertial object) whose mass is in contact with the source and an appendage system consisting of a series of oscillators each of which is attached to the base mass. The driving point response is a result of the coupling between the two subsystems. (C) 2014 Elsevier Ltd. All rights reserved.

Random regression models using different functions to model milk flow in dairy cows

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Optimized delta expansion for relativistic nuclear models

The optimized delta-expansion is a nonperturbative approach for field theoretic models which combines the techniques of perturbation theory and the variational principle. This technique is discussed in the lambda phi(4) model and then implemented in the Walecka model for the equation of state of nuclear matter. The results obtained with the delta expansion are compared with those obtained with the traditional mean field, relativistic Hartree and Hartree-Fock approximations.

Three-dimensional finite element analysis of anterior single implant-supported prostheses with different bone anchorages

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Representation of transmission lines: comparison of models and parameters distributed discrete parameters

A transmission line is characterized by the fact that its parameters are distributed along its length. This fact makes the voltages and currents along the line to behave like waves and these are described by differential equations. In general, the differential equations mentioned are difficult to solve in the time domain, due to the convolution integral, but in the frequency domain these equations become simpler and their solutions are known. The transmission line can be represented by a cascade of π circuits. This model has the advantage of being developed directly in the time domain, but there is a need to apply numerical integration methods. In this work a comparison of the model that considers the fact that the parameters are distributed (Universal Line Model) and the fact that the parameters considered concentrated along the line (π circuit model) using the trapezoidal integration method, and Simpson's rule Runge-Kutta in a single-phase transmission line length of 100 km subjected to an operation power. © 2003-2012 IEEE.

Influence of crown-to-implant ratio, retention system, restorative material, and occlusal loading on stress concentrations in single short implants

Purpose: The aim of this study was to assess the contributions of some prosthetic parameters such as crown-to-implant (C/I) ratio, retention system, restorative material, and occlusal loading on stress concentrations within a single posterior crown supported by a short implant. Materials and Methods: Computer-aided design software was used to create 32 finite element models of an atrophic posterior partially edentulous mandible with a single external-hexagon implant (5 mm wide × 7 mm long) in the first molar region. Finite element analysis software with a convergence analysis of 5% to mesh refinement was used to evaluate the effects of C/I ratio (1:1; 1.5:1; 2:1, or 2.5:1), prosthetic retention system (cemented or screwed), and restorative material (metal-ceramic or all ceramic). The crowns were loaded with simulated normal or traumatic occlusal forces. The maximum principal stress (σmax) for cortical and cancellous bone and von Mises stress (σvM) for the implant and abutment screw were computed and analyzed. The percent contribution of each variable to the stress concentration was calculated from the sum of squares analysis. Results: Traumatic occlusion and a high C/I ratio increased stress concentrations. The C/I ratio was responsible for 11.45% of the total stress in the cortical bone, whereas occlusal loading contributed 70.92% to the total stress in the implant. The retention system contributed 0.91% of the total stress in the cortical bone. The restorative material was responsible for only 0.09% of the total stress in the cancellous bone. Conclusion: Occlusal loading was the most important stress concentration factor in the finite element model of a single posterior crown supported by a short implant.

Predicting the behaviour of structures under impact loads using geometrically distorted scaled models

When a scaled structure (model or replica) is used to predict the response of a full-size compound (prototype), the model geometric dimensions should relate to the corresponding prototype dimensions by a single scaling factor. However, owing to manufacturing technical restrictions, this condition cannot be accomplished for some of the dimensions in real structures. Accordingly, the distorted geometry will not comply with the overall geometric scaling factor, infringing the Pi theorem requirements for complete dynamic similarity. In the present study, a method which takes geometrical distortions into account is introduced, leading to a model similar to the prototype. As a means to infer the performance of this method, three analytical problems of structures subjected to dynamic loads are analysed. It is shown that the replica developed applying this technique is able to accurately predict the full-size structure behaviour even when the studied models have some of their dimensions severely distorted. (C) 2012 Elsevier Ltd. All rights reserved.

Equation of Motion Governing the Dynamics of Vertically Collapsing Buildings

The present paper aims at contributing to a discussion, opened by several authors, on the proper equation of motion that governs the vertical collapse of buildings. The most striking and tragic example is that of the World Trade Center Twin Towers, in New York City, about 10 years ago. This is a very complex problem and, besides dynamics, the analysis involves several areas of knowledge in mechanics, such as structural engineering, materials sciences, and thermodynamics, among others. Therefore, the goal of this work is far from claiming to deal with the problem in its completeness, leaving aside discussions about the modeling of the resistive load to collapse, for example. However, the following analysis, restricted to the study of motion, shows that the problem in question holds great similarity to the classic falling-chain problem, very much addressed in a number of different versions as the pioneering one, by von Buquoy or the one by Cayley. Following previous works, a simple single-degree-of-freedom model was readdressed and conceptually discussed. The form of Lagrange's equation, which leads to a proper equation of motion for the collapsing building, is a general and extended dissipative form, which is proper for systems with mass varying explicitly with position. The additional dissipative generalized force term, which was present in the extended form of the Lagrange equation, was shown to be derivable from a Rayleigh-like energy function. DOI: 10.1061/(ASCE)EM.1943-7889.0000453. (C) 2012 American Society of Civil Engineers.

Constraining H-0 in general dark energy models from Sunyaev-Zeldovich/X-ray technique and complementary probes

In accelerating dark energy models, the estimates of the Hubble constant, Ho, from Sunyaev-Zerdovich effect (SZE) and X-ray surface brightness of galaxy clusters may depend on the matter content (Omega(M)), the curvature (Omega(K)) and the equation of state parameter GO. In this article, by using a sample of 25 angular diameter distances of galaxy clusters described by the elliptical beta model obtained through the SZE/X-ray technique, we constrain Ho in the framework of a general ACDM model (arbitrary curvature) and a flat XCDM model with a constant equation of state parameter omega = p(x)/rho(x). In order to avoid the use of priors in the cosmological parameters, we apply a joint analysis involving the baryon acoustic oscillations (BA()) and the (MB Shift Parameter signature. By taking into account the statistical and systematic errors of the SZE/X-ray technique we obtain for nonflat ACDM model H-0 = 74(-4.0)(+5.0) km s(-1) Mpc(-1) (1 sigma) whereas for a fiat universe with constant equation of state parameter we find H-0 = 72(-4.0)(+5.5) km s(-1) Mpc(-1)(1 sigma). By assuming that galaxy clusters are described by a spherical beta model these results change to H-0 = 6(-7.0)(+8.0) and H-0 = 59(-6.0)(+9.0) km s(-1) Mpc(-1)(1 sigma), respectively. The results from elliptical description are in good agreement with independent studies from the Hubble Space Telescope key project and recent estimates based on the Wilkinson Microwave Anisotropy Probe, thereby suggesting that the combination of these three independent phenomena provides an interesting method to constrain the Bubble constant. As an extra bonus, the adoption of the elliptical description is revealed to be a quite realistic assumption. Finally, by comparing these results with a recent determination for a, flat ACDM model using only the SZE/X-ray technique and BAO, we see that the geometry has a very weak influence on H-0 estimates for this combination of data.

A method to investigate the shrinkage stress developed by resin-composites bonded to a single flat surface

Objectives. To purpose a method for predicting the shrinkage stress development in the adhesive layer of resin-composite cylinders that shrink bonded to a single flat surface, by measuring the deflection of a glass coverslip caused by the shrinkage of the bonded cylinders. The correlation between the volume of the bonded resin-composite and the stress-peak was also investigated. Methods. A glass coverslip deflection caused by the shrinkage of a bonded resin-composite cylinder (diameter: d = 8 mm, 4 mm, or 2 mm, height: h = 4 mm, 2 mm, 1 mm, or 0.5 mm) was measured, and the same set-up was simulated by finite element analysis (3D-FEA). Stresses generated in the adhesive layer were plotted versus two geometric variables of the resin-composite cylinder (C-Factor and volume) to verify the existence of correlations between them and stresses. Results. The FEA models were validated. A significant correlation (p < 0.01, Pearson's test) between the stress-peak and the coverslip deflection when the resin-composites were grouped by diameter was found for diameters of 2 and 4 mm. The stress-peak of the whole set of data showed a logarithmic correlation with the bonded resin-composite volume (p < 0.001, Pearson's test), but did not correlate with the C-Factor. Significance. The described method should be considered for standardizing the stress generated by the shrinkage of resin-composite blocks bonded to a single flat surface. (C) 2012 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Twinlike models for kinks and compactons in flat and warped spacetime

This work deals with the presence of twinlike models in scalar field theories. We show how to build distinct scalar field theories having the same extended solution, with the same energy density and linear stability. Here, however, we start from a given but generalized scalar field theory, and we construct the corresponding twin model, which also engenders generalized dynamics. We investigate how the twinlike models arise in both flat and curved spacetimes. In the curved spacetime, we consider a braneworld model with the warp factor controlling the spacetime geometry with a single extra dimension of infinite extent. In particular, we study linear stability in both flat and curved spacetimes, and in the case of curved spacetime-in both the gravity and the scalar field sectors-for the two braneworld models. DOI: 10.1103/PhysRevD.86.125021

Influence of the magnetization damping on dynamic hysteresis loops in single domain particles

This article reports on the influence of the magnetization damping on dynamic hysteresis loops in single-domain particles with uniaxial anisotropy. The approach is based on the Neel-Brown theory and the hierarchy of differential recurrence relations, which follow from averaging over the realizations of the stochastic Landau-Lifshitz equation. A new method of solution is proposed, where the resulting system of differential equations is solved directly using optimized algorithms to explore its sparsity. All parameters involved in uniaxial systems are treated in detail, with particular attention given to the frequency dependence. It is shown that in the ferromagnetic resonance region, novel phenomena are observed for even moderately low values of the damping. The hysteresis loops assume remarkably unusual shapes, which are also followed by a pronounced reduction of their heights. Also demonstrated is that these features remain for randomly oriented ensembles and, moreover, are approximately independent of temperature and particle size. (C) 2012 American Institute of Physics. [doi:10.1063/1.3684629]