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

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

DYNAMIC ANALYSIS OF SPINNING DEPLOYMENT OF A SOLAR SAIL COMPOSED OF VISCOELASTIC MEMBRANES

近年来, 可用于航天器推进的太阳帆自旋展开技术引起人们广泛关注. 这类太阳帆可视为由中心旋转毂轮、若干柔性绳索、太阳帆薄膜和集中质量等组成的刚柔耦合多体系统.为了对系统中的太阳帆薄膜进行建模, 提出了基于绝对节点坐标方法描述的黏弹性薄板单元, 并对其有效性进行了验证.针对简化的"IKAROS"自旋展开太阳帆系统, 采用结合自然坐标方法与绝对节点坐标方法的绝对坐标方法对其进行建模, 采用广义-α方法对大规模系统动力学方程进行求解.研究了黏弹性太阳帆薄膜自旋展开过程的动力学特性, 讨论了薄膜的黏弹性阻尼对自旋展开过程的影响规律.     

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

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

THE APPLICATION OF VARIATIONAL MULTISCALE METHOD ON A ONE-DIMENSIONAL MECHANICAL MODEL

将变分多尺度方法应用于一维缆索模型,导出受力缆索的宏观有限元模型并求得细观位移解析解,总结出变分多尺度方法应用于具体模型的关键点和缺陷. 假定刚度为常值,数值模拟一定边界和受力下的缆索,得到宏观和细观位移. 将细观与宏观位移叠加,相比于精确位移得出：细观位移可视为常规有限元模型的后验误差. 变分多尺度方法在一维力学模型中的成功应用,推进了其实用性,为其在更多力学及工程问题中的运用和发展提供了参考.     

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 2D ELASTICITY

针对边界元法中高阶单元中几乎奇异积分计算难题,解剖了二维边界元法高阶单元的几何特征,定义源点相对高阶单元的接近度。将高阶单元上奇异积分核函数用近似奇异函数逼近,从而分离出积分核中主导的奇异函数部分,其奇异积分核分解为规则核函 数和奇异核函数两项积分之和。规则核函数用常规高斯数值积分,再对奇异核函数积分导出解析公式,从而建立了一种新的半解析法,用于高阶边界单元上几乎强奇异和超奇异积分计算。给出3个算例,采用边界元法高阶单元的半解析法计算了弹性力学薄体结构和近边界点位移/应力,并与线性边界元正则化算法结果作了比较,结果表明提出的二次元的半解析算法更加有效。特别是分析薄体结构,采用正则化算法的线性边界元分析比有限元有显著优势,而用提出的二次边界元半解析算法分析比其线性元的有效接近度又减小了4个量级。     

NONLINEAR LARGE DEFORMATION BENDING PROPERTIES OF THE TRACHEA1)

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

PREDICTIONS OF EFFECTIVE OUT-PLANE SHEAR MODULUS AND SIZE EFFECT OF HEXAGONAL HONEYCOMB

给出了预测六边形蜂窝材料等效剪切模量及其尺寸效应的圆筒扭转力学模型和扭转能量法,建立了等效面外剪切模量G₁₃相对于材料体分比ν、周向单胞数n、圆筒半径r和单胞层数参数m变化的解析表达式;同时将扭转能量法、有限元数值模拟计算和G-A经典细观力学方法进行了比较,从理论上揭示并验证了尺寸效应的存在性. 结果表明,当蜂窝体胞尺寸相对结构尺寸无穷小时,预测结果趋近于细观力学方法的结果. 此外,利用周期性蜂窝材料的结构对称特性,使用体胞子结构有限元计算模型进行等效面外剪切模量及其尺寸效应的预测,在不影响计算结果的前提下极大地提高了计算效率.     

BAR SIMULATION METHOD FOR BENDING NORMAL STRESS OF I-BEAM IN CONSIDERATION OF SHEAR EFFECT 1)

随着工字形短深梁和宽翼缘梁结构的发展,截面非线性剪切变形对弯曲应力的影响愈加突出,导致传统设计中所采用的初等梁理论计算结果误差较大,不再适用。本文基于比拟杆法综合考虑剪切效应,推导出工字形梁横力弯曲应力解析计算公式,并与有限元及现有解析计算方法进行对比分析。结果表明：当跨高比较小,翼缘腹板面积比较大时剪切效应对弯曲变形有显著影响。同时相比于现有解析方法,本文计算结果精度较高且适用范围更广,可用于梁结构设计。     

