共查询到20条相似文献，搜索用时 171 毫秒
1.
框架结构损伤识别的改进直接解析法 总被引：1，自引：1，他引：0
在一般直接解析法的基础上进行模态截尾误差的改进,从而大大降低了测量模态数目,而且识别效果也很好,但迭代次数多且收敛比较慢,本文通过改进方程解法来提高其速度。作者将约束最小二乘法迭代迅速的特点和极小最小二乘法搜索寻优功能相结合,提出了一种混合迭代求解方法,改进的方法充分组合了一般方法的搜索能力强和约束最小二乘法求解精度高的优点,从而解决了改进模态截尾误差后迭代收敛慢的问题。通过对五层两跨的平面框架结构进行损伤识别数值模拟分析后表明改进的方程解法是有效的,与一般的方法相比,算法收敛速度明显加快,求解精度明显提高,验证了改进后方法的可行性与准确性。 相似文献
2.
利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性. 相似文献
3.
《世界地震工程》2017,(1)
提出了一种基于曲线拟合的改进型FDD算法Im CFDD,以加速度自功率谱矩阵的SVD分解得到的主奇异值为对象,推导了主奇异值与模态参数间的关系。推导表明:主奇异值与加速度自功率谱函数是等价的,奇异值可以分解为关于频率ω的分子和分母均为4阶的比值形式。据此建立了求解模态参数的表达式,并给出了迭代拟合算法。所提算法同频域分解算法FDD、增强型的频域分解算法EFDD、曲线拟合频域分解算法CFDD、频域空间分解算法FSDD有着本质的不同,在公式的推导上做了最小的简化,同时考虑了负频率部分和共轭极点部分的贡献,得到的模态参数求解方程是非线性的,需要迭代求解。最后,采用数值模型验证了新算法的有效性。 相似文献
4.
5.
土层参数的随机性对场地传递函数的影响 总被引：6，自引：0，他引：6
本文采用秦勒级数展开估计土层参数随机性对场地传递函数的影响。对单层土模型在土层参数服从独立高斯分布假定下进行了公式推导和实例计算,比较了不同的泰勒级数阶次、土层参数变异系数以及阻尼比对传递系数随机性的综合效应。结果表明,即使在小参数摄动范围内采用一阶泰勒级数也是不够的;一般情况下采用二阶泰勒级数即可给出满意的估计。 相似文献
6.
传统有限差分系数是通过泰勒级数展开求取的,这样导致所计算的频散曲线在大波数区域会产生较强的数值误差.针对二阶空间偏导数的显式有限差分离散,本文发展了一种新的优化差分系数方法：首先将泰勒级数展开与多点采样方法结合应用于空间频散关系,基于最大范数建立直观有效的优化目标函数,采用Remez算法求解该目标函数,从而获得最优化差分系数.利用优化有限差分方法求解三维垂直对称轴横向各向同性（VTI）介质中的声波和弹性波方程.另外,本文将二维混合吸收边界条件推广到三维VTI介质中,用于吸收人工截断边界反射;基于各向异性特征,合理调整了边界区域的速度值来提高吸收效果.考虑到三维情况下计算效率的问题,本文波场外推过程中采用图形处理器（GPU）取代传统的中央处理器（CPU）.数值精度分析表明,相比较于传统的泰勒级数展开方法,优化有限差分方法在大波数区域对频散误差的压制效果更明显.在三维均匀和修改的Hess VTI模型中的数值模拟实验证明了本文方法具有更高的精度与效率,混合吸收边界条件在三维VTI介质中具有良好的边界吸收效果. 相似文献
7.
质点的轨迹计算是半拉格朗日模式的重要基础,传统的数值计算方法由于采用时间差分代替微分,只能得到质点运动轨迹终点的速度,因此质点的移动轨迹（位移）只能靠风速外推的方法计算,导致了模式计算不稳定等问题.借鉴精细积分法中使用半解析解的思路,利用正压原始方程研究了用运动方程的半解析解构建数值模式的可能性.求解了运动方程的一阶和二阶微分方程组的半解析解,通过时间积分半解析解计算质点运动轨迹.数值试验表明,一阶微分方程组的半解析解比差分解略有优势.二阶微分方程组的半解析解在时间步长增大时优势非常明显,而且在保证计算精度的前提下,节省计算时间,这对提高模式性能有重要作用. 相似文献
8.
在前人工作的基础上,本文推导了电导率任意各向异性介质的海洋可控源电磁三维谱元法正演方程.采用一次场/二次场分离算法结合混合阶矢量基函数,可以有效避免源点的奇异性的影响,从而提高数值解的精度.采用任意六面体单元离散研究区域,有利于模拟复杂地形和地电结构.利用不完全LU分解的Induced Dimension Reduction(IDR(s))迭代算法求解线性方程组,有效地提高了求解的效率.设计典型的地电模型进行正演计算,并将计算结果与有限元解进行对比,对比结果表明本文提出的基于混合阶矢量基函数的海洋可控源电磁三维谱元数值模拟算法是正确的、有效的.本文算法具有良好的通用性,可推广用于电导率呈任意各向异性的陆地电磁、井中电磁等数值模拟研究. 相似文献
9.
本文以波动理论为基础, 半解析化求解地震勘探中常用的SH波方程. 获得的主要结果包括: 给出了二维均匀介质中SH波方程的解析解; 利用Cagniard-de Hoop方法详细推导了二维双层介质中SH波方程的解析解, 获得了透射波的解析解表达式. 同时, 基于SH波方程的解析表达式, 给出了包含各种波(如直达波、反射波、首波以及透射波)的解析解和波形图. 对于比较复杂的积分型解析解, 利用数值积分方法给出了数值结果, 并与优化的近似解析离散化方法(ONADM)和4阶Lax-Wendroff修正方法(LWC)的数值结果进行了比较, 以验证解析解的正确性. 本文的研究成果有望在检验波动方程数值新方法的有效性、波传播理论分析等方面得到应用. 相似文献
10.
利用泰勒级数展开和弹性力学的几何方程,推导了球坐标系中由GPS位移数据解算图形单元应变和旋转张量的解析公式. 通过线性化处理,并利用误差传播定律,详细推导了图形单元应变和旋转张量的误差公式. 利用川滇地区最新的GPS测站位移速率数据,采用图形单元法解算了该区的面应变率及最大剪切应变率分布,并对其进行了初步的分析. 同时阐述了图形单元应变数学模型的局限性,讨论了图形单元法计算应变的意义以及图形单元的选择问题,进一步分析了应变计算的定权方法,还讨论了GPS观测网图形单元的地心半径差异与应变的关系问题. 相似文献
11.
A step‐by‐step approximate procedure taking into consideration high‐frequency modes, usually neglected in the modal analysis of both classically and non‐classically damped structures, is presented. This procedure can be considered as an extension of traditional modal correction methods, like the mode‐acceleration method and the dynamic correction method, which are very effective for structural systems subjected to forcing functions described by analytical laws. The proposed procedure, herein called improved dynamic correction method, requires two steps. In the first step, the number of differential equations of motion are reduced and consequently solved by using the first few undamped mode‐shapes. In the second step, the errors due to modal truncation are reduced by correcting the dynamic response and solving a new set of differential equations, formally similar to the original differential equations of motion. The difference between the two groups of differential equations lies in the forcing vector, which is evaluated in such a way as to correct the effects of modal truncation on applied loads. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
12.
A Taylor basis for kinematic nonlinear real-time simulations. Part I: The complete modal derivatives
Despite todays computational power, only small nonlinear numerical substructures may be simulated in real time. The size restriction on the substructures in nonlinear finite element analysis is primarily due to the time-consuming evaluation of the internal restoring forces, which is performed element-by-element in every iteration step. The present work constitutes the first of two papers presenting a method to simulate kinematic nonlinear structures more efficiently. It involves applying a reduced basis with modal derivatives representing the nonlinearities of the system in an efficient way. Previously, the modal derivatives have been determined from a set of approximate governing equations. In the present paper, a novel set of equations governing the complete modal derivatives is derived. This is done by introducing a Taylor series into the free undamped kinematic nonlinear equations of motion. Also, the approximate governing equations are improved by introducing a novel geometric restriction. By way of an example, it is shown that only the modal derivatives determined from the complete set of equations are consistent with the Taylor series. In the second paper, it is shown that the novel modal derivatives may be used in a so-called Taylor basis and that they improve the computational time and stability significantly. 相似文献
13.
《Journal of Hydrology》2003,281(4):251-264
Practical application of geostatistical inversion to coupled problems is hampered by a number of difficulties. In this paper, we address two of them: first, the computational cost of sensitivity (Jacobian) matrices and, second, the evaluation of the relative weights of different types of data. Regarding the first, we revise the adjoint state equations to propose a form whose cost is independent of the number of unknown parameters and only grows with the number of observation wells. Regarding the second, we derive expressions for the relative weights of different types of data. These expressions are based on minimizing the expected likelihood, rather than the likelihood itself. The efficiency of both improvements is tested on a synthetic example. The example analyzes a wide range of groundwater flow and solute transport conditions. Yet, the expected likelihood consistently yields the optimal weights. The proposed form of the adjoint state equations leads to one order of magnitude reduction in CPU time with respect to the conventional sensitivity equations. 相似文献
14.
Dynamic equilibrium equations of structural systems with non‐proportional damping are coupled through the damping terms. Such coupling invalidates application of the classical modal superposition method. In this paper, a mode‐superposition pseudo‐force method is proposed. The coupled equilibrium equations are solved by an iterative process in which the coupling terms are treated as pseudo‐forces. A scale factor for each mode of the system is obtained by optimizing the iteration convergence. Through these uniquely solved scale factors, the modified modal equations not only converge much faster but also yield results with higher accuracy. A proof of the convergence of the iterative process is also presented. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
15.
16.
R. Parashar Rao S. Govindaraju 《Stochastic Environmental Research and Risk Assessment (SERRA)》2006,20(1-2):95-105
Temporal and spatial moment analysis of one-dimensional equations governing fate and transport of parent compounds along with
their transformation products is useful for parameter estimation of model parameters, and for understanding the average attributes
of contaminant behavior. The objective of this paper is to present analytical expressions for the lower order moments of members
in a sequential chain reaction, where members undergo first-order decay to produce the next member in the chain. Specifically,
moments up to second order are discussed for the first two members. For the case of purely advective transport (Peclet number
tending to infinity), temporal moment expressions are provided for more members of the chain. The sensitivity of temporal
moments is examined with respect to Peclet number and transformation rates. Spatial moments are derived by two methods—one
using Fourier transforms and another using moment generating differential equations. The behavior of lower order moments for
the first few members of a chain can be elucidated from their mathematical expressions. However, expressions for higher order
moments tend to be very complicated especially for members further down the chain. 相似文献
17.
利用混凝土拱坝地震记录识别的模态参数，可以揭示结构在地震过程中实际动力特性的变化情况，为结构地震反应分析和震后损伤评估提供重要信息。首先对模态识别常用方法的基本原理进行介绍；然后利用龙羊峡拱坝两次地震观测数据，采用"输入—输出"型和"仅考虑输出"型两类方法对大坝模态参数进行识别；最后将获得的模态识别结果与工程经验值和其他学者的研究成果进行对比，以分析评价识别结果的合理性。相关研究成果可为基于地震记录的混凝土拱坝模态参数识别的工程应用提供经验借鉴。 相似文献
18.
基于Hilbert-Huang变换和随机子空间识别技术提出了两种土木工程结构的模态参数识别方法。方法一是基于Hilbert-Huang变换和自然激励技术,通过经验模态分解和Hilbert变换提取信号的瞬时特性,进而利用自然激励技术和模态分析的基本理论识别结构的模态参数;方法二是基于经验模态分解和随机子空间识别技术,通过经验模态分解对信号进行预处理,进而运用随机子空间识别方法处理得到的结构单阶模态响应以识别结构的模态参数。利用这两种方法,通过对一12层钢筋混凝土框架模型振动台试验测点加速度记录的处理,识别了该模型结构的模态参数。识别结果与传统的基于傅里叶变换的识别结果及有限元分析结果的对比验证了这两种方法的可行性和实用性。 相似文献
19.
A modal procedure for non-linear analysis of multistorey structures with high-damping base-isolation systems was proposed. Two different isolation devices were considered in the analysis: an high-damping laminated rubber bearing and a lead-rubber bearing. Starting from deformational properties verified by tests, the isolation systems were characterized using three different analytical models (an Elastic Viscous, a Bilinear Hysteretic and a Wen's Model) with parameters depending from maximum lateral strain. After non-linear modelling of isolation and lateral-force-resisting systems, the effects of material non-linearities were considered as pseudo-forces applied to the equivalent linear system (Pseudo-Force Method) and the formally linearized equations of motion were uncoupled by the transformation defined by the complex mode shapes. The modal responses were finally obtained with an extension of Nigam–Jennings technique to non-linear and non-classically damped systems, in conjunction with an iterative technique searching for non-linear contributions satisfying equations of motion and constitutive laws. Since the properties of the isolated structure usually change with maximun lateral strain of isolation bearings, the integration of a new set of governing equations was required for each design-displacement value. The procedure proposed was described in detail and then applied for the determination of modal and total seismic responses in some real cases. At first, a very good agreement between non-linear responses obtained with the proposed mode superposition and with a direct integration method was observed. Then a comparison of results obtained with the three different analytical models of the isolation bearings was carried out. At last, the exact modal response obtained with analytical models depending from the design displacement of the isolation bearings was compared with two different approximated solutions, evaluated using mode shapes and isolation properties, respectively, calculated under simplified hypothesis.© 1998 John Wiley & Sons, Ltd. 相似文献
20.
A method for seismic design of plane steel moment resisting frames based on the use of equivalent modal damping ratios is developed. The method determines the design base shear of the structure through spectrum analysis using rationally obtained equivalent modal damping ratios instead of the crude strength reduction (behavior) factor. An equivalent linear structure, which retains the mass and initial stiffness of the original non-linear structure and takes into account geometrical non-linearity and inelasticity in the form of equivalent, time-invariant, modal damping ratios is established. The equivalent damping ratios for the first few significant modes are numerically computed by first iteratively forming a frequency response transfer function modulus until it satisfies certain smoothness criteria and then by solving a set of non-linear algebraic equations. Thus, design equations providing equivalent damping ratios as functions of period and allowable deformation and damage are constructed using extensive numerical data coming from plane steel moment resisting frames excited by various seismic motions. These equations can be used in conjunction with a design spectrum, appropriately constructed for high damping values, and modal synthesis tools to calculate the seismic design forces of the structure. The proposed method is illustrated by numerical examples. It is concluded that unlike the usual approach of seismic codes employing a single common value of the strength reduction factor value for all modes, the proposed approach working with deformation and damage dependent equivalent modal damping ratios leads to more accurate and deformation and damage controlled results. 相似文献