炮口冲击式火炮后坐模拟试验的数值分析

为了获得炮口冲击试验装置的冲击参数,保证火炮后坐模拟试验的精度和有效性,以某型地面火 炮为对象,基于MSC.ADAMS软件,开发了虚拟的火炮后坐模拟试验平台,实现了火炮动态后坐、复进过程 的动力学模拟,开展了冲击参数的均匀试验设计。通过试验与数值模拟结果的比较,验证了冲击参数组合的 合理性和模拟试验的有效性。     

磁流体特性对磁流变火炮后坐阻尼器性能的影响

针对火炮后坐磁流变阻尼器的特点,以某单管25 mm火炮实验用磁流变阻尼器为研究对象,基于Herschel-Bulkley本构模型,建立了该炮磁流变后坐阻尼器的轴对称一维层流模型,获得了不同磁场作用下阻尼力随活塞速度的变化规律。运用火炮的后坐运动方程,计算了不同磁流体特性指数下,火炮的后坐位移和后坐速度。计算结果表明,磁流体特性指数的变化对磁流变后坐阻尼器的性能影响显著。     

降落伞膜索非线性有限元模拟

用膜单元和索单元模拟降落伞织物绳索系统,基于完全拉格朗日格式的非线性有限元方法编程计算降落伞的结构动力学特性。采用增量与迭代混合方法改善非线性计算的收敛特性并结合HHT隐式时间推进方法减小整体迭代计算量。使用修正应力应变张量导数的方法模拟膜单元单向应力状态并针对膜单元和索单元分别进行了非线性有限元计算验证。最后针对C-9型降落伞建立三维有限元模型,根据设定流速对伞衣施加均匀压强载荷,将模拟展开的结果与使用相同模型、不同方法商业软件的文献进行对比,显示了隐式非线性有限元方法模拟降落伞膜索系统大变形动力学的能力。     

阻尼夹层筋板结构有限元动力分析

在小变形线弹性理论基础上，导出了粘弹性阻尼夹层筋板壳结构的有限元动力方程，构造了壳单元、阻尼夹层壳单元和开口薄壁梁单元。采用这些单元可以很好地模拟具有不连续阻尼夹层处理的筋板壳结构。数值分析采用［３］中提出的两次渐近法。通过一块阻尼加筋板的特征分析和振动测量，证实了所建立的有限元两次渐近法的良好精度和广泛的工程适用性。     

A COMPLEX PARAMETER MODEL UPDATING METHOD FOR DAMPING CHARACTERISTICS OF FINITE ELEMENT MODELS

阻尼对于结构动力学响应具有重要的影响,但有限元模型一般很难对阻尼特性进行精确建模.基于实测频响函数,研究了一种有限元模型阻尼特性的复参数修正方法.以待修正区域各单元质量、刚度矩阵的比例修正系数为复修正参数,建立了单元矩阵比例修正的灵敏度方程直接算法,并对比分析了复修正参数与不同阻尼特性之间的数学关系.以六自由度集中参数模型和25杆平面桁架模型为例,验证了复参数修正方法在阻尼特性修正中的有效性.     

关于粘弹性结构和复合结构有限元动力学方程的探讨

对粘弹性结构有限元动力学方程和复合结构有限元动力学方程进行了探讨。对于粘弹性结构有限元动力学方程，研究了表示成MCK形式的可能性，并对几种典型模型进行了分析研究。根据对阻尼的假设，给出了考虑弹性材料阻尼及结构外部粘性阻尼的拉氏域内弹性－－粘弹性复合结构有限元动力学方程。     

有限元模型阻尼特性的复参数修正方法研究

阻尼对于结构动力学响应具有重要的影响,但有限元模型一般很难对阻尼特性进行精确建模.基于实测频响函数,研究了一种有限元模型阻尼特性的复参数修正方法.以待修正区域各单元质量、刚度矩阵的比例修正系数为复修正参数,建立了单元矩阵比例修正的灵敏度方程直接算法,并对比分析了复修正参数与不同阻尼特性之间的数学关系.以六自由度集中参数模型和25杆平面桁架模型为例,验证了复参数修正方法在阻尼特性修正中的有效性.     

考虑空气耦合的薄膜结构动力分析

