关于“用迭代法解位移法和力法典型方程时的收敛性”的讨论

李作骏同志通过线代方法证明了位移法之反力矩阵(刚度矩阵)及力法之位移矩阵(柔度矩阵)均为对称正定矩阵。因而断定,用 Seidel 迭代法总是收敛的。其实通过结构力学方法,亦可予以证明。И.М.拉宾诺维奇在文献中,已试图证明力法之位移矩阵的各级主子式大于零。但由于他没有给出与各级子矩阵相应之结构关系,因此在逻辑上不够明确。现本文建立了与反力矩阵及位移矩阵之各级子矩阵相应的结构关系,并分别证明它们的各级主子式大于零,从而判定它们均为对称正定矩阵。     

有阻尼振动系统的自由界面子结构综合法

研究了阻尼振动系统的自由界面子结构综合法.为了在最终给出的系统综合方程中考虑剩余柔度矩阵的影响,本文利用质量矩阵、阻尼矩阵和结构保留模态矩阵构造了一种与保留模态矩阵加权正交的矩阵.在计算自由-自由结构的剩余柔度矩阵时,通过这种加权正交矩阵能够避免直接对子结构刚度矩阵的求逆运算.在此基础上提出了一种新的有阻尼振动系统的自由界面子结构综合法.算例分析表明,本文的方法具有很高的计算精度.     

随机减量特征矩阵和相关函数矩阵的关系

本文推导出多路信号的随机减量特征矩阵和相关函数矩阵之间的一般关系。对于承受平稳随机过程激励的线性时不变系统这种特殊情况,输出的随机减量特征矩阵正比于输出的协方差函数矩阵。作为例子,计算了运载火箭加速度响应时间历程的随机减量特征矩阵与相关函数矩阵。     

多自由度系统动力响应灵敏度分析与振动控制

采用脉冲响应矩阵将多自由度动力系统响应及其设计变量的灵敏度表示成积分形式,它适用于质量矩阵,阻尼矩阵和刚度和矩阵非对称及非比例阻尼情况。     

组合结构总体节点参数自动标定的基本定理

在有限元位移法中,单元矩阵(刚度矩阵、几何矩阵、一致质量矩阵、载荷列阵等)往往是在单元局部坐标下计算的(以下简称局部方位);结构整体矩阵的形成,单元应力和变形的计算,都要求确定总体节点参数(自由度)和单 ...     

斜交连续梁桥数值分析研究

在分析了矩阵位移法和传递矩阵法的理论基础上，将二者组合而成矩阵混合法。矩阵混合法的实质就是将结构物的一端未知向量通过矩阵的简单乘法传递到另一端，再将已知条件代入求解的过程。在矩阵传递的过程中引入位移法中的前进代入的概念，即将未知向量不断地以新的未知向量替换，可以防止误差的传递和积累。将矩阵混合法运用于斜交连续梁桥内力计算，完成了矩阵混合法在斜交连续梁桥内力计算中的理论推导和程序的设计，给出了实际算例并与有限元法进行了比较。     

埃尔米特梁单元的分块对角与高阶质量矩阵

埃尔米特梁单元常用的集中质量矩阵,是由挠度自由度对应的一致质量矩阵元素通过行求和或节点积分构造。然而,数值结果表明该集中质量矩阵在求解包含自由端的梁振动问题时,会出现频率精度掉阶现象。本文首先从保障质量矩阵最优收敛性的数值积分精度出发,分别针对三次和五次梁单元,发展了质量矩阵的梯度增强节点积分方案。利用梯度增强节点积分方案,可以得到具有分块对角形式的单元质量矩阵,而其组装的整体质量矩阵除边界节点外仍然呈现对角形式。对于两种单元,其分块对角质量矩阵分别具有4阶最优精度和6阶次优精度。再者,将标准一致质量矩阵和具有同阶精度的梯度增强节点积分质量矩阵进行优化组合,建立了具有超收敛特性的高阶质量矩阵。最后,通过数值算例系统验证了三次和五次单元的分块对角与高阶质量矩阵的频率计算精度。     

一个统一单元矩阵计算的方法

1.引言有限元应用软件系统中,普遍采用有限元位移法。但单元矩阵(刚度矩阵、质量矩阵、阻尼矩阵、以至载荷列阵)等的计算,多数是逐一给出,部分是以元素簇给出的。本文将讨论不同类型元素的单元矩阵的统一算法和程序实施技巧,给出计算精密单元和几何非线性单元方便、有效的途径,达到以积累插值函数积累单元矩阵库的目的。     

滚筒式洗衣机减隔振系统刚度阵的计算

本文针对滚筒式洗衣机弹簧减隔振系统的刚度阵的计算，提出了一种方便而概念清晰的方法。该法根据弹簧的主刚度方向，引入一个方向余弦矩阵和一个座标变换矩阵，通过两次矩阵变换，即可得到所求矩阵。     

