求解含缺陷一维周期结构动力响应的高效算法

基于周期结构的动力特性和群理论,建立了一种高效求解含缺陷一维周期结构动力响应的数值方法。在求解结构动力响应时,高效求解结构对应的线性代数方程组最为关键。采用凝聚技术,可减小结构对应线性代数方程组的规模。基于周期结构动力系统中线性代数方程组的特性,通过一个小规模含缺陷结构和一维周期结构的响应分析,可得到含缺陷一维周期结构的动力响应。同理,一维周期结构的动力响应可通过一系列小规模结构的响应分析得到,且小规模结构的动力响应可基于群理论高效求解。数值算例表明,本文算法有较高的求解效率。     

周期结构动力响应的高效数值方法

提出一种计算周期结构动力响应的高效率算法. 以精细积分方法为基础, 利用周期结构的对称性和动力问题的物理特性, 分析了周期结构对应矩阵指数的特殊结构, 并基于此给出一种计算周期结构对应矩阵指数的高效率方法. 在高效和精确计算周期结构对应矩阵指数的基础上, 得到了周期结构动力响应的高效率和高精度算法. 数值算例表明, 该方法效率高且节省存储要求.      

应用同伦分析法研究Mathieu-Duffing振子的周期解

应用同伦分析法研究了Mathieu-Duffing振子的周期解,展示了Mathieu-Duffing振子的周期1和周期2解的求解过程,通过求解构造的非线性代数方程组而获得周期解,应用Floquet理论判别了周期解的稳定性。比较了同伦分析方法得到的周期解和数值方法得到的周期解,结果表明两者具有一致性。     

一维离散结构能带结构与表面态的辛分析方法

本文基于辛几何方法推导了一维离散周期结构、半无穷周期结构和含杂质半无穷周期结构的本征方程,力求建立一个完整的辛分析体系。通过辛分析,将一维离散半无限周期结构转化到一个元胞上求解,大大简化了计算量。对于含杂质半无穷周期结构,结合辛分析和W-W算法,给出求解含杂质半无穷周期结构本征值问题的精确、稳定和高效算法。数值算例说明了本文算法的有效性。     

有限元瞬态响应分析软件的一种并行开发方法

瞬态响应分析是有限元动力分析的重要内容之一，而串行计算机上运行的有限元软件在解题规模和速度上都受到很大限制。为此，基于系统集成思想对串行有限元软件进行并行开发。分析得出了瞬态响应分析并行开发的重点——位移方程组的瞬态响应，给出了线性方程组并行求解的思路和实现方法。用一个实例系统的实现验证了上述开发思路，从而也为并行应用软件开发探索了一条新的途径；最后对并行求解程序进行了算例验证。     

加性非平稳激励下结构速度响应的FPK方程降维

结构在随机激励下的非线性响应分析是具有高度挑战性的困难问题. 对于白噪声或过滤白噪声激励,求解FPK方程将获得结构响应 的精确解. 遗憾的是,对于非线性多自由度系统,FPK方程难以直接求解. 事实上,其数值解法严重受限于方程维度,而解析求解 则仅适用于少数特定的系统,且多是稳态解. 因此,将FPK方程进行降维,是求解高维随机动力响应分析问题的重要途径. 本文针 对幅值调制的加性白噪声激励下多自由度非线性结构的非平稳随机响应分析问题,将联合概率密度函数满足的高维FPK方程进行降 维. 针对结构速度响应概率密度函数求解,通过引入等价漂移系数,原FPK方程可转化为一维FPK型方程. 建议了构造等价漂移系数 的条件均值函数方法. 进而,采用路径积分方法求解降维FPK型方程,得到速度概率密度函数的数值解答. 结合单自由度Rayleigh 振子、十层线性剪切型框架和非线性剪切型框架结构在幅值调制的加性白噪声激励下的非平稳速度响应求解,讨论了本文方法的精 度和效率,验证了其有效性.      

加性非平稳激励下结构速度响应的FPK方程降维

结构在随机激励下的非线性响应分析是具有高度挑战性的困难问题.对于白噪声或过滤白噪声激励,求解FPK方程将获得结构响应的精确解.遗憾的是,对于非线性多自由度系统,FPK方程难以直接求解.事实上,其数值解法严重受限于方程维度,而解析求解则仅适用于少数特定的系统,且多是稳态解.因此,将FPK方程进行降维,是求解高维随机动力响应分析问题的重要途径.本文针对幅值调制的加性白噪声激励下多自由度非线性结构的非平稳随机响应分析问题,将联合概率密度函数满足的高维FPK方程进行降维.针对结构速度响应概率密度函数求解,通过引入等价漂移系数,原FPK方程可转化为一维FPK型方程.建议了构造等价漂移系数的条件均值函数方法.进而,采用路径积分方法求解降维FPK型方程,得到速度概率密度函数的数值解答.结合单自由度Rayleigh振子、十层线性剪切型框架和非线性剪切型框架结构在幅值调制的加性白噪声激励下的非平稳速度响应求解,讨论了本文方法的精度和效率,验证了其有效性.     

