Legendre Orthogonal Polynomials Research Articles

Guided thermoelastic wave propagation in layered plates without energy dissipation

In this paper, the propagation of guided thermoelastic waves in laminated orthotropic plates subjected to stress-free, isothermal boundary conditions is investigated in the context of the Green-Naghdi (GN) generalized thermoelastic theory (without energy dissipation). The coupled wave equations and heat conduction equation are solved by the Legendre orthogonal polynomial series expansion approach. The validity of the method is confirmed through a comparison. The dispersion curves of thermal modes and elastic modes are illustrated simultaneously. Dispersion curves of the corresponding pure elastic plate are also shown to analyze the influence of the thermoelasticity on elastic modes. The displacement and temperature distributions are shown to discuss the differences between the elastic modes and thermal modes.

Estimation of (co)variance components and genetic parameters for weights of red-winged tinamou using random regression models

The objective of this work was to determine genetic parameters for body weight of tinamou in captivity. It was used random regression models in analyses of data by considering the direct additive genetic (DA) and permanent environmental effects of the animal (PE) as random effects. Residual variances were modeled by using a fifth-order variance function. The mean population growth curve was fitted by sixth-order Legendre orthogonal polynomials. Direct additive genetic effects and animal environmental permanent effect were modeled by using Legendre polynomials of order two to nine. The best results were obtained by models with orders of fit of 6 for direct additive genetic effect and of order 3 for permanent effect by Akaike information criterion and of order 3 for both additive genetic effect and permanent effect by Schwarz Bayesian information criterion and likelihood ratio test. Heritability estimates ranged from 0.02 to 0.57. The first eigenvalue explained 94% and 90% of the variation from additive direct and permant environmental effects, respectively. Selection of tinamou for body weight is more effective after 112 days of age.

Vibration analysis of a cracked rotating tapered beam using the p-version finite element method

Vibration characteristics of cracked rotating tapered beam are investigated by the p-version finite element method. In the present formulation, the shape functions enriched with the shifted Legendre orthogonal polynomials are employed to represent the transverse displacement field within the rotating tapered beam element. The crack element stiffness matrix and the p-version finite element model of the structural system are obtained by using fracture mechanics and the Lagrange equation, respectively. The common integral terms appeared in the rotating tapered beam element matrices are derived in analytical forms so as to speed up the generation of these matrices during the implementation of the present method. Comparisons have been made for several configurations between results obtained by the present method with those available in literature, and the proposed method shows excellent accuracy and convergence. The effects of crack location, crack size, rotating speed and hub radius on vibration characteristics of a cracked rotating tapered beam are investigated. The mode shapes of the cracked rotating tapered beam are also obtained and analyzed with the spatial wavelet transform approach to detect the slight perturbation at the crack position. The present method is of general applicability and can be applied to investigate vibration of arbitrarily continuous non-uniform rotating beams.

Genome-Wide Association Studies for Bivariate Sparse Longitudinal Data

Objective: Longitudinal measurements with bivariate response have been analyzed by several authors using two separate models for each response. However, for most of the biological or medical experiments, the two responses are highly correlated and hence a separate model for each response might not be a desirable way to analyze such data. A single model considering a bivariate response provides a more powerful inference as the correlation between the responses is modeled appropriately. In this article, we propose a dynamic statistical model to detect the genes controlling human blood pressure (systolic and diastolic). Methods: By modeling the mean function with orthogonal Legendre polynomials and the covariance matrix with a stationary parametric structure, we incorporate the statistical ideas in functional genome-wide association studies to detect SNPs which have significant control on human blood pressure. The traditional false discovery rate is used for multiple comparisons. Results: We analyze the data from the Framingham Heart Study to detect such SNPs by appropriately considering gender-gene interaction. We detect 8 SNPs for males and 7 for females which are most significant in controlling blood pressure. The genotype-specific mean curves and additive and dominant effects over time are shown for each significant SNP for both genders. Simulation studies are performed to examine the statistical properties of our model. The current model will be extremely useful in detecting genes controlling different traits and diseases for humans or non-human subjects.

