基于mbS模式的有限单元法

mbS模式及其有限元法是在固体和结构分析模型中引入薄膜、弯曲和剪切理论,且采用纯拉压、纯弯和纯剪单元进行分析的数值方法。在时空系中剖分物质单元和时间单元上构造以指数函数和贝塞尔函数为插入函数且按Lagrange插值条件的薄膜、弯曲和剪切等基本位移函数,由此得到更加完备和耦合的固体和结构实体单元的变形模式,根据能量泛函变分原理得到静动力有限元基本方程的一致格式。研究表明,mbS模式及其有限元法可用于梁柱和板壳等结构的静动力分析及屈曲分析。     

跳跃振子平衡及其分岔的能量分析

用能量方法研究了跳跃振子的平衡与分岔.用势能驻值条件确定了平衡位置所满足的方程,通过势能极值判断平衡的稳定性.在不同的弹簧构型下,数值计算了平衡随系统弹性刚度和质量比变化的分岔图.结果表明,弹性刚度和质量比较小时,系统只有一个稳定平衡点和一个不稳定平衡点;刚度和质量比充分大时,系统分岔出一个新的稳定平衡点和一个新的不稳定平衡点.     

薄膜单元和六面体单元的耦合计算模型

本文提出薄膜单元和六面体二次等参单元的耦合计算模型.由所编制的计算程序进行实例计算,计算结果与解析解一致.     

螺旋压榨机的受力分析及平衡问题

对螺旋压榨机的受力分析问题进行了详细的讨论;用平衡方程、虚位移原理和类比的方法推导出螺杆旋进与旋出时驱动力偶矩与压力的关系及自锁条件;并将有关结论推广到其它螺旋机构.     

脉冲液体射流泵能量平衡试验研究

根据脉冲液体射流泵主要流动部件的能量损失压力比,对其能量平衡进行试验,研究了主要过流部件的能量转化与脉冲频率之间的关系。与恒定液体射流泵试验数据进行对比,证明了脉冲频率、脉冲装置对脉冲液体射流泵的主要流动部件的能量损失压力比有一定的影响;且脉冲射流是提高射流泵传能及传质效率的有效途径。     

直接基于单元平衡的变截面梁反应分析方法

固定形状的单元位移插值函数不能合理地近似变截面梁内部的位移变化,从而影响了传统梁单元用于计算变截面梁的精度.采用直接基于单元平衡的思想给出了计算变截面梁反应的有限元方法,解决了单元位移插值函数局限性所带来的问题.导出了变截面梁单元的单元刚度矩阵、单元等效节点荷载和单元一致质量矩阵.在此基础上,利用编制的程序进行了算例验证与分析.算例验证了本文理论的正确性,表明本文方法具有很高的计算精度.     

有限单元质量质心集中方法的探讨

提出利用单元质量质心集中方法来优化单元惯性效应表征、改善单元力学效应匹配的思想: 将单元质量质心集中能够建立精确描述单元质点系惯性效应的动力学方程, 进而达到优化单元惯性效应表征以及改善单元惯性和变形两种力学效应匹配的目的. 本文通过一维问题的频散效应分析等对该方法的有效性进行了详细的分析和验证.     

准零刚度驱动式压电低频振动能量采集方法

如何高效且经济环保地为数以万计的传感器网络节点供电,是物联网快速发展和大范围应用的瓶颈性难题.将振动能转换为电能以实现传感器自供电是物联网传感器网络节点供能的潜在方案.但是,环境振动中低频成分占比较大,而传统振动能量采集方法较难实现低频(<10 Hz)振动能量的高效转化,这限制了振动能量采集技术在物联网领域的大范围应用.文章提出了一种准零刚度驱动式压电振动能量采集装置,可将环境、人体及部分机械设备的低频振动能量高效地转化为电能.首先利用能量法得到压电俘能单元的机电耦合方程,并用谐波平衡法获得系统动力学及电学响应的解析表达式,同时对比了数值解与解析表达式的结果;进一步探究了阻尼比、激励幅值等参数对动力学响应及电学输出的影响.最后加工制备了准零刚度驱动式压电振动能量采集装置样机,搭建了实验平台,测试了系统的动力学响应与电学输出,验证了理论结果的正确性.研究结果显示当频率为2.5 Hz时,准零刚度驱动式压电振动能量采集装置中单个能量转化单元的最大峰值电压达到25 V.文章提出的准零刚度压电能量采集装置有望克服传统共振型压电能量采集装置能量采集频带依赖于能量采集系统固有频率以及多稳态能量...     

