Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Adaptive Algorithm for Constrained Least-Squares Problems

This paper is concerned with the implementation and testing of an algorithm for solving constrained least-squares problems. The algorithm is an adaptation to the least-squares case of sequential quadratic programming (SQP) trust-region methods for solving general constrained optimization problems. At each iteration, our local quadratic subproblem includes the use of the Gauss–Newton approximation but also encompasses a structured secant approximation along with tests of when to use this approximation. This method has been tested on a selection of standard problems. The results indicate that, for least-squares problems, the approach taken here is a viable alternative to standard general optimization methods such as the Byrd–Omojokun trust-region method and the Powell damped BFGS line search method.     

拆边法求最短路径

本文利用局部比较法 ,在图中定义子图、无效路径、以及可去边 .利用推导的有关定理 ,拆去可去边 ,利用最短路径相同的等价性 ,达到化简图 ,从而求出最短路径     

最速下降规则不失为一可用规则

本从统计意义上平均迭代次数的观点出发，指出求解线性规划的最速下降规则是可取的，好用的。     

关于AOR迭代法的研究

本文论证了严格对角占优矩阵之AOR法的误差估计式中的误差估计常数hγ,ω（0≤γ≤ω0）的最小值是h1,1.     

线性分式规划最优解集的求法

本文使用多面集的表示定理，导出了线性分式规划最优解集的结构，并给出确定全部最优解的计算步骤。     

Truncated Newton method for sparse unconstrained optimization using automatic differentiation

When solving large complex optimization problems, the user is faced with three major problems. These are (i) the cost in human time in obtaining accurate expressions for the derivatives involved; (ii) the need to store second derivative information; and (iii), of lessening importance, the time taken to solve the problem on the computer. For many problems, a significant part of the latter can be attributed to solving Newton-like equations. In the algorithm described, the equations are solved using a conjugate direction method that only needs the Hessian at the current point when it is multiplied by a trial vector. In this paper, we present a method that finds this product using automatic differentiation while only requiring vector storage. The method takes advantage of any sparsity in the Hessian matrix and computes exact derivatives. It avoids the complexity of symbolic differentiation, the inaccuracy of numerical differentiation, the labor of finding analytic derivatives, and the need for matrix store. When far from a minimum, an accurate solution to the Newton equations is not justified, so an approximate solution is obtained by using a version of Dembo and Steihaug's truncated Newton algorithm (Ref. 1).This paper was presented at the SIAM National Meeting, Boston, Massachusetts, 1986.     

混沌实验数据处理及仿真

大学物理实验中混沌实验大多采用观察现象的方法进行,本实验采用蔡氏电路(Chua s circuit)产生混沌行为.在观察不同初始值条件下出现的倍周期分岔、阵发混沌、奇异吸引子等相图及现象的基础上,通过对采集数据进行处理,对负电阻伏安特性进行分段线性拟合,用功率频谱法、计算机仿真方法(龙格-库塔数值积分法)对混沌现象进行描绘,将实验数据与非线性方程组的数值解相结合,呈现出混沌现象的本质.     

脉冲电场屏蔽效能测试系统及测试方法

 自主研制了一套小型脉冲电场屏蔽效能测试系统，该系统由脉冲电场发射设备和脉冲电场测试设备构成。脉冲电场发射设备的天线口面前60 cm处可产生峰值达7.5 kV/m的脉冲电场，测试设备的动态测试范围达97 dB。系统采用光纤测量设备，在测量工作中能有效抑制强电磁场干扰。采用Matlab编写自动测量及数据处理程序，实现了数据采集与处理自动化。实验测量了金属桥架、控制柜、屏蔽帐篷、导电水泥混凝土房的脉冲电场屏蔽效能，其脉冲电场峰值衰减量分别为52,64,66,30 dB。实验表明可用脉冲电场的峰值衰减量来评估屏蔽体的脉冲电场屏蔽效能。     

The influence of pulse duration on the stress levels in ablation of ceramics: A finite element study

We present a finite element model to investigate the dynamic thermal and mechanical response of ceramic materials to pulsed infrared radiation. The model was applied to the specific problem of determining the influence of the pulse duration on the stress levels reached in human dental enamel irradiated by a CO₂ laser at 10.6 μm with pulse durations between 0.1 and 100 μs and sub-ablative fluence. Our results indicate that short pulses with durations much larger than the characteristic acoustic relaxation time of the material can still cause high stress transients at the irradiated site, and indicate that pulse durations of the order of 10 μs may be more adequate both for enamel surface modification and for ablation than pulse durations up to 1 μs. The model presented here can easily be modified to investigate the dynamic response of ceramic materials to mid-infrared radiation and help determine optimal pulse durations for specific procedures.