Nonstationary probability densities of system response of strongly nonlinear single-degree-of-freedom system subject to modulated white noise excitation

The nonstationary probability densities of system response of a single-degree-of -freedom system with lightly nonlinear damping and strongly nonlinear stiffness subject to modulated white noise excitation are studied.Using the stochastic averaging method based on the generalized harmonic functions,the averaged Fokker-Planck-Kolmogorov equation governing the nonstationary probability density of the amplitude is derived. The solution of the equation is approximated by the series expansion in terms of a set...     

基于Wilson-θ和Newmark-β法的非线性动力学方程改进算法

Wilson-θ法和Newmark-β法是非线性动力学方程求解的常用方法。它们的一个基本步骤是,将方程改写为增量平衡的形式,在每一个积分步长内用状态参量修正平衡方程的系数矩阵,其本质是在单个步长内对系统的非线性环节进行了线性化处理。本文基于增量思想分别改进了Wilson-θ法和Newmark-β法,根据即时解给出下一步的猜测解,然后对猜测解进行迭代校正,最终得到收敛的近似解。算例表明,改进算法的精度更高,且收敛准则简单。更为重要的是,本文方法无须对非线性项进行线性化处理,因而计算效率更高,适应范围更广。     

Analysis of efficient preconditioned defect correction methods for nonlinear water waves

Robust computational procedures for the solution of non‐hydrostatic, free surface, irrotational and inviscid free‐surface water waves in three space dimensions can be based on iterative preconditioned defect correction (PDC) methods. Such methods can be made efficient and scalable to enable prediction of free‐surface wave transformation and accurate wave kinematics in both deep and shallow waters in large marine areas or for predicting the outcome of experiments in large numerical wave tanks. We revisit the classical governing equations are fully nonlinear and dispersive potential flow equations. We present new detailed fundamental analysis using finite‐amplitude wave solutions for iterative solvers. We demonstrate that the PDC method in combination with a high‐order discretization method enables efficient and scalable solution of the linear system of equations arising in potential flow models. Our study is particularly relevant for fast and efficient simulation of non‐breaking fully nonlinear water waves over varying bottom topography that may be limited by computational resources or requirements. To gain insight into algorithmic properties and proper choices of discretization parameters for different PDC strategies, we study systematically limits of accuracy, convergence rate, algorithmic and numerical efficiency and scalability of the most efficient known PDC methods. These strategies are of interest, because they enable generalization of geometric multigrid methods to high‐order accurate discretizations and enable significant improvement in numerical efficiency while incuring minimal storage requirements. We demonstrate robustness using such PDC methods for practical ranges of interest for coastal and maritime engineering, that is, from shallow to deep water, and report details of numerical experiments that can be used for benchmarking purposes. Copyright © 2014 John Wiley & Sons, Ltd.     

Nonlinear transient transfinite element thermal analysis of thick-walled FGM cylinders with temperature-dependent material properties

An algorithm for investigation of nonlinear systems by the transfinite element method is presented. Basically, the transformation techniques have been developed for linear systems. Nonlinear transient heat transfer of a thick FGM cylinder with temperature-dependent material properties is investigated in the present paper to clarify the proposed algorithm. Two main novelties of the present research are: (1) incorporating the temperature-dependency of the material properties in the thermal analysis which lead to highly non-linear governing equations and (2) proposing an updating numerical transfinite element procedure to solve the resulted highly nonlinear governing equations. To reduce the effect of the artificial local heat source generation at the mutual boundaries of the elements, second order elements are used. Influences of various boundary conditions, geometric parameters, and volume fraction indices on the temperature distribution are investigated. Results of the proposed transfinite element technique show a good agreement with those obtained using the iterative time integration or analytical method. Furthermore, results reveal the significant effect of the temperature-dependency of the material properties. The present solution algorithm prevents numerical oscillations and damping, and accumulated time integration errors. The present technique may be used to obtain relatively accurate and stable results in a less computational time.     

窄带随机噪声作用下单自由度非线性碰撞系统的响应

