组合结构总体节点参数自动标定的基本定理

在有限元位移法中,单元矩阵(刚度矩阵、几何矩阵、一致质量矩阵、载荷列阵等)往往是在单元局部坐标下计算的(以下简称局部方位);结构整体矩阵的形成,单元应力和变形的计算,都要求确定总体节点参数(自由度)和单 ...     

充满液体弹塑性圆柱壳轴向冲击屈曲的数值模拟

基于对充满液体的弹塑性圆柱壳轴向冲击屈曲的实验观察 ,建立了相应的计算模型 ,借用 LS-DYNA软件包 ,在 HP-C36 0工作站上完成了数值模拟 ,得到了动画显示 ,给出的屈曲模态及一些特征量的时程曲线与实验结果相当一致。文中对壳的几何、物理参数对屈曲性能的影响进行了较详细的讨论     

基于因素敏感性分析的人工排土场高边坡稳定性研究

在考虑人工排土场松散高边坡基本特征、变形破坏方式等的基础上,采用因素敏感性分析法,计算讨论了影响排土场高边坡稳定的主要因素--坡体几何形态特征、日常爆破震动、水等因素的变化对其稳定性的影响程度,确定了坡体几何形态特征因素为影响其稳定性的关键因素以及在确保排土场高边坡稳定条件下的单台阶合理坡角、坡高、预留宽度等人工堆积参数及最大堆积高度。     

基于区间分析的结构非概率可靠性模型

本文用非概率的凸集模型模拟结构的不确定性,将结构的不确定参数描述为区间变量,基于区间分析,提出了一种新的非概率可靠性度量体系分析方法,从物理、几何意义等方面解释了文中理论的合理性,其计算方法简便,衫,给出了算例分析。     

结点位移计算的一种简单方法

 推导了一种结点位移计算的简单方法,在平面桁架结点位移的计算和静不定杆系结构的变形几何关系的确定中应用方便.     

简支梁大变形的优化算法

提出一种求解几何非线性问题的优化算法,并研究了简支梁的几何非线性大变形问题。首先取简支梁大变形后的平衡状态为研究对象,分别创建它的变分模型和微分模型;然后基于微分模型,通过动坐标的迭代关系式求得微段端点坐标,构建微段端点未知坐标的目标函数;最后确定简支梁几何非线性大变形的最优化问题,并编制相应优化程序进行求解。通过分析典型算例,并同有限元方法的计算结果相比较,表明提出的优化算法在求解强几何非线性大变形问题中的正确性,为处理几何非线性大变形问题提供了一种新方法。     

可重构力学超材料的设计与波动特性研究

力学超材料中的弯曲梁双稳态结构由于其主动调控性强且调控精度高等优点近年来受到广泛关注. 文章利用中心受压弯曲梁的不稳定性设计了六角型双稳态结构, 首先建立了等效弯曲梁模型, 基于梁变形微分方程及能量最低原理探明了结构双稳态特性的产生基理, 之后利用有限元数值计算研究了结构几何参数对其整体力学性能的影响, 分别得到了具备自恢复及双稳态性能的结构几何参数范围, 绘制了几何参数与力学性能之间的相图. 同时, 可重构结构的可控变形能力有助于调整整体的色散特性, 利用数值仿真研究了具备双稳态特性的结构在拉伸和压缩两种构型下的色散关系, 对比分析了不同结构几何参数及构型转变对结构产生的带隙位置及范围的影响, 之后对由不同构型单胞组成的周期性结构进行了频响分析来验证带隙计算的准确性. 通过六角型可重构结构的力学特性、色散特性研究及频响分析表明可以通过结构几何参数的设计实现对结构整体性能的主动调控, 为可逆向设计的弹性波超材料结构研究分析提供了一条可靠路径.      

井孔约束下水平井段钻柱的浑沌运动

考虑了钻柱即时构形中轴向力对钻柱弯曲的影响,计及由于弯曲变形而产生的轴向附加力,得到了水平井段钻柱在周期性波动钻压激励下的非线性参数激励振动系统,导出了根据钻柱的物理和几何性质以及井孔尺寸计算系统参数的公式,选通过Calerkin法对控制方程在空间域内加权消残,然后在时间时间域内用Melnikov-Holmes方法得到了钻柱可能发生浑沌振动的参数激励的阈值,文中结果可能对钻井轨迹的控制具有一定的指导意义。     

基于一阶剪切变形理论的变曲率曲梁的几何非线性方程

基于一阶剪切变形理论和轴线可伸长的精确几何非线性理论,推导了变曲率曲梁在热机载荷共同作用下的几何非线性控制方程。通过引入轴线伸长率,变形后的轴线弧长被当作基本未知量之一,基本未知量均被表示为变形前的轴线坐标的函数,使问题的求解区间仍为未变形时的曲梁轴线长度;给出了在给定曲梁轴线参数方程时,利用本文控制方程进行几何分析所需的初始曲率、变形前曲梁几何关系的数学表达式;介绍了几种常见的曲梁边界条件。所给数学模型可为轴线可伸长的变曲率曲梁的几何非线性分析和计算提供理论参考。     

