覆冰三分裂导线舞动的数值模拟方法

以覆冰三分裂导线为研究对象,提出了一种模拟覆冰分裂导线舞动的数值方法。用Hamilton变分原理建立系统的动力学平衡方程,利用罚函数法引入子导线上间隔棒连接点的运动约束条件;采用Newmark法进行时间积分、Newton-Raphson法迭代求解非线性方程。通过数值算例验证了方法的正确性。在该方法的应用中由于尾流影响而导致迎风侧子导线和背风侧子导线所受空气动力载荷不同,数值模拟结果反映了这一因素对各子导线舞动轨迹的影响。此方法为分裂导线舞动的深入研究提供了一种有效途径。     

温度效应对覆冰导线舞动特征的影响研究

以某耐张段单跨导线作为研究对象,基于哈密顿原理和增量热场理论,建立了考虑温度效应的覆冰导线非线性舞动模型并推导出其舞动方程。分别使用近似模态和复杂模态得到离散后的面内-面外耦合的有限维度舞动方程,接着采用多尺度法得到该方程的近似解。数值算例研究表明,温度对频率的影响与Irvine参数有关,对张力的影响与初始张力的范围有关,而且降温会增大导线舞动幅值;采用精确模态计算得到的舞动幅值显著小于使用近似模态计算得到的幅值。     

两种离散方法对覆冰四分裂导线舞动特征的影响

首先基于Hamilton变分准则推导了覆冰导线三自由度耦合动力学方程组,接着采用两种离散方法离散动张力应变,一种是直接使用Galerkin法离散,另一种是先将动张力应变等效处理,然后再使用Galerkin法离散.由风洞试验测得新月形覆冰四分裂导线的气动力系数,并将各子导线的气动力系数等效,接着利用泰勒展开式拟合气动力系...     

一种新的覆冰导线舞动非线性有限元分析方法

覆冰导线舞动严重威胁输电线路的安全稳定运行,论文针对覆冰导线的大幅舞动问题,基于Hamilton原理建立了耦合三个平动自由度和一个扭转自由度以及偏心覆冰作用下导线的非线性舞动方程,提出了一种新的覆冰导线舞动分析非线性有限元模型,将相邻跨导线和绝缘子等效为线性弹簧单元,同时考虑覆冰导线非线性气动力和大幅舞动的几何非线性,采用Newmark-β时间积分结合修正的Newton-Raphson非线性迭代求解舞动有限元方程.论文数值解与D形覆冰导线舞动实测值无论在振幅还是频率方面均吻合得非常好,证明了论文新提出的导线舞动有限元模型的准确性.论文研究表明:舞动是一种以上下运动为主的低频振动,通常发生在一阶上下自振频率附近,其振幅和频率由输电线路的物理参数和风载荷唯一确定,与初始运动状态无关.通过与现有方法对比,论文在覆冰导线的舞动模型和气动力处理方面有较大的创新和改进,计算精度更高,效率较传统有限元法有显著提升,能为输电线路舞动研究和防舞抑舞控制提供一种新的方法与思路.     

双分裂导线相间间隔棒防舞效果数值模拟研究

采用ABAQUS有限元软件建立覆冰双分裂导线特征段线路模型,针对三角形、垂直和水平三种导线排布方式,数值模拟特征段线路的舞动过程,以及在安装相间间隔棒后的动力响应过程.基于线路在安装相间间隔棒前后的动力响应结果,对三种导线排布方式下相间间隔棒的防舞效果进行对比分析,得到相间间隔棒对垂直和三角形排布线路有较明显的防舞效果,而对水平排布线路几乎无效果的结论.     

气动力⁃安培力⁃谐波激励作用下载流导线的非线性振动特征

研究了气动力-安培力-受迫激励联合作用下载流新月形覆冰导线的振动特征,在传统只考虑安培力的载 流导线振动方程中引入了气动力和受迫激励项,使得模型更加符合实际工程．首先建立了载流新月形覆冰输电 导线模型,接着推导了载流覆冰导线模型的振动方程,利用Galerkin方法将该振动方程转变为有限维度的常微 分方程,采用多尺度法进行摄动求解,求得幅-频响应函数,并对比分析主共振情况下考虑气动力与不考虑气 动力时载流、间距、风速、张力、激励幅值对振动幅值的影响．研究结果表明：考虑气动力的振动幅值较未考 虑气动力时小．各线路参数对导线的幅值都有一定程度的影响,并且使系统表现出明显的软弹簧特征和复杂的 动力学行为．本文研究成果能为实际工程提供一定参考价值     

