Un algorithme efficace d'intégration plastique pour un matériau obéissant au critère anisotrope de HillAn efficient algorithm for plastic integration of material satisfying the Hill's anisotropic yield criterion.

This Note deals with an efficient algorithm to carry out the plastic integration and compute the stresses due to large strains for materials satisfying the Hill's anisotropic yield criterion. The classical algorithm of plastic integration such as ‘Return Mapping Method’ is largely used for nonlinear analyses of structures and numerical simulations of forming processes, but it requires an iterative schema and may have convergence problems. A new direct algorithm based on a scalar method is developed which allows us to directly obtain the plastic multiplier without an iteration procedure; thus the computation time is largely reduced and the numerical problems are avoided. To cite this article: I. Titeux et al., C. R. Mecanique 332 (2004).     

Analysis of HIV/AIDS Epidemic and Socioeconomic Factors in Sub-Saharan Africa

Sub-Saharan Africa has been the epicenter of the outbreak since the spread of acquired immunodeficiency syndrome (AIDS) began to be prevalent. This article proposes several regression models to investigate the relationships between the HIV/AIDS epidemic and socioeconomic factors (the gross domestic product per capita, and population density) in ten countries of Sub-Saharan Africa, for 2011–2016. The maximum likelihood method was used to estimate the unknown parameters of these models along with the Newton–Raphson procedure and Fisher scoring algorithm. Comparing these regression models, there exist significant spatiotemporal non-stationarity and auto-correlations between the HIV/AIDS epidemic and two socioeconomic factors. Based on the empirical results, we suggest that the geographically and temporally weighted Poisson autoregressive (GTWPAR) model is more suitable than other models, and has the better fitting results.     

Optimum replacement policies for systems subject to shocks

Summary: A system experiences a shock of magnitude γ _i at time τ _i . Each shock deteriorate system to some extent and due to the decrease in efficiency, the system becomes more expensive to run. Assuming that the shock process is a general birth process and the cost structure depends on the magnitude of the shocks and time, the optimum replacement period of the system has been derived Optimum replacement periods for particular cases of general birth process are discussed in detail with suitable examples.     

A numerical study of fluids with pressure‐dependent viscosity flowing through a rigid porous medium

Much of the work on flow through porous media, especially with regard to studies on the flow of oil, are based on ‘Darcy's law’ or modifications to it, such as Darcy–Forchheimer or Brinkman models. While many theoretical and numerical studies concerning flow through porous media have taken into account the inhomogeneity and anisotropy of the porous solid, they have not taken into account the fact that the viscosity of the fluid and drag coefficient could depend on the pressure in applications, such as enhanced oil recovery (EOR). Experiments clearly indicate that the viscosity varies exponentially with respect to the pressure and the viscosity can change, in some applications, by several orders of magnitude. The fact that the viscosity depends on pressure immediately implies that the ‘drag coefficient’ would also depend on the pressure. In this paper we consider modifications to Darcy's equation wherein the drag coefficient is a function of pressure, which is a realistic model for technological applications, such as EOR and geological carbon sequestration. We first outline the approximations behind Darcy's equation and the modifications that we propose to Darcy's equation, and derive the governing equations through a systematic approach using mixture theory. We then propose a stabilized mixed finite element formulation for the modified Darcy's equation. To solve the resulting nonlinear equations we present a solution procedure based on the consistent Newton–Raphson method. We solve representative test problems to illustrate the performance of the proposed stabilized formulation. One of the objectives of this paper is to show that the dependence of viscosity on the pressure can have a significant effect both on the qualitative and quantitative nature of the solution. Copyright © 2010 John Wiley & Sons, Ltd.     

Remarks on the updated Hessian matrix methods

Optimizing a function with respect to a set of variables using the quasi‐Newton–Raphson method implies updating the Hessian matrix at each iteration. The Broyden–Fletcher–Goldfarb–Shanno update formula is used for minimization and the Murtagh–Sargent–Powell update formula for optimization of first‐order saddle points. Two new formulae are proposed to update the Hessian matrix. One of these formulae is derived using exponential weights and should be used to locate first‐order saddle points. The second formula is a modification of the TS–Broyden–Fletcher–Goldfarb–Shanno update and could used for both minimum and first‐order saddle point optimizations. These two update Hessian matrix formulae present a performance that is the same and in many cases better that the Broyden–Fletcher–Goldfarb–Shanno and Murtagh–Sargent–Powell formulae. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 324–332, 2003     