求解对称矩阵特征值及特征向量的初等变换法

 在力学中有一类量的求解可归结为矩阵特征值和特征向量的求解，而求解 矩阵的特征值将要求解高次方程的根，这在数学上将遇到难以克服的困难. 应用初等变 换方法，将对称矩阵的特征矩阵对角化. 利用这种方法，同时可求出该矩阵所有的特征值和 正交的特征向量，避免了求解高次方程根的困难与把各特征向量正交化的麻烦.     

Analysis of stress concentrations in unidirectional composites taking into account the tensile load in the matrix

A modified shear-lag model accounting for the effect of the tensile stiffness of the matrix is proposed for solving the stress redistribution due to the failure of fibers and matrix in unidirectionally fibre-reinforced composites. The advantages of this model are simple, reasonable and accurate by comparison with the other similar modified shear-lag models. It can be further extended to study the stress redistribution with interfacial damage between fibres and matrix. This paper quantitatively discusses the influence of the tensile stiffness ratio of matrix to fibre and of the fibre volume fraction on the stress concentration in the fibres and matrix adjacent to cut fibres and matrix, and suggests that the influence of the matrix stiffness on the stress concentration can be neglected when the matrix stiffness is low, such as polymer matrix composites, and the fibre volume fraction is high. For other cases such as ceramic and metal matrix composites, the tensile load of the matrix cannot be neglected in the shearlag analysis. The project supported by the Guangdong Provincial Natural Science Foundation of China.     

Oil Displacement for One-Dimensional Three-Phase Flow with Phase Change in Fractured Media

In this article, the numerical simulations for one-dimensional three-phase flows in fractured porous media are implemented. The simulation results show that oil displacement in matrix is dominated by oil–water capillary pressure only under certain conditions. When conditions are changed to decrease the amount of water entering into the fractured media from the boundary of the flow field, water in fracture may be vaporized to superheated steam. In these cases, the appearance of superheated steam in fracture rather than in matrix will decrease the fracture pressure and generate the pressure difference between matrix and fracture, which results in oil flowing from matrix to fracture. Assuming that oil is wetting to steam, the matrix steam–oil capillary pressure will decrease the matrix oil-phase pressure as the matrix steam saturation increases. After the steam–oil capillary pressure finally exceeds the pressure difference due to the appearance of superheated steam in fracture, the oil displacement in matrix will stop. It is also shown that variations of the water relative permeability curve in matrix do not result in different mechanisms for oil displacement in matrix. The simulation results suggest that the amount of liquid water supply from the boundary of flow field fundamentally influence the mechanisms for oil displacement in matrix.     

多体系统传递矩阵法研究进展

作为一种多体系统动力学新方法, 多体系统传递矩阵法由于其无需系统总体动力学方程和快速计算的特点, 已被广泛用于各种多管火箭、自行火炮、舰炮等复杂大型机械系统动力学分析与设计. 本文介绍了该方法的研究进展, 包括: 线性多体系统传递矩阵法、多体系统离散时间传递矩阵法、二维系统传递矩阵法、受控多体系统传递矩阵法、多体系统传递矩阵法和通常动力学方法的混合方法等, 给出了该方法解决自行火炮、多管火箭武器多体系统动力学的重大工程应用实例.      

Evidence of Multi-Process Matrix Diffusion in a Single Fracture from a Field Tracer Test

Compared to values inferred from laboratory tests on matrix cores, many field tracer tests in fractured rock have shown enhanced matrix diffusion coefficient values (obtained using a single-process matrix-diffusion model with a homogeneous matrix diffusion coefficient). To investigate this phenomenon, a conceptual model of multi-process matrix diffusion in a single-fracture system was developed. In this model, three matrix diffusion processes of different diffusion rates were assumed to coexist: (1) diffusion into stagnant water and infilling materials within fractures, (2) diffusion into a degraded matrix zone, and (3) further diffusion into an intact matrix zone. The validity of the conceptual model was then demonstrated by analyzing a unique tracer test conducted using a long-time constant-concentration injection. The tracer-test analysis was conducted using a numerical model capable of tracking the multiple matrix-diffusion processes. The analysis showed that in the degraded zone, a diffusion process with an enhanced diffusion rate controlled the steep rising limb and decay-like falling limb in the observed breakthrough curve, whereas in the intact matrix zone, a process involving a lower diffusion rate affected the long-term middle platform of slowly increasing tracer concentration. The different matrix-diffusion-coefficient values revealed from the field tracer test are consistent with the variability of matrix diffusion coefficient measured for rock cores with different degrees of fracture coating at the same site. By comparing to the matrix diffusion coefficient calibrated using single-process matrix diffusion, we demonstrated that this multi-process matrix diffusion may contribute to the enhanced matrix-diffusion-coefficient values for single-fracture systems at the field scale.     

