索网天线的参变量变分及非线性有限元方法

针对大型周边桁架式索网天线由拉索拉压模量不同引起的本构非线性和结构大变形引起的几何非线性问题,给出了基于参变量变分原理的几何非线性有限元方法. 首先针对含预应力索单元拉压模量不同分段描述的本构关系,通过引入参变量,导出了基于参变量及其互补方程的统一描述形式,避免了传统算法需要根据当前变形对索单元张紧/松弛状态的预测,提高了算法收敛性. 然后利用拉格朗日应变描述索网天线结构大变形问题,结合几何非线性有限元法,建立了基于参变量的非线性平衡方程和线性互补方程;并给出了牛顿-拉斐逊迭代法与莱姆算法相结合的求解算法. 数值算例验证了本文提出的算法比传统算法具有更稳定的收敛性和更高的求解精度,特别适合于大型索网天线结构的高精度变形分析和预测.      

浅曲梁大挠度分析的余能有限元方法

以高玉臣提出的弹性大变形余能原理为基础,利用Lagrange乘子,放松平衡方程和力边界条件对余能泛函的约束,推导出广义的余能原理.根据极分解定理,将变形分为刚性转动和纯变形两部分,则余能也包含相应的两部分,一部分与刚性转动有关,而另一部分与纯变形有关.使用线弹性本构关系,建立了可用于几何非线性计算的有限元模型.应用更新的Lagrange列式法,给出了增量形式的有限元公式.数值计算结果表明,该方法可用于浅曲粱的几何大变形计算.     

NONLINEAR LARGE DEFORMATION BENDING PROPERTIES OF THE TRACHEA1)

研究了人工气管和生物气管拉伸和弯曲的力学性能及仿真方法。根据气管拉伸实验数据,拟合建立了人工气管的线弹性和生物气管的非线性幂函数本构模型,其与实验的最大误差为18.8%。根据非线性弯曲理论,建立了气管大变形的弯曲变形方程,并得到了其解析解,计算的弯曲载荷与实验的最大误差为8.7%。基于显式有限元法提出了气管弯曲变形的有限元仿真方法,犬气管的弯曲角度的有限元计算结果与实验的误差为3%。     

基面力概念在几何非线性余能有限元中的应用

以基面力为基本未知量描述一个弹性系统的应力状态并表征单元的余能，将大变形的余能分 解为变形余能部分和转动余能部分，采用Lagrange乘子法放松单元的平衡方程，利 用已有的弹性大变形余能原理建立了一种几何非线性显式有限元模型，编制了相应的几何非线 性余能原理有限元程序. 数值算例表明：该方法具有较好的收敛性和计算精度，可进行大 载荷步的大位移、大转动计算.     

基于余能原理的有限变形问题有限元列式

利用基面力概念,推导了一种基于余能原理的有限变形问题显式有限元列式,可应用于结构的大位移、大转动问题。以基面力为状态变量来表达单元的余能,将有限变形情况下的单元余能分解为变形余能部分和转动余能部分,利用Lagrange乘子法推导出余能原理有限元的控制方程,编制了相应的非线性有限元程序。通过算例分析,说明该列式和程序的可靠性和精确性。     

NONLINEAR VIBRATIONS OF AXIALLY ACCELERATING VISCOELASTIC TIMOSHENKO BEAMS

研究了轴向加速黏弹性Timoshenko梁的非线性参数振动。参数激励是由径向变化张力和轴向速度波动引起的。引入了取决于轴向加速度的径向变化张力,同时还考虑了有限支撑刚度对张力的影响。应用广义哈密尔顿原理建立了Timoshenko梁耦合平面运动的控制方程和相关的边界条件。黏弹性本构关系采用Kelvin模型并引入物质时间导数。耦合方程简化为具有随时间和空间变化系数的积分-偏微分型非线性方程。采用直接多尺度法分析了Timoshenko梁的组合参数共振。根据可解性条件得到了Timoshenko梁的稳态响应,并应用Routh-Hurvitz判据确定了稳态响应的稳定性。最后通过一系列数值例子描述了黏弹性系数、平均轴向速度、剪切变形系数、转动惯量系数、速度脉动幅值、有限支撑刚度参数以及非线性系数对稳态响应的影响。     

Nonlinear numerical simulation method for galloping of iced conductor