研究了单自由度非线性单边碰撞系统在窄带随机噪声激励下的次共振响应问题。用Zhuravlev变换将碰撞系统转化为速度连续的非碰撞系统,然后用随机平均法得到了关于慢变量的随机微分方程。在没有随机扰动情形,得到了系统响应幅值满足的代数方程;在有随机扰动的情形下,给出了系统响应稳态矩计算的迭代公式。讨论了系统阻尼项、非线性项、随机扰动项和碰撞恢复系数等参数对于系统响应的影响。理论计算和数值模拟表明,系统响应幅值将在激励频率接近于次共振频率时达到最大。而当激励频率逐渐偏离次共振频率时,系统响应迅速衰减。     

A thermo-mechanically coupled finite strain model considering inelastic heat generation

The procedure for reuse of finite element method (FEM) programs for heat transfer and structure analysis to solve advanced thermo-mechanical problems is presented as powerful algorithm applicable for coupling of other physical fields (magnetic, fluid flow, etc.). In this case, nonlinear Block-Gauss–Seidel partitioned algorithm strongly couples the heat transfer and structural FEM programs by a component-based software engineering. Component template library provides possibility to exchange the data between the components which solve the corresponding subproblems. The structural component evaluates the dissipative energy induced by inelastic strain. The heat transfer component computes the temperature change due to the dissipation. The convergence is guaranteed by posing the global convergence criterion on the previously locally converged coupled variables. This enables reuse of software and allows the numerical simulation of thermo-sensitive problems.     

基于特征分裂有限元准隐格式的共轭传热整体耦合数值模拟方法

共轭传热现象在科学和工程领域中大量存在. 随着计算能力的发展, 对共轭传热现象进行准确有效的数值模拟, 成为科学研究和工程设计上的重要挑战.共轭传热数值模拟的方法可以分为两大类: 分区耦合和整体耦合.本文采用有限元法对共轭传热问题进行整体耦合模拟. 固体传热求解采用标准的伽辽金有限元方法.流动求解采用基于特征分裂的有限元方法. 该方法是一种重要的求解流动问题的有限元方法, 可以使用等阶有限元. 该方法的准隐格式与其他格式相比, 具有时间步长大的特点. 将稳定项中的时间步长与全局时间步长分开, 改进了准隐格式的稳定性. 基于改进的特征分裂有限元方法的准隐格式, 发展了一种层流共轭传热数值模拟的整体耦合方法. 采用这种方法可以将流体计算域和固体计算域作为一个整体划分有限元网格, 并且所有变量都可以采用相同的插值函数, 从而有利于程序的实现. 通过对典型问题的模拟, 验证了这种方法的准确性. 本工作还研究了固体区域时间步长对定常共轭传热问题数值模拟收敛性的影响.      

Modeling of simultaneous heat and mass transfer processes in ammonia–water absorption systems from general correlations

This paper presents a general differential mathematical model to analyze the simultaneous heat and mass transfer processes that occur in different components of an ammonia–water absorption system: absorber, desorber, rectifier, distillation column, condenser and evaporator. Heat and mass transfer equations are considered, taking into account the heat and mass transfer resistances in the liquid and vapour phases. The model considers the different regions: vapour phase, liquid phase and an external heating or cooling medium. A finite difference numerical method has been considered to solve the resulting set of nonlinear differential equations and an iterative algorithm is proposed for its solution. A map of possible solutions of the mass transferred composition z is presented when varying the interface temperature, which enables to establish a robust implementation code. The analysis is focused on the processes presented in ammonia–water absorption systems. The model is applied to analyze the ammonia purification process in an adiabatic packed rectification column and the numerical results show good agreement with experimental data.     

Effects of thermophoresis particle deposition and of the thermal conductivity in a porous plate with dissipative heat and mass transfer

