非线性结构时程分析的高阶单步法

本文提出了一种新的求解运动方程的直接积分法－高阶单步β法。     

海相沉积软土蠕变BP神经网络本构模型

海相沉积软土具有很强的蠕变特性,传统分级构建蠕变本构模型方法不实用,且很难真正反映岩土流变的非线性特性。为此,本文引入具有超强非线性映射和容错能力的BP神经网络模型,通过改进BP算法,根据江门软土的室内直剪蠕变试验结果,建立了海相沉积软土BP神经网络蠕变本构模型,避免了传统方法为满足试验曲线变化规律和蠕变特性而需要建立复杂的本构数学表达式。最后,利用上海地区软土蠕变实验结果对本文提出的方法进行了验证,并对BP神经网络蠕变模型在描述软土流变方面的特点进行了讨论。结果表明,本文建模方法简单,并能很好地描述软土的非线性蠕变问题。     

VTI介质中非线性对地震波频率和频宽变化的约束

采用改进型FCT有限差分法分别模拟各向同性、具有垂直对称轴的横向各向同性(VTI)以及非线性VTI介质中地震波场,随后采用小波分析探讨VTI介质中非线性对地震波频率和频宽变化的约束.数值模拟及分析结果表明：介质各向异性会引发地震记录中各波型的能量随炮检距的变化,而介质各向异性与非线性的耦合作用则会进一步加剧这种变化.在遵守能量守恒原理的前提下,介质的非线性会引发地震波场中各波型间能量的再分配.同时还会引发地震波速度的波动,且随着入射角的增大,波动幅度变小.这些现象在频率域内的响应表现为主频和频宽的降低或增大、地震记录中与各波型相对应的瞬时主频和频宽峰值出现时刻的偏移以及它们随炮检距变化而表现出的幅值波动.在本文给定的弱非线性介质中（非线性炽较薄且远离震源）,在震源能量较低的情况下,介质的非线性效应会造成地震波瞬时主频和频宽相对于线性介质情况分别为1/10～1/4和1/10～1/3左右的变化量.此外,随着入射波能量的下降,介质非线性所造成的地震波瞬时特征（振幅、主频、频宽）的变化量减小,但随着地震波在穿越非线性介质时路径的加长,非线性作用的累积效应使得波形畸变和频率分布特征的变化（尤其是频宽）等现象得以维持.     

声波全波形反演目标函数性态

地震波传播的复杂性所引起的地震反演中强烈的非线性问题是目前全波形反演在技术上遇到的最大难题,了解全波形反演中不同的目标函数随不同物性参数的不同摄动尺度的变化性态,对选择合理的反演方法和反演策略具有重要意义.本文参照Jannane等对波形反演目标函数性态的分析方法,通过变密度声波方程,分析了多种地震数据子集的不同目标函数随物性参数的摄动尺度的变化关系,重点分析了它们的非线性程度,为进行分步骤、分尺度全波形反演方法和反演策略的选择提供了理论指导.     

汶川地震中宝珠寺重力坝地震响应的三维有限元模拟

汶川地震中宝珠寺水电站遭受的地震烈度为8°(相当于水平峰值加速度0.2 g),远超过大坝的设计地震水平(0.1 g),震后大坝未见明显震损.为解释大坝在地震中的抗震现象,构建了坝址区三维模型.考虑坝体横缝非线性以及三个方向地震作用的不同组合方式,对汶川地震中大坝的动力响应进行有限元模拟.在此基础上,针对震后提高的抗震设防标准,进一步选取典型坝段,采用二维弹塑性方法对大坝进行抗震复核并分析可能的破坏模式.模拟结果表明:横河向地震分量起主导作用而顺河向地震作用相对较弱是宝珠寺重力坝在汶川地震中免于发生损坏的主要原因.坝顶混凝土发生挤压破碎缘于永久横缝在地震中高频渐开渐合行为引起的剧烈碰撞.宝珠寺重力坝对设计地震0.27 g的强震可以保持整体的安全性,对校核地震0.32 g的强震整体安全性降低,水库正常运行及抵抗余震的能力将受到影响.     

Transmission of obliquely incident waves at low-crested breakwaters: Theoretical interpretations of experimental observations

