求解非线性方程组Newton—PCG并行算法

基于向量机（Ｖｅｃｔｏｒ Ｃｏｍｐｕｔｅｒ），应用二次迭代法和不完全的Ｃｈｏｌｅｓｋｉ不塞入分解原理，建立了求解非线性方程组的Ｎｅｗｔｏｎ－ＰＣＧ并行算法。     

强非线性振动系统的渐近解法

本文推广文[1]中方法于求解广泛一类强非线性振动系统,导出了适于近似定性分析和定量计算的简便公式。作为例子研究了修正的Vander Pol振子,最后给出数字Poincaré映射结果,从近似定性分析和定量计算两方面证明了本文求得的方法的有效性。     

求解非线性方程组Newton－PCG并行算法

基于向量机（ＶｅｃｔｏｒＣｏｍｐｕｔｅｒ），应用一次迭代法和不完全的Ｃｈｏｌｅｓｋｉ不塞入分解原理，建立了求解非线性方程组的Ｎｅｗｔｏｎ－ＰＣＧ并行算法。     

弹塑性接触问题数值计算的二次规划算法

针对弹塑性接触问题所推得的数值求解式子,运用二次规划法具体设计了算法,该算法采用有限的基底交换运算就可得到收敛的数值解,具有较好的收敛性及较小的计算工作量.工程计算算例结果表明文章所提出的接触问题的求解方法是有效的,针对结构中接触问题所建立的数值计算模型能真实反映实际工作状况是可靠的.文章中还详细给出了接触单元相关矩阵和向量的具体推导形式.     

结构优化的内点序列二次规划方法

本文给出一种对偶内点二次规划算法，通过利用代理约束的思想，构造了约束为单纯形的对偶问题，并通过Ｋａｒｍａｒｋａｒ的投影悄度变换来求解，同时以无约束极小点来判断所参与运算的约束，提高了算法的效率。目的在于提高序列二次规划的求解效率，并使之用于结构优化问题。     

二次分叉数值分析方法及在多孔介质热对流上的应用

本文给出了数值计算具有Ｚ２对称性和双参数的非线性动力学方程二次分叉问题的具体方法,并用此方法计算了不同长宽比矩形区域内多孔介质热对流的二次分叉点和相应的流场和温度场分布。     

求解非线性方程组的混沌分形方法

混沌分形是动力系统普遍出现的一种现象,牛顿-拉夫森NR(Newton-Raphson)方法是重要的一维及多维迭代技术,其迭代本身对初始点非常敏感,该敏感区是牛顿-拉夫森法所构成的非线性离散动力系统Julia集,在Julia集中迭代函数会呈现出混沌分形现象,提出了一种寻找牛顿-拉夫森函数的Julia点的求解方法,利用非线性离散动力系统在其Julia集出现混沌分形现象的特点,提出了一种基于牛顿-拉夫森法的非线性方程组求解的新方法,计算实例表明了该方法的有效性和正确性.     

非线性方程组的解法:局部弧长法

描述了一个新的非线性方程组的求解方法——局部弧长法．该方法是在弧长法的基础上发展起来的适合于材料非线性有限元分析的数值解法．其约束方程充分利用了结构中破坏区域内的非线性变形信息，有效地解决了材料非线性分析中的稳定性与收敛性问题．数值计算表明，该方法不仅适合于求解结构的极限承载能力，也适合于求解结构达到极限承载参力以后的荷载－变形的全过程     

弹塑性问题参数二次规划分析的隐式算法

传统的二次规划算法求解弹塑性问题时一般要经过对问题的线性化,如对屈服条件的一阶近似展开等,这在一定程度上会造成数值解的误差。为此,本文提出一种改进的策略,引入迭代与规划算法相结合的技术对问题进行处理,算法收敛平稳迅速,在大步长荷载增量下使算法的精度大大提高。由于本文的算法属于隐式算法,因而也就弥补了原二次规划算法求解弹塑性问题时只有显式算法的不足,从而达到了对原算法的进一步完善。     

Parametric quadratic programming method for viscoplasticity

Perzyna model in viscoplasticity has been studied by the parametric variationalprinciple,which could be transformed into solving the parametric quadratic programmingproblem.The FEM form of this problem and its implementation have also been discussed inthe paper.     

Mixed energy method for solution of quadratic programming problems and elastic-plastic analysis of truss structures

A new algorithm for the solution of quadratic programming problems is put forward in terms of the mixed energy theory and is further used for the incremental solution of elastic-plastic truss structures. The method proposed is different from the traditional one, for which the unknown variables are selected just in one class such as displacements or stresses. The present method selects the variables in the mixed form with both displacement and stress. As the method is established in the hybrid space, the information found in the previous incremental step can be used for the solution of the present step, making the algorithm highly efficient in the numerical solution process of quadratic programming problems. The results obained in the examples of the elastic-plastic solution of the truss structures verify what has been predicted in the theoretical analysis. Project supported by the National Natural Science Foundation of China (No. 50178916, No. 19732020 and No. 19872016), the National Key Basic Research Special Foundation (No. G1999032805), the Special Funds for Major State Basic Research Projects and the Foundation for University Key Teachers by the Ministry of Education of China.     

斜拉桥索力优化的强次可行序列二次规划法

