陀螺系统随机振动分析的辛本征展开方法

探讨了受随机载荷作用下陀螺阻尼系统随机动力响应问题.虚拟激励法作为随机振动分析的一种高效、精确方法已经广泛应用于结构抗震、抗风等工程领域.在以单类物理变量描述的Lagrange(拉格朗日)体系框架下,振型分解方法已被有效应用于上述随机振动问题的模型自由度缩减.然而,对于陀螺系统的随机振动问题,由于陀螺效应的存在,基于Rayleigh商本征值的振型分解方法受到很大限制.对此,首先给出了陀螺系统辛本征值问题的一般形式.然后对于受平稳随机载荷激励的陀螺系统(无阻尼或有阻尼)引入虚拟激励法,基于辛本征空间展开推导了系统随机振动响应功率谱的求解列式;对于仅考虑陀螺效应的保守系统(无阻尼),该求解列式可以表述为一个显式表达式.在数值算例中,应用该文提出的方法分析了平稳随机载荷作用下一类阻尼陀螺系统的随机振动响应问题,通过与其它方法进行对比,验证了该方法的精确性和有效性.     

车桥耦合系统随机振动的时域显式解法

在桥面和轨道随机不平顺作用下,车桥耦合系统振动是一个典型的非平稳随机振动问题.笔者分别建立表征物理演变机制的车辆系统和桥梁系统的动力响应显式表达式,然后利用车桥之间的运动相容条件,建立车桥之间接触力关于桥面不平顺的显式表达式.在此基础上,即可直接利用统计矩运算法则,获得车桥接触力的统计矩演化规律,并进一步计算车辆系统和桥梁系统关键响应的演变统计矩.此外,也可以基于车桥接触力关于桥面不平顺的显式表达式,高效地进行随机模拟（即Monte Carlo模拟, MCS）,以获得车桥耦合系统关键响应的演变统计矩及其他统计信息.在上述过程中,由于实现了车桥耦合系统物理演变机制和概率演化规律的相对分离,在响应统计矩计算中,无需反复求解车桥耦合系统的运动微分方程,且可以仅针对车桥接触力及其他所关注的关键响应开展降维计算,大幅提高了车桥耦合系统随机振动的计算效率.数值算例表明,所提出的方法具有理想的计算精度和计算效率.     

三维粘弹性分层介质中平稳随机波的传播

研究平稳随机波在粘弹性分层横观各向同性介质中的传播问题．将岩层考虑为分层介质,各层性质不同,岩层位于基岩上面,并且认为基岩比岩层刚很多,在基岩处给出随机激励．在频率和波数域中将控制方程化为常微分方程求解．对常微分方程,应用两点边值问题的精细积分法进行求解．因此,近年来发展的应用于结构随机振动的虚拟激励法可推广于当前分层岩层响应的计算．     

大跨度叠合梁斜拉桥多维多点非平稳随机地震响应分析

为研究地震动非平稳特性对大跨度叠合梁斜拉桥随机地震响应的影响,采用绝对位移求解的虚拟激励法在通用有限元软件中对某座叠合梁斜拉桥进行多维多点非平稳随机地震响应分析.研究结果表明:采用绝对位移求解的虚拟激励法,能在通用有限元软件中高效地实现对大跨度叠合梁斜拉桥多维多点非平稳随机地震响应分析;设计中假定地震动平稳性,得到的结果偏于保守;行波效应对大跨度叠合梁斜拉桥结构响应有显著影响,考虑行波效应对主塔塔顶位移和塔底内力是有利的,但应注意其对主梁跨中位移和内力的不利影响.     

非平稳随机激励下系统首次穿越概率的近似解法

提出了非平稳Gausss白噪声激励下线性系统条件首次穿越概率的近似解析解.该近似解基于VanMarcke近似,但是,因为引进了随机过程和界限水平的标准化,VanMarcke公式中的期望衰减率可由响应的二阶统计矩获得,而不需要知道响应的相关函数或谱密度函数.给出了非平稳激励下线性系统响应的显式二阶统计矩.调制白噪声激励下单自由度线性系统的首次穿越概率分析说明了该方法的精度、效率和应用过程.     

地震作用下管道-土耦合结构非平稳随机振动分析

