首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 107 毫秒
1.
基于分析动力学与分析结构力学在数学理论上的一致性,在有限元分析方面,同时对时间、空间坐标离散组成混和的时空有限元网格.然后利用Hamilton变分原理,取一次变分为零,导出三维混和元列式;混和元列式矩阵的对称性,保证了混和有限元保辛的性质.数值例题表明,时-空混和有限元能灵活地处理多尺度波动问题和变动边界问题.  相似文献   

2.
飞行器跨声速气动弹性数值分析   总被引:5,自引:1,他引:4  
杨国伟  钱卫 《力学学报》2005,37(6):769-776
将流体和结构运动方程分别构造为含子迭代的计算格式,发展了一种紧耦合气动弹性分析方法.其中流体计算的空间离散采用改进的HLLEW(Harten—Lax-van Leer-Einfeldt-Wada)格式. TFI(transfinite inter- polation)方法用于生成随结构变形的自适应多块动网格.利用所发展的方法,对-翼-身-尾气动外形,数值预测了马赫数在0.3-1.3范围内的气动颤振边界.并详细研究了时间步长、子迭代步数、初始流场、耦合方法、疏密网格对颤振计算结果的影响.  相似文献   

3.
以圆柱绕流为研究对象,针对圆形边界,采用O型网格对流场进行离散,用二阶精度的中心差分有限体积法作空间离散,用二阶精度的中心差分处理时间问题,用双时间方法求解了二维非定常Navier-Stokes方程,系统研究了计算方法对收敛精度、时间步长和网格数量的依赖性.计算结果表明,对于长时间历程的非定常问题,虽然双时间方法收敛性很好,但对于分离流而言,时间步长的选取并非没有限制;每一步伪时间的推进中,收敛精度也有要求;而要模拟圆柱分离流的非线性气动力现象,计算网格至少要达到260×80的数量.  相似文献   

4.
计算含动边界非定常流动的无网格算法   总被引:1,自引:0,他引:1  
在无网格算法中考虑了含动边界的流动问题,研究了可以计算处理包含一定位移及扭转动边界非定常流动的算法.创建了无网格算法的动点法则,并引入抗扭方法对弹簧方法进行改进来处理离散点运动,提高了方法的可用度及精度.发展了求解基于无网格的ALE方程组的算法,在点云离散的基础上采用曲面逼近计算空间导数及HLLC格式计算数值通量,运用四步龙格-库塔法进行时间推进.在跨、超音速条件下,计算模拟了典型翼型简谐振动流场,计算结果与实验结果及文献对比吻合,验证了该算法的正确性.  相似文献   

5.
为了在三角形和四边形网格上采用水平集方法模拟轴对称爆轰波阵面与曲率相关的运动,假设爆 轰波阵面的法向速度是曲率的线性函数,通过坐标变换得到了轴对称坐标系下的水平集方程。水平集方程的 曲率无关项采用正格式离散,曲率项采用伽辽金等参有限元方法空间离散,时间离散采用半隐格式。算例表 明,在轴对称三角形网格和四边形网格上,含曲率的水平集方程的离散格式为强一阶精度。给出了三角形和 四边形混合网格上不光滑界面以曲率收缩的运动过程,收缩过程未出现不稳定现象。多个爆轰波阵面相互作 用的算例说明本文的格式可有效地模拟曲率相关的爆轰波的相互作用问题  相似文献   

6.
采用HLL格式,在三角形非结构网格下采用有限体积离散,建立了求解二维浅水方程的高精度的数值模型.本文采用多维重构和多维限制器的方法来获得高精度的空间格式以及防止非物理振荡的产生,时间离散采用三阶Runge-Kutta法以获得高阶的时间精度.基于三角形网格,底坡源项采用简单的斜底模型离散,为保证计算格式的和谐性,对经典的HLL格式计算的数值通量中的静水压力项进行了修正.算例证明本文提出的方法的和谐性并具有高精度的间断捕捉能力和稳定性.  相似文献   