Random regression models using different functions to model test-day milk yield of Brazilian Holstein cows

We analyzed 152,145 test-day records from 7317 first lactations of Holstein cows recorded from 1995 to 2003. Our objective was to model variations in test-day milk yield during the first lactation of Holstein cows by random regression model (RRM), using various functions in order to obtain adequate and parsimonious models for the estimation of genetic parameters. Test-day milk yields were grouped into weekly classes of days in milk, ranging from 1 to 44 weeks. The contemporary groups were defined as herd-test-day. The analyses were performed using a single-trait RRM, including the direct additive, permanent environmental and residual random effects. In addition, contemporary group and linear and quadratic effects of the age of cow at calving were included as fixed effects. The mean trend of milk yield was modeled with a fourth-order orthogonal Legendre polynomial. The additive genetic and permanent environmental covariance functions were estimated by random regression on two parametric functions, Ali and Schaeffer and Wilmink, and on B-spline functions of days in milk. The covariance components and the genetic parameters were estimated by the restricted maximum likelihood method. Results from RRM parametric and B-spline functions were compared to RRM on Legendre polynomials and with a multi-trait analysis, using the same data set. Heritability estimates presented similar trends during mid-lactation (13 to 31 weeks) and between week 37 and the end of lactation, for all RRM. Heritabilities obtained by multi-trait analysis were of a lower magnitude than those estimated by RRM. The RRMs with a higher number of parameters were more useful to describe the genetic variation of test-day milk yield throughout the lactation. RRM using B-spline and Legendre polynomials as base functions appears to be the most adequate to describe the covariance structure of the data.

Statistical modelling of growth using a mixed model with orthogonal polynomials

In statistical modelling, the effects of single-nucleotide polymorphisms (SNPs) are often regarded as time-independent. However, for traits recorded repeatedly, it is very interesting to investigate the behaviour of gene effects over time. In the analysis, simulated data from the 13th QTL-MAS Workshop (Wageningen, The Netherlands, April 2009) was used and the major goal was the modelling of genetic effects as time-dependent. For this purpose, a mixed model which describes each effect using the third-order Legendre orthogonal polynomials, in order to account for the correlation between consecutive measurements, is fitted. In this model, SNPs are modelled as fixed, while the environment is modelled as random effects. The maximum likelihood estimates of model parameters are obtained by the expectation–maximisation (EM) algorithm and the significance of the additive SNP effects is based on the likelihood ratio test, with p-values corrected for multiple testing. For each significant SNP, the percentage of the total variance contributed by this SNP is calculated. Moreover, by using a model which simultaneously incorporates effects of all of the SNPs, the prediction of future yields is conducted. As a result, 179 from the total of 453 SNPs covering 16 out of 18 true quantitative trait loci (QTL) were selected. The correlation between predicted and true breeding values was 0.73 for the data set with all SNPs and 0.84 for the data set with selected SNPs. In conclusion, we showed that a longitudinal approach allows for estimating changes of the variance contributed by each SNP over time and demonstrated that, for prediction, the pre-selection of SNPs plays an important role.

A Discrete Generalized Multigroup Energy Expansion Theory

This study describes the generalized multigroup energy treatment for the neutron transport equation. Discrete Legendre orthogonal polynomials (DLOPs) are used to expand the energy dependence of the angular flux into a set of flux moments. The leading (zeroth)-order equation is identical to a standard multigroup solution, while the higher-order equations are decoupled from each other and only depend on the leading-order solution because of the orthogonality property of the DLOPs. This decoupling leads to computational times comparable to the coarse-group calculation but provides an accurate fine-group energy spectrum. One-dimensional single-assembly and core calculations were performed to demonstrate the potential of the discrete generalized multigroup method. Computational results show that the discrete generalized multigroup method can produce an accurate fine-group whole-core solution for less computational time. A source update process is also introduced that provides improvement of integral quantities such as eigenvalue and reaction rates over the coarse-group solution.

Orthogonal polynomial expansions for the matrix exponential