考虑地震载荷的非平稳特性和相干特性,建立了地下管道-土耦合结构随机振动分析的一种频域方法.应用Fourier-Stieltjes积分描述非平稳随机过程,由核函数刻画地震载荷的幅值和频率成份等时变特性;同时引入指数型衰减函数描述地震载荷的空间分布特性.在频域上结合虚拟激励法和Fourier(傅里叶)分析技术,推导了地下管道-土耦合结构受具有相干特征非平稳随机载荷作用下的响应演变功率谱的闭合解,建立了随机振动输入和输出的频域关系.数值算例与传统的分析方法进行对比验证,说明了建立方法的正确性和有效性;进一步研究了不同土参数变化以及不同端部约束条件下管道结构随机振动行为规律.     

A Complex Mode Method for Wind-Induced Responses of 6-Parameter Practical Viscoelastic Damping Energy Dissipation Structures Based on the Davenport Wind Speed Spectrum北大核心CSCD

针对六参数实用黏弹性阻尼耗能结构,基于Davenport风速谱系列响应问题进行了系统的研究.首先,利用六参数黏弹性阻尼器的微分型本构关系,建立了耗能结构基于Davenport风速谱激励下的运动方程;然后,运用复模态法将耗能结构的运动方程由二阶微分方程转化为一阶方程,获得了耗能结构系统对风振激励响应的频域解和功率谱密度函数表达式;最后,利用数学恒等式,基于随机振动理论获得了耗能结构系统在Davenport风速谱激励下的响应和阻尼器受力的解析解.该文方法不仅考虑了结构系统在风振激励作用下全振型展开的结果,表达式较现有结果更为简便,效率及精度更高,且适用于非经典阻尼结构.     

车辆随机振动功率谱分析的虚拟激励法概述

车辆受到路面或轨道不平度激励会产生随机振动,对这类随机振动进行快速有效的分析对于提高车辆性能具有非常重要的意义．而这类随机振动的分析,历来都只能将车辆模拟为自由度很低的弹簧质量体系,其计算精度难于保证．特别在需要计算疲劳应力集中的情况,需要使用网格很密的有限元模型,计算困难就更为突出．此外,在考虑多个车轮受到路面随机激励时,就更难应用效率相对较高的频域分析来计算车辆的随机振动;而只能借助于一条或几条路面不平度样本,让车轮在上面移动,借助数值积分工具而粗略地得到随机响应的统计特性．效率很低而且精度不高．虚拟激励法对于克服上述困难具有很好的效果;目前已在我国的汽车、火车、磁浮列车等领域获得日益广泛的应用．该文对此作一概略的叙述．     

广义例外华罗庚域的Bergman核函数

计算出两类广义例外华罗庚域的Bergman核函数的显式表达式. 同时利用Appell的多变量的超几何函数给出了例外华罗庚域的Bergman核函数在边界点附近的渐近性质.     

非齐次小周期复合材料稳态热问题的多尺度渐近展开和收敛性分析

利用多尺度渐近展开和均匀化思想讨论了小周期复合材料的稳态热问题,得到了非齐次边界条件下二阶椭圆型方程的渐近解,并给出了原始解与渐近解之间的误差估计,数值结果表明了结论的正确性.     

Multiscale damage simulation

Numerous materials show a softening behaviour at dynamic loading. The decrease of stress is caused by the evolution in the microscale in terms of areas where the local stiffness is reduced, e.g. due to micro-void growth. For a numerical treatment of this material behaviour, phenomenological damage approaches are used in daily engineering practice. For a better understanding of the micromechanical process of such phenomenological models, multiscale methods are becoming increasingly important. The physical quantities that are responsible for the microstructural evolution associated with the damage process are transferred into the numerical model. In this context, the method of configurational forces will be used to describe the geometrical changes of damaged areas. With the help of homogenization, macro- and microscale will be coupled. In consequence, each Gaussian point of the macroscale is modelled by an own microstructure (RVE), where the microscale evolves during the loading process according to observable damage phenomena. Hereby, we present the general case of hyperelastic materials at finite strains. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The heterogeneous multiscale method as a framework for coupling the finite element method and molecular dynamics

Hierarchical two-scale methods are computationally very powerful as there is no direct coupling between the macro- and microscale. Such schemes develop first a microscale model under macroscopic constraints, then the macroscopic constitutive laws are found by averaging over the microscale. The heterogeneous multiscale method (HMM) is a general top-down approach for the design of multiscale algorithms. While this method is mainly used for concurrent coupling schemes in the literature, the proposed methodology also applies to a hierarchical coupling. This contribution discusses a hierarchical two-scale setting based on the heterogeneous multi-scale method for quasi-static problems: the macroscale is treated by continuum mechanics and the finite element method and the microscale is treated by statistical mechanics and molecular dynamics. Our investigation focuses on an optimised coupling of solvers on the macro- and microscale which yields a significant decrease in computational time with no associated loss in accuracy. In particular, the number of time steps used for the molecular dynamics simulation is adjusted at each iteration of the macroscopic solver. A numerical example demonstrates the performance of the model. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A mixed multiscale finite element method for elliptic problems with oscillating coefficients