Based on the principle of virtual work, an updated Lagrangian finite element formulation for the geometrical large deformation analysis of galloping of the iced conductor in an overhead transmission line is developed. In numerical simulation, a three-node isoparametric cable element with three translational and one torsional degrees-of-freedom at each node is used to discretize the transmission line. The nonlinear dynamic system equation is solved with the Newmark time integration method and the Newton-Raphson nonlinear iteration. Numerical examples demonstrate the efficiency of the presented method and the developed finite element program. A new possible galloping mode, which may reflect the saturation phenomenon of a nonlinear dynamic system, is discovered under the condition that the lowest order of vertical natural frequency of the transmission line is approximately two times of the horizontal one.     

含随机参数的非线性黏弹性有限元计算

进行了粗粒土与结构接触面单调和循环加载试验，基于宏细观测量结果, 扩展了 损伤概念以 描述该类接触面在受载过程中的物态演化, 及由于物态演化导致的力学特性从初始状态到最终 稳定状态的连续变化过程. 揭示了接触面损伤的细观物理基础主要是接触面内土的颗粒破碎 和剪切压密这两种物态演化；指出接触面的剪胀体应变可以划分为可逆性和不可逆性剪胀体 应变两部分，其中不可逆性剪胀体应变可作为接触面损伤发展的宏观量度，因此其归一化 形式可作为一种损伤因子的定义；提出了建立粗粒土与结构接触面一种损伤本构关系的基本思路.     

索网天线找形分析的参变量变分有限元法

针对柔性空间索网天线的非线性特点,建立了基于参变量变分描述索网拉压非线性和共旋列式描述几何非线性的有限元控制方程,应用Lernke与改进牛顿法进行求解.通过对索网预张力平衡计算,证明改进牛顿法比Newton-Raphson法具有更强的收敛能力.进一步将力密度法迭代原理与有限元法结合应用于索网天线的非线性找形分析中,获得了理想的索网构型.本文的索网找形方法可广泛应用于空间索网天线结构的设计.     

Semiclosed-form solution for static nonlinear analysis of extensible cables

The analysis of a cable deforming under statically applied body forces leads to a series of differential equations that are inherently nonlinear due to coupling of geometry and loading. These equations have been solved by nonlinear finite element technique, finite difference formulation, or a two-step shooting method. However, a semiclosed-form solution is presented herein in which the governing differential equations are integrated to yield exact expressions for the unknown cable forces and the parameters defining the geometry of the loaded cable. The formulation presented herein is capable of taking care of the large sag and deformation in the cable without employing any numerical solution procedure.     

大型索网天线机电耦合动力学建模与主动变形控制

索网天线因其折叠收纳比高、质量轻以及口径大等优点,成为众多空间任务的首选天线形式,但反射面精度较低限制了大型索网天线的在轨性能指标。在索网天线竖向张力索中加入压电陶瓷堆(PZT)作动器,通过主动控制提高反射面的形面精度。首先,将张紧的索网天线分为含有作动器的主动索单元和不含作动器的被动索单元,通过压电本构方程、几何方程和哈密顿原理推导出主动索单元的动力学模型,并与被动索单元一起组装得到索网天线的机电耦合动力学模型。然后,根据天线反射面的初始形面误差和所建立的天线机电耦合动力学模型,基于模型预测控制(MPC)方法计算出各作动器的控制输入电压曲线。最后,以一个10m口径的索网天线为算例,验证本文所建立的模型和控制方法的有效性。     

含芯拧绞绳非线性弯曲动力学特性分析与研究

建立细长缆索大柔性多体动力学模型时,现实存在的复杂捻制几何构型多不予考虑,而是将柔索简化为材料均匀梁进行描述,致使运动仿真模型与物理实际存在一定差距. 为此,研究一种典型非线性拧绞绳股的大变形等效动力学建模方法,考虑准静态与大范围运动情况下绳股内的线接触,计算了受摩擦力及弯曲曲率影响的绳股可变弯曲刚度,通过等效梁模型避免了绳股精细建模时的大规模计算消耗. 基于连续介质力学与绝对节点坐标方法,建立了拧绞绳惯性广义坐标下的多柔体动力学模型. 为了验证等效模型的可行性,与基于有限段方法建立的精细模型进行对比仿真分析,通过位形验证了等效模型的精度. 进一步地,根据力载作用下的准静态构型,研究了特定构型绳股弯曲刚度沿轴向的分布规律;通过自重力下一端固定柔性绳摆自由运动仿真并与传统均匀梁模型相比,研究了模型弯曲特性的差异. 最后,根据能量守恒原理分析了摩擦耗散系统内各种能量间的相互转化. 拧绞绳大变形等效动力学模型能够提高绳索动力系统运动预测的仿真计算效率,还能为钢丝绳参数与构型设计提供依据.      

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