基于能量原理的连续体结构安定分析

基于安定分析的静力定理,建立了一种以结构吸收的塑性能为安定判断准则,应用一种特殊的加载方式求解连续体结构安定问题的数值方法.应用这种能量方法对两个经典安定分析算例进行了计算,与相关结果进行对比表明能量法的有效性和可靠性.该方法摆脱了传统分析方法中的数学规划运算,有较高的求解效率,可方便地应用于实际工程中.     

基于等效单元概念的梁结构损伤定位方法

研究目的是为梁结构提供一种基于等效单元概念的新损伤定位方法.区分刚度损伤和质量损伤两种情况, 分别对连续梁和简支梁给出了具体的损伤定位方法. 利用的具体指标是支座反力和跨中位移. 文中的新方法都可以得到以小区间表示的损伤位置结果, 有效性得到了数值模拟算例和模型实验的验证.     

Formulation of thermo-hydro-mechanical coupling behavior of unsaturated soils based on hybrid mixture theory

Thermo-Hydro-Mechanical （THM） coupling pro- cesses in unsaturated soils are very important in both theoretical researches and engineering applications. A coupled formulation based on hybrid mixture theory is derived to model the THM coupling behavior of unsaturated soils. The free-energy and dissipative functions for different phases are derived from Taylor＇s series expansions. Constitutive relations for THM coupled behaviors of unsaturated soils, which include deformation, entropy change, fluid flow, heat conduction, and dynamic compatibility conditions on the interfaces, are then established. The number of field equations is shown to be equal to the number of unknown variables; thus, a closure of this coupling problem is established. In addition to modifications of the physical conservation equations with coupling effect terms, the constitutive equations, which consider the coupling between elastoplastic deformation of the soil skeleton, fluid flow, and heat transfer, are also derived.     

连续胶体材料化学-力学耦合有限元理论及其在ABAQUS中UEL程序实现

基于哈密顿原理,得到水凝胶的化学-力学耦合控制方程的等效积分形式和有限元形式。在整体坐标系下推导出用形函数表示的化学-力学耦合应变矩阵和单元刚度矩阵,并且得到在局部坐标系下的离散化形式。结合ABAQUS软件,编制了用户单元子程序UEL,通过数值算例验证了所开发单元的正确性,为在ABAQUS软件中实现各种耦合问题的有限元UEL编程提供了参考依据。     

Finite volume element methods for nonequilibrium radiation diffusion equations

Nonequilibrium radiation diffusion problems are described by the coupled radiation diffusion and material conduction equations. Because of the highly nonlinear, strong discontinuous, and tightly coupled phenomena, solving this kind of problems is a challenge. We construct two finite volume element schemes for the equations. One of them is monotone on many kinds of meshes, which is proved theoretically and verified by numerical tests. The other one is hard to satisfy the monotonicity, but this defect can be corrected by different repair techniques. Numerical results show that these new methods are practical and efficient on distorted meshes.Copyright © 2013 John Wiley & Sons, Ltd.     

Perturbational solutions for fuzzy-stochastic finite element equilibrium equations (FSFEEE)

I.IntroductionInthefieldofengineeringstructureanalysis,suchparametersasmodulusofmaterial-Poissonratio,geometricalsizesofstructure,boundaryconditionsandexternalloadingsprobablyhavefuzzinessor/andrandomness.Thestochasticfiniteelementmethod(SFEM)1'1disposeso…     