Many different algorithms have been suggested for computing the matrix exponential. In this paper, we put forward the idea of expanding in either Chebyshev, Legendre or Laguerre orthogonal polynomials. In order for these expansions to converge quickly, we cluster the eigenvalues into diagonal blocks and accelerate using shifting and scaling.

Generalized thermoelastic waves in functionally graded plates without energy dissipation

Abstract In this paper, a dynamic solution of the propagating thermoelastic waves in functionally graded material (FGM) plate subjected to stress-free, isothermal boundary conditions is presented in the context of the Green–Naghdi (GN) generalized thermoelastic theory. The FGM plate is composed of two orthotropic materials. The materials properties are assumed to vary in the direction of the thickness according to a known variation law. The coupled wave equation and heat conduction equation are solved by the Legendre orthogonal polynomial series expansion approach. The convergency of the method is discussed through a numerical example. The dispersion curves of the inhomogeneous thermoelastic plate and the corresponding pure elastic plate are compared to show the characteristics of thermal modes and the influence of the thermoelasticity on elastic modes. The displacement, temperature and stress distributions of elastic modes and thermal modes are shown to discuss their differences. A plate with a different gradient variation is calculated to illustrate the influence of the gradient field on the wave characteristics.

Wave Characteristics in Magneto-electro-elastic Functionally Graded Spherical Curved Plates

Piezoelectric-piezomagnetic functionally graded materials (FGM), with a gradual change of the mechanical, electric and magnetic properties, have great promise for different applications. Based on the Legendre orthogonal polynomial series expansion approach, a dynamic solution is presented for the propagation of harmonic waves in piezoelectric-piezomagnetic FGM spherical curved plates. The material properties are assumed to vary in the direction of the thickness according to a known variation law. The dispersion curves of the piezoelectric-piezomagnetic FGM spherical curved plate and the corresponding non-piezoelectric, non-piezomagnetic plates are calculated to show the influences of the piezoelectricity and piezomagnetism. Electric potential and magnetic potential distributions are also obtained to illustrate the different influences of the piezoelectricity and piezomagnetism. Finally, a spherical curved plate at different ratio of radius to thickness is calculated to show the influence f the ratio on the piezoelectric effect and piezomagnetic effect.

Estimação de funções de covariância para características de crescimento da raça Tabapuã utilizando modelos de regressão aleatória

Foram utilizados 35.732 registros de peso do nascimento aos 660 dias de idade de 8.458 animais da raça Tabapuã para estimar funções de covariância utilizando modelos de regressão aleatória sobre polinômios de Legendre. Os modelos incluíram: como aleatórios, os efeitos genético aditivo direto, materno, de ambiente permanente de animal e materno; como fixos, os efeitos de grupo de contemporâneo; como covariáveis, a idade do animal à pesagem e a idade da vaca ao parto (linear e quadrática); e sobre a idade à pesagem, polinômio ortogonal de Legendre (regressão cúbica) foi considerado para modelar a curva média da população. O resíduo foi modelado considerando sete classes de variância e os modelos foram comparados pelos critérios de informação Bayesiano de Schwarz e Akaike. O melhor modelo apresentou ordens 4, 3, 6, 3 para os efeitos genético aditivo direto e materno, de ambiente permanente de animal e materno, respectivamente. As estimativas de covariância e herdabilidades, obtidas utilizando modelo bicaracter, e de regressão aleatória foram semelhantes. As estimativas de herdabilidade para o efeito genético aditivo direto, obtidas com o modelo de regressão aleatória, aumentaram do nascimento (0,15) aos 660 dias de idade (0,45). Maiores estimativas de herdabilidade materna foram obtidas para pesos medidos logo após o nascimento. As correlações genéticas variaram de moderadas a altas e diminuíram com o aumento da distância entre as pesagens. A seleção para maiores pesos em qualquer idade promove maior ganho de peso do nascimento aos 660 dias de idade.

Orthogonal Polynomials and Related Special Functions Applied in Geosciences and Engineering Computations