Shape-preserving interpolation by G1 and G2 PH quintic splines

The interpolation of a planar sequence of points p₀, ..., p_Nby shape-preserving G¹ or G² PH quintic splines with specifiedend conditions is considered. The shape-preservation propertyis secured by adjusting ‘tension’ parameters thatarise upon relaxing parametric continuity to geometric continuity.In the G² case, the PH spline construction is based on applyingNewton–Raphson iterations to a global system of equations,commencing with a suitable initialization strategy—thisgeneralizes the construction described previously in NumericalAlgorithms 27, 35–60 (2001). As a simpler and cheaperalternative, a shape-preserving G¹ PH quintic spline schemeis also introduced. Although the order of continuity is lower,this has the advantage of allowing construction through purelylocal equations.     

A Geometric Approach to Maximum Likelihood Estimation of the Functional Principal Components From Sparse Longitudinal Data

In this article, we consider the problem of estimating the eigenvalues and eigenfunctions of the covariance kernel (i.e., the functional principal components) from sparse and irregularly observed longitudinal data. We exploit the smoothness of the eigenfunctions to reduce dimensionality by restricting them to a lower dimensional space of smooth functions. We then approach this problem through a restricted maximum likelihood method. The estimation scheme is based on a Newton–Raphson procedure on the Stiefel manifold using the fact that the basis coefficient matrix for representing the eigenfunctions has orthonormal columns. We also address the selection of the number of basis functions, as well as that of the dimension of the covariance kernel by a second-order approximation to the leave-one-curve-out cross-validation score that is computationally very efficient. The effectiveness of our procedure is demonstrated by simulation studies and an application to a CD4+ counts dataset. In the simulation studies, our method performs well on both estimation and model selection. It also outperforms two existing approaches: one based on a local polynomial smoothing, and another using an EM algorithm. Supplementary materials including technical details, the R package fpca, and data analyzed by this article are available online.     

前馈神经网络和Logit回归的比较研究

近年来 ,前馈神经网络广泛地应用在 Logit回归作为标准统计方法的分析领域 .但却很少作它们之间的直接比较 ,本文是 Logit回归和前馈神经网络“比较研究”的一个尝试 ,说明了一些理论结果和特性 ,讨论了在它们的应用中碰到的一些实际问题 ,还进一步用分析的和模拟的两种方法研究了一些重要的渐近概念、过分拟合以及模型选择等问题 ,最后讨论并给出一些结论     

Globally convergent computation of chemical equilibrium composition

We report the Newton-Raphson based globally convergent computational method for determination of chemical equilibrium composition. In the computation of chemical equilibrium composition, an appearance of nonpositive value of number of moles of any component leads to discrepancy. The process of conditional backtracking and adaptive set of refining factors for Newton-Raphson steps are employed to resolve the problem. The mathematical formulation proposed by Heuze et al. (J Chem Phys 1985, 83, 4734) has been solved using proposed computational method, instead of empirical iterative formulation, as proposed by them. Results for the same numerical example, used by Heuze et al. (J Chem Phys 1985, 83, 4734) and White et al. (J Chem Phys 1958, 28, 751) are presented in addition to decomposition of Cyclotrimethylenetrinitramine for fixed temperature and pressure. It is observed that the proposed method is efficient and globally convergent. An even noteworthy finding is that the set of refining factors can be chosen from the range 0.1 to eta, where eta may be greater than one depending on how smoothly system of nonlinear equations is dependant on corresponding variable. Related analysis and results are discussed.     

Load-Sharing Model under Lindley Distribution and Its Parameter Estimation Using the Expectation-Maximization Algorithm

A load-sharing system is defined as a parallel system whose load will be redistributed to its surviving components as each of the components fails in the system. Our focus is on making statistical inference of the parameters associated with the lifetime distribution of each component in the system. In this paper, we introduce a methodology which integrates the conventional procedure under the assumption of the load-sharing system being made up of fundamental hypothetical latent random variables. We then develop an expectation maximization algorithm for performing the maximum likelihood estimation of the system with Lindley-distributed component lifetimes. We adopt several standard simulation techniques to compare the performance of the proposed methodology with the Newton–Raphson-type algorithm for the maximum likelihood estimate of the parameter. Numerical results indicate that the proposed method is more effective by consistently reaching a global maximum.