Nonlinear dynamic analysis of dielectric elastomer membrane with electrostriction

The dielectric elastomer (DE) is an important intelligent soft material widely used in soft actuators, and the dynamic response of the DE is highly nonlinear due to the material properties. In the DE, electrostriction denotes the deformation-dependent permittivity. In the present study, we formulate the nonlinear dynamic governing equations of the DE membrane considering the electrostriction effect. The free vibration and parametric excitation of the DE membrane with different geometric sizes are calculated. The free vibration bifurcations induced by the initial location and the voltage are both discussed according to an energy-based approach. The amplitude-frequency characteristics and bifurcation diagrams of parametric excitation are also given. The results show that electrostriction decreases the free vibration amplitude and increases the frequency, but it has less influence on the parametric excitation oscillation frequency and decreases the parametric excitation amplitude except when the membrane resonates. The initial location and the applied voltage can induce the snap-through instability of the free vibration. A large geometric size will lead to a much lower resonance frequency. The resonance amplitudes increase while the resonance frequencies decrease with the increase in the applied voltage. The critical voltage of snap-through instability for the parametric excitation is larger than that for the free vibration one.

     

Nonlinear dynamics and performance analysis of modified snap-through vibration energy harvester with time-varying potential function

Vibration energy harvesting has emerged as a promising method to harvest energy for small-scale applications. Enhancing the performance of a vibration energy harvester(VEH) incorporating nonlinear techniques, for example, the snap-through VEH with geometric non-linearity, has gained attention in recent years. A conventional snap-through VEH is a bi-stable system with a time-invariant potential function, which was investigated extensively in the past. In this work, a modified snap-through VEH wit...     

MPS-FEM数值分析带自由面的流固耦合问题

随着计算科学的发展,研究人员为探索流固耦合问题的物理机理而提出了众多的数值方法。其中,耦合的移动粒子半隐式方法 MPS(Moving Particle Semi-Implicit method)和有限单元法FEM(Finite Element method)为流固耦合问题的数值仿真工作提供了新的途径。本文所有流场的数值模拟工作均采用课题组自主开发的无网格法求解器MLParticle-SJTU来完成。该求解器在原始的MPS法基础上,对核函数、压力梯度模型、压力泊松方程的求解和自由面判断方式等方面进行了改进。此外,在该求解器框架内,基于FEM法拓展了针对结构场进行求解的功能。首先,对MPS和FEM方法的理论模型及其耦合策略进行了介绍。然后,采用该自研MPS-FEM耦合求解器,数值模拟了溃坝流动对弹性结构的冲击及其相互作用的标准问题。通过将结构变形及自由面波型变化等结果与已发表结果进行对比,验证了该求解器在处理带自由面剧烈变化的粘性流体和柔性变形结构的耦合作用问题上的可行性。     

无单元法（EFGM）——在岩土工程上有限元法的有力补充

作为一种新的岩土工程数值计算方法，无网格伽辽金法具有局部化技术的无单元特性。该法与传统FEM相比，在精度和后处理等方面有着明显的计算优势。本文全面介绍了EFGM，并结合岩土工程领域中的土体固结变形问题，推导了EFGM的一维、二维固结问题的离散方程；同时对算例进行计算，将其结果与精确解和传统FEM解进行比较后，得出结论：精度是较高的，处理边界是准确的。尤其当与FEM联合使用时，既能发挥FEM的经验优质，又能体现EFGM的无网格特点，计算效果十分显著，因此，EFGM是FEM的有力补充。     

Geometric and Material Nonlinear Static and Dynamic Analysis of Space Truss Structures

A formulation of the nonlinear FEM truss element that takes into account both the geometric nonlinearity and material nonlinearity of the structure is derived. The associated iterative algorithm for the application is described. Based on two inelastic material models (isotropic hardening and dynamic hardening), three space truss structures are subjected to static and to dynamic loads in order to illustrate the application of nonlinear FEM. The nonlinear FEM described can accurately trace the complex structural behavior of space truss structures, including snap-through, and buckling. The FEM results match well with the theoretical results.     