D形覆冰六分裂导线空气动力系数数值模拟

利用FLUENT流体动力学分析软件计算D形覆冰六分裂导线在特定风速下的绕流问题,得到各覆冰子导线空气动力系数随风攻角的变化曲线.分析了迎风侧子导线尾流对背风侧子导线气动特性的影响.结果表明,D形覆冰六分裂导线各子导线的空气动力系数存在明显差异,在研究舞动及其防舞技术时有必要考虑这种差异.     

A numerical method for simulation of attached cavitation flows

A new numerical algorithm for attached cavitation flows is developed. A cavitation model is implemented in a viscous Navier–Stokes solver. The liquid–vapour interface is assumed as a free surface boundary of the computation domain. Its shape is determined with an iterative procedure to match the cavity surface to a constant pressure boundary. The pressure distribution, as well as its gradient along the wall, is taken into account in updating the cavity shape iteratively. A series of computations are performed for the cavitating flows across three kinds of headform/cylinder bodies: conic, ogival and hemispheric heads. A range of cavitation numbers is investigated for each headform/cylinder body. The obtained results are reasonable and the iterative procedure of cavity shape updating is quite stable. The superiority of the developed cavitation model and algorithm is demonstrated. Copyright © 2006 John Wiley & Sons, Ltd.     

弹性曲梁在机械和热载荷共同作用下的几何非线性模型及其数值解

基于Euler-Bernoulli梁的几何非线性理论,建立了弹性曲梁在任意分布机械载荷和热载荷共同作用下的几何非线性静平衡控制方程。该模型不仅计及了轴线伸长,同时也精确地考虑了梁的初始曲率对变形的影响以及轴向变形与弯曲变形之间的相互耦合效应。应用打靶法数值求解了半圆形曲梁在横向均匀升温作用下的非线性弯曲问题,数值比较了轴向伸长对曲梁变形的影响。     

混凝土中爆炸数值仿真算法研究

用AUTODYN动力学软件对装药在混凝土介质靶中爆炸作数值仿真,研究了不同的网格描述方法,以及网格划分方法对计算结果的影响。通过实际算例比较了Lagrange算法和SPH算法以及不同的网格划分细度对靶径向、轴向不同距离的最大冲击压力、靶背过载以及损伤区、炸坑大小的模拟计算结果。分析得到了混凝土介质靶中爆炸问题数值仿真的合适算法,即拉氏算法对炸点近处的较高的最大冲击压力有较好的敏感能力,SPH算法对仿真计算爆坑形态以及抛掷体速度等具有明显的优势。对于算例问题,拉氏算法10 mm网格和SPH算法3 mm粒子即可满足一定的计算精度。     

不同自由度下覆冰四分裂导线舞动特征分析

基于哈密顿变分准则,推导了单档导线的振动控制方程,通过建立的气动荷载模型得到考虑气动荷载影响的覆冰导线舞动控制方程．进行风洞试验并得到覆冰四分裂导线的三分力气动系数,利用泰勒定律将气动系数在攻角α = 0处展开,得到了气动系数的三次拟合曲线．采用摄动法求解了二自由度下覆冰导线的位移响应,接着采用数值法求解了三自由度下覆冰导线的位移响应．最后,通过具体的数值算例分析了两者位移响应的区别．算例分析表明：二自由度下覆冰导线的位移响应与三自由度下覆冰导线的位移响应存在一定的区别,三自由度下覆冰导线的舞动更符合导线真实状态下的舞动,其结果更具有参考价值．     

海冰动力学数值模拟中改进的PIC方法