大变形下初始斜交异性本构方程

采用材料主轴法,建立了初始斜交异性材料在变形构形（Euler描述）下的斜交异性本构方程,以及在初始构形（Lagrange描述）下的形式。具体给出了斜交异性线弹性材料方程的显式,它在Lagrange描述下形式简洁,可方便地用于有限元计算。文中指出,在变形构形下是线弹性的材料,在Lagrange描述下其本构方程一般已成为非线性,我们称之为本构转换非线性。这种非线性在实际的有限元计算中还未引起重视。为理论简明,本构方程是对二维给出的。     

NONLINEAR DYNAMIC ANALYSIS OF A LAMINATED HYBRID COMPOSITE PLATE SUBJECTED TO TIME-DEPENDENT EXTERNAL PULSES

Nonlinear dynamic responses of a laminated hybrid composite plate subjected to time-dependent pulses are investigated. Dynamic equations of the plate are derived by the use of the virtual work principle. The geometric nonlinearity effects are taken into account with the von Kármán large deflection theory of thin plates. Approximate solutions for a clamped plate are assumed for the space domain. The single term approximation functions are selected by considering the nonlinear static deformation of plate obtained using the finite element method. The Galerkin Method is used to obtain the nonlinear differential equations in the time domain and a MATLAB software code is written to solve nonlinear coupled equations by using the Newmark Method. The results of approximate-numerical analysis are obtained and compared with the finite element results. Transient loading conditions considered include blast, sine, rectangular, and triangular pulses. A parametric study is conducted considering the effects of peak pressure, aspect ratio, fiber orientation and thicknesses.     

AN APPLICATION OF THE CBS SCHEME IN THE FLUID-MEMBRANE INTERACTION

提出了一种不可压缩流体与弹性薄膜耦合问题的特征线分裂有限元解法. 首先, 给出了流场和结构的控制方程. 然后, 对流场、结构以及流固耦合的具体求解过程进行了描述. 其中, 流场求解采用改进特征线分裂方法和双时间步方法相结合的隐式求解方式, 并利用艾特肯加速法对每个时间步的迭代收敛过程进行了加速处理;结构部分的空间离散和时间积分分别采用伽辽金有限元方法和广义方法, 并通过牛顿迭代法对所得非线性代数方程组进行了求解;流场网格的更新采用弹簧近似法;流场、结构两求解模块之间采用松耦合方式.最后, 采用该方法对具有弹性底面的方腔顶盖驱动流问题进行了求解, 验证了算法的准确性和稳定性.此外, 计算结果表明艾特肯加速法可以显著地提高双时间步方法迭代求解过程的收敛速度.     

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

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

Base force element method of complementary energy principle for large rotation problems

Using the concept of the base forces, a new finite element method （base force element method, BFEM） based on the complementary energy principle is presented for accurate modeling of structures with large displacements and large rotations. First, the complementary energy of an element is described by taking the base forces as state variables, and is then separated into deformation and rotation parts for the case of large deformation. Second, the control equations of the BFEM based on the complementary energy principle are derived using the Lagrange multiplier method. Nonlinear procedure of the BFEM is then developed. Finally, several examples are analyzed to illustrate the reliability and accuracy of the BFEM.     

A combined parametric quadratic programming and precise integration method based dynamic analysis of elastic-plastic hardening/softening problems

The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems. The project supported by the National Key Basic Research Special Foundation (G1999032805), the National Natural Science Foundation of China (19872016, 50178016, 19832010) and the Foundation for University Key Teacher by the Ministry of Education of China     

有初始几何缺陷的壳体弹塑性大变形分析

本文在等温小变形弹塑性内时本构方程偏量形式的基础上，导出了适用于大位移、大转动、小应变分析的弹塑性内时本构方程，进一步推导出了带有初始几何缺陷的非线性弹塑性问题的有限元方程，可用于分析缺陷对结构非线性弹塑性反应的影响，也可用于带缺陷的非线性问题求解及稳定性分析．