Electro-mechanical coupling bifurcation and bulging propagation in a cylindrical dielectric elastomer tube

This paper explores the critical and post-bulging bifurcation of a cylindrical dielectric elastomer (DE) tube undergoing finite deformation under electro-mechanical coupling loading. Explicit expressions for the critical conditions of electro-mechanical bifurcation are derived by using a simplified mathematical method. The post-bifurcation path is comprehensively investigated by specifying the material model as ideal dielectric elastomer. In the post-bifurcation analysis, we analytically establish conditions for the phase coexistence of steady propagation and analyze the physical implications. We demonstrate a global instability under force or voltage control and a localized instability under volume or charge control. Cylindrical tube experiments have been carried out under electro-mechanical coupling loading to verify the theoretical predictions. Good agreements on the critical conditions as well as the post-bifurcation path are obtained. This work characterizes the bifurcation mechanism of rubber-like materials under complex coupling loading.     

竖井开挖的不稳定性和岩爆的力学机制

基于稳定性理论和弹塑性厚壁筒的解析解,研究了竖井开挖过程的平衡路径,并将屈服后的井壁压力--位移平衡路径曲线应用到围岩的稳定性分析.研究表明,竖井围岩的不稳定及其岩爆的发生不仅取决于岩石材料应力应变全过程曲线的峰后特性,还取决于竖井开挖过程对应的平衡路径曲线的类型.平衡路径曲线能够描述围岩结构加载屈服后的力学行为.平衡路径的曲线形状只可能有两种类型,一种曲线上的点都处于稳定的平衡状态,虽有塑性变形的累积,但不会失稳崩落;另一种曲线则发生极值点型井壁失稳并伴有位移突跳.竖井开挖过程的平衡路径曲线是判定围岩岩爆发生的关键因素.极值点型失稳和位移突跳是均质围岩中岩爆可能发生的一种力学机制.      

INSTABILITY OF SHAFT EXCAVATION AND MECHANISM OF ROCKBURST

基于稳定性理论和弹塑性厚壁筒的解析解,研究了竖井开挖过程的平衡路径,并将屈服后的井壁压力--位移平衡路径曲线应用到围岩的稳定性分析.研究表明,竖井围岩的不稳定及其岩爆的发生不仅取决于岩石材料应力应变全过程曲线的峰后特性,还取决于竖井开挖过程对应的平衡路径曲线的类型.平衡路径曲线能够描述围岩结构加载屈服后的力学行为.平衡路径的曲线形状只可能有两种类型,一种曲线上的点都处于稳定的平衡状态,虽有塑性变形的累积,但不会失稳崩落;另一种曲线则发生极值点型井壁失稳并伴有位移突跳.竖井开挖过程的平衡路径曲线是判定围岩岩爆发生的关键因素.极值点型失稳和位移突跳是均质围岩中岩爆可能发生的一种力学机制.     

Mixed finite element and differential quadrature method for free and forced vibration and buckling analysis of rectangular plates

This paper presents a combined application of the finite element method (FEM) and the differential quadrature method (DQM) to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.     

RESEARCH DEVELOPMENTS OF SMOOTHED PARTICLE HYDRODYNAMICS METHOD AND ITS COUPLING WITH FINITE ELEMENT METHOD

拉格朗日型的有限元法和光滑粒子法在模拟材料大变形问题时各存优缺点, 而有限元与光滑粒子耦合算法实现了在小变形区域采用有限元法计算, 在局部的大变形区域采用光滑粒子法计算, 有效地综合了有限元法计算效率高和光滑粒子法能够自然地模拟材料大变形问题的特点.重点论述了有限元法、光滑粒子法以及有限元与光滑粒子耦合算法的研究现状及应用进展, 并讨论了各方法中需要进一步解决的问题.最后通过算例对3种方法的计算精度和计算效率进行了分析, 供研究人员参考.     

Self-adaptive one-dimensional nonlinear finite element method based on element energy projection method

The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.     

Axisymmetric static and dynamics snap-through phenomena in a thermally postbuckled temperature-dependent FGM circular plate

