求解Rayleigh阻尼系数的加权最小二乘法

在地震反应分析过程中,提出了一种优化方法以解决Rayleigh阻尼系数计算时选择两阶合理参考频率的难题。该方法是以反应谱理论为基础,以结构位移峰值误差最小为目标函数。将位移反应谱用一阶Taylor级数近似计算,从而将目标函数简化为加权最小二乘法的方程。随后以框架结构为例,讨论了模态个数和阻尼比模型对Rayleigh阻尼系数计算的影响,并与传统方法、最小二乘法及基于多参考振型的加权最小二乘法进行比较。计算结果表明,最小二乘法由于忽略了模态贡献的影响,不是计算Rayleigh阻尼系数的合理方法。当模态个数所包含的累积振型参与质量达90%以上,本文方法所得Rayleigh阻尼系数计算结果稳定,结构动力反应的计算精度高。     

基于微分求积法求解非粘滞阻尼系统的时程响应

研究一种含有指数型非粘滞阻尼线性多自由度振动系统的时程分析问题。该非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表述为质点速度与核函数的卷积。由于阻尼模型的改变,常用的数值积分方法（如Newmark-β法、Wilson-θ法）不能直接应用于这种非粘滞阻尼系统。基于一种无条件稳定的微分求积方法,给出了这种非粘滞阻...     

用虚拟激励法求解非比例阻尼线性体系的非平稳随机地震响应

应用复振型分解方法,将非比例阻尼线性体系在地震作用下的动力方程求解问题转化为若干个广义复振子的求解与叠加问题。通过假定地震地面运动为一零均值的非平稳随机激励,应用虚拟激励法原理,推导得到了广义复振子动力坐标计算的一般公式,进而得到了非比例阻尼线性体系非平稳随机地震响应计算的一般解答。由于可以选择少量共轭复振型的影响进行计算,对于大型复杂非比例阻尼结构,其随机地震响应计算工作量可以大幅度减小。算例证实了这种方法的可靠性及可行性。     

基于数值分析的模态阻尼系数获取方法

结构模态阻尼系数是影响振动疲劳特性的主要因素,获取模态阻尼系数对于结构振动疲劳的分析和仿真计算有重要作用,对于揭示金属材料振动疲劳损伤形成机理有直接意义.本研究针对2024-O铝合金,进行了大量的元件级试件的振动疲劳仿真分析、试验研究以及数值分析计算,并提出了一种基于数值分析的快速获取结构模态阻尼系数的方法,适合于获取试验件振动疲劳试验过程中的模态阻尼系数变化趋势.研究表明:在不中断试件疲劳试验的情况下,本文方法可以快速准确地得到试验件在整个振动疲劳历程中的模态阻尼系数,固有频率的相对偏差小于0.06％,模态阻尼系数的相对偏差小于1％.为进一步揭示金属材料振动疲劳损伤的形成机理奠定了基础.     

建立P-S-N曲线的双加权最小二乘法

为充分利用成组疲劳试验数据和疲劳强度试验数据,获得较高精度的P-S-N曲线,提出了一种建立P-S-N曲线的双加权最小二乘法。该方法充分考虑了疲劳试验方法、疲劳试验数据分散性和试件样本容量对P-S-N曲线建立的影响,构建了不同来源误差的分配模型。它首先利用疲劳试验的最小样本容量判据处理成组试验数据,导出试验数据分布p分位点的置信区间长度,并进行第一次加权拟合。再结合疲劳强度试验数据,进行第二次加权拟合,推导出了任意可靠度下的P-S-N曲线。作者对17组金属材料和PMMA(聚甲基丙烯酸甲酯)材料进行了统计分析以验证所建立的P-S-N曲线的拟合能力、合理性和预测精度,结果表明建立的P-S-N曲线具有较高的精度和可信度。     

非比例阻尼线性体系平稳随机地震响应计算的虚拟激励法

