A self-optimizing inverse analysis method for estimation of cyclic elasto-plasticity model parameters

In this paper, a novel inverse analysis methodology call a Self-Optimizing Inverse Method (Self-OPTIM) has been presented, which inversely estimates cyclic elasto-plastic constitutive model parameters using global forces and displacement on the same partial boundaries and full-(or partial-) field displacement data. A novelty of the methodology is that it automatically self-estimates material parameters by updating “full-field” reference stresses and strains through two parallel nonlinear finite element simulations. Although a well-known classical cyclic plasticity model is chosen in this paper, it must be emphasized that the proposed Self-OPTIM method is a model-independent method, which means that any advanced model can be naturally integrated with the proposed methodology. Thus, using numerically generated test data of low-carbon steel specimens (AISI 1010), the proposed Self-OPTIM method has been verified showing its successful performance to estimate nonlinear isotropic and kinematic hardening parameters, yield stress, Young’s modulus and Poisson ratio. The effects of experimental noises from CCD camera and measurement errors of the boundary forces are also investigated for the Self-OPTIM method.     

关于求解带线性互补约束的数学规划问题正则方法的一个注记

本文主要研究在某些较弱条件下求解带线性互补约束的数学规划问题(MPLCC)正则方法的收敛性.若衡约束规划线性独立约束规范条件(MPEC-LICQ)在由正则方法产生的点列的聚点处成立,且迭代点列满足二阶必要条件,同时,若比在文[7]中渐近弱非退化条件Ⅰ更弱的渐近弱非退化条件Ⅱ在聚点处也成立,那么所有这些聚点都是B-稳定点.此外,在弱MPEC-LICQ,二阶必要条件及渐近弱退化条件Ⅱ下,我们仍能证明通过正则方法所得的聚点都是B-稳定点.     

On a nonlocal Cahn-Hilliard/Navier-Stokes system with degenerate mobility and singular potential for incompressible fluids with different densities

We consider a diffuse interface model describing flow and phase separation of a binary isothermal mixture of (partially) immiscible viscous incompressible Newtonian fluids having different densities. The model is the nonlocal version of the one derived by Abels, Garcke and Grün and consists in a inhomogeneous Navier-Stokes type system coupled with a convective nonlocal Cahn-Hilliard equation. This model was already analyzed in a paper by the same author, for the case of singular potential and non-degenerate mobility. Here, we address the physically more relevant situation of degenerate mobility and we prove existence of global weak solutions satisfying an energy inequality. The proof relies on a regularization technique based on a careful approximation of the singular potential. Existence and regularity of the pressure field is also discussed. Moreover, in two dimensions and for slightly more regular solutions, we establish the validity of the energy identity. We point out that in none of the existing contributions dealing with the original (local) Abels, Garcke Grün model, an energy identity in two dimensions is derived (only existence of weak solutions has been proven so far).     

关于两类污染数据回归分析的参数估计

研究简单回归模型：ｙｊ＝ａ＋βｘｊ＋εｊ，ｊ＝１，２，…，ｎ，其中Ｅεｊ＝０，Ｅε＾２ｊ＝σ＾２ｊ；但ｙ１，ｙ２，…，ｙｎ受到另一独立同分布随机变量序列ｔ，ｔ２，…，ｔｎ两种不同方式的污染，ｔｊ与ｙｊ独立。本文给出了两种污染方式下的ａ、β和污染参数的估计。     

Approximation of the steepest descent direction for the O-D matrix adjustment problem

In this paper, a method to approximate the directions of Clarke's generalized gradient of the upper level function for the demand adjustment problem on traffic networks is presented. Its consistency is analyzed in detail. The theoretical background on which this method relies is the known property of proximal subgradients of approximating subgradients of proximal bounded and lower semicountinuous functions using the Moreau envelopes. A double penalty approach is employed to approximate the proximal subgradients provided by these envelopes. An algorithm based on partial linearization is used to solve the resulting nonconvex problem that approximates the Moreau envelopes, and a method to verify the accuracy of the approximation to the steepest descent direction at points of differentiability is developed, so it may be used as a suitable stopping criterion. Finally, a set of experiments with test problems are presented, illustrating the approximation of the solutions to a steepest descent direction evaluated numerically. Research supported under Spanish CICYT project TRA99-1156-C02-02.     

Logistic模型参数估计及预测实例

本文提出了对Ｌｏｇｉｓｔｉｃ模型中的参数进行迭代估计的新算法，通过比较分析，说明了本文算法的有效性。     

Implicit Regularization and Renormalization of QCD

We apply the Implicit Regularization (IR) technique in a non-abelian gauge theory. We show that IR preserves gauge symmetry as encoded in relations between the renormalization constants required by the Slavnov–Taylor identities at the one-loop level of QCD. Moreover, we show that the technique handles divergencies in massive and massless QFT on equal footing.     

Identification of shallow sea models

In this paper we consider a parameter estimation procedure for shallow sea models. The method is formulated as a minimization problem. An adjoint model is used to calculate the gradient of the criterion which is to be minimized. In order to obtain a robust estimation method, the uncertainty of the open boundary conditions can be taken into acoount by allowing random noise inputs to act on the open boundaries. This method avoids the possibility that boundary errors are interpreted by the estimation procedure as parameter fluctuations. We apply the parameter estimation method to identify a shallow sea model of the entire European continental shelf. First, a space-varying bottom friction coefficient is estimated simultaneously with the depth. The second application is the estimation of the parameterization of the wind stress coefficient as a function of the wind velocity. Finally, an uncertain open boundary condition is included. It is shown that in this case the parameter estimation procedure does become more robust and produces more realistic estimates. Furthermore, an estimate of the open boundary conditions is also obtained.     

关于迭代Tikhonov正则化的最优正则参数选取

本文讨论了算子和右端都近似给定的第一类算子方程的迭代Ｔｉｋｈｏｎｏｖ正则化，给出了不依赖于准确解的任何信息但能得到最优收敛阶的正则参数选取法。     

金、铜、铂、银表面光学二次谐波参数的研究

以1.06μm激光为基频光,测量了金、铜、铂、银膜的反射二次谐波随基频光偏振角的变化,在b=-1的条件下拟合求得了这几种样品面谐波电流的唯象参数“a”值。