Efficient approximation of Sparse Jacobians for time‐implicit reduced order models

This paper introduces a sparse matrix discrete interpolation method to effectively compute matrix approximations in the reduced order modeling framework. The sparse algorithm developed herein relies on the discrete empirical interpolation method and uses only samples of the nonzero entries of the matrix series. The proposed approach can approximate very large matrices, unlike the current matrix discrete empirical interpolation method, which is limited by its large computational memory requirements. The empirical interpolation indices obtained by the sparse algorithm slightly differ from the ones computed by the matrix discrete empirical interpolation method as a consequence of the singular vectors round‐off errors introduced by the economy or full singular value decomposition (SVD) algorithms when applied to the full matrix snapshots. When appropriately padded with zeros, the economy SVD factorization of the nonzero elements of the snapshots matrix is a valid economy SVD for the full snapshots matrix. Numerical experiments are performed with the 1D Burgers and 2D shallow water equations test problems where the quadratic reduced nonlinearities are computed via tensorial calculus. The sparse matrix approximation strategy is compared against five existing methods for computing reduced Jacobians: (i) matrix discrete empirical interpolation method, (ii) discrete empirical interpolation method, (iii) tensorial calculus, (iv) full Jacobian projection onto the reduced basis subspace, and (v) directional derivatives of the model along the reduced basis functions. The sparse matrix method outperforms all other algorithms. The use of traditional matrix discrete empirical interpolation method is not possible for very large dimensions because of its excessive memory requirements. Copyright © 2016 John Wiley & Sons, Ltd.     

大型辛矩阵本征问题的逆迭代法

基于共轭辛子空间迭代法,求解了大型辛矩阵的主要本征解。随着迭代的进行,可以无限地逼近其精确解。     

Study on application of random walk method to calculate water exchange in large-scale bay

The water exchange matrix is an efficient tool to study the water exchange among the sub-areas in large-scale bays. The application of the random walk method to calculate the water exchange matrix is studied. Compared with the advection-diffusion model, the random walk model is more flexible to calculate the water exchange matrix. The forecast matrix suggested by Thompson et al. is used to evaluate the water exchange characteristics among the sub-areas fast. According to the theoretic analysis, it is found that the precision of the predicted results is mainly affected by three factors, namely, the particle number, the generated time of the forecast matrix, and the number of the sub-areas. The impact of the above factors is analyzed based on the results of a series of numerical tests. The results show that the precision of the forecast matrix increases with the increase of the generated time of the forecast matrix and the number of the particles. If there are enough particles in each sub-area, the precision of the forecast matrix will increase with the number of the sub-areas. Moreover, if the particles in each sub-area are not enough, the excessive number of the sub-areas can result in the decrease of the precision of the forecast matrix.     

Micromechanical prediction of the two-way shape memory effect in shape memory alloy composites

The two-way shape memory effect in monolithic shape memory alloys has been widely investigated both theoretically and experimentally. In the present study, this effect is analyzed for shape memory alloy composites by employing a micromechanical model. To this end, the responses of polymeric matrix and metal matrix unidirectional composites with embedded shape memory alloy fibers are determined. For the polymeric matrix composite, the effect of axial, transverse and shear loadings as well as the fiber volume fraction on the resulting two-way shape memory behavior are studied. The local distributions of stresses among the shape memory alloy fiber and epoxy matrix in the low- and high-temperature shapes of the composite are also investigated. Two training procedures that generate the two-way shape memory effect in the metal matrix composite are offered. The present analysis shows that the two-way shape memory effect in the chosen type of metal matrix composite is not as useful as in the polymeric matrix one. Finally, for a polymeric matrix composite that is subjected to a transverse normal loading, the effect of imperfect bonding between the shape memory alloy fibers and the neighboring matrix is investigated.     

Some new conditions for generalized H-matrices

By using a continuous transition method of a matrix and the estimate for spectral radius of a sub-matrix etc.,decision methods for a generalized H-matrix under positive definite matrix conditions are researched.Some new sufficient conditions for generalized H-matrices are obtained.When a block matrix degenerates a point matrix, these conditions namely become sufficient conditions of H-matrix.     

Solutions of matrix equations in problems of mechanics and control

The paper considers algorithms for solving linear matrix equations related to problems of mechanics and control, namely, the Lyapunov and Sylvester matrix equations and Riccati-type nonlinear matrix equations. These algorithms are capable of solving both linear equations and linear matrix inequalities. Algorithms based on the Bass relations are used to solve Riccati-type nonlinear matrix equations in so-called special cases where some eigenvalues of the matrix pencil are on a unit circle. These algorithms are compared with those of other authors by way of examples. It is shown that the algorithms can be implemented in symbolic computing routings, which allows solving these equations with high accuracy