广义特征值问题的并行块Jacobi-Davidson方法及应用

给出了对称矩阵广义特征值问题AX=λBX的并行块Jacobi-Davidson方法.该方法使用投影技术将大型矩阵特征值问题转变成低维子空间中矩阵特征值问题,并利用Neumann级数展开对校正方程进行预处理.该方法可同时并行计算广义特征值问题的几个极端特征对,具有良好的并行性.将这一方法应用于某型号机翼及挂架的结构动力分析并行计算,在IBM-P650并行计算机上的数值试验结果表明,在相同迭代精确度的条件下,Jacobi-Davidson方法比子空间迭代法使用较少的迭代次数和运算时间,并具有更高的加速比和并行效率.     

不确定结构动力特征值区间分析的一种算法

将结构系统中的不确定性参数用区间数来表示,对获得的广义区间特征值方程的求解方法进行了讨论,提出了一种区间逐步离散的方法。此方法通过独立的不确定性参数取区间离散点的值,将广义区间特征值方程的求解转化为相应的确定性问题,再搜索方程解中的最大最小值来确定各阶特征值边界。用数学算例对此算法的正确性和有效性进行了验证,然后应用于工程算例的特征值区间分析,并与其它算法结果进行了比较。计算结果表明该算法的计算效率较高,准确性较好。     

广义特征值问题的Collatz包含定理

工程中已发展了许多矩阵特征值问题的近似求解方法,由Duncley法给出固有频率基频的下界,Rayleigh-Ritz近似法建立的方程,给出基频的上界,以及通常的矩阵迭代法给出的矩阵的固有频率程序中是以某一元素迭代前后比值确定的,这样实际上很难说是上界或下界。Collatz包含定理仅适用于对称标准特征值问题,可以给出特征值上、下界。采用矩阵Cholesky三角分解的原理,把Collatz包含定理推广到适用于具有对称矩阵的一般结构系统的广义征值问题,对于分解刚度矩阵或质量矩阵可给出基频,或最高因有频率。为了验证理论的正确性,给出了算例。     

求解特征值问题的广义最小二乘法

本文研究线性特征值问题Aφ-λBφ=0,其中A、B为线性对称运算子且B正定,提出了能够求得特征近似值,同时求得其上下限的广义最小二乘法,证明了特征值上下限定理,导出了特征值的第一上限、第一下限、第二上限、第二下限。本文方法可以求解微分方程的特征值及微分方程组的特征值、矩阵特征值等问题,后两类问题的两个算例结果精度很高。     

线性广义特征值问题的摄动迭代法

本文给出了一种适用于迭代计算的矩阵摄动法,它是进行广义特征值问题Ax=λBx的摄动重分析的一种高精度算法,同时也可用于改进由其它矩阵摄动分析方法提供的近似解的精度。实际算例表明,当结构参数修改量不太大时,采用这种摄动迭代法进行特征值问题的精确重分析是十分有效的。     

大型陀螺特征值问题的广义Arnoldi减缩算法

基于Ａｒｎｏｌｄｉ法，建立陀螺特征值问题的广义Ａｒｎｏｌｄｉ格式，并利用系统矩阵的反对称特性，得到极其简洁的甚至比对称矩阵Ｌａｎｃｚｏｓ法更为简单的递推格式，可称为陀螺Ａｒｎｏｌｄｉ减缩算法。     

空间刚架结构的广义逆力法

在结构计算中，根据算法中所采用的基本未知量的不同．结构分析方法可以分为力法、位移法和混合法。其中位移法由于适宜计算机处理而在结构计算领域得到了广泛的应用．经典力法相比之下应用就远不如位移法普遍，虽然力法本身在力学上有其独特的优势。广义逆矩阵做为一种较新的数学工具，自二十世纪五十年代诞生以来正日益表现出越来越旺盛的生命力。广义逆力法就是一种基于力法和广义逆矩阵理论的新的迭代算法。这种算法是一种完全适合计算机处理的力法方法。该算法的思路以及对于求解线弹性空间刚架结构问题的具体公式均在文中给出并给出了算例。从算例计算结果可以看到广义逆力法有着较好的计算效率和计算精度。该算法的提出为力法在计算机计算领域的应用开拓了新的发展空间。该算法在材料非线性问题和结构并行计算方面也有着较好的应用前景。     

