基于神经网络的大型空间网架结构的有限元模型修正

本文提出了依据神经网络方法完成大型空间网架结构节点固结系数的识别，实现空间网架结构有限元模型修正的方法。文中，通过有效独立法优化配置有限数量传感器的安放位置。在此基础上，利用从优化位置处的结构时程响应中获取的结构自振频率和局部位移模态信息，本文建立了基于径向基神经网络的空间网架结构有限元模型修正的计算方法。数值仿真的计算结果表明，该方法能有效地用于空间网架结构的有限元模型修正。     

预条件共轭梯度法在地震数据重建方法中的应用

基于最小平方的Fourier地震数据重建方法最终转化为求解一个线性方程组, 其系数矩阵是Toeplitz矩阵,可以用共轭梯度法求解该线性方程组.共轭梯度法的迭代次数受系数矩阵病态程度的影响,地震数据的非规则采样程度越高,所形成的系数矩阵病态程度越高,就越难收敛和得到合理的计算结果.本文研究了基于Toeplitz矩阵的不同预条件的构造方法,以及对共轭梯度法收敛性的影响.通过预条件的使用,加快了共轭梯度法的迭代速度, 改进了共轭梯度算法的收敛性,提高了计算的效率.数值算例和实际地震数据重建试验证明了预条件共轭梯度法对计算效率有很大的提高.     

用改进的光滑约束最小二乘正交分解法实现电阻率三维反演

对三维电阻率反演问题进行了深入研究,提供了一种利用地表观测数据实现三维反演的实用算法.该方法应用有限差分求正演解,并通过对粗糙度矩阵元素进行适当改进,使之适用于各种情况下粗糙度矩阵的求取,进而建立在模型的总粗糙度极小条件下的反演方程.对反演方程采用收敛速度快且稳定的最小二乘正交分解（LSQR）法进行迭代求解,在迭代求解过程中只需利用偏导数矩阵和其转置矩阵乘以一个向量的结果,回避了直接求偏导数矩阵的繁琐计算,节省了内存,加快了反演的计算速度.不同的计算实例表明上述方法是求解大规模三维电阻率反演问题的有效方法.     

煤矿井下微震震源高精度定位研究

煤矿井下微震震源准确定位,对于动力灾害监测预警具有重要意义.由于微震震源需要通过井下传感器接收信息反演确定,传感器的安装位置限制于煤矿井下巷道周围,传感器沿巷道近平面的不合理布置将大大降低震源定位精度.针对由传感器信息反演震源位置引起的病态问题,本文提出了基于微震监测测点优化布置的震源高精度定位算法.首先通过计算系数矩阵条件数,判定病态问题;然后利用中心化法和行平衡法联合进行病态矩阵预处理.对预处理后的矩阵A、b利用L曲线法计算正则参数,结合Tikhonov正则化算法计算得到震源坐标正则解.研究结果表明,中心化法有效降低了矩阵数量级,行平衡预处理降低了病态条件数,预处理后Tikhonov正则解的震源坐标误差最小可以达到3.09m,与预处理前的高斯消去解相比误差大大降低.通过上述优化处理,实现了井下受限空间微震监测震源高精度定位.     

基于特征矩阵联合近似对角化和经验模态分解的储层识别

地震信号是非线性和非平稳的信号.本文提出了基于特征矩阵联合近似对角化和经验模态分解相结合的算法,将地震信号分解为多个相互独立的固有模态函数分量.然后利用已知井储层发育情况,选取对储层识别有效的IMF分量进行储层预测.仿真和实际地震数据应用表明了该方法的有效性.     

基于正交时频原子的地震信号快速匹配追踪

匹配追踪能够根据地震信号自身的特点进行自适应分解,但因计算量巨大使其不能被广泛应用.本文提出了一种正交时频原子快速匹配追踪方法,该方法在迭代分解过程中,充分利用地震信号的局部性特征作为先验信息点,在先验信息点附近采用动态搜索策略寻找最佳匹配时频原子,同时对时频原子进行正交变换,消除了时频原子库中的冗余分量,最终将信号分解为一系列正交时频原子的线性叠加.测试结果表明,该方法不仅保持了匹配追踪的分解精度,而且使计算效率有了质的提高.     

基于不等式约束的最小二乘法三维电阻率反演及其算法优化