连续油管井下力学行为模拟实验研究

连续油管作业系统在井下受力复杂,借助相关实验研究可为连续油管实际入井深度、连续油管屈曲、锁死及解除锁死等问题提供解决思路. 根据相关理论计算结果确定了合适的实验系统几何参数. 设计了实验台架,进行了不同规格模拟连续油管和不同弯曲段长度实验. 分析了屈曲、螺旋屈曲锁死、首端屈服等实验现象,提出了螺旋屈曲锁死的条件. 实验结果与理论计算结果吻合较好. 该实验为连续油管井下力学行为分析提供了依据.     

高速碰撞中的有限元方法及其应用

讨论了连续介质力学守恒方程的有限元离散方法,针对高速碰撞的物理特点,提出了数值计算中所应重视的问题及其有效的解决办法。文中给出的典型算例说明本方法对各类高速碰撞问题的数值模拟计算是十分有用和成功的。     

THE PARAMETRIC VARATIONAL PRINCIPLE AND NON-LINEAR FINITE ELEMENT METHOD FOR ANALYSIS OF ASTROMESH ANTENNA STRUCTURES

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

黏弹流体流动的数值模拟研究进展

综述了黏弹流体流动数值模拟的研究进展,突出介绍近十年来有限元法在黏弹流体流动数值模拟研究中取得的成果,通过动量方程的适当变形和本构方程离散权函数的合理选择,可以显著增强数值计算的稳定性。得到较高Weissenberg数下的解,同时文中对黏弹流体流动数值模拟中本构方程的应用、非等温情况和三维空间下的研究进行了介绍。     

修正金属本构模型在超高速撞击模拟中的应用

为更加准确地计算93钨合金弹超高速撞击Q345钢板问题,构建了修正的金属本构模型。引入GRAY三相物态方程描述材料相态变化,采用Johnson-Cook强度模型描述撞击后期材料的力学行为。结合封加波损伤演化模型以及Johnson-Cook失效模型描述不同应力三轴度下材料的拉伸、剪切失效行为;引入曹祥提出的断裂演化模型,描述材料失效后应力归零的过程。通过对比超高速撞击数值模拟结果与实验结果,验证了本构模型的适用性,并进一步分析了典型弹靶撞击条件下破片群的空间分布特征。研究结果表明:基于修正金属本构模型获得的超高速撞击靶板穿孔直径、弹体侵蚀长度、破片群扩展速度结果与实验结果一致;GRAY三相物态方程能够相对准确地给出弹体撞击首层靶板以及剩余弹体、破片群撞击第2层靶板时弹靶材料的熔化情况;封加波损伤演化模型能够准确判断超高速撞击过程中靶板是否产生层裂破坏;综合封加波损伤演化模型、Johnson-Cook失效模型以及曹祥提出的断裂演化模型后,数值模拟获得的破片群撞击后效靶板的穿孔面积与累积数量的统计曲线结果与实验结果一致;获得了典型条件下的柱形93钨弹体超高速撞击Q345靶板破片群空间分布结果,破片群的前端具有较高的质量、轴向动量以及横向动量（绝对值）。     

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

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

Hybrid analysis of plate problems

For the problem of plates on yielding subgrades, with deflections larger than allowed according to the classical theory of thin plates, a geometric nonlinear theory has been introduced, which led to two simultaneous differential equations depending on one another. These constitutive equations cannot be solved without further information about the response of the subgrade. Therefore, and hybrid technique has been applied. By means of moiré and photoelastic methods, experimental data have been obtained. These are introduced into the constitutive equations, which then can be solved without any mathematical difficulties.     

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

在等温小变形弹塑性内时本构方程偏量形式的基础上，导出了适用于大位移、小应变分析的弹塑性内时本构方程。并导出了带有初始几何缺陷的非线性弹塑性问题的有限元方程。文中给出的算例表明本方法是可行有效的。     

光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用

扼要讨论了光滑粒子法的离散思想,充分利用光滑粒子法和有限元方法各自的优点,提出了一种初始时刻用有限元建模,计算过程中大变形单元自动转换为光滑粒子的耦合算法。高速碰撞的系列算例说明,耦合算法不但适宜于计算大变形冲击动力学问题,而且由于集两种方法的优点于一身,可以更高效地模拟一些高速碰撞问题,提高计算效率。     

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.     

形状记忆合金管接头空间轴对称有限元分析

本文采用形状记忆合金（ＳＭＡ）的三维本构方程和有限变形理论，考虑拉、压不同应力状态对相变点移动的规律，编制了ＳＭＡ轴对称大变形的有限元程序，与单向拉伸下解析所得的应力、应变曲线相比，证实程序的正确性．文末计算一ＳＭＡ管接头，并指出按空间轴对称计算的必要性．