论文部分内容阅读
摘 要：为了更好地研究预制基础桩的非线性动力特征及其它特性，考虑桩-土之间的摩擦影响，利用kelvin模型，建立在轴向荷载下粘弹性桩基的非线性振动动力学模型，通过力学模型分析得到结构的动力偏微分方程；运用微分求积法（DQM）将偏微分方程在空间域进行网格划分并离散化，进而导出粘弹性桩基的常微分动力方程；最后用matlab数值模拟得到不同桩基弹性模量下的相平面图、功率谱图、庞加莱截面图和时程曲线图.结果表明：在轴向荷载作用下粘弹性桩基会发生混沌运动，且桩基弹性模量越大桩基越容易发生混沌运动，同时更加直观形象地验证了混沌效应的基本特征.
关键词：粘弹性桩；kelvin模型；微分求积法；混沌
中图分类号：TU311.3 文献标志码：A
0 引言
随着科技的发展，大量新型材料被广泛运用到实际工程中，使得传统结构呈现出粘弹性的性能，并引起了大量的学者对其进行研究.桩基作为工程中常见的结构，建立其非线性动力学模型对研究桩的动力学分析有重要的意义.胡春林等[1-2]研究非线性粘弹性嵌岩桩纵向振动中运用Galerkin方法将非线性偏微分方程进行简化得到简化的偏微分-积分动力学方程.刘明等[3]分析桩基轴向振动动力学时运用多时间尺度法对动力学偏微分方程进行近似简化.骆文和等[4]运用分离变量法对动力学方程简化分析了桩-土之间的竖向振动.随后，学者运用复模态法[5]、分数导数微分[6]等方法研究分析了桩基非线性动力学振动.
为了更好地分析桩基的非线性动力学行为，如何建立粘弹性桩的动力学方程就显得尤为重要.本文采用微分求积法[7-9]对粘弹性桩基的动力学偏微分方程在空间域进行网格划分离散简化，得到粘弹性预制基础桩基的非线性动力学常微分方程模型，最后采用四阶Runge-Kutta法运用matlab数值模拟分析桩基在轴向荷载下发生的复杂的非线性动力学行为.
1 力学模型的建立
一个等截面圆形基础桩.设桩长为l，直径为d，截面面积为A.材料的密度为ρ，桩在轴向荷载的作用下，截面始终保持为平面.假定桩基向上为正半轴，桩基底部的圆心为坐标系的原点处，如图1所示.考虑桩-土之间的相互作用力，假定周围土对桩基的反力大小与位移成正比，且周围土对桩身的阻力大小与其位移速率也成正比，则预制基础桩的动力学模型为：
3 数值模拟
采用四阶Runge-Kutta法对非线性常微分动力学模型运用matlab数值模拟分析.初始参数为桩的直径A=0.1 m2，材料的密度ρ=2.5 g/cm3，土的刚度系数k=2.2×105 Pa，阻尼系数C=2.3×103 Pa·s，桩身的粘弹性系数？浊=1×105 Pa·s.外部轴向荷载作用力为F=500 N.
假设材料弹性模量E的变化范围为2.1×108 Pa～2.1×1010 Pa，以步长为10增加，进而得到3组所对应的相平面图、功率谱图、庞加莱截面图、时程曲线图，如图2～图4所示.
分别对比分析图2～图4中所对应的相平面图、功率谱图、庞加莱截面图、时程曲线图，可发现随着桩基弹性模量的增加， 桩基将从周期运动进入混沌运动， 弹性模量越大越容易发生混沌运动， 同时图3和图4也符合混沌效应的基本特征，具体分析如下：
1）对于相平面图分析，图2中的相平面图形中明显是个封闭的运动轨线，图3和图4的图形中明显扭结在一起形成一个不封闭的运动轨迹，具有无限层次的自相似性；
2）对于功率谱图分析，图2中功率谱图中对应尖峰，而图3和图4中图形出现“噪声背景”和宽峰；
3）对于庞加莱截面图分析，图2中的庞加莱截面图上是个闭合的曲线图形，而图3和图4中的图形渐渐呈现出成片的密集点，越来越具有层次结构性；
4）对于时程曲线图分析，图2中时程曲线图在一段时间后出现了一定的规律，没有了跳动性变化，非常稳定，而图3和图4中的功率谱图形随着时间的推移，一直有跳跃性变动.
4 结论
本文考虑了桩-土之间摩擦影响，建立了在轴向荷载下粘弹性预制基础桩的动力学偏微分运动学方程.利用微分求积法在空间域进行网格划分离散对其简化得到常微分方程，这种方法更加简洁，在matlab的分析中更为准确.数值分析表明：受轴向荷载作用下的粘弹性桩基可以呈现出周期运动和混沌运动；桩身弹性模量对桩基的混沌运动有影响，弹性模量越大桩基越容易发生混沌运动；同时图形也更加直观形象地反映出了混沌运动的基本特征.
参考文献
[1] 胡春林，程昌钧.非线性粘弹性基桩纵向混沌运动[J].岩土力学，2005，26（10）：1541-1544.
[2] 任九生，程昌钧.非线性粘弹性桩耦合运动中的混沌分析[J].振動与冲击， 2006， 25（5）：21-23，189.
[3] 刘民，陈中学，胡春林.基于多时间尺度法的桩基轴向振动分析研究[J].长江科学院院报， 2009，26（12）：102-105.
[4] 骆文和，闫启方.考虑桩-土相互作用的粘弹性土中管桩的纵向动力阻抗分析[J].昆明理工大学学报（理工版），2010， 35（5）：28-32，36.
[5] 刘民，陈中学.基于复模态法的粘弹性基桩横向振动研究[J].长江科学院院报，2010，27（12）：76-78，82.
[6] 刘林超， 闫启方，杨骁.分数导数粘弹性土层模型中桩基竖向振动特性研究[J].工程力学，2011，28（8）：177-182.
[7] 刘静，朱媛媛，胡育佳.基础-饱和土地基耦合系统分析的微分求积单元法[J].力学季刊，2014，35（1）：54-65.
[8] 陈姗，琚宏昌.简谐荷载作用下粘弹性梁振动的非线性动力学模型及其简化[J].广西科技大学学报，2014， 25（4）：30-33，39. [9] 王冬梅，张伟，李慕荣.用微分求积法分析轴向移动粘弹性梁的非平面非线性振动[J].动力学与控制学报， 2015，13（1）：23-27.
Abstract： To study the nonlinear dynamic characteristics and other characteristics of prefabricated foundation piles better， considering the influence of friction between pile and soil and based on the Kelvin model， nonlinear vibration dynamic model of viscoelastic pile under axial load is established. The dynamic partial differential equation of structure is obtained by the analysis of the mechanical model， the partial differential equations are partitioned and discretized in space domain by differential quadrature method （DQM）， and the ordinary differential dynamic equations of viscoelastic foundation pile are derived. Finally， phase diagram， power spectrum， poincare section and time history curve are gained with numerical simulation of matlab. The results show that the chaotic motion can occur in the viscoelastic foundation pile under axial load； the bigger the elastic modulus of pile foundation， the more prone to chaotic motion； meanwhile， the basic features of chaotic effects are more intuitively and vividly verified.
Key words： viscoelastic foundation pile； Kelvin model； differential quadrature method； chaos
（學科编辑：黎 娅）
关键词：粘弹性桩；kelvin模型；微分求积法；混沌
中图分类号：TU311.3 文献标志码：A
0 引言
随着科技的发展，大量新型材料被广泛运用到实际工程中，使得传统结构呈现出粘弹性的性能，并引起了大量的学者对其进行研究.桩基作为工程中常见的结构，建立其非线性动力学模型对研究桩的动力学分析有重要的意义.胡春林等[1-2]研究非线性粘弹性嵌岩桩纵向振动中运用Galerkin方法将非线性偏微分方程进行简化得到简化的偏微分-积分动力学方程.刘明等[3]分析桩基轴向振动动力学时运用多时间尺度法对动力学偏微分方程进行近似简化.骆文和等[4]运用分离变量法对动力学方程简化分析了桩-土之间的竖向振动.随后，学者运用复模态法[5]、分数导数微分[6]等方法研究分析了桩基非线性动力学振动.
为了更好地分析桩基的非线性动力学行为，如何建立粘弹性桩的动力学方程就显得尤为重要.本文采用微分求积法[7-9]对粘弹性桩基的动力学偏微分方程在空间域进行网格划分离散简化，得到粘弹性预制基础桩基的非线性动力学常微分方程模型，最后采用四阶Runge-Kutta法运用matlab数值模拟分析桩基在轴向荷载下发生的复杂的非线性动力学行为.
1 力学模型的建立
一个等截面圆形基础桩.设桩长为l，直径为d，截面面积为A.材料的密度为ρ，桩在轴向荷载的作用下，截面始终保持为平面.假定桩基向上为正半轴，桩基底部的圆心为坐标系的原点处，如图1所示.考虑桩-土之间的相互作用力，假定周围土对桩基的反力大小与位移成正比，且周围土对桩身的阻力大小与其位移速率也成正比，则预制基础桩的动力学模型为：
3 数值模拟
采用四阶Runge-Kutta法对非线性常微分动力学模型运用matlab数值模拟分析.初始参数为桩的直径A=0.1 m2，材料的密度ρ=2.5 g/cm3，土的刚度系数k=2.2×105 Pa，阻尼系数C=2.3×103 Pa·s，桩身的粘弹性系数？浊=1×105 Pa·s.外部轴向荷载作用力为F=500 N.
假设材料弹性模量E的变化范围为2.1×108 Pa～2.1×1010 Pa，以步长为10增加，进而得到3组所对应的相平面图、功率谱图、庞加莱截面图、时程曲线图，如图2～图4所示.
分别对比分析图2～图4中所对应的相平面图、功率谱图、庞加莱截面图、时程曲线图，可发现随着桩基弹性模量的增加， 桩基将从周期运动进入混沌运动， 弹性模量越大越容易发生混沌运动， 同时图3和图4也符合混沌效应的基本特征，具体分析如下：
1）对于相平面图分析，图2中的相平面图形中明显是个封闭的运动轨线，图3和图4的图形中明显扭结在一起形成一个不封闭的运动轨迹，具有无限层次的自相似性；
2）对于功率谱图分析，图2中功率谱图中对应尖峰，而图3和图4中图形出现“噪声背景”和宽峰；
3）对于庞加莱截面图分析，图2中的庞加莱截面图上是个闭合的曲线图形，而图3和图4中的图形渐渐呈现出成片的密集点，越来越具有层次结构性；
4）对于时程曲线图分析，图2中时程曲线图在一段时间后出现了一定的规律，没有了跳动性变化，非常稳定，而图3和图4中的功率谱图形随着时间的推移，一直有跳跃性变动.
4 结论
本文考虑了桩-土之间摩擦影响，建立了在轴向荷载下粘弹性预制基础桩的动力学偏微分运动学方程.利用微分求积法在空间域进行网格划分离散对其简化得到常微分方程，这种方法更加简洁，在matlab的分析中更为准确.数值分析表明：受轴向荷载作用下的粘弹性桩基可以呈现出周期运动和混沌运动；桩身弹性模量对桩基的混沌运动有影响，弹性模量越大桩基越容易发生混沌运动；同时图形也更加直观形象地反映出了混沌运动的基本特征.
参考文献
[1] 胡春林，程昌钧.非线性粘弹性基桩纵向混沌运动[J].岩土力学，2005，26（10）：1541-1544.
[2] 任九生，程昌钧.非线性粘弹性桩耦合运动中的混沌分析[J].振動与冲击， 2006， 25（5）：21-23，189.
[3] 刘民，陈中学，胡春林.基于多时间尺度法的桩基轴向振动分析研究[J].长江科学院院报， 2009，26（12）：102-105.
[4] 骆文和，闫启方.考虑桩-土相互作用的粘弹性土中管桩的纵向动力阻抗分析[J].昆明理工大学学报（理工版），2010， 35（5）：28-32，36.
[5] 刘民，陈中学.基于复模态法的粘弹性基桩横向振动研究[J].长江科学院院报，2010，27（12）：76-78，82.
[6] 刘林超， 闫启方，杨骁.分数导数粘弹性土层模型中桩基竖向振动特性研究[J].工程力学，2011，28（8）：177-182.
[7] 刘静，朱媛媛，胡育佳.基础-饱和土地基耦合系统分析的微分求积单元法[J].力学季刊，2014，35（1）：54-65.
[8] 陈姗，琚宏昌.简谐荷载作用下粘弹性梁振动的非线性动力学模型及其简化[J].广西科技大学学报，2014， 25（4）：30-33，39. [9] 王冬梅，张伟，李慕荣.用微分求积法分析轴向移动粘弹性梁的非平面非线性振动[J].动力学与控制学报， 2015，13（1）：23-27.
Abstract： To study the nonlinear dynamic characteristics and other characteristics of prefabricated foundation piles better， considering the influence of friction between pile and soil and based on the Kelvin model， nonlinear vibration dynamic model of viscoelastic pile under axial load is established. The dynamic partial differential equation of structure is obtained by the analysis of the mechanical model， the partial differential equations are partitioned and discretized in space domain by differential quadrature method （DQM）， and the ordinary differential dynamic equations of viscoelastic foundation pile are derived. Finally， phase diagram， power spectrum， poincare section and time history curve are gained with numerical simulation of matlab. The results show that the chaotic motion can occur in the viscoelastic foundation pile under axial load； the bigger the elastic modulus of pile foundation， the more prone to chaotic motion； meanwhile， the basic features of chaotic effects are more intuitively and vividly verified.
Key words： viscoelastic foundation pile； Kelvin model； differential quadrature method； chaos
（學科编辑：黎 娅）