基于光滑约束的最小二乘法是三维电阻率反演的主要方法,但该方法在某些情况下存在着多解性较强的问题,且普遍耗时较长,严重制约了三维反演方法的推广与发展.为改善上述问题,将表征模型参数变化范围的不等式约束作为先验信息引入最小二乘线性反演方法中,有效地改善了反演结果的精度,降低了反演的多解性问题.为了解决耗时较长的问题,基于预条件共轭梯度(PCG)算法和Cholesky分解法的特点提出了一套优化三维电阻率反演计算效率的计算方案.在该方案中,Cholesky分解法被用来求解敏感度矩阵计算中的多个点源场的正演问题,Cholesky分解法只需对总体系数矩阵进行一次分解,然后对不同的右端向量进行回代即可.将预条件共轭梯度法引入到三维电阻率反演方程的求解中,将雅可比迭代中的对角阵作为预处理矩阵,其具有求逆方便、无需内存空间的特点,有效地加快了收敛速度.对合成数据以及实测数据的反演算例表明,借助不等式约束和反演效率优化方案,最小二乘反演方法可得到较为精确的反演结果,有效地提高了反演计算效率,具有良好的推广前景.     

求解弹性波动方程的频率域近似解析离散化波场模拟方法

为提高频率域弹性波动方程数值求解的计算效率,本文引入近似解析离散化(NAD)方法将其进行数值离散并得到大型线性代数方程组.在详细分析了相应系数矩阵的稀疏分块结构与数学性质之后,本文提出采用不精确旋转分块三角预处理子加速Krylov子空间迭代方法来快速求解该线性方程组,并利用数值试验证实这种方法在弹性波场模拟方面的数值效率.通过与另外两种经典数值方法(常规有限差分方法和交错网格有限差分方法)对多种介质模型进行波场模拟、数值频散分析以及与解析解的波形对比,NAD方法显示了其在压制数值频散和提高计算效率方面的优势以及对复杂介质模型弹性波场数值模拟的有效性.     

钢筋混凝土框架结构非线性全过程分析的一种新方法

在极限承载力附近的计算发散问题一直是困扰混凝土结构非线性全过程分析的关键问题,现有的迭代方法中只有弧长法被认为是最为有效的迭代方法,但由于混凝土结构具有高度的材料非线性,弧长法的稳定性难以保证。为此本文提出了一种新的方法———虚拟单元法,该方法在原结构上添加虚拟单元,虚拟单元的刚度系数在每个加载步中均保持为常数,该法在保证结构比例加载的前提下,可以有效改善结构刚度矩阵的性能,为得到结构非线性全过程分析曲线,特别是下降段部分提供了可能。本文采用虚拟单元法,通过自编的程序进行了实例分析,计算中,选取一榀混凝土框架进行模拟计算,杆件采用纤维模型,分析结果与ANSYS计算结果进行了对比,表明了该法的有效性,并对建立在纤维模型基础上的结构极限荷载进行了一定的探讨。     

基于预处理广义极小残量法的大地电磁正演计算

有限单元法求解大地电磁正演问题,会形成大型稀疏的带状复系数矩阵方程组,其条件数大,易造成病态矩阵,求解极困难.为此,采用每步迭代都具有最优性的广义极小残量法(Generalized Minimal Residual)求解;但是,每步迭代的计算量和存储量都线性增长,造成迭代收敛速度慢,因此,引入不完全LU分解预处理技术以降低矩阵的条件数,加快广义极小残量法的收敛速度.通过二维模型和层状介质模型的电磁响应的计算表明,基于预处理广义最小残量法精确而稳定.     

Optimal transducer placement for health monitoring of long span bridge

In experimental modal testing, the measurement locations and the number of measurements have a major influence on the quality of the results. In general, there are several alternative schemes for sensor placement, and the accuracy of the data increases as the number of sensors utilized increases. However, the number of transducers that can be attached to a real structure is limited by economic constraints. Therefore, algorithms that address the issue of limited instrumentation and its effects on resolution and accuracy are important from the standpoint of experimental modal analysis. The authors are particularly interested in structural dynamics based damage evaluation of large structures, and the development and implementation of suitable sensor location algorithms are critical for such a problem. A kinetic energy optimization technique (EOT) has been derived, and numerical issues are addressed and applied to real experimental data obtained from a model of an asymmetric long span bridge. Using experimental data from the bridge model, the algorithm proposed in this paper is compared to Kammer's EIM algorithm, which optimizes the transducer placement for identification and control purposes.     

地震作用下多层剪切型组合结构的合理综合阻尼比取值探讨