能量有限元预示复合材料结构动响应研究进展

复合材料结构高频动响应预示是飞行器等结构设计中的重要研究内容之一.为了探讨精确预示复合材料结构高频动响应方法,分析比较了目前较为通用的3种动力学响应预示方法,指出能量有限元法最适合求解具有各向异性特征的复合材料结构高频动响应问题.紧接着概述了国内外关于能量有限元方法和该方法在复合材料结构高频动响应预示方面的研究进展.在此基础上分析了能量有限元法在预示复合材料结构高频动响应问题中尚待深入研究的问题.     

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

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

一种XFEM断裂分析的裂尖单元新型改进函数

提出了一种适用于裂尖改进单元的新型改进函数, 基于三角变换的方法, 保留裂纹尖端场的应力奇异性和裂纹上、下表面的位移不连续性, 将常规扩展有限元法裂尖改进单元的4 项改进函数缩减为2 项, 裂尖改进单元的结点由常规的8 个改进自由度减少为4 个. 采用2 个正交的水平集函数表征材料内部裂纹面, 详细阐述了改进单元类型的判别方法, 给出一种改进单元的分区域积分方案. 最后, 若干断裂力学问题经典算例的数值计算结果表明:建议的裂尖改进函数具有较高的数值精度, 该方法是十分有效的.     

斜圈弹簧力学性能的有限元分析

斜圈弹簧在医疗器械、高压开关、汽车工业等领域应用广泛. 但其结构特殊, 受力复杂, 目前尚无有效的理论计算方法得出其载荷位移曲线, 只能通过数值仿真的方法研究其力学性能. 利用有限元软件ANSYS参数化建立环形斜圈弹簧有限元模型, 通过对模型的分析求解, 得出了其主要参数线径、节距、倾斜角与斜圈弹簧力学性能之间的关系. 研究了斜圈弹簧工作状态下的受力变形情况, 计算分析结果为斜圈弹簧的优化设计和合理应用提供了理论依据.     

基于应变梯度中厚板单元的石墨烯振动研究

基于应变梯度理论建立了单层石墨烯等效明德林(Mindlin) 板动力学方程,推导了四边简支明德林中厚板自由振动固有频率的解析解. 提出了一种考虑应变梯度的4 节点36 自由度明德林板单元,利用虚功原理建立了单层石墨烯的等效非局部板有限元模型. 通过对石墨烯振动问题的研究,验证了应变梯度有限元计算结果的收敛性. 运用该有限元法研究了尺寸、振动模态阶数以及非局部参数对石墨烯振动特性的影响. 研究表明,这种单元能够较好地适用于研究考虑复杂边界条件石墨烯的尺度效应问题. 基于应变梯度理论的明德林板所获得石墨烯的固有频率小于基于经典明德林板理论得到的结果. 尺寸较小、模态阶数较高的石墨烯振动尺度效应更加明显. 无论采用应变梯度理论还是经典弹性本构关系,考虑一阶剪切变形的明德林板模型预测的固有频率低于基尔霍夫(Kirchho) 板所预测的固有频率.      

NUMERICAL SIMULATION FOR NAVIER-STOKES EQUATIONS BY PROJECTION/CHARACTERISTIC-BASED OPERATOR-SPLITTING FINITE ELEMENT METHOD

提出了一种求解非定常不可压缩纳维-斯托克斯方程（N-S方程）的新型有限元法：基于投影法的特征线算子分裂有限元法.在每一个时间层上将N-S方程分裂成扩散项、对流项、压力修正项.对流项采用多步显式格式,且在每一个对流子时间步内采用更加精确的显式特征线-伽辽金法进行时间离散,空间离散采用标准伽辽金法.应用此算法对平面泊肃叶流、方腔流和圆柱绕流进行数值模拟,所得结果与基准解符合良好.尤其对于Re=10000的方腔流,给出了方腔中分离涡发展和运动的计算结果,并发现在该雷诺数下存在周期解,表明该算法能较好地模拟流体流动中的小尺度物理量以及流场中分离涡的运动.     

ANALYTICAL SOLUTION OF THE FREE VIBRATION OF VISCOELASTIC SANDWICH ANNULAR PLATE