Finite element implementation of poroelasticity theory for swelling dynamics of hydrogels

Hydrogel can swell to many times of its dry volume, resulting in large deformation which is vital for its function. The swelling process is regulated by many physical and chemical mechanisms, and can, to some extent, be fairly described by the poroelasticity theory. Implementation of the poroelasticity theory in the framework of finite element method would aid the design and optimization of hydrogel-based soft devices. Choosing chemical potential and displacement as two field variables, we present the implementation of poroelasticity tailored for hydrogel swelling dynamics, detail the normalization of physical parameters and the treatment of boundary conditions. Several examples are presented to demonstrate the feasibility and correctness of the proposed strategy.     

脆性固体损伤和局部软化的有限元分析

本文研究混凝土、岩石一类工程中常用的应变软化材料的有限元分析方法。在作者以往有关粘塑性损伤本构模型的工作基础上，给出了一组便于有限元计算的本构方程表达式。包括损伤弹性矩阵和局部损伤软化矩阵，分别运用于计算硬化和软化阶段的有限元刚度矩阵；对所提出的本构方程的实验验证计算和一些算例的有限元数值分析，表明文中给出的本构方程是可行的，相应的有限元算法能较好地对损伤固体的局部软化效应进行数值分析，并可成功地追踪应力应变响应的软化曲线     

基于Cosserat理论的广义协调元法

首次将广义协调元理论与Cosserat理论相结合,构造了一个基于Cosserat理论的平面四节点广义协调等参元。依据常应力与线性应力下的广义协调条件,推导了广义协调位移,进而得到有限元列式。利用广义协调元,一方面克服了协调元过于刚硬的缺点,另一方面消除了非协调元不一定收敛的弱点。分析了带有圆孔的应力集中问题,可以看出,...     

Parallel finite element algorithm based on full domain partition for stationary Stokes equations

Based on the full domain partition, a parallel finite element algorithm for the stationary Stokes equations is proposed and analyzed. In this algorithm, each subproblem is defined in the entire domain. Majority of the degrees of freedom are associated with the relevant subdomain. Therefore, it can be solved in parallel with other subproblems using an existing sequential solver without extensive recoding. This allows the algorithm to be implemented easily with low communication costs. Numerical results are given showing the high efficiency of the parallel algorithm.     

Finite element method for coupled diffusion-deformation theory in polymeric gel based on slip-link model

A polymeric gel is an aggregate of polymers and solvent molecules, which can retain its shape after a large deformation. The deformation behavior of polymeric gels was often described based on the Flory-Rehner free energy function without considering the influence of chain entanglements on the mechanical behavior of gels. In this paper, a new hybrid free energy function for gels is formulated by combining the Edwards-Vilgis slip-link model and the Flory-Huggins mixing model to quantify the time-dependent concurrent process of large deformation and mass transport. The finite element method is developed to analyze examples of swelling-induced deformation. Simulation results are compared with available experimental data and show good agreement. The influence of entanglements on the time-dependent deformation behavior of gels is also demonstrated. The study of large deformation kinetics of polymeric gel is useful for diverse applications.     

Local and parallel finite element algorithms based on two-grid discretization for steady Navier-Stokes equations

Local and parallel finite element algorithms based on two-grid discretization for Navier-Stokes equations in two dimension are presented. Its basis is a coarse finite element space on the global domain and a fine finite element space on the subdomain. The local algorithm consists of finding a solution for a given nonlinear problem in the coarse finite element space and a solution for a linear problem in the fine finite element space, then droping the coarse solution of the region near the boundary. By overlapping domain decomposition, the parallel algorithms are obtained. This paper analyzes the error of these algorithms and gets some error estimates which are better than those of the standard finite element method. The numerical experiments are given too. By analyzing and comparing these results, it is shown that these algorithms are correct and high efficient.