In applications of mathematics involving either the Laplace or the Helmholtz equation in spherical coordinates the associated Legendre equation occurs. Its solutions are called associated Legendre functions. They have some relations to classical Legendre orthogonal polynomials and classical Jacobi orthogonal polynomials. The classical Legendre polynomials and the associated Legendre functions have been used by some authors at the approximation of length of curves and cylindrical surfaces and the approximation of Earth shape in geosciences and engineering computations. In the present paper we investigate some relations between these functions pointing out possible applications in geosciences.

On the efficiency of treating singularities in triatomic variational vibrational computations. The vibrational states of H+3 up to dissociation

Several techniques of varying efficiency are investigated, which treat all singularities present in the triatomic vibrational kinetic energy operator given in orthogonal internal coordinates of the two distances-one angle type. The strategies are based on the use of a direct-product basis built from one-dimensional discrete variable representation (DVR) bases corresponding to the two distances and orthogonal Legendre polynomials, or the corresponding Legendre-DVR basis, corresponding to the angle. The use of Legendre functions ensures the efficient treatment of the angular singularity. Matrix elements of the singular radial operators are calculated employing DVRs using the quadrature approximation as well as special DVRs satisfying the boundary conditions and thus allowing for the use of exact DVR expressions. Potential optimized (PO) radial DVRs, based on one-dimensional Hamiltonians with potentials obtained by fixing or relaxing the two non-active coordinates, are also studied. The numerical calculations employed Hermite-DVR, spherical-oscillator-DVR, and Bessel-DVR bases as the primitive radial functions. A new analytical formula is given for the determination of the matrix elements of the singular radial operator using the Bessel-DVR basis. The usually claimed failure of the quadrature approximation in certain singular integrals is revisited in one and three dimensions. It is shown that as long as no potential optimization is carried out the quadrature approximation works almost as well as the exact DVR expressions. If wave functions with finite amplitude at the boundary are to be computed, the basis sets need to meet the required boundary conditions. The present numerical results also confirm that PO-DVRs should be constructed employing relaxed potentials and PO-DVRs can be useful for optimizing quadrature points for calculations applying large coordinate intervals and describing large-amplitude motions. The utility and efficiency of the different algorithms is demonstrated by the computation of converged near-dissociation vibrational energy levels for the H molecular ion.

Non-linear free vibration of elliptic sector plates by a curved triangular p-element

A curved triangular p-element is developed and applied to geometrically non-linear free vibration of isotropic elliptic sector plates. The formulation takes into account shear deformation, rotary inertia, and geometric non-linearity. The element can describe the geometry of the elliptic sector plate exactly and is therefore suitable for this type of plate. The shifted Legendre orthogonal polynomials are used as enriching hierarchic shape functions. The element stiffness and mass matrices are integrated numerically using the Gauss–Legendre quadrature. The non-linear equations of free motion are obtained using the harmonic balance method and solved iteratively by the linearized updated mode method. Numerical results for the linear and non-linear frequencies of clamped elliptic sector plates are obtained. The efficiency of the proposed p-element is demonstrated through convergence studies and comparisons with published results. The effects of the ellipse eccentricity, thickness ratio, and sector angle on the backbone curves are examined. It is shown that these parameters influence the hardening behavior.

Analysis of Feed Intake of Individually Fed Beef Bulls Using Random Regression Models

Selapa, N.W., Nephawe, K.A., Maiwashe, A., Norris, D. and Ngambi, J.W. 2009. Analysis of feed intake of individually fed beef bulls using random regression models. J. Appl. Anim. Res., 36: 223–226. To estimate genetic parameters for weekly feed intake (FI) of beef bulls records of FI on Bonsmara bulls (n = 3 012) for 5 years were analyzed using random regression model (RRM). The model included the fixed linear regression of cumulative feed intake on third order orthogonal Legendre polynomials (OLP) of the actual days on test (d 7 to 84 weekly) for starting age group and contemporary group effects. Random regressions on third order OLP were included for additional random effect of weaning-herd-year (WHY) effect and additive genetic effect of the animal and on fourth order for the additional random permanent environmental effect of the animal. The residual effects were assumed to be independently distributed with heterogeneous variance for each measurement period. The first additive genetic eigenfunction (with 95.6% additive genetic variance contribution) revealed that selection for decreased cumulative FI at any test day would tend to decrease cumulative feed intake throughout the test and vice versa. Random regression model could be used in genetic evaluation of cumulative feed intakes of individually fed beef bulls in South Africa.