分数阶广义热弹性理论下中空柱热弹性分析

基于Laplace变换及特征值法,推导并给出了分数阶广义热弹性理论下中空柱内表面作用有热冲击情况的解析解,通过Laplace数值逆变换法求解得到了位移场、温度场、应力场的分布规律。结果表明:特征值法能准确给出Laplace域内方程组的解;分数阶参数对温度场和应力场有较大影响,对位移场影响较小。作为广义热弹性理论的一种推广,在处理热传导问题时,通过分数阶广义热弹性理论进行研究更科学、全面。     

计算具有区间参数结构特征值范围的一种新方法

基于区间效学的包含单调性和区间函效所表述的实际物理意义，把广义区间特征值问题转化为两个以非确定参效为优化变量，以关心的特征值为目标函效的全局优化问题，并采用遗传算法对优化问题求解，计算得到结构特征值的区间范围。通过效值算例对本文方法的有效性进行了验证，并和区间摄动法的计算结果进行了比较。     

0rr－Sommerfeld方程数值解法中的复广义矩阵特征值问题

Ｏｒｒ－Ｓｏｍｍｅｒｆｅｌｄ方程的求解通常可以化为一个复广义矩阵特征值问题ＡＸ＝ωＢＸ。本文用酉变换分别约化Ａ，Ｂ为上Ｈｅｓｓｅｎｂｅｒｇ阵和上三角阵，然后利用Ｍｕｌｌｅｒ求根方法可以求出其全部特征值，其中特征多项式的值由Ｈｙｍａｎ方法给出。当仅需要判断有无不稳定模态时，利用一个简单的矩阵变换将其化为强特征值的求解问题，从而可使用最简单的幂迭代，Ｃｈｅｂｙｓｈｅｖ配置点法算例表明两种算法均快速有效。     

Existence and algorithm of solutions for a system of generalized mixed implicit equilibrium problems in Banach spaces

A new system of generalized mixed implicit equilibrium problems is introduced and studied in Banach spaces. First, the notion of the Yosida proximal mapping for generalized mixed implicit equilibrium problems is introduced. By using the notion, a system of generalized equation problems is considered, and its equivalence with the system of generalized mixed implicit equilibrium problems is also proved. Next, by applying the system of generalized equation problems, we suggest and analyze an iterative algorithm to compute the approximate solutions of the system of generalized mixed implicit equilibrium problems. The strong convergence of the iterative sequences generated by the algorithm is proved under quite mild conditions. The results are new and unify and generalize some recent results in this field.     

哈密顿系统正则变换在时变最优控制中的应用

利用哈密顿系统正则变换和生成函数理论求解线性时变最优控制问题,构造了新的最优控制律形式并提出了控制增益计算的保结构算法. 利用生成函数求解最优控制导出的哈密顿系统两端边值问题,并构造线性时变系统的最优控制律,由第2类生成函数所构造的最优控制律避免了末端时刻出现无穷大反馈增益. 控制系统设计中需求解生成函数满足的时变矩阵微分方程组. 根据生成函数与哈密顿系统状态转移矩阵之间的关系,从正则变换的辛矩阵描述出发,导出了求解这组微分方程组的保结构递推算法.为了保持递推计算中的辛矩阵结构,哈密顿系统状态转移矩阵的计算中利用了Magnus级数.      

Generalized complementarity problems for fuzzy mappings

In this paper we introduce a new class of generalized complementarity problems for the fuzzy mappings and construct a new iterative algorithm. We also discuss the existence of solutions for the generalized complementarity problems and the convergence of tterative sequence.     

Auxiliary principle and approximation solvability for system of new generalized mixed equilibrium problems in reflexive Banach spaces

A new system of generalized mixed equilibrium problems (SGMEPs) involving generalized mixed variational-like inequality problems is introduced and studied in reflexive Banach spaces. A system of auxiliary generalized mixed equilibrium problems (SAGMEPs) for solving the SGMEPs is first introduced. Existence and uniqueness of the solutions to the SAGMEPs is proved under quite mild assumptions without any coercive conditions in reflexive Banach spaces. Using the auxiliary principle technique, a new iterative algorithm for solving the SGMEPs is proposed and analyzed. Strong convergence of the iterative sequences generated by the algorithm is also proved under quite mild assumptions without any coercive conditions. These results improve, unify, and generalize some recent results in this field.     