为了准确地模拟海冰的动力过程,需要建立精确有效的数值方法。本文结合质点网格法(PIC)和光滑质点流体动力学方法(SPH)发展了一种改进的PIC方法。该方法在欧拉坐标下对海冰动量方程进行差分计算,在拉格朗日坐标下进行海冰质点位移、厚度和密集度计算,并采用Gauss函数进行欧拉网格点与拉格朗日质点间海冰参数的交互插值。采用改进的PIC方法对规则区域内的海冰堆积过程进行了数值试验,对渤海海冰的动力过程进行了72小时数值模拟。计算结果均表明改进的PIC方法具有计算量小,计算结果平稳精确的优点,可很好地适用于海冰动力作用过程的数值模拟。     

Nonlinear three-dimension analysis for axially symmetrical circular plates and multilayered plates

Analytic nonlinear three-dimension solutions are presented for axially symmetrical homogeneous isotropic circular plates and multilayered plates with rigidly clamped boundary conditions and under transverse load.The geometric nonlinearily from a moderately large deflection is considered.A developmental perturbation method is used to solve the complicated nonlinear three-dimension differential equations of equilibrium.The basic idea of this perturbation method is using the two-dimension solutions as a basic form of the corresponding three-dimension solutions,and then processing the perturbation procedure to obtain the three-dimension perturbation solutions.The nonlinear three-dimension results in analytic expressions and in numerical forms for ordinary plates and multilayered plates are presented.All of the plate stresses are shown in figures.The results show that this perturbation method used to analyse nonlinear three-dimension problems of plates is effective.     

New method for geometric nonlinear analysis of large displacement drill strings

Introduction Drillstringsina3Dlargedisplacementwellboreareusuallyregardedas3Dcurved beamswithlargedisplacementandlargerotationbutwithsmallstrain,thus,therelationship betweenstressandstrainislinear.Thedeformationofdrillstringscanbecomposedoftwo parts:oneiscalledtheinitialdeformationcausedwhentheaxisofthedrillstringisdeformed fromitsstraightconfigurationtothepositioncoincidentwiththewellaxis,involvinglarge rotationsandwiththeimplicationoffiniterotationalproblem.Theotheristhedeformation fromthew…     

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.     

岩石双孔爆破过程的流形元法模拟

在流形元基本理论基础上,引入断裂力学的裂纹产生及扩展判据,应用二阶流形元对岩石双孔爆破过程进行了研究,模拟了双孔同时起爆和毫秒延时起爆两种条件下裂纹的产生和扩展、块体的形成以及漏斗形成过程,分析了爆炸载荷作用下岩石的破坏规律以及起爆条件对该过程的影响,验证了流形元法在研究冲击载荷作用下岩石从连续体到不连续体破坏过程的准确性和有效性,为研究类似问题提供一个新的思路和方法。     

计算力学非线性分析的现状与展望

本文对计算力学非线性分析近年来的发展做了简单回顾。文中说明了非线性数值分析日益受到重视并且还需要走很长的路,进而介绍了近廿年来非线性问题算法的两个最要方面:延续算法和单纯形算法。文末对数值方法未来十年的可能发展做了展望。     

模拟长杆弹侵彻混凝土靶的MCA方法

采用改进的移动元胞自动机法数值软件,对长杆弹侵彻混凝土靶进行了二维数值模拟,给出了弹靶的破坏变形过程,得到了侵彻深度与冲击速度的关系曲线。数值模拟结果与已有的实验现象吻合较好,说明该计算方法可以有效地计算和模拟高速侵彻问题。     

动力松弛法在板为成形过程数值模拟中的应用

论文阐述了动力松驰法的基本原理,推导出它在中心差分动力显式算法中的应用列式,并将其用于锥形凸缘的圆形杯拉深成形过程的数值模拟研究,获得了满意的结果。     

入水问题的Euler-Lagrange流固耦合数值模以技术研究

针对结构入水过程的特点,建立了一种把Eulerian计算过程和Lagrangian计算过程耦合起来的数值技术,编写了计算程序,把结构对水作用的过程和结构动力学响应过程的计算结合起来,实现了结构与水之间流固耦合作用过程的模拟.其中水、空气的运动利用Eulerian过程的流体弹塑性有限差分计算程序进行模拟,结构的运动变形利用Lagrangian过程的结构动力学有限元计算程序模拟结构响应.通过板入水问题的计算结果与商业软件计算结果进行比较,验证了计算方法对流固间相互作用处理的正确性.本文计算程序可用于一般入水问题的结构响应-流场运动分析中.