The paper presents a series of analytical and numerical investigations of oblique wave transmission at low-crested breakwaters. For a smooth breakwater, two important features of wave height and direction are analyzed to establish the generic nature of the wave transmission process at oblique incidence. The proposed framework of research is validated against laboratory data from the EU-sponsored project, DELOS. The numerical simulations exhibit a significant decrease of the transmission coefficient with increasingly oblique incidence at a smooth breakwater. The roles of wave-breaking, nonlinearity, wave-induced currents and set-up in determining the characteristics of oblique wave transmission are demonstrated in the paper. It is found that both the amplitude-dependent phase velocity and the decrease of mean wave period contribute to the change of mean wave direction on the transmission side. An attempt has also been made to qualitatively explain the different behaviour of oblique wave transmission at a rubble-mound breakwater.     

Occurrence probability of a freak wave in a nonlinear wave field

The Edgeworth’s form of a cumulative expansion of the probability density function (PDF) of surface elevation expands the maximum wave height distribution to predict the occurrence probability of freak waves. This study investigated the enhancement of the occurrence probability of freak waves due to the fourth order moment of surface elevation, kurtosis, change and found that the nonlinear effects on the occurrence probability of a freak wave linearly depends on kurtosis for a small number of waves N=250. The statistical theory was compared with field data, and freak waves sometimes appear when not expected by the Rayleigh theory, but they were predicted by the proposed theory.     

Evolution of wave shape over a low-crested structure

This paper investigates the evolution of wave shape over a low-crested structure (LCS) using a 2-D RANS-VOF model. The model predictions of surface elevation and wave skewness and asymmetry are in good agreement with the recent measurements collected in a small scale wave channel at the University of Cantabria (UCA). The empirical formulae relating wave skewness and asymmetry to local Ursell number by Peng et al. (2009) have been extended to include the effect of wave reflection and the ramp in front of LCS and a wider range of Ursell number in the present study. In the presence of LCS, wave skewness decreases slightly above the seaward slope, then increases rapidly up to a maximum value above the structure crest, and decreases drastically above the leeward slope. Wave asymmetry decreases sharply above the seaward slope to a negative minimum value at the structure crest, and then increases rapidly to a positive value above the leeward slope. Our bispectral analysis indicates that sum interactions increase skewness and decrease asymmetry while difference interactions have opposite effects and that the former dominate above the seaward slope and on the structure crest but the latter dominate above the leeward slope of LCS. The observed wave shape evolution over a LCS can be attributed to the changes in the interplay of sum and difference interactions. We found that incident wave height and wave period, relative structure freeboard, structure crest width and structure porosity are the controlling factors for wave shape evolution over LCS. This study provides new insights on the role of wave skewness and asymmetry in the breakwaters stability and sediment transport around the structure and on the beaches behind it.     

基于变形P-L模型的矩阵方程迭代精细横波预测

横波速度预测问题的关键有两个,一是如何建立合理的岩石物理模型,二是针对建立的横波预测目标函数,如何准确高效地求解.针对第一个问题,对Pride模型和Lee模型(P-L模型)进行变形,提出拟固结指数的概念,将干岩石模量和岩石基质模量相联系,变形后的P-L模型在没有降低P-L模型准确度的情况下简化了问题的复杂度,建立起了饱和流体岩石弹性模量与干岩石模量、岩石基质模量、混合流体模量之间的关系,进而计算理论上的纵波速度,并通过比较实测纵波速度与计算的理论纵波速度大小,最终建立了横波预测的目标函数.针对第二个问题,借鉴了地震反演的思路,将该目标函数的最优化问题转化为线性矩阵方程组迭代求解问题,通过几步迭代就可以求解出合适的拟固结指数,进而得到预测横波速度.实际验证和应用表明,该横波预测方法具有很好的稳定性和准确性,并且岩石物理模型的构建和目标函数的求解思路可用于其他储集类型地层的横波预测.     

变密度声波方程多参数全波形反演策略

多参数全波形反演中各参数之间的相互耦合增加了反演的非线性程度.通过分析各参数之间的相互影响,提出合理的多参数反演策略是解决该问题的有效途径.本文从变密度声波方程出发,首先研究了密度在速度反演中的重要作用,然后分析了速度对密度反演的影响程度,进而提出了一种有利于速度、密度分步联合反演的策略.第一步,利用给定的初始模型对速度、密度进行同时反演,得到比较可靠的速度反演结果;第二步,利用第一步反演得到的速度和给定的初始密度作为初始模型,继续进行双参数同时反演,这样可以同时得到比较可靠的速度、密度反演结果.为了进一步提高反演精度,将第二步反演得到的速度、密度作为初始模型,再进行下一轮双参数联合反演.二维理论模型实验结果充分说明了本文提出的这种反演策略的有效性.