对于组合结构,提出了一种新的基于单元瑞利阻尼模型的应变能振型阻尼比,并证明了其在特定整体阻尼矩阵下与强迫解耦法的等价性;推导了具有明确理论依据的综合阻尼比计算公式,并基于相应的应变能振型阻尼比得出了结构的刚度综合阻尼比和瑞利综合阻尼比。分别采用复振型分解法和振型分解法对算例结构进行了地震荷载作用下的弹性时程分析,结果表明,瑞利综合阻尼比对于以剪切变形为主的多层组合结构具有良好的计算精度和适用性。     

DPCM-based vibration sensor data compression and its effect on structural system identification

Due to the large scale and complexity of civil infrastructures, structural health monitoring typically requires a substantial number of sensors, which consequently generate huge volumes of sensor data. Innovative sensor data compression techniques are highly desired to facilitate efficient data storage and remote retrieval of sensor data. This paper presents a vibration sensor data compression algorithm based on the Differential Pulse Code Modulation (DPCM) method and the consideration of effects of signal distortion due to lossy data compression on structural system identification. The DPCM system concerned consists of two primary components: linear predictor and quantizer. For the DPCM system considered in this study, the Least Square method is used to derive the linear predictor coefficients and Jayant quantizer is used for scalar quantization. A 5-DOF model structure is used as the prototype structure in numerical study. Numerical simulation was carried out to study the performance of the proposed DPCM-based data compression algorithm as well as its effect on the accuracy of structural identification including modal parameters and second order structural parameters such as stiffness and damping coefficients. It is found that the DPCM-based sensor data compression method is capable of reducing the raw sensor data size to a significant extent while having a minor effect on the modal parameters as well as second order structural parameters identified from reconstructed sensor data.     

Structural physical parameter identification based on evolutionary-simplex algorithm and structural dynamic response

Evolutionary computation based on the idea of biologic evolution is one type of global optimization algorithm that uses self-adaptation, self-organization and random searching to solve optimization problems. The evolutionary-simplex algorithm is introduced in this paper. It contains floating encoding which combines the evolutionary computation and the simplex algorithm to ovcrcomc the problems encountered in the genetic algorithm and evolutionary strategy methods.Numerical cxpcrimcnts arc performed using seven typical functions to verify the algorithm. An inverse analysis method to identify structural physical parameters based on incomplete dynamic responses obtained from the analysis in the time domain is prcscntcd by using the evolutionary-simplex algorithm. The modal evolutionary-simplex algorithm converted from the time domain to the modal domain is proposed to improve the inverse efficiency. Numerical calculations for a 50-DOF system show that whcn compared with other methods, the evolutionary-simplex algorithm offers advantages of high precision,cfficient searching ability, strong ability to resist noise, independence of initial value, and good adaptation to incomplete information conditions.     

Estimation of rainfall intensity–duration–frequency curves at ungauged locations using quantile regression methods

Intensity–duration–frequency (IDF) curves of extreme rainfall are used extensively in infrastructure design and water resources management. In this study, a novel regional framework based on quantile regression (QR) is used to estimate rainfall IDF curves at ungauged locations. Unlike standard regional approaches, such as index-storm and at-site ordinary least-squares regression, which are dependent on parametric distributional assumptions, the non-parametric QR approach directly estimates rainfall quantiles as a function of physiographic characteristics. Linear and nonlinear methods are evaluated for both the regional delineation and IDF curve estimation steps. Specifically, delineation by canonical correlation analysis (CCA) and nonlinear CCA (NLCCA) is combined, in turn, with linear QR and nonlinear QR estimation in a regional modelling framework. An exhaustive comparative study is conducted between standard regional methods and the proposed QR framework at sites across Canada. Overall, the fully nonlinear QR framework, which uses NLCCA for delineation and nonlinear QR for estimation of IDF curves at ungauged sites, leads to the best results.     

大地电磁三维数据空间反演并行算法研究

目前大地电磁三维反演实际应用的主要问题是计算效率低.在对大地电磁三维数据空间反演算法进行深入分析的基础上,本文提出了基于频点和矩阵划分的大粒度并行反演方案和具体实现步骤,并在曙光TC5000A高性能计算平台上实现了基于MPI的大地电磁三维数据空间反演并行算法.该算法实现了包括三维正演、灵敏度矩阵、叉积矩阵以及模型改正量的并行执行,不仅计算效率高,而且每个节点机上灵敏度矩阵的存储空间只需原来微机上的2/N(N是参加并行计算的节点机个数),大大地减少了内存开销.通过两个理论模型合成的数据对实现的三维数据空间反演并行算法进行试算,对比分析了多个节点机下程序的执行效率.测试结果表明,所实现的三维数据空间反演并行算法是可行的、高效的,与单机相比,不仅可以提高运行速度,缩短计算时间,而且还可以扩大计算规模,极大地推动了大地电磁三维反演的实用化.     

