论文部分内容阅读
结构与复杂层状无限地基相互作用是无限域静力和动力分析中重要环节,准确模拟层状无限域辐射阻尼是解决结构与无限地基相互作用分析的重要基础。大量研究表明,复杂层状无限地基对上部结构的影响不容忽略,其中各向异性无限地基的影响尤为明显。在实际工程中地基常呈现层状分布,而大多学者为简化计算选取均质地基为研究对象,因此研究复杂层状地基静力和动力问题求解更具有普遍性和实际意义。针对层状无限域与结构相互作用问题,广大学者开展了大量的研究工作,实验和多种数值方法被用于求解层状无限域问题。根据结构与无限地基相互作用研究现状可知,现有大多方法求解大型模型时存在计算过程较繁琐,计算效率不高等问题,使很多研究成果不能大范围推广。因此,提出一种新型高效的数值方法分析结构与复杂层状介质相互作用问题是十分必要的,同时也具有一定的实际意义和科学研究价值。本文基于相似中心轴概念演化的改进比例边界有限元法,提出了一种适于分析复杂无限层状介质的新型数值方法,推导了层状无限域的新型连分式解和传递边界,建立了改进比例边界有限元时域方程,对复杂层状地基的静力和动力特性进行广泛深入研究。由于比例边界有限元法仅需离散问题域边界,大大降低了计算量,且与边界元法相比不需要基本解,适用于模拟复杂形状模型。将相似中心用相似轴线代替,可直接模拟水平层状介质,避免了将相似中心置于无穷远处近似求解,同时基于子结构法,可对倾斜层状无限域问题进行求解。本文研究主要内容如下:1.基于改进比例边界有限元法,建立了复杂层状无限时域计算模型,求解了复杂层状无限域的动力响应。本模型仅需离散问题域边界,不需要基本解,且自动满足无限域辐射边界条件。由于采取相似中心轴概念,可直接求解层状无限域问题。推导过程中分别采用传统加权余量技术和哈密顿变分原理,建立了各向同性和各向异性无限域改进比例边界有限元控制方程,给出了改进比例边界有限元动刚度的新型连分式解。通过引入辅助变量,将二阶偏微分方程转换为一阶偏微分方程,简化计算难度。采用移谱法消除系统内虚假模态导致的奇异性,提高数值结果的稳定性,并且减小计算规模。将精细积分法和改进比例边界有限元法耦合求解各向异性全域运动方程,提高了计算效率和精度。提出了改进比例边界有限元子结构法分析倾斜层状模型,扩展本文算法的适用范围。利用本文算法求解各向同性和各向异性层状无限域问题,对含有软弱层、隧洞、倾斜层等模型和各向异性参数等进行详细分析,并与参考解对比验证了本文算法的准确性。2.建立层状无限域改进比例边界有限元静力模型,给出侧边荷载和侧边约束位移作用下的静力解。建立新型侧边荷载和侧边约束位移表达式,由于两者均展成径向坐标的指数形式,使其适合分析层状介质静力问题。为便于分析含孔洞有限域问题,利用子结构法求解孔洞结构,扩展了改进比例边界有限元法在分析有限域问题中的应用范围。通过讨论和分析数值算例的静力响应,验证了本文解法的准确性,并对各向异性参数、侧边荷载形式、荷载组合形式、含洞模型等进行了详细分析,为实际工程设计提供技术依据。3.基于加速度单位脉冲响应分析,建立一种适于水平层状介质的新型波动传递边界。此方法适用于求解标量和矢量无限域波动问题,扩展了改进比例边界有限元法的应用范围。在时域内直接进行求解,避免时域和频域转换导致的误差,采用分段线性假设和外推参数法离散加速度单位脉冲响应,即当作用时间超过线性截止时间,加速度单位脉冲响应函数呈线性变化,避免求解时域分析中的卷积积分,提高计算效率。同时给出结构与无限地基相互作用力的新型积分离散格式,建立了层状无限地基时域控制方程,求解时域运动方程得到复杂无限介质的动力响应。利用该算法分别计算标量和矢量无限域问题,分析截止时间和外推参数对加速度单位脉冲响应函数的影响,验证了本文算法的准确性,并进行了详细的参数研究,分析了复杂隧洞结构中隔震沟的隔震效应,说明该算法具有广泛地适用性。4.推导出三维层状半空间波动解耦方程,采用傅里叶-贝塞尔变换建立横观各向同性层状半空间在频域-波数域内波动方程,将复杂偏微分波动方程解耦成二阶常微分方程进行求解,采用精细积分法进行分层求解,根据相邻层间位移与应力关系进行组装,利用逆变换得到不同边界条件下层状地基的动刚度或动柔度阵,分析层状半空间土体的动力特性。该算法可直接准确模拟三维层状半空间无限域,表明本文算法具有广泛的实际意义。各向同性和各向异性数值算例结果表明本文算法具有较高精度,并讨论了各向异性参数对模型动力特性的影响。