Thermal postbuckling analysis and the axisymmetric static and dynamic snap-through phenomena due to static/sudden uniform lateral pressure in a thermally postbuckled functionally graded material circular plate are performed in this research. Plate is formulated using the first order shear deformation plate theory. Thermo-mechanical properties of the plate are assumed to be temperature dependent where dependency is described according to the higher order Touloukian representation. Two types of temperature loading are considered. Uniform temperature rise and heat conduction across the thickness direction. The one dimensional heat conduction equation in the thickness direction is obtained and discreted via the central finite difference method. The obtained system of equations is nonlinear since the thermal conductivity itself is a function of the unknown nodal temperatures. Using the von-Kármán assumptions, the governing equations of the plate are obtained in a matrix representation with the aid of the conventional Ritz method whose shape functions are developed using the Gram-Schmidt process. At first thermal postbuckling analysis is performed which is a nonlinear problem with respect to both temperature and displacements. Afterwards, response of the bulged thermally postbuckled plate is obtained under the static and dynamic uniform pressure. Snap-through phenomenon may be observed in both static and dynamic loading cases, due to the immovability of the edge of the plate and the initial deflection caused by postbuckling deflection. To capture the snapping phenomenon and trace the path beyond the limit loads, cylindrical arch-length technique is used. In dynamic snap-through analysis, the effect of structural damping is also included. Numerical results of this study reveal that the structure is sensitive to the initial deflection caused by thermal postbuckling load. Increasing the temperature prior to mechanical loads enhances the snap-through intensity and also increases both the upper and lower limit loads. As shown, dynamic snap-through loads are lower than the static ones, however dynamic snap-through intensity is more than the static snap-though intensity. Furthermore, structural damping enhances the dynamic buckling loads of the plate and decreases the dynamic postbuckling deflection of the plate.     

铁路道床有砟-无砟过渡段动力特性的DEM-FEM耦合分析

针对铁路道床有砟-无砟过渡段的结构特点,采用离散元-有限元耦合模型分析散体道砟和无砟道床间过渡段的动力特性。散体道砟道床和无砟道床分别采用离散元方法 DEM和有限元方法 FEM模拟,而在过渡段将道砟颗粒嵌入无砟道床以增加道砟颗粒与无砟道床间的咬合力,并在离散元和有限元耦合区域实现了力学参数的传递。采用以上DEM-FEM耦合方法对有砟-无砟道床及其过渡段在列车荷载作用下的沉降过程进行了数值分析。计算结果表明,离散元方法中道砟颗粒间的力链呈现非对称梯形分布,其与有限元方法中的应力分布趋势一致;采用嵌入式道砟颗粒的方法可以增加有砟-无砟过渡段道砟间的咬合力,有效约束道砟颗粒的位移,减少有砟-无砟道床间的沉降差异。本文计算模型可以合理地分析有砟道床的力链分布以及无砟道床的应力分布,确定列车荷载下道床有砟-无砟过渡段的动力学行为。     

用均匀化理论分析蜂窝结构的等效弹性参数

在线弹性范围内,根据均匀化理论,并结合有限元方法推导出适用于二维周期性结构的均匀化的有限元格式（Homo FEM）,计算出不同相对密度下的规则蜂窝结构的等效弹性模量Ee和泊松比νe.同时,利用蜂窝结构的代表胞元模型,用常规的有限元方法（FEM）模拟计算出相应的等效弹性参数.最后将两种数值计算结果与己有的理论公式进行了比较和分析讨论.结果表明:在考察的相对密度全场范围内（0～0.4）,HOmO FEM得到的蜂窝结构的 Ee和νe 与 FEM使用平面实体单元模拟计算得到的结果一致吻合,反映出 Homo FEM数值方法的客观准确性和可行性.而 Gibson公式和 W-K得到的等效弹性模量值 Ee只是在较小相对密度的情况下（小于0.15）与数值计算结果吻合.当结构的相对密度较大时,必须考虑胞棱附近区域由应力集中导致的复杂的应力和应变分布的影响.     

固体介质中SPH方法的拉伸不稳定性问题研究进展

