在海洋工程中,波浪荷载及其它循环荷载作用下海底沉积物变形与强度特性直接影响结构地基的稳定性。海床失稳或地基承载力失效则对海洋结构产生破坏,轻则造成一定量的经济损失,重则引发海洋生态灾害。在低应力条件下,砂土具有可液化特性,使得砂土海床地基在外界动力作用下更容易引发灾害问题。数值方法是目前评价海洋工程场地动力效应最为有效的技术手段,但需要结合能够描述复杂加载的土体本构模型。本文基于砂土基本力学特性开发了旋转硬化本构模型并结合数值平台Abaqus将其应用于海洋工程中。开展的工作及取得的成果如下:1.构建砂土等向硬化本构模型:通过分析砂土压缩特性及其在0)-lg平面内非线性特点,提出更适合于描述砂土在等向压缩路径下的参考压缩线,并使用对数-幂函数进行数学描述;在此基础上引入统一硬化(UH)参量构建了砂土体积硬化参量。提出适用于描述砂土剪胀特性的屈服面表达式。通过分析屈服面形状参数对屈服面形状的影响,选取了合适的屈服面形状参数并在等有效压力()路径下给出了另外一个屈服面形状参数a的确定方法。利用当前应力比所对应的压缩线作为状态参量参考线,以调控硬化参数与剪胀方程,进而反映砂土的压缩及剪胀特性。这种做法实现了砂土由压缩至剪切过程中状态参量的统一表述。利用Toyoura砂的等向压缩试验、不排水剪切试验及排水剪切试验对等向硬化本构模型进行验证,并分析了屈服面形状参数及参考压缩线斜率对本构模型预测结果的影响。研究结果表明等向硬化本构模型能够很好描述各向同性条件下Toyoura砂在不同孔隙比和不同压力下的力学特性。2.建立旋转硬化本构模型:在等向硬化本构模型基础上引入旋转硬化规律,以描述应力诱导各向异性对砂土力学特性影响。建议的屈服面随着各向异性的发展而逐渐变为扁平,并能够退回至等向硬化本构模型屈服面。采用插值方法获取各向异性压缩线,使屈服面旋转时参考压缩线随之改变。在体积硬化参量中引入统一硬化(UH)参量的同时,也考虑应力诱导各向异性影响,使屈服面的大小随应力诱导各向异性的发展而改变。引入考虑各向异性影响的状态参量来调控统一硬化(UH)参量和剪胀方程,以反映砂土的剪胀特性;利用Toyoura砂与Neveda砂的三轴单调加载试验与循环加载试验验证旋转硬化本构模型预测能力,结果表明本文建立的旋转硬化本构模型很好地描述了砂土单调加载与循环加载下的力学特性;此外,对比分析旋转硬化参数b_l、b_r与剪胀修正参数R_k对模型预测的影响,并针对Toyoura砂给出了合理的b_l、b_r、R_k取值范围。3.开发基于旋转硬化本构模型的Abaqus子程序:根据UH模型加卸载特点,修正了次阶应力点法中屈服面与当前应力点之间的相关关系,从而使次阶应力点法更适用于旋转硬化本构模型。基于优化的次阶应力点法将旋转硬化本构模型嵌入至Abaqus子程序UMAT中,实现了有限元软件的二次开发并完成了程序的调试。通过对比不排水剪切试验、排水剪切试验及循环加载试验的单元预测与Fortran程序预测结果,证明了基于次阶应力点法的有限元二次开发程序完全正确并且修正的次阶应力点法也较为合理。4.基于旋转硬化本构模型分析波浪-砂质海床相互作用:利用嵌入旋转硬化本构模型的Abaqus研究波浪作用下砂土海床力学特性,在此过程中利用弹性理论验证波浪荷载与模型边界条件的正确性,采用离心机水槽试验结果验证有限元模型在分析波土相互作用的适用性。研究极端海况下孔隙水压力响应、剪切应力变化特性及波浪引起海床瞬态液化过程,结果表明弹性海床与弹塑性海床对波浪的响应过程完全不同,波浪作用下弹塑性海床存在孔隙水压力累积现象,引起海床有效应力严重下降并部分达到液化状态,而弹性海床则不存在这一问题。