假定空气是可压缩、无粘性和有势的,推导了薄膜结构的非线性刚度矩阵和空气与薄膜结构耦合作用的气动力表达式。根据薄膜结构的受力特点,建立了薄膜结构与空气耦合系统的运动方程。采用有限元方法对无限域的空气进行简化,建立了空气与薄膜结构的耦合有限元模型,模拟了空气对薄膜结构的附加质量和声致阻尼影响。考虑薄膜结构的几何非线性,对结构在空气中的自由振动和受迫简谐振动进行了系统分析,同时探讨了流体环境对薄膜结构动力响应的影响因素。     

一个统一单元矩阵计算的方法

1.引言有限元应用软件系统中,普遍采用有限元位移法。但单元矩阵(刚度矩阵、质量矩阵、阻尼矩阵、以至载荷列阵)等的计算,多数是逐一给出,部分是以元素簇给出的。本文将讨论不同类型元素的单元矩阵的统一算法和程序实施技巧,给出计算精密单元和几何非线性单元方便、有效的途径,达到以积累插值函数积累单元矩阵库的目的。     

张拉结构非线性分析的五节点等参单元

本文针对张拉结构的特点，提出了一种五节点等参数单元有限元模型，采用四次多项式作为位移插值函数及单元初始形状函数，并假定索是理想柔性的且满足虎克定律，基于修正的Ｌａｇｒａｎｇｉａｎ坐标描述法，建立了非线性有限元基本方程和切线刚度矩阵，利用Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ法进行了实例计算。结果表明：本文方法精度极高，可供大跨度索网，索穹顶，拉线塔等张拉结构分析，设计时采用。     

Post-divergence and post-flutter analysis of sub-tangentially loaded and damped Beck column using a nonlinear FE procedure

Both post-divergence and post-flutter behaviors of damped Beck columns subjected to a sub-tangentially follower force are rigorously explored using an exact co-rotational frame element. First a linear stability theory of the damped Beck column is summarized using a stability map. A geometrically nonlinear frame element based on the co-rotational formulation is then formulated including mass matrix, Rayleigh damping matrix, and load-correction stiffness matrix due to circulatory forces. The dynamic FE analysis is performed using Newmark integration method. Finally a Beck column model is parametrically analyzed in order to investigate non-linear stability characteristics of the internally and externally damped non-conservative system. In particular, interesting nonlinear behaviors of Beck column quite different from those predicted by the linear theory are reported through static and dynamic nonlinear analysis with variation of sub-tangentiality of the follower force.     

半平面问题远域辐射阻尼的时域边界元法模拟研究

辐射阻尼在岩石基坑爆破开挖、边坡稳定、结构抗震以及结构-地基动力相互作用等实际工程问题中具有重要意义.为了模拟半平面问题的远域辐射阻尼,以时域边界元法(TD-BEM)理论为基础,根据应力波在弹性介质中的传播特性,在时域内提出了一种新的单元,即自适应半无限边界单元,专门用于离散远域半无限边界.该单元外侧节点是一个始终处于...     

沉管隧道横向动力响应分析的多体动力学方法

沉管隧道横向地震响应分析是沉管隧道结构抗震设计的重点.根据沉管隧道的结构特点和接头剪力键的节点构造,建立了沉管隧道的多刚体-剪切铰-黏性阻尼铰多体动力学模型,并应用多体动力学中的离散时间传递矩阵法(multibody system discrete time transfer matriX method,MS—DT-TMM)建立了数学模型及表达式,同时编写Matlab程序进行求解实现.设定算例分别采用MS—DT—TMM和传统有限元法(finite element method,FEM)对沉管隧道横向地震响应进行对比分析,两者吻合较好,表明MS—DT-TMM应用于沉管隧道横向动力响应分析的可行性和有效性.     

简谐激励下黏弹性非线性能量阱的研究

黏弹性材料作为一种良好的减振材料,广泛应用于机械、航空和土木等领域.本文用黏弹性Maxwell器件代替传统非线性能量阱中的阻尼元件,提出一种新型的黏弹性非线性能量阱,并对该模型在简谐激励下的减振性能进行分析.首先,根据牛顿第二定律建立系统的动力学方程,采用谐波平衡法求解系统的幅频响应曲线,并利用MATLAB中的Runge-Kutta数值方法验证解析解的正确性,结果吻合良好.然后,分析黏弹性非线性能量阱的减振性能和参数的影响.最后,分析了不同质量比下非线性刚度比和阻尼比同时变化时减振效果的变化趋势,并讨论了黏弹性非线性能量阱的最佳取值范围.研究结果表明:主系统的最大振幅随着非线性刚度的增加先减小后增大;当参数选取恰当时,黏弹性非线性能量阱比传统非线性能量阱的减振效果更优;另外,随着质量比的增加,主系统最大振幅的最小值出现先减小后趋于不变的现象,且非线性刚度比和阻尼比的最佳取值范围有所增大.以上结论对黏弹性非线性能量阱的实际应用提供了一定的理论依据.     