基于薄板小挠度理论和Kelvin-Voigt 黏弹性本构方程,建立了黏弹性夹层环形薄板振动控制方程.采用分离变量法计算了内边固支、外边自由黏弹性夹层环形薄板的固有频率和振型,并与有限元计算结果进行比较. 分别讨论了夹心层比和内外半径比对固有频率及衰减系数的影响. 研究表明:系统频率随夹心层厚度增大,先增大后减小,而衰减系数一直增大;系统频率和衰减系数随内外半径比增大而增大.     

ANALYSIS OF SHEAR LAG EFFECT OF TWIN-CELL BOX GIRDERS

针对单箱双室箱梁,考虑各翼板间剪力滞翘曲的差异,并结合全截面轴力自平衡条件,定义了箱梁各翼板的剪滞翘曲位移函数. 利用最小势能原理,建立了双室箱梁考虑剪力滞效应的控制微分方程. 对一典型的单箱双室简支箱梁,利用空间板壳数值方法和本文解析解方法,研究了满跨均布载荷和跨中集中力作用下截面的剪力滞分布规律. 结果表明,本文提出的剪力滞翘曲位移模式能够反映双室箱梁各翼板间剪力滞翘曲的差异,本文解析解与有限元数值解吻合良好. 双室箱梁中腹板部位顶、底板处的剪力滞效应与边腹板部位有一定差异,对算例结构,中腹板部位的顶、底板应力小于边腹板部位的应力.     

EVALUATION ON STRESS INTENSITY FACTOR OF CRACK UNDER DYNAMIC LOAD USING NUMERICAL MANIFOLD METHOD

相较于传统有限元,数值流形方法（numerical manifold method, NMM） 的一个显著优点是在处理裂纹问题时网格无需与裂纹重合,这就方便了岩体破坏过程的模拟. 基于包含裂尖增强函数的NMM,采用Newmark 隐式动力学算法进行时间积分,重点研究了动力载荷条件下裂纹动态应力强度因子（dynamic stress intensity factor,DSIF） 的求解方法. 针对典型的线弹性动力裂纹问题,给出了NMM 的数值算例. 结果表明NMM 能够准确计算动载荷条件下裂纹的DSIF,并且具有较好的收敛性.     

A p-TYPE SUPERCONVERGENT RECOVERY METHOD FOR FE ELASTIC STABILITY ANALYSIS OF EULER BEAMS1)

本文将有限元p型超收敛算法应用于欧拉梁弹性稳定分析。该法基于有限元解答中失稳载荷和失稳模态结点位移的超收敛特性,建立了单元上失稳模态近似满足的线性常微分方程边值问题,在每个单元上,对该边值问题采用一个高次元进行求解,获得失稳模态的超收敛解,再将失稳模态的超收敛解代入瑞利商的解析表达式,最终获得失稳载荷的超收敛解。该法思路简明,通过少量计算即可显著提高失稳载荷和失稳模态的精度与收敛阶。数值算例表明,该法高效、可靠,值得进一步研究和推广到各类杆系结构。     

RELIABILITY SENSITIVITY STUDY OF THE BENDING VIBRATION OF THE COMPRESSOR ROTOR

以转子系统H858 型离心式压缩机为研究对象,利用MSC.Patran&Nastran 软件对该转子系统弯曲振动情况下的固有频率进行分析,定义了系统振动问题的频率可靠性模式,并以此为基础,应用随机摄动技术、可靠性理论和灵敏度方法,推导出参数为正态随机变量时频率可靠性灵敏度的计算公式,并应用于转子系统弯曲振动情况下的频率可靠性灵敏度研究中,同时将理论分析结果与基于Nessus软件的有限元分析结果进行比较,进一步验证了所提方法的有效性.     

钢--混凝土组合结构抗震性能研究进展简

钢-混凝土组合结构的抗震性能受到了工程研究领域的广泛关注. 本文总结了国内外研究者在组合梁、组合柱、组合节点及组合框架结构抗震性能方面的试验研究概况,分析了组合构件的作用机理及地震损伤演化累积效应,讨论了钢-混凝土组合结构地震反应分析的数值模型,包括微观模型、宏观模型以及近年来迅速发展的多尺度模型,并阐述了组合结构抗震性能评价指标及抗震性能水平的确立方法.