MODAL SYNTHESIS METHOD FOR NORM COMPUTATION OF H∞ DECENTRALIZED CONTROL SYSTEMS (Ⅱ)

When using H ∞ techniques to design decentralized controllers for large systems,the whole system is divided into subsystems, which are analysed using H∞ control theory before being recombined. An analogy was established with substructural analysis in structural mechanics, in which H ∞ decentralized control theory corresponds to substructural modal synthesis theory so that the optimal H ∞ norm of the whole system corresponds to the fundamental vibration frequency of the whole structure. Hence, modal synthesis methodology and the extended Wittrick-Williams algorithm were transplanted from structural mechanics to compute the optimal H∞ norm of the control system. The orthogonality and the expansion theorem of eigenfunctions of the subsystems H∞ control are presented in part (Ⅰ)of the paper. The modal synthesis method for computation of the optimal H ∞ norm of decentralized control systems and numerical examples are presented in part (Ⅱ).     

MODAL SYNTHESIS METHOD FOR NORM COMPUTATION OF H ∞ DECENTRALIZED CONTROL SYSTEMS ( Ⅰ )

When using H ∞ techniques to design decentralized controllers for large systems,the whole system is divided into subsystems, which are analysed using H ∞ control theory before being recombined. An analogy was established with substructural analysis in structural mechanics, in which H∞ decentralized control theory corresponds to substructural modal synthesis theory so that the optimal H∞ norm of the whole system corresponds to the fundamental vibration frequency of the whole structure. Hence, modal synthesis methodology and the extended Wittrick-Williams algorithm were transplanted from structural mechanics to compute the optimal H ∞ norm of the control system. The orthogonality and the expansion theorem of eigenfunctions of the subsystems H∞ control are presented in part ( Ⅰ ) of the paper. The modal synthesis method for computation of the optimal H∞ norm of decentralized control systems and numerical examples are presented in part ( Ⅱ ).     

PREDICTOR-CORRECTOR ALGORITHMS FOR SOLVING GENERALIZED MIXED IMPLICIT QUASI-EQUILIBRIUM PROBLEMS

A new class of generalized mixed implicit quasi-equilibrium problems （GMIQEP） with four-functions is introduced and studied. The new class of equilibrium problems includes many known generalized equilibrium problems and generalized mixed implicit quasi-variational inequality problems as many special cases. By employing the auxiliary principle technique, some predictor-corrector iterative algorithms for solving the GMIQEP are suggested and analyzed. The convergence of the suggested algorithm only requires the continuity and the partially relaxed implicit strong monotonicity of the mappings     

Multi-symplectic method for generalized Boussinesq equation

The generalized Boussinesq equation that represents a group of important nonlinear equations possesses many interesting properties. Multi-symplectic formulations of the generalized Boussinesq equation in the Hamilton space are introduced in this paper. And then an implicit multi-symplectic scheme equivalent to the multi-symplectic Box scheme is constructed to solve the partial differential equations （PDEs） derived from the generalized Boussinesq equation. Finally, the numerical experiments on the soliton solutions of the generalized Boussinesq equation are reported. The results show that the multi-symplectic method is an efficient algorithm with excellent long-time numerical behaviors for nonlinear partial differential equations.     

An operative algebraic formulation for the unilaterally-constrained mechanical problem of smart tensegrities

Kinematics and statics of tensegrities are addressed by means of a novel algebraic formulation. The inequality constraints, associated to cable-type unilateral structural members, are explicitly enforced in the equilibrium and compatibility problems. Fundamental tensegrity properties (rigidity, pre-stressability, and stability) are focused by a novel structural perspective and algebraic criteria for their assessment are established. Some classical results are generalized to the case of tensegrity models involving both deformable and non-deformable structural members. An operative algorithm for the analysis of the large-displacement elastic tensegrity response is proposed, not limited by special requirements in terms of structural symmetries or member connectivity, and therefore resulting a general design tool. Exemplary applications highlight the effectiveness of the proposed approach for designing tensegrity structures endowed with smart global behavior related to the optimal tuning of structural stiffness.