粘弹性饱和土体中半封闭圆形隧洞的动力响应分析

基于Biot波动方程,研究分析了粘弹性饱和土体中半封闭圆形隧洞的动力响应问题.假定衬砌材料为多孔介质,引入了更符合工程实际的半封闭边界条件.通过引入势函数,在Laplace变换域中得到隧洞边界上作用轴对称荷载和流体压力条件下应力、位移和超孔隙水压力的解答.利用Laplace数值逆变换得到时域中的解,分析了隧洞边界的半透水特性对隧洞动力响应问题的影响,结果表明:隧洞边界的半透水特性对应力、位移场的变化和超孔隙水压力的消散有很大的影响,透水和不透水下条件的解仅是本文的两个特例.     

基于分数导数本构的粘弹性土层—隧洞衬砌系统的稳态动力响应

将混凝土衬砌材料视为具有分数导数本构的粘弹性体,在频率域研究了深埋圆形隧洞粘弹性土的稳态动力响应。利用衬砌的内边界条件和混凝土衬砌与土体界面处的连续性条件,得到了粘弹性土体和衬砌稳态振动的应力和位移解析解。进行了隧洞衬砌厚度、土体阻尼比、本构阶数、材料参数比等对土体的位移和应力幅值影响的算例分析,结果表明:材料参数比对系统动力特性有较大影响,随着材料参数比的增加,位移和应力幅值减小;随着土体的阻尼比、衬砌厚度的增加,位移和应力幅值减小;当材料参数比Tσ/Tε=3时,随着分数导数阶数增大,响应幅值减小。     

分数阶导数粘弹性模型的有限元法

在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。     

考虑老化的混凝土粘弹性分数导数模型

混凝土是一种具有分形结构的材料。采用分数微积分模型来研究具有分形结构材料的老化规律目前尚未见到。本文的目的是采用含分数阶导数的类标准线性体来模拟考虑老化的混凝土的蠕变和松弛规律。给出了分数导数与Abel核之间的关系。讨论了类标准线性体的蠕变柔量和松弛模量及其在考虑老化的混凝土中的应用。与传统的混凝土流变模型相比较表明，类标准线性体可以更好地同时拟合混凝土在不同龄期的蠕变和松弛曲线。而且其形式简单、统一，在计算过程中需要调整的参数很少。可以预见，类标准线性体在混凝土的结构设计和计算中将有着广泛的应用前景。     

上覆分数导数粘弹性场地土地震放大效应

基于一维波动模型和分数导数粘弹性本构关系,分析了在竖直方向上传播的剪切地震波作用下,基岩上分数导数粘弹性模型描述的场地土的横向振动问题,用直接刚度矩阵法求得了场地土的地震放大效应系数,并用数值算例讨论了相关参量对分数导数粘弹性场地土地震放大效应系数的影响。研究结果表明:在简谐剪切地震波作用下,分数导数粘弹性场地土存在共振现象;分数导数的阶数、模型参数和基岩与上覆场地土层底部之间的阻抗比对场地土的地震放大效应系数有较大的影响。     

粘弹性基支粘弹板轴对称问题的动力响应

本文讨论粘弹性半空间地基上粘弹性圆板受轴对称载荷时的动力响应,将问题化为在Laplace变换空间中的第一类Fr积分方程,通过数值解和进行数值逆变换求得问题的解答。     

具有分数导数本构关系的非线性粘弹性Timoshenko梁动力学行为分析

本文利用分数导数型本构关系建立了在有限变形情况下Timoshenko梁的控制方程并利用Galerkin方法进行简化。然后利用一种存储部分历史数据的分数积分的计算方法对梁的控制方程进行求解。考察了载荷参数和分数导数参数对梁振动的影响,并采用非线性动力学中的各种数值方法,如时程曲线、功率谱、相图、Poincare截面等,揭示了非线性粘弹性Timoshenko梁丰富的动力学行为。     

分数导数粘弹性模型描述的土中管桩竖向振动的频域解

在分数导数粘弹性本构模型的基础上综合考虑桩周土和桩芯土的平衡方程和几何方程建立了桩周土和桩芯土的竖向运动的控制方程．在频率域内利用分离变量法和分数导数的性质求解了桩周土和桩芯土竖向振动控制方程．考虑管桩与桩周土、管桩与桩芯土的边界连续性条件以及三角函数的正交性得到了分数导数粘弹性模型描述的土中管桩的竖向振动,通过数值分析研究了管桩和土体模型参数和几何参数对管桩的桩顶复刚度的影响规律．结果显示：桩芯土本构模型的分数导数的阶数对管桩竖向振动的影响较桩周土本构模型的阶数要小,且与频率有一定关系;桩芯土与桩周土的模型参数比τ1 和τ2 对等效阻尼的影响较对刚度因子的影响要大;管桩桩周和桩芯的直径比d 越小,管桩复刚度的实部和虚部就越大;土体力学性能对管桩竖向振动的影响要比管桩桩身力学性能的影响小．     

纤维加强复合材料的总体粘弹性响应

在某些纤维增强复合材料(FRC)中使用金属或高分子聚合物作为基体材料。在高温等情况下,这类材料具有明显的粘弹性特性。本文采用Riemann—Liouville形式的分数阶导数模型描述基体的粘弹性特性。通过渐近均匀化方法给出了预测FRC整体三维本构关系的解析表达式。给出了应用于基体具有Makris粘弹性关系的具体形式。以圆截面纤维正方形排列的情形为例,给出了等效模量随纤维体积比的变化曲线。结果说明,这类复合材料仍具有粘弹性特性,其整体粘弹性本构关系的弹性部分综合了纤维弹性和基体弹性的贡献,粘性部分来自基体粘性的贡献,复合材料具有和基体相同的粘性系数和分数阶。为分析微结构特征对整体特性的贡献,须求解两类局部问题。可以看出,在整体的等效模量中包含了局部变形的贡献,局部变形增加了复合材料的耦合刚度。     