应用复振型分解方法,将非比例阻尼线性体系在地震作用下的动力方程求解问题转化为若干个广义复振子的求解与叠加问题。通过假定地震地面运动为一零均值的平稳随机激励,应用虚拟激励法原理,推导得到了广义复振子动力坐标的解析计算公式,进而得到了以复振型为基础的非比例阻尼线性体系随机地震响应计算的一般实数解析解答。算例证实了这种方法的可靠性及高效率。     

卷积型最小二乘法求解梁的动力学问题

提出了卷积型最小二乘法,并用其计算了不同初始条件和不同边界条件下梁的动力学问题,算例表明,方法概念简单,计算方便,精度高,计算工作量少,卷积型最小二乘法是计算结构动力学问题的一种简单、高效的方法,试函数可以不满足边界条件.     

求解非比例阻尼体系复模态的实模态摄动法

根据工程结构的实际情况,建立了非比例阻尼结构体系复模态特性的近似求解方法------实模态摄动法. 这一方法以复Ritz向量展开原理为基础,把非比例阻尼结构体系复模态特性的分析过程分解为两个基本步骤,首先以结构体系的实模态向量构建复Ritz向量的求解子空间,然后通过非线性复代数方程组的求解代替扩阶后的复特征值方程的求解,从而简化了计算过程.通过两个算例表明:这一方法不仅计算简便,而且具有较高的计算精度和执行效率,对于复杂的非比例阻尼系统是很适用的,具有一定的工程应用价值.     

门字型输电塔风振响应分析及风振系数的确定

门字型塔是一种综合性能比较好的钻越塔,风荷载对这种高柔结构的影响很大,但实际设计中风振系数按规范统一取值不是很合理,且取值较大,偏于保守.故本文针对门字型输电塔在脉动风作用下的响应、风振系数的求法及合理取值进行专门研究.主要建立了门字型塔ANSYS有限元模型,进行了模态分析;以Davenport风速谱为脉动风功率自谱,应用AR法编制程序顺风模拟了风荷载时程曲线,将由模拟风计算出的风载加载至门字型塔各塔段上对其进行动力时程分析,计算了风振系数并与高耸结构设计规范值进行了比较.研究结果表明:在风载作用下,塔身和横担均会产生较为明显的动力响应;塔身、横担风振系数可取均值,分别为1.39、1.59;规范计算得到的风振系数和本文时程分析得到的风振系数变化规律基本接近,但规范取值较大.     

基于最小二乘法的反应流拉格朗日分析在冲击起爆中的研究

现有反应流拉格朗日分析方法在已知粒子速度情况下,求解方法仍有不足。本文针对这一情况,将反解法与自洽检验法相结合,提出了基于最小二乘法的反应流拉格朗日反解法,该方法的理论精度能够实现应力沿路径线的 M阶导数为零（M为迹线数）,并且该方法能够满足自洽检验法。为了验证该方法的有效性,利用该方法对一组小隔板冲击起爆试验数据进行了处理,对比了本文方法、试验处理以及传统的反解法处理结果,表明该方法不仅可以适当减小迹线上的偶然误差,还能够使得迹线函数更好地反映各物理量沿迹线的变化性态。     

非粘滞阻尼系统时程响应分析的精细积分方法

考虑一个具有非粘滞阻尼特性的多自由度系统响应的时程分析问题.该非粘滞阻尼模型假设阻尼力与质点速度的时间历程相关,数学表达式体现为阻尼力等于质点速度与某一核函数的卷积.在利用状态空间方法将系统运动方程转换成一阶的状态方程的基础上,采用精细积分方法对状态方程进行数值求解,得到一种求解该阻尼系统时程响应的精确、高效的计算方法.通过两个数值算例表明,采用该方法得到几乎精确的数值计算结果,而且计算效率有成数量级的提高.     

The Rayleigh quotient method for computing eigenvalue bounds of vibrational systems with interval parameters