7.
航空、航天和兵器技术等领域的研究中存在大量包含运动边界的流场。非结构重叠网格方法是一种高效的处理动边界问题的新方法。围绕相对运动的每个物体单独生成非结构网格,在网格重叠区域通过搜索和插值完成网格系之间的信息传递,提出了动态八叉树搜索算法,发展了绝对坐标系和相对坐标系相结合的流场求解方式,采用二阶精度Van Leer/Hanel格式和四阶Runge-Kutta法分别进行空间和时间离散,形成了一种新的非结构重叠网格算法。对三维Riemann问题的求解结果与精确解能很好吻合,证明了本文的重叠网格算法具有较好的时空离散精度和插值精度。对7.62mm步枪射击过程进行了数值模拟,描述了弹丸离开膛口后膛口流场的发展过程,与实验结果体现的发展过程较为吻合,验证了本文提出的非结构网格算法体系具有较好的计算性能,是研究含动边界复杂流场的一种有效手段。  相似文献   

8.
无网格法是基于散点信息求解偏微分方程问题的数值方法,无网格法可减少或完全消除对网格的依赖,数值实施更加灵活.因此,考虑采用基于径向基函数的无网格插值法求解一类分段连续型延迟偏微分方程.首先,利用θ-加权有限差分法得到方程时间上的离散格式,利用基于径向基函数的无网格插值法近似空间导数,得到了全离散数值格式.采用的基函数是Multiquadric (MQ)径向基函数,MQ径向基函数在精度及稳定性等方面都优于其他径向基函数.其次,采用傅里叶分析方法对该方法进行稳定性分析,得到了该方法稳定的条件,且该条件只与时间步长有关.最后,通过数值算例验证了方法的收敛性和稳定性,从而说明了方法的有效性和适用性.  相似文献   

9.
在空间域上采用只与结点有关的无网格方法离散,在时间域上采用精细积分方法求 解. 无网格离散过程中,利用伽辽金积分等效弱形式代替微分形式的控制方程,并 用修正变分原理满足位移边界条件,采用移动最小二乘法求解离散的形函数,把形 函数代入等效积分弱形式得到离散的二阶方程;精细积分过程中非齐次项采 用Romberg积分. 同时给出了两种不同边界条件的谐响 应求解的两个数值算例,得到了精确的数值结果.  相似文献   

10.
段庆林  李锡夔 《力学学报》2007,39(6):749-759
在有限增量微积分(finite increment calculus, FIC)的理论框架下,通过引入一个附加变量,发展了压力稳定型分步算法,有效改善了经典 分步算法的压力稳定性,同时还避免了标准FIC方法中存在的空间高阶导数的计算. 为保证 数值方法同时具有较快的计算速度和较好的健壮性,发展了有限元与无网格的耦合空间离散 方法. 该方案可在网格发生扭曲的区域采用无网格法空间离散以保证求解的精度和稳定性, 而在网格质量较好的区域以及本质边界上保留使用有限元法空间离散以提高计算效率和便于 施加本质边界条件. 方腔流考题的数值模拟结果突出地显示了所发展的压力稳定型分步算 法比经典分步算法具有更好的压力稳定性,能够有效消除速度-压力插值空间违反LBB条件而 导致的压力场的虚假数值振荡. 平面Poisseuille流动和一个典型型腔充填过程的数值模拟 结果, 表明了发展的耦合离散方案相对于单一的有限元法和单一的无网格法在综合考虑计 算效率和算法健壮性方面的突出优点.  相似文献   

11.
We present an extended finite element method (XFEM) for the direct numerical simulation of the flow of viscoelastic fluids with suspended particles. For moving particle problems, we devise a temporary arbitrary Lagrangian–Eulerian (ALE) scheme which defines the mapping of field variables at previous time levels onto the computational mesh at the current time level. In this method, a regular mesh is used for the whole computational domain including both fluid and particles. A temporary ALE mesh is constructed separately and the computational mesh is kept unchanged throughout the whole computations. Particles are moving on a fixed Eulerian mesh without any need of re-meshing. For mesh refinements around the interface, we combine XFEM with the grid deformation method, in which nodal points are redistributed close to the interface while preserving the mesh topology. Our method is verified by comparing with the results of boundary fitted mesh problems combined with the conventional ALE scheme. The proposed method shows similar accuracy compared with boundary fitted mesh problems and superior accuracy compared with the fictitious domain method. If the grid deformation method is combined with XFEM, the required computational time is reduced significantly compared to uniform mesh refinements, while providing mesh convergent solutions. We apply the proposed method to the particle migration in rotating Couette flow of a Giesekus fluid. We investigate the effect of initial particle positions, the Weissenberg number, the mobility parameter of the Giesekus model and the particle size on the particle migration. We also show two-particle interactions in confined shear flow of a viscoelastic fluid. We find three different regimes of particle motions according to initial separations of particles.  相似文献   