Network simulation method(NSM) is used to solve the laminar heat and mass transfer of an electricallyconducting,heat generating/absorbing fluid past a perforated horizontal surface in the presence of viscous and Joule heating problem. The governing partial differential equations are non-dimensionalized and transformed into a system of nonlinear ordinary differential similarity equations,in a single independent variable,η. The resulting coupled,nonlinear equations are solved under appropriate transformed boundary conditions. Computations are performed for a wide range of the governing flow parameters,viz Prandtl number,thermophoretic coeffcient(a function of Knudsen number),thermal conductivity parameter,wall transpiration parameter and Schmidt number. The numerical details are discussed with relevant applications. The present problem finds applications in optical fiber fabrication,aerosol filter precipitators,particle deposition on hydronautical blades,semiconductor wafer design,thermo-electronics and problems including nuclear reactor safety.     

处理一种非线性问题迭代法的计算机方法

本文依非线性问题的积分方程组,按迭代法得出一般迭代表达式.由此,任意次迭代解的所有待定系数可用计算机计算,使研究复杂的非路性问题成为可能.     

半空间弹性波双参数强非线性逆散射

半空间弹性波强散射源逆问题是一个实用而又难度很大的强非线性问题，而经典的Born迭代方法只适用于弱非线性问题．基于DBI方法（Distorted-Born Iterative方法），提出了一种求解弹性波强非线性逆散射问题的迭代方法．在数值模拟运算时利用矩量法进行离散处理，并采用正则化原则避免求解病态矩阵方程。计算结果表明该方法具有较快的收敛速度及较高的精度。     

一类非线性周期系统响应的精细积分法

对于一类非线性周期／变系数微分方程,提出基于精细积分法的数值解法,处理非线性周期／变系数微分方程系统的响应问题,其积分策略是：采用精细积分格式处理常系数部分;采用线性插值格式处理非线性周期／变系数部分,既继承精细积分方程高度准确的特点,又保证足够的精度与较小的计算量。通过数值算例,与以往与用的微分方程直接数值积分法（如预估－校正哈明法）求得的解加以比较表明,对于给定的精度要求,精细积分法更经济有效     

一种求解瞬态热传导方程的无条件稳定方法

瞬态热传导问题普遍存在于航空航天、土木和冶金等领域中, 对这类问题精确、高效的数值求解方法一直备受关注. 为此, 本文提出了一种无条件稳定的单步时间积分方法. 在所提出的方法中, 拉格朗日插值函数被用来近似真实的温度场及其一次导数场, 之后, 加权残量法被用来建立二者之间的关系. 通过对算法参数的巧妙设计, 本文提出的方法具有二阶精度和L型数值耗散, 其中, L型耗散使得本文方法能够快速过滤掉虚假的高频振荡. 多数现有时间积分方法仅对线性系统具有无条件稳定性, 对非线性系统则是条件稳定的. 为此, 本文改进了Hughes针对一阶非线性热传导问题提出的时间积分方法稳定性评估理论, 并将改进的理论用于方法的参数设计中. 理论分析的结果表明本文方法对线性和非线性热传导系统均是无条件稳定的. 即使对于著名的Crank-Nicolson方法失稳的非线性热传导问题, 本文方法仍能给出稳定且精确的预测. 数值测试结果显示, 所提出的方法相较于当前流行的方法具有明显的精度、耗散和稳定性优势.      

Non-Darcian effects on conjugate natural convection between horizontal concentric cylinders filled with a porous medium

An analysis is presented of steady conjugate free convection between two horizontal concentric cylinders filled with a fluid-saturated porous medium; the innermost cylinder surface is maintained at a high temperature and the outermost cylinder surface at a lower one. The velocity-pressure-gradient relation is taken to be nonlinear, with departure from the linear Darcy situation measured by a parameter F₀. The investigation is based on the numerical solution, by a finite-difference method, of the full momentum and energy equations. The streamline and isotherm patterns as well as the local and mean Nusselt numbers are plotted for several physical parameters to show some of the flow and heat transfer characteristics. It is found that all parameters play an important role in the flow and heat transfer characteristics. The model can be applied to a variety of engineering problems.     

Mixed convection flow of a micropolar fluid with radiation and chemical reaction