光滑粒子流体动力学法(smoothed particle hydrodynamics, SPH)是一种基于核估计的无网格Lagrange数值方法.它用粒子方程离散流体动力学的连续方程, 既可以处理有限元难于处理的大变形和严重扭曲问题, 又可以处理有限差分法不易处理的自由边界和材料界面的问题, 在固体力学中的冲击、爆炸和裂纹模拟中具有广阔的发展前景.但是, 该算法的拉伸不稳定性(tensile instability)问题是它在固体力学领域中应用的最大障碍.对SPH稳定性分析表明, 算法不稳定性的条件仅与应力状态和核函数的2阶导数有关.目前, 应力点法(stress points)、Lagrange核函数法、人工应力法(artificialstress)、修正光滑粒子法(corrective smoothed particle method, CSPM)和守恒光滑法(conservativesmoothing)以及其他一些方法成功地改善了SPH的拉伸不稳定性, 但是每一种方法都不能彻底解决SPH的拉伸不稳定性问题.本文介绍了SPH法的方程和Von Neumann稳定性分析的思想, 以及国内外在这几个方面的研究成果及其最新进展, 同时指出目前研究中存在的问题和研究的方向.      

Reduced-order finite element method based on POD for fractional Tricomi-type equation

The reduced-order finite element method (FEM) based on a proper orthogonal decomposition (POD) theory is applied to the time fractional Tricomi-type equation. The present method is an improvement on the general FEM. It can significantly save memory space and effectively relieve the computing load due to its reconstruction of POD basis functions. Furthermore, the reduced-order finite element (FE) scheme is shown to be unconditionally stable, and error estimation is derived in detail. Two numerical examples are presented to show the feasibility and effectiveness of the method for time fractional differential equations (FDEs).     

小净距公路隧道稳定性数值模拟分析

常(德)-吉(首)高速公路上的阿娜隧道是因隧道净距小,双洞间影响较大,故应用二维弹塑性有限元分析方法对其稳定性进行分析研究,计算了衬砌和围岩的内力与位移,进行了变形性状的分析,为高速公路类似的小净距隧道工程设计、施工提供一定的科学依据和参考价值。     

Coupling between finite element method and material point method for problems with extreme deformation

As a Lagrangian meshless method, the material point method (MPM) is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method (FEM) for small deformation problems. Therefore, an algorithm for the coupling of FEM and MPM is proposed to take advantages of both methods. Furthermore, a conversion scheme of elements to particles is developed. Hence, the material domain is firstly discretized by finite elements, and then the distorted elements are automatically converted into MPM particles to avoid element entanglement. The interaction between finite elements and MPM particles is implemented based on the background grid in MPM framework. Numerical results are in good agreement with experimental data and the efficiency of this method is higher than that of both FEM and MPM.     

Coupled MPS-modal superposition method for 2D nonlinear fluid-structure interaction problems with free surface

In this paper, a coupled MPS-modal superposition method is developed for 2D nonlinear fluid-structure interaction problems. In this method, the rigid-body and relatively small elastic deformation are coupled together, which considers the mutual effect between them. The elastic deformation of the structure is represented by a mode superposition formulation, which is more efficient compared with FEM, regardless of the size of the structure. For 2D cases, if the first three modes are chosen to represent the flexible deformation of the structure, it only results in a 6×6 matrix equation to be solved. For the fluid motion, the modified Moving Particle Semi-implicit (MPS) method, which significantly reduces the fluctuation of pressure calculation of the original MPS method, is used.Two nonlinear problems, i.e. breaking-water-dam impacting a floating beam and flexible wedge slamming into the water are simulated to demonstrate the performance of the developed method. The numerical simulations show that this coupling model is capable of providing stable results that are generally in good agreement with the available experimental data. For the highly nonlinear case with very large rigid motions, the mutual effect between elastic deformation and rigid motions could accumulate to a relatively remarkable level shown by the curves of trajectories or acceleration history of the body mass centre. This also indicates the importance of mutual effect to analyse highly nonlinear FSI problems with large rigid-body motions and relatively small flexible deformation.