In this paper, an effective method, the Rayleigh quotient iteration, for computing the bounds of structures with interval parameters is presented. When the uncertainty of the structural parameters is described by using the interval, i.e. an unknown-but-bounded version, the structural vibration eigenvalues will become the interval. It is desirable to give these interval eigenvalues. An application of this method is illustrated by a numerical example. The results show that the Rayleigh quotient iteration method is very effective for constructing the upper and lower bounds on the eigenvalues of structures with interval parameters.The Project Supported by National Natural Science Foundation of China     

场地土的随机地震反应分析

研究了剪切模量随深度呈线性变化时场地土的随机地震反应分析方法, 分别建立了基岩输入地震动加速度功率谱函数为白噪声和过滤白噪声时场地土的随机地震反应分析理论。     

基于状态观测器的带分布参数体系地震响应控制

针对带分布参数体系地震响应控制问题,建立带分布参数体系的运动控制方程,推导出应用于带分布参数体系的序列最优模态控制算法。根据序列最优模态控制算法振动控制时对结构全状态反馈的需要,提出了针对序列最优模态控制算法的全维状态观测器,形成基于全维状态观测器的序列最优模态控制算法,然后通过测量容易测量的结构地震响应,由全维状态观...     

变截面梁动力响应的频域传递矩阵计算方法

为提高变截面梁振动分析的计算效率,提出了基于频域传递矩阵法的动力计算算法．首先,选择线速度、角速度、弯矩和剪力作为求解变量,通过Laplace变换将变截面梁的动力响应时域偏微分方程转换为频域常微分方程;然后,通过求解频域方程并结合协调和边界条件建立变截面梁的频域传递矩阵;通过构造傅里叶级数展开形式的时域响应函数,对变截面梁传递矩阵方法求解的频响函数进行Laplace逆变换,建立了变截面梁的固有特性计算和时域瞬态响应计算方法,最后,借助数值仿真软件,开发了变截面梁动力响应分析的计算程序．完成对算例的仿真计算和分析,并与有限元计算结果进行对比,数值结果验证了该方法的正确性和有效性．     

有限元模型阻尼特性的复参数修正方法研究

阻尼对于结构动力学响应具有重要的影响,但有限元模型一般很难对阻尼特性进行精确建模.基于实测频响函数,研究了一种有限元模型阻尼特性的复参数修正方法.以待修正区域各单元质量、刚度矩阵的比例修正系数为复修正参数,建立了单元矩阵比例修正的灵敏度方程直接算法,并对比分析了复修正参数与不同阻尼特性之间的数学关系.以六自由度集中参数模型和25杆平面桁架模型为例,验证了复参数修正方法在阻尼特性修正中的有效性.     

爆破地震波三要素对多层砌体结构弹塑性地震响应的影响

利用Matlab编写了一套Newmark时程分析法与刚度退化二线型恢复力模型相结合求解多层砌体结构体系弹塑性地震响应的程序.以一个具体的4层砌体结构为分析对象,输入荷载用实测及人工模拟爆破地震波作用的形式,并按某原则变换输入爆破地震波的特性参数,根据所求出的结构各弹塑性地震响应幅值结果,分别讨论了爆破地震波幅值、主频和...     

沉管隧道横向动力响应分析的多体动力学方法

沉管隧道横向地震响应分析是沉管隧道结构抗震设计的重点.根据沉管隧道的结构特点和接头剪力键的节点构造,建立了沉管隧道的多刚体-剪切铰-黏性阻尼铰多体动力学模型,并应用多体动力学中的离散时间传递矩阵法(multibody system discrete time transfer matrix method,MS-DT-TMM)建立了数学模型及表达式,同时编写Matlab程序进行求解实现.设定算例分别采用MS-DT-TMM和传统有限元法(finite element method,FEM)对沉管隧道横向地震响应进行对比分析,两者吻合较好,表明MS-DT-TMM应用于沉管隧道横向动力响应分析的可行性和有效性.