The recently introduced multiscale finite element method for solving elliptic equations with oscillating coefficients is designed to capture the large-scale structure of the solutions without resolving all the fine-scale structures. Motivated by the numerical simulation of flow transport in highly heterogeneous porous media, we propose a mixed multiscale finite element method with an over-sampling technique for solving second order elliptic equations with rapidly oscillating coefficients. The multiscale finite element bases are constructed by locally solving Neumann boundary value problems. We provide a detailed convergence analysis of the method under the assumption that the oscillating coefficients are locally periodic. While such a simplifying assumption is not required by our method, it allows us to use homogenization theory to obtain the asymptotic structure of the solutions. Numerical experiments are carried out for flow transport in a porous medium with a random log-normal relative permeability to demonstrate the efficiency and accuracy of the proposed method.     

Nonlinear Frequency Analysis of FGM Pipes Based on the Homotopy Method北大核心CSCD

Based on the homotopy analysis method, the nonlinear vibration of porous functionally graded material (FGM) conveying pipes under generalized boundary conditions was studied. Based on the power-law distribution of the FGM and the Voigt model, the physical properties of the porous pipe material were described. Under the Euler-Bernoulli beam theory and the von Kármán nonlinear theory, and by means of Hamilton’s variational principle, the dynamic control equations and generalized boundary conditions for porous FGM conveying pipes were established. The homotopy analysis method was used to solve the nonlinear vibration characteristics of the porous FGM conveying pipe under generalized boundary conditions. The numerical results show that, the translation spring has little effect on the critical velocity of instability, while the rotation spring increases the critical velocity of instability, making the system more stable; in the nonlinear system, the viscoelastic coefficient does not change the critical velocity; the pipe length, the power-law exponent and the porosity all influence the nonlinear free vibration of the porous FGM conveying pipe. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Macroscale modelling of multilayer diffusion: Using volume averaging to correct the boundary conditions

This paper investigates the form of the boundary conditions (BCs) used in macroscale models of PDEs with coefficients that vary over a small length-scale (microscale). Specifically, we focus on the one-dimensional multilayer diffusion problem, a simple prototype problem where an analytical solution is available. For a given microscale BC (e.g., Dirichlet, Neumann, Robin, etc.) we derive a corrected macroscale BC using the method of volume averaging. For example, our analysis confirms that a Robin BC should be applied on the macroscale if a Dirichlet BC is specified on the microscale. The macroscale field computed using the corrected BCs more accurately captures the averaged microscale field and leads to a reconstructed microscale field that is in excellent agreement with the true microscale field. While the analysis and results are presented for one-dimensional multilayer diffusion only, the methodology can be extended to and has implications on a broader class of problems.     

模糊随机桁架结构的动力特性双因子分析法

针对模糊随机桁架结构的动力特性分析,提出了一种新的模糊随机有限元方法．当结构的物理参数和几何尺寸同时具有模糊随机性时,利用模糊因子法和随机因子法建立了结构刚度矩阵和质量矩阵;从结构振动的Rayleigh商表达式出发,利用区间运算推导出结构动力特性模糊随机变量的计算表达式;之后利用随机变量的矩法和代数综合法,推导出结构特征值的数字特征的计算式．通过算例分析了模糊随机桁架结构参数的模糊随机性对其动力特性的影响．该方法的优点是能准确反映结构某一参数的模糊随机性对结构特征值及其数字特征的影响．     

Derivation of Macroscale Equations for a Class of Physical Applications

A macroscale formulation is constructed from a system of partialdifferential equations which govern the microscale dependent variables. Theconstruction is based upon transformations of equations on second ordercontact manifolds under conditions of integrability. Necessary conditions onthe structure of the macroscale equations are obtained under the requirementthat the solutions of the macroscale equations satisfy, in some approximatesense, the equations associated with the microscale. This approach offers amethod whereby one can construct only those macroscale equations that can bevalidated by a condition of consistency based on the model error. Themethodology is employed to construct a turbulence closure model forincompressible flow. It is shown that the large eddy viscosity, whichsatifies contemporary tests based on Galilean invariance, fails theconsistency condition defined here.