Covariance functions for body weight from birth to maturity in Nellore cows1

The objective of this study was to estimate (co)variance functions using random regression models on Legendre polynomials for the analysis of repeated measures of BW from birth to adult age. A total of 82,064 records from 8,145 females were analyzed. Different models were compared. The models included additive direct and maternal effects, and animal and maternal permanent environmental effects as random terms. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of animal age (cubic regression) were considered as random covariables. Eight models with polynomials of third to sixth order were used to describe additive direct and maternal effects, and animal and maternal permanent environmental effects. Residual effects were modeled using 1 (i.e., assuming homogeneity of variances across all ages) or 5 age classes. The model with 5 classes was the best to describe the trajectory of residuals along the growth curve. The model including fourth- and sixth-order polynomials for additive direct and animal permanent environmental effects, respectively, and third-order polynomials for maternal genetic and maternal permanent environmental effects were the best. Estimates of (co)variance obtained with the multi-trait and random regression models were similar. Direct heritability estimates obtained with the random regression models followed a trend similar to that obtained with the multi-trait model. The largest estimates of maternal heritability were those of BW taken close to 240 d of age. In general, estimates of correlation between BW from birth to 8 yr of age decreased with increasing distance between ages.

Generalized thermoelastic waves in spherical curved plates without energy dissipation

In this article, the propagation of thermoelastic waves in orthotropic spherical curved plates subjected to stress-free, isothermal boundary conditions is investigated in the context of the Green–Naghdi (GN) generalized thermoelastic theory (without energy dissipation). The theoretical formulation is based on the linear GN thermoelastic theory. The coupled wave equation and heat conduction equation expressed by the displacement and temperature are obtained. By the Legendre orthogonal polynomial series expansion approach, the coupled controlling equations are solved. The convergence of the method is demonstrated through a numerical example. The dispersion curves of thermal modes and elastic modes are illustrated simultaneously. Dispersion curves of the corresponding purely elastic spherical plate are also shown to analyze the influence of thermoelasticity on elastic modes. The displacement, temperature and stress distributions of both elastic modes and thermal modes are calculated to show their differences. A thermoelastic spherical plate with a different ratio of radius to thickness is considered to show the influence of the ratio on the characteristics of thermoelastic waves.

Circumferential thermoelastic waves in orthotropic cylindrical curved plates without energy dissipation

In this article, the propagation of guided thermoelastic waves in the circumferential direction of orthotropic cylindrical curved plates subjected to stress-free, isothermal boundary conditions is investigated in the context of the Green–Naghdi (GN) generalized thermoelastic theory (without energy dissipation). The coupled wave equations and heat conduction equation are solved by the Legendre orthogonal polynomial series expansion approach. The convergency of the method is discussed through a numerical example. The dispersion curves of thermal modes and elastic modes are illustrated simultaneously. Dispersion curves of the corresponding pure elastic cylindrical plate are also shown to analyze the influence of the thermoelasticity on elastic modes. The displacement, temperature and stress distributions are shown to discuss the differences between the elastic modes and thermal modes. A thermoelastic cylindrical plate with a different ratio of radius to thickness is considered to discuss the influence of the ratio on the characteristics of circumferential thermoelastic waves.

Indentation of a round punch into a rough elastic half-space

We consider the axisymmetric contact problem on the indentation of a round punch with plane base into an elastic half-space. Roughness is taken into account in the framework of the well-known model in which the local surface displacements are proportional to some power of the contact pressures at the same point. To study the problem, we use the theory of nonlinear integral equations of Hammerstein type together with the mathematical technique of orthogonal Legendre polynomials and the contraction mapping principle. We obtain an approximate analytic solution of the problem, numerically analyze the results, and reveal typical laws of variation of the main mechanical variables.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on 1/6/2010orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.