提出了一种斜拉桥索力优化的实用方法-强次可行序列二次规划法. 该方法通过建立斜拉 桥索力优化的非线性规划模型，以主梁和索塔的弯曲应变能为目标函数，斜拉索的索力为设 计变量，结构应力和索力为约束条件，计入大跨度斜拉桥各种几何非线性因素的影响，采用 强次可行序列二次规划算法进行优化求解，确定斜拉桥成桥合理状态的索力. 运用该方法对 某斜拉桥进行索力优化，结果表明该方法简单、有效.     

Stochastic level-value approximation for quadratic integer convex programming

We propose a stochastic level value approximation method for a quadratic integer convex minimizing problem in this paper. This method applies an importance sampling technique, and make use of the cross-entropy method to update the sample density functions. We also prove the asymptotic convergence of this algorithm, and report some numerical results to illuminate its effectiveness.     

A combined parametric quadratic programming and precise integration method based dynamic analysis of elastic-plastic hardening/softening problems

The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems. The project supported by the National Key Basic Research Special Foundation (G1999032805), the National Natural Science Foundation of China (19872016, 50178016, 19832010) and the Foundation for University Key Teacher by the Ministry of Education of China     

Numerical simulation of contact problems in vane machinery by a parametric quadratic programming method

Interference fits are widely used for connecting impeller and shaft assembly that are forced together slowly by pressing. The interference fit design ensures stable balance behavior and allows for positive contact between the impeller and shaft assembly throughout the range of operating speeds. In addition to maintaining radial contact, sufficient net radial interface pressure must remain in order to transmit torque when the rotational speed is very high. Therefore, the interference fit between the impeller and the shaft assembly is one of the most important factors influencing the performance of the turbo unit in the design of turbocharger compressor. A suitable fit tolerance needs to be considered in the structural design. A locomotive-type turbocharger compressor with 24 blades under combined centrifugal and interference fit loading is used for the analysis. The finite-element (FE) parametric quadratic programming (PQP) method developed based on the parametric variational principle (PVP) is used for the analysis of the stress distribution in the three-dimensional (3D) contact problem of impeller. The advantages of the parametric programming method compared with conventional approaches are that the penalty factors can be canceled and that solutions can be obtained directly without tedious iterative procedures such as the general incremental iterative method. To save time in the computation, a~multi-substructure technique is adopted for structural modeling. This not only simplifies the calculation, but also provides a convenient service for process computer-aided design (CAD) by means of FE simulation. The effects of the fit tolerance, coefficient of friction and rotational speed (centrifugal force), wall thickness of the shaft sleeve and the contact stress on the interference-fitting surfaces are studied in detail in the numerical computation. It is found that a nonuniform initial amount of interference in the structural design avoids the relative displacement generated and ensures uniformity of the contact stress. To assure quality of press-fitting, the amount of interference between the shaft sleeve and shaft should be strictly controlled to avoid the rapid increase of the contact stress. The numerical results demonstrate the high accuracy and good convergence of the algorithm presented here, which provides an effective approach that achieves more-reliable interference-fitted connections and more-precise assembly accuracy with lower manufacturing cost in the structural design.     

NEWTON METHOD FOR SOLVING A CLASS OF SMOOTH CONVEX PROGRAMMING

1ConvexProgrammingandExactPenaltyFunction Weconsiderthefollowingconvexprogramming:(P)minf(x)s.t.x∈S={x∈Rn:gi(x)≤0,i=1,…,m}.SupposethatSisacompactset.ThusthereexistsalargeboundedboxX,suchthatS={x∈Rn:gi(x)≤0,i=1,…,m}intX.Assumethatf(x),gi(x),i=1,…,m     

A method of parametric quadrature programming for the problem of osmosis-solidification in elastoplastic rock-soil medium

The process of deformation in a medium of the kind of rock-soil is a dynamic process, because there is osmosis-solidification resulting in elastic-plastic deformation. In this paper, the parametric variational method and corresponding FEM are established for dealing with this problem. The problem is reduced to solving a quadratic programming problem with restrictive conditions (constitutive state equations). The choice of element form and the specific process of implementation are discussed. Two examples are given.     

A mixture differential quadrature method for solving two-dimensional incompressible Navier-Stokes equations

Ｔｈｅｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ（ＤＱＭ）ｐｒｏｐｏｓｅｄｂｙＲ．Ｂｅｌｌｍａｎ［１,２］ｈａｓｂｅｅｎｓｕｃｃｅｓｓｆｕｌｌｙｅｍｐｌｏｙｅｄｉｎｎｕｍｅｒｉｃａｌｃｏｍｐｕｔａｔｉｏｎｓｏｆｐｒｏｂｌｅｍｓｉｎｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｈｙｓｉｃａｌｓｃｉｅｎｃｅ．ＢｅｃａｕｓｅｔｈｅｉｎｆｏｒｍａｔｉｏｎｏｎａｌｌｇｒｉｄｐｏｉｎｔｓｉｓｕｓｅｄｔｏｆｉｔｔｈｅｄｅｒｉｖａｔｉｖｅｓａｔｇｒｉｄｐｏｉｎｔｓｉｎｔｈｅＤＱＭ,ｉｔｉｓｅｎｏｕｇｈｔｏｏｂｔａ…