Effects of heat and mass transfer on the mixed convection flow of a magnetohydrodynamic (MHD) micropolar fluid bounded by a stretching surface have been investigated. Homotopy analysis procedure is adopted for computations of a set of coupled nonlinear ordinary differential equations. Numerical values of skin friction coefficient and Nusselt and Sherwood numbers are worked out. A comparative study is provided with the limiting available numerical solution. Copyright © 2010 John Wiley & Sons, Ltd.     

On a Rayleigh Wave Equation with Boundary Damping

In this paper an initial-boundary value problem for a weakly nonlinear string (or wave) equation with non-classical boundary conditions is considered. One end of the string is assumed to be fixed and the other end of the string is attached to a dashpot system, where the damping generated by thedashpot is assumed to be small. This problem can be regarded as a simple model describing oscillations of flexible structures such as overhead transmission lines in a windfield. An asymptotic theory for a class ofinitial-boundary value problems for nonlinear wave equations is presented. Itwill be shown that the problems considered are well-posed for all time t. A multiple time-scales perturbation method incombination with the method of characteristics will be used to construct asymptotic approximations of the solution. It will also be shown that all solutions tend to zero for a sufficiently large value of the damping parameter. For smaller values of the damping parameter it will be shown how the string-system eventually will oscillate. Some numerical results are alsopresented in this paper.     

基于FPK方程的减摆器的非线性阻尼测试方法

受高斯白噪声外激的一阶非线性动力学方程能通过求解对应的FPK方程得到精确稳态解．本文基于这一结果导出减摆器非线性阻尼力与系统速度输出的概率结构的关系,将动力学系统中非线性阻尼力参数的测试问题转化测量系统的概率结构,并通过仿真进行了验证．     

Monolithic approach to thermal fluid-structure interaction with nonconforming interfaces

This paper presents a monolithic approach to the thermal fluid-structureinteraction(FSI) with nonconforming interfaces.The thermal viscous flow is governedby the Boussinesq approximation and the incompressible Navier-Stokes equations.Themotion of the fluid domain is accounted for by an arbitrary Lagrangian-Eulerian(ALE)strategy.A pseudo-solid formulation is used to manage the deformation of the fluid do-main.The structure is described by the geometrically nonlinear thermoelastic dynamics.An efficient data transfer strategy based on the Gauss points is proposed to guarantee theequilibrium of the stresses and heat along the interface.The resulting strongly coupledset of nonlinear equations for the fluid,structure,and heat is solved by a monolithicsolution procedure.A numerical example is presented to demonstrate the robustness andefficiency of the methodology.     

任意四边形薄板大挠度弯曲问题研究

本文用pb-2 Ritz能量法求解任意四边形薄板的几何非线性弯曲问题.首先通过坐标变换将任意四边形区域转换到一个2×2单位正方形求解区域,并建立求解域内的能量泛函,然后取一个完备的二元多项式级数(p-2)与描述边界形状的一个基本函数(b)的乘积作为Ritz函数,由能量最小原理建立结构的刚度方程,非线性代数方程采用拟牛顿法求解.文中给出了详细的数学公式,并对四边简支方板和四边固定菱形板进行了数值分析,算例表明,本方法具有公式简单、精度较高的优点,用以分析大挠度问题是非常有效的.     

Transient multiboiling in a pin fin with temperature dependent thermal conductivity

Transient heat conduction in fins undergoing different kinds of convective processes (film, transition, nucleated boiling and natural convection) as in multiboiling processes take place, is a strongly non-linear problem because of the abrupt changes in the heat transfer coefficient that occur at certain temperatures. Transient equations for the thermal fields and fluxes are solved simultaneously, giving the time constant of the process; the stationary solution is compared with the numerical or experimental values of other authors. Temperature dependencies of the heat transfer coefficient and the thermal conductivity is assumed due to the large interval of temperatures occurring. Network Simulation Method is used for the numerical solution, which gives simultaneously thermal field of temperatures and heat fluxes.