含阻尼效应的黏弹性杆稳态响应的解析解

基于格林函数法和叠加原理,给出了简谐激励下开尔文黏弹性杆稳态响应的解析解。运用拉普拉斯变换法得到了任意边界条件下黏弹性杆强迫振动的格林函数解。设计了一个统一方案,将杆的各种边界表示为矩阵形式,便于确定格林函数中的待定常数。算例表明,材料黏性阻尼和外部阻尼导致杆位移在空间和时间上不可分离。

     

一种工程实用的分数指数黏弹性固体模型

从热力学对模型的限制出发,提出了一种工程实用的、新固体流变模型-分数指数模型,它可与实验结果很好吻合,具有简单实用、计算速度快、自相谐调的优点,克服了目前固体模型收敛速度太慢、计算量太大的严重缺点,可以广泛应用于线性或非线性黏弹性被子示问题的计算之中。     

含黏弹性夹芯FG-GRC后屈曲梁的低速冲击响应

研究了含黏弹性夹芯的功能梯度石墨烯增强复合材料(functionally graded graphene reinforced composite, FG-GRC)后屈曲梁在低速跌落冲击下的跳跃振荡行为.采用修正Halpin-Tsai细观模型预测FG-GRC的材料宏观属性.使用赫兹点接触模型确定冲击器和梁之间的接触力.提出了考虑轴向预应力的复合材料层本构关系和阻尼层的Kelvin型黏弹性本构.通过一种广义高阶剪切变形锯齿梁模型建立夹芯梁的非线性位移场. 基于Hamilton 能量变分原理, 推导了动力学控制方程组. 通过两步分析,首先获得弹性后屈曲平衡路径作为冲击问题的初始状态. 随后, 结合四阶龙格库塔法,拓展了两步摄动-伽辽金法计算接触力的时程曲线以及后屈曲梁的位移时程曲线.研究了后屈曲梁在单次和两次撞击下双稳态大幅振荡过程的动力学特征.讨论了轴向载荷、冲击速度、黏弹性阻尼特性、冲击器材料等因素对于碰撞接触力以及后屈曲梁动力响应的影响规律.结果表明, 接触力仅对冲击速度较为敏感,一定的结构碰撞参数设计可以在接触力变化不大的情况下,使得后屈曲梁由单势能阱运动转变为双阱大幅振荡.      

广义Maxwell黏弹性流体在两平板间的非定常流动

将分数阶微积分运算引入Maxwell黏弹性流体的本构方程，研究了黏弹性流体在两平板问的非定常流动．对于广义Maxwell黏弹性流体的分数阶导数模型，导出了对时间具有分数阶导数的特殊运动方程，利用分数阶微积分的Laplace变换理论，得到了流动的解析解．     

三种分形和分数阶导数阻尼振动模型的比较研究

标准的整数阶导数方程不能准确描述粘弹性材料的记忆性参考文献[1]和阻尼的分数次幂频率依赖[2],因此分形导数、分数阶导数及正定分数阶导数被用于描述粘弹性介质中的阻尼振动.该文通过分析模型和数值模拟,比较了三种模型描述的振动过程.结果显示,当p小于约O.75或大于约1.9时(p为非整数阶导数的阶数),分形导数模型衰减最快;当P大于约0.75且小于约1.9时,正定分数阶导数模型衰减最快,衰减最慢的分别为分数阶导数模型(p1).且正定分数阶导数模型衰减快于分数阶导数模型,当p接近2时,两种模型较为相近.     

公路隧道新型衬砌结构双孔开挖与支护方案的研究

论文以公路隧道为研究对象,取Ⅳ级围岩,利用有限元程序对台阶开挖法和环形开挖法进行了数值模拟,分析了其支护后衬砌结构的应力和变形特征.在此基础上,针对目前隧道衬砌结构在修建和运营过程中的不足,提出一种双孔支护方案,并进行了较系统的数值模拟.研究结果表明:环形开挖法的最大应力和位移明显小于台阶开挖法,双孔支护方案能很好地减小衬砌的受力和变形,为隧道的设计和施工提供一种新的思路.     

丙烯酸弹性体的率相关分数阶黏弹性模型研究

丙烯酸弹性体VHB 4910作为一种重要的介电弹性体,在软体机器人、致动器、俘能器和智能隔振器等领域有很好的应用前景.但材料的非线性黏弹性对其力学行为有显著影响.近来分数阶模型在复杂材料的建模中取得了成功.本文基于分数阶有限变形Kelvin-Voigt流变学模型建立弹性体的三维张量本构,并进一步推导单向拉伸情况下的本构...     

变截面黏弹性直杆纵振时复频率分析的差分解法

从一维黏弹性本构方程出发,导出了黏弹性变截面直杆纵向振动微分方程的一般形式,采用了有限差分法,并以二阶矩阵表示的递推形式,建立了该问题的复特征值方程组。两种Maxwell黏弹性变截面（指数指数、线性函数）直杆的数值计算表明,该方法运算简单,计算精度高,能适用于求解任意变截面黏弹性直属的纵向自由振动问题。