基于转换矩阵的FEM/MLPG耦合算法

首次基于有限元的转换矩阵(TMF)和无网格的转换矩阵(TMM),提出有限单元法(FEM)和无网格局部彼得罗夫-伽辽金法(MLPG)的耦合算法。编制了相应算法的三维程序,计算分析了三维柱体的拉伸和弯曲问题,并将计算结果与ABAQUS软件计算结果以及理论解进行了比较。结果表明,本文给出的耦合算法计算精度高,收敛性好,可以用以模拟裂纹扩展等问题。     

Time-dependent hydroelastic analysis of cavitating propulsors

A 3-D potential-based boundary element method (BEM) is coupled with a 3-D finite element method (FEM) for the time-dependent hydroelastic analysis of cavitating propulsors. The BEM is applied to evaluate the moving cavity boundaries and fluctuating pressures, as well as the added mass and hydrodynamic damping matrices. The FEM is applied to analyze the dynamic blade deformations and stresses due to pressure fluctuations and centrifugal forces. The added mass and hydrodynamic damping matrices are superimposed onto the structural mass and damping matrices, respectively, to account for the effect of fluid–structure interaction. The problem is solved in the time-domain using an implicit time integration scheme. An overview of the formulation for both the BEM and FEM is presented, as well as the BEM/FEM coupling algorithm. Numerical and experiment validation studies are shown. The effects of fluid–structure interaction on the propeller performance are discussed.     

Dynamics of a spacecraft with large flexible appendage constrained by multi-strut passive damper

This paper is concerned with the dynamics of a spacecraft with multi-strut passive damper for large flexible appendage.The damper platform is connected to the spacecraft by a spheric hinge,multiple damping struts and a rigid strut.The damping struts provide damping forces while the rigid strut produces a motion constraint of the multibody system.The exact nonlinear dynamical equations in reducedorder form are firstly derived by using Kane’s equation in matrix form.Based on the assumptions of small velocity and small displacement,the nonlinear equations are reduced to a set of linear second-order differential equations in terms of independent generalized displacements with constant stiffness matrix and damping matrix related to the damping strut parameters.Numerical simulation results demonstrate the damping effectiveness of the damper for both the motion of the spacecraft and the vibration of the flexible appendage,and verify the accuracy of the linear equations against the exact nonlinear ones.     

Robust tracking control method based on composite nonlinear feedback technique for linear systems with time-varying uncertain parameters and disturbances

In this paper, the composite nonlinear feedback control method is considered for robust tracking and model following of uncertain linear systems. The control law guarantees that the tracking error decreases asymptotically to zero in the presence of time varying uncertain parameters and disturbances. For performance improvement of the dynamical system, the proposed robust tracking controller consists of linear and nonlinear feedback parts without any switching element. The linear feedback law is designed to allow the closed loop system have a small damping ratio and a quick response while the nonlinear feedback law increases the damping ratio of the system as the system output approaches the output of the reference model. A new collection of different nonlinear functions used in the control law are offered to improve the reference tracking performance of the system. The proposed robust tracking controller improves the transient performance and steady state accuracy simultaneously. Finally, the simulations are provided to verify the theoretical results.     

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.     

Time domain analysis of second-order wave diffraction by an array of vertical cylinders

A time domain method is employed to analyze interactions of water waves and a group or an array of cylinders. The nonlinear free-surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a finite element method (FEM). The mesh required is generated based on a 2-D unstructured grid on a horizontal plane and its extension in the vertical direction. The matrix equation of the FEM is solved through an iteration and the initial solution is obtained from the result at the previous time step. The radiation condition is imposed through a combination of the damping zone method and the Sommerfeld–Orlanski equation. Results for various configurations, including two cylinders, four cylinders, an array of 10 cylinders and two arrays of eight cylinders are provided to show the effect of the interaction and their behaviour near the trapped mode.