不可压饱和粘弹性多孔介质的变分原理及其无网格模拟

基于饱和多孔介质理论,在固相和液相微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,建立了流体饱和粘弹性多孔介质动力响应的若干Gurtin型变分原理,包括Hu-Washizu变分原理.利用所建立的变分原理,导出了流体饱和粘弹性多孔介质动力响应无网格数值模拟的离散控制方程,此方程是一个关于时间的对称微分方程组,便于分析计算.作为数值例子,研究了流体饱和粘弹性多孔柱体的一维动力响应,数值结果揭示了流体饱和粘弹性多孔柱体中波的传播特性以及固相粘性的影响.     

一种求解局部非线性结构稳态响应的方法

为了降低求解局部非线性结构稳态响应的计算量,基于子结构和阻抗缩聚提出了一种用于求解局部非线性结构稳态响应的计算方法.将局部非线性结构分解为线性子结构和非线性子结构,利用谐波平衡构造各个子结构的阻抗方程,对线性子结构进行缩聚,将局部非线性动力学方程转化为求解一组非线性代数方程组问题,通过迭代求解非线性代数方程组,求解系统的稳态响应.     

阶梯式Timoshenko梁自由振动的DCE解

本文基于微分容积法和区域叠加技术提出了微分容积单元法（Differential Cubature Element method，以下简称DCE方法），并用之求解阶梯式变截面Timoshenko梁的自由振动问题。根据梁的变截面情况将其划分为几个单元，在每个单元内应用微分容积法将梁的控制微分方程和边界约束方程离散成为一组关于该单元内配点位移的线性代数方程组，将这些方程组写在一起并在各单元之间应用连续性条件和平衡条件得到一组关于整个域内各点位移的齐次线性代数方程组，这是一广义特征值问题，由子空间迭代法求解该特征问题便可求得系统的自振动频率。数值算例表明，本方法能稳定收敛、并有较高的数值精度和计算效率。     

一种结构动力响应分析的快速高阶算法

为了提高基于高阶格式的结构动力响应微分求积分析方法的计算效率,发展了一种求解动力方程的快速算法．利用微分求积原理将结构动力方程转化为标准Sylvester方程的形式,通过对系数矩阵进行矩阵分解,进而将动力响应Sylvester方程化为一系列标准线性方程组,采用相关成熟算法求解这些线性方程组后即可获得结构动力时程响应的全部解答．结构动力响应微分求积分析方法为高阶数值方法,一步计算可以获得多个时点处的动力响应．基于本文快速算法,不必直接对矩阵方程进行求解．数值算例表明,本文快速算法能够准确地计算出结构动力响应,具有数值精度高、收敛性好的优点．     

滑动轴承油膜力Jacobi矩阵的一种快速算法

基于变分不等方程理论,把滑动轴承油膜力及其Jacobi矩阵的求解转换为求解一组三对角矩阵代数方程,采用一种修正的追赶法同步快速求解。同时将系数矩阵分解为轴颈运动相关量和常矩阵乘积叠加的形式,常矩阵一次求得反复调用,大大减少冗余运算。算例表明,采用本文算法,结果精度得到保证,运算时间大幅减少,能很好地揭示滑动轴承-转子系统的非线性动力特性。     

Steady state response analysis for fractional dynamic systems based on memory-free principle and harmonic balancing

A semi-analytic approach is proposed to analyze steady state responses of dynamic systems containing fractional derivatives. A major purpose is to efficiently combine the harmonic balancing (HB) technique and Yuan–Agrawal (YA) memory-free principle. As steady solutions being expressed by truncated Fourier series, a simple yet efficient way is suggested based on the YA principle to explicitly separate the Caputo fractional derivative as periodic and decaying non-periodic parts. Neglecting the decaying terms and applying HB procedures result into a set of algebraic equations in the Fourier coefficients. The linear algebraic equations are solved exactly for linear systems, and the non-linear ones are solved by Newton–Raphson plus arc-length continuation algorithm for non-linear problems. Both periodic and triple-periodic solutions obtained by the presented method are in excellent agreement with those by either predictor–corrector (PC) or YA method. Importantly, the presented method is capable of detecting both stable and unstable periodic solutions, whereas time-stepping integration techniques such as YA and PC can only track stable ones. Together with the Floquet theory, therefore, the presented method allows us to address the bifurcations in detail of the steady responses of fractional Duffing oscillator. Symmetry breakings and cyclic-fold bifurcations are found and discussed for both periodic and triple-periodic solutions.     