12.
This study carries out a complete analysis of time-space solution of hydrodynamics of pentagonal/decagonal quasicrystals. The behaviors of wave propagation for phonons and diffusion for phasons and coupling between phonon-phason fields are explored explicitly. Comprehensive discussion on physical time-space variations of all hydrodynamic field variables of the alloy quasicrystals is given. The computational specimen is simple, convenient in testing computational results, and provides a possibility that is easy to test experimentally. The quantitative results of mass density, viscosity velocities, phonon displacements, phason displacements, phonon stresses, phason stresses, viscosity stresses, and their time-space variations help us understand the motion of solid quasicrystals in a hydrodynamic condition (long-wavelength and low-frequency). The analysis presented in this paper can be used for octagonal and dodecagonal quasicrystals and is easily extended to other two-dimensional quasicrystals and three-dimensional icosahedral quasicrystals. Some problems explored by the computational results are also discussed.  相似文献   

13.
K型管节点的应力集中系数研究   总被引:1,自引:0,他引:1  
采用有限元法分析了K型管节点在承受轴力作用下的应力集中系数. 用三维固体单元模 拟结构,模拟中对不同区域分别划分网格,每个区域的网格独立产生后,再形成整个结构的 有限元网格. 使用ABAQUS软件分析了200组在轴向循环载荷作用下的节点数值模型,获得 了K节点在轴向循环载荷作用下沿焊缝的应力集中系数(SCF),并对主管和支管的SCF进行 研究,得出了节点几何参数对主管和支管的SCF及极值位置的影响规律.  相似文献   

14.
利用非线性有限元方法分析了轴向力作用下多平面TT节点的极限强度。在数值分析中,采用三维20结点固体单元模拟管道结构和焊缝形状,将结构有限元网格划分为不同区域,每个区域的网格独立产生,通过合并形成整个结构的有限元网格。通过控制位移增量法得到了加载过程中载荷和位移之间的关系曲线。使用ABAQUS软件分析了TT节点在支管端部承受轴向载荷的变形及与外部载荷之间的关系,得到了不同参数影响下的TT节点极限强度。  相似文献   

15.
混凝土细观随机骨料结构与有限元网格剖分   总被引:23,自引:1,他引:22  
在细观层次上,混凝土被认为是一种由粗骨料、水泥砂浆及二者间的粘结带所组成三相非均质复合材料。本文首先基于蒙特卡罗随机抽样原理,用“取和放”方法在计算机上产生形状、尺寸和骨料颗粒分布与真实混凝土相似的随机骨料结构,再使用有限元分析软件ANSYS对骨料区域及砂浆区域分别划分网格,并编程在骨料和砂浆之间生成三角形三结点可控制厚度粘结单元,从而使三相网格缝合为一个整体,为混凝土非线性有限元分析提供可靠的细观计算模型。最后利用建立的模型进行混凝土轴心受拉和轴心受压的仿真模拟,在细观层次研究的基础上揭示出混凝土的宏观力学性能。  相似文献   

16.
Simulating transient compressible flows involving shock waves presents challenges to the CFD practitioner in terms of the mesh quality required to resolve discontinuities and prevent smearing. This paper discusses a novel two‐dimensional Cartesian anisotropic mesh adaptation technique implemented for transient compressible flow. This technique, originally developed for laminar incompressible flow, is efficient because it refines and coarsens cells using criteria that consider the solution in each of the cardinal directions separately. In this paper, the method will be applied to compressible flow. The procedure shows promise in its ability to deliver good quality solutions while achieving computational savings. Transient shock wave diffraction over a backward step and shock reflection over a forward step are considered as test cases because they demonstrate that the quality of the solution can be maintained as the mesh is refined and coarsened in time. The data structure is explained in relation to the computational mesh, and the object‐oriented design and implementation of the code is presented. Refinement and coarsening algorithms are outlined. Computational savings over uniform and isotropic mesh approaches are shown to be significant. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