Structural identification with incomplete output‐only data and independence of measured information for shear‐type systems

Structural identification is the inverse problem of estimating physical parameters of a structural system from its vibration response measurements. Incomplete instrumentation and ambient vibration testing generally result in incomplete and arbitrarily normalized measured modal information, often leading to an ill‐conditioned inverse problem and non‐unique identification results. The identifiability of any parameter set of interest depends on the amount of independent available information. In this paper, we consider the identifiability of the mass and stiffness parameters of shear‐type systems in output‐only situations with incomplete instrumentation. A mode shape expansion‐cum‐mass normalization approach is presented to obtain the complete mass normalized mode shape matrix, starting from the incomplete non‐normalized modes identified using any operational modal analysis technique. An analysis is presented to determine the minimum independent information carried by any given sensor set‐up. This is used to determine the minimum necessary number and location of sensors from the point of view of minimum necessary information for identification. The different theoretical discussions are illustrated using numerical simulations and shake table experiments. It is shown that the proposed identification algorithm is able to obtain reliably accurate physical parameter estimates under the constraints of minimal instrumentation, minimal a priori information, and unmeasured input. The sensor placement rules can be used in experiment design to determine the necessary number and location of sensors on the monitored system. John Wiley & Sons, Ltd.     

基于损伤部位向量法的钢框架损伤识别研究

对损伤部位向量(DLV)法作了简单介绍,并用该方法对钢框架进行了损伤识别和损伤定位。该方法假定结构损伤前后为线性,对结构损伤前后柔度矩阵差进行奇异值分解,将奇异值为零所对应的向量,作为静荷载施加在无损结构的测点位置,则应力为零的单元为可能损伤的单元。对3种不同工况的钢框架进行了振动模态试验,用前3阶模态参数构造框架的柔度矩阵,按照DLV法对其进行了损伤识别,识别结果与已知损伤情况相一致。从测试自由度不完备、噪声和振型质量归一化系数这3个方面对识别效果进行了分析,结果表明:当损伤使结构动力特性有微小改变时,使用该方法不易定位损伤,应结合局部损伤识别方法进行判定;当损伤使结构动力特性有较大改变时,该方法能有效识别损伤的单元。DLV方法概念简单,理论明确,不受结构类型的限制,不需要结构的数学模型和模型缩聚或扩展技术,只需获得结构损伤前后的前几个低阶模态参数,即可识别结构一处或多处损伤,实际应用时可操作性强。     

Blind modal identification of output‐only structures in time‐domain based on complexity pursuit

Output‐only modal identification is needed when only structural responses are available. As a powerful unsupervised learning algorithm, blind source separation (BSS) technique is able to recover the hidden sources and the unknown mixing process using only the observed mixtures. This paper proposes a new time‐domain output‐only modal identification method based on a novel BSS learning algorithm, complexity pursuit (CP). The proposed concept—independent ‘physical systems’ living on the modal coordinates—connects the targeted constituent sources (and their mixing process) targeted by the CP learning rule and the modal responses (and the mode matrix), which can then be directly extracted by the CP algorithm from the measured free or ambient system responses. Numerical simulation results show that the CP method realizes accurate and robust modal identification even in the closely spaced mode and the highly damped mode cases subject to non‐stationary ambient excitation and provides excellent approximation to the non‐diagonalizable highly damped (complex) modes. Experimental and real‐world seismic‐excited structure examples are also presented to demonstrate its capability of blindly extracting modal information from system responses. The proposed CP is shown to yield clear physical interpretation in modal identification; it is computational efficient, user‐friendly, and automatic, requiring little expertise interactions for implementations. Copyright © 2013 John Wiley & Sons, Ltd.     

Modal solutions in stratified multi-layered fluid-solid half-space

A new systematic and efficient algorithm to calculate the modal solutions of multi-layered ocean-Earth model is presented. Our algorithm distinguishes itself as terseness of formulation, efficiency in numerical computation, and stableness at high frequencies, thus, thoroughly solving the problem of loss-of-precision at high frequencies. This new algorithm is expected to be very useful in studying the surface waves propagating in ocean-Earth model and other related topics.