Construction of approximate analytical solutions to strongly nonlinear damped oscillators

An analytical approximate method for strongly nonlinear damped oscillators is proposed. By introducing phase and amplitude of oscillation as well as a bookkeeping parameter, we rewrite the governing equation into a partial differential equation with solution being a periodic function of the phase. Based on combination of the Newton’s method with the harmonic balance method, the partial differential equation is transformed into a set of linear ordinary differential equations in terms of harmonic coefficients, which can further be converted into systems of linear algebraic equations by using the bookkeeping parameter expansion. Only a few iterations can provide very accurate approximate analytical solutions even if the nonlinearity and damping are significant. The method can be applied to general oscillators with odd nonlinearities as well as even ones even without linear restoring force. Three examples are presented to illustrate the usefulness and effectiveness of the proposed method.     

Simulation of elastic wave diffraction by multiple strip-like cracks in a layered periodic composite

The problem of numerical simulation of the steady-state harmonic vibrations of a layered phononic crystal (elastic periodic composite) with a set of strip-like cracks parallel to the layer boundaries is solved, and the accompanying wave phenomena are considered. The transfer matrix method (propagator matrix method) is used to describe the incident wave field. It allows one not only to construct the wave fields but also to calculate the pass bands and band gaps and to find the localization factor. The wave field scattered by multiple defects is represented by means of an integral approach as a superposition of the fields scattered by all cracks. An integral representation in the form of a convolution of the Fourier symbols of Green’s matrices for the corresponding layered structures and a Fourier transform of the crack opening displacement vector is constructed for each of the scattered fields. The crack opening displacements are determined by the boundary integral equation method using the Bubnov-Galerkin scheme, where Chebyshev polynomials of the second kind, which take into account the behavior of the solution near the crack edges, are chosen as the projection and basis systems. The system of linear algebraic equations with a diagonal predominance of components arising when the system of integral equations is discretized has a block structure. The characteristics describing qualitatively and quantitatively the wave processes that take place under the diffraction of plane elastic waves by multiple cracks in a phononic crystal are analyzed. The resonant properties of a system of defects and the influence of the relative positions and sizes of defects in a layered phononic crystal on the resonant properties are studied. To obtain clearer results and to explain them, the energy flux vector is calculated and the energy surfaces and streamlines corresponding to them are constructed.     

Theory of dynamic processes in mechanical systems and materials of regular structure

The theory of vibrations and waves in natural and synthesized materials of regular structure is analyzed. Models based on different averaging and continualization methods are outlined. Emphasis is on periodically inhomogeneous structures. The exact solutions are obtained and analyzed using the closed-form solution of infinite algebraic systems, representing equations in Hamiltonian operator form and solving them based on the theory of differential equations with periodic coefficients, mode selection rule, and methods of drawing wave shapes at limit and arbitrary frequencies     

非饱和土-结构动力响应的多耦合周期性有限元法

真实的地基土体-隧道系统中土体及结构性质往往沿线路纵向变化.为考虑土体与结构沿纵向的变化特性,提出了一种非饱和土-结构系统动力响应分析的多耦合周期性有限元法.首先基于非饱和土的实用波动方程,采用Galerkin法推导了单节点5个自由度的非饱和土ub-pl-pg格式有限元表达式,相比于单节点9个自由度的ub-v-w格式有...     

非线性转子系统周期响应的分析

基于时域的时间有限元法,将描述转子系统动力学特征的非线性微分方程组离散成一组非线性代数方程,然后应用吴消去法的特征列思维对所得到的非线性代数方程组进行降维求解,进而得到待求节点位移响应的解形式,并据此对一具有非线性支撑的柔性Jeffcott转子模型响应的性质进行了分析。     

Stress concentrations in the particulate composite with transversely isotropic phases

The accurate series solution have been obtained of the elasticity theory problem for a transversely isotropic solid containing a finite or infinite periodic array of anisotropic spherical inclusions. The method of solution has been developed based on the multipole expansion technique. The basic idea of method consists in expansion the displacement vector into a series over the set of vectorial functions satisfying the governing equations of elastic equilibrium. The re-expansion formulae derived for these functions provide exact satisfaction of the interfacial boundary conditions. As a result, the primary spatial boundary-value problem is reduced to an infinite set of linear algebraic equations. The method has been applied systematically to solve for three models of composite, namely a single inclusion, a finite array of inclusions and an infinite periodic array of inclusions, respectively, embedded in a transversely isotropic solid. The numerical results are presented demonstrating that elastic properties mismatch, anisotropy degree, orientation of the anisotropy axes and interactions between the inclusions can produce significant local stress concentration and, thus, affect greatly the overall elastic behavior of composite.