17.
欧阳武  袁小阳  杨培基  纪峰  陈红斌 《应用力学学报》2012,29(3):325-329,357,358
在转子轴承系统振动信号处理中,针对平稳信号的传统傅里叶变换精度较低、快变启动过程的非平稳信号频谱分析方法较复杂的问题,本文仿真构造了两类响应信号。通过对比给定信号参数与信号识别参数的误差研究了几种谱分析方法或过程的简便性和准确性。对转子系统振动平稳信号离散频谱分析时存在的误差进行了定量分析,利用比例插值法对误差进行校正,开发了高精度谱分析测试软件;分析了转子轴承系统快变过程非平稳振动信号的特征,探索了一种将t时空域非平稳信号转变为tn时空间域平稳信号的办法或过程,然后结合比例插值校正法对其进行频谱分析,再返回到t时空域获得某时刻的谱特征参数;构造了转子系统振动仿真信号检验了上述过程的准确性。研究结果表明:比例插值法提取的谱特征数据近乎与仿真信号设定值相等;针对本文构造的快变过程非平稳仿真信号,利用本文给出的谱分析过程产生的频率误差最大值为0.47%,幅值误差最大值为0.2%。本文的仿真研究为提出和考证新的谱分析方法提供了手段。  相似文献   

18.
覆盖件冲压仿真计算模型中网格密度分布的合理性与网格单元形态的优劣,对仿真结果的准确性有很大的影响.提出了一种改进的全四边形网格细分方法,使网格的密度分布适于覆盖件冲压分析计算要求,可保证细分后网格的协调性,并将算法推广以处理非结构化四边形网格和三角形四边形混合网格的细分.提出的网格细分策略,有助于提高细分后网格的质量.提出了适用于细分后四边形网格和非结构四边形网格的拓扑形态优化操作,可有效的提高网格模型的形态质量.  相似文献   

19.
王年华  鲁鹏  常兴华  张来平  邓小刚 《力学学报》2021,53(10):2682-2691
网格自动化生成和自适应是制约计算流体力学发展的瓶颈问题之一, 网格生成质量、效率、灵活性、自动化程度和鲁棒性是非结构网格生成的关键问题. 在非结构网格生成中, 网格空间尺度分布控制至关重要, 直接影响网格生成质量、效率和求解精度. 采用传统的背景网格法进行空间尺度分布控制需要在背景网格上求解微分方程得到背景网格上的尺度分布, 再将网格尺度从背景网格插值到真实空间点, 过程十分繁琐且耗时. 本文从效率和自动化角度提出两种网格尺度控制方法, 首先发展了基于径向基函数(RBF)插值的网格尺度控制方法, 通过贪婪算法实现边界参考点序列的精简, 提高了RBF插值的效率. 同时, 还采用人工神经网络进行网格尺度控制, 初步引入相对壁面距离和相对网格尺度作为神经网络输入输出参数, 建立人工神经网络训练模型, 采用商业软件生成二维圆柱和二维翼型非结构三角形网格作为训练样本, 通过训练和学习建立起相对壁面距离和相对网格尺度的神经网络关系. 进一步实现了二维圆柱、不同的二维翼型的尺度预测, RBF方法和神经网络方法的效率与传统背景网格法相比提高了5~10倍, 有助于提高网格生成的效率. 最后, 将方法推广应用于各向异性混合网格尺度预测, 得到的网格质量满足要求.   相似文献   

20.
A multi-velocity formulation is proposed for the solution of an Eulerian representation of an inert, disperse, and dilute particle-phase of a gas-particle flow. Single-velocity formulations are capable of predicting regions of zero particle concentration but are problematic with crossing particle trajectories or compression waves. The multi-velocity formulation described here can account for crossing particle trajectories by splitting the particle-phase into distinct velocity families which are transported separately in the flow. Switching of the particle families at solid boundaries and due to momentum transfer with the gas-phase is conducted in a manner that enforces conservation of mass, momentum, and energy. This numerical method is combined with a parallel block-based adaptive mesh refinement algorithm that is very effective in treating problems with disparate length scales. The block-based data structure lends itself naturally to domain decomposition and thereby enables efficient and scalable implementations of the algorithm on distributed-memory multi-processor architectures. Numerical results are described to demonstrate the capabilities of the approach for predicting gas-particle flows.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号