非饱和土的应力状态变量研究

以热力学和连续介质力学的应力理论为基础,从非饱和土的土骨架平衡方程出发,导出了一系列应力状态变量,目前文献中所使用的应力状态变量都在其中。研究发现：非饱和土的应力状态变量有多种组合形式（5组以上）,不是惟一的;大多数应力状态变量与土的孔隙率、饱和度等物性指标孪生相伴;所有应力状态变量都具有应力的量纲,单纯孔隙率或饱和度及其组合都不构成应力状态变量;简化的Bishop有效应力公式和Fredlund倡导的两个应力状态变量都是简化了的应力状态变量,前者忽略了孔隙率的影响,而后者忽略了与孔隙率及饱和度的关联性。在选择应力状态变量时应考虑理论合理、应用方便、逻辑关系正确的原则,还应做一定的验证工作。     

深隧道围岩分区破裂的数学模拟

该研究为先前深隧道围岩分区破裂现象内变量梯度塑性模型的进一步发展。利用应变梯度模型研究了深隧道围岩的分区破裂现象。作为额外的状态变量,在此引入应变梯度这一新变量。利用虚功原理得到了深隧道围岩的平衡方程、边界条件和流动准则,利用Clausius-Duhem不等式获得了岩体的本构方程。对于圆形深隧道,由上述模型的一般方程得到了弹性变形情况下、具有下降段的弹塑性变形情况下和不考虑弹性变形的塑性变形情况下圆形深隧道围岩的支配方程,得到了解析解,并讨论了解析解的性质。这一模型不仅扩展了隧道围岩的经典弹塑性模型,也为下一步数值研究深隧道围岩的分区破裂现象奠定了理论基础。     

基于流固耦合理论某尾矿坝失稳特性及稳定性分析

选用典型的尾矿库工程实例,采用流固耦合和强度折减法相结合对其尾矿坝进行稳定性分析,确定尾矿库渗流场分布及浸润线的位置,该浸润线位置与实测浸润线位置吻合较好。建议强度折减过程中尾矿坝的失稳准则,计算安全系数,确定潜在滑裂面的位置,并与极限平衡法计算出的安全系数及临界滑裂面位置进行比较,表明强度折减法得出的结果与Bishop法结果接近,流固耦合-强度折减法在尾矿库稳定性分析的可行性。在该基础上研究尾矿坝潜在的失稳模式,将尾矿坝潜在的失稳模式分为局部失稳和整体性失稳两种,局部失稳为尾矿坝的部分坝面发生滑移,整体失稳为整个尾矿坝坝体发生失稳,并分析导致该尾矿坝两种失稳模式的主要因素,认为浸润线埋深过浅是导致局部失稳的主要原因。     

隧道水平冻结施工引起地表冻胀的历时预测模型

隧道水平冻结施工过程中,土体冻结引起体积膨胀,进而会在地表产生冻胀现象。实际工程一般采用多根冻结管形成冻结壁。冻结壁交圈前,地表冻胀由多个冻土柱的叠加膨胀变形引起;冻结壁交圈后,地表冻胀则由整个冻结壁的膨胀变形引起。鉴于此,考虑冻结壁的形成过程,基于随机介质理论,建立了隧道水平冻结施工引起地表冻胀位移的历时预测模型。同时对冻结外锋面半径和冻胀区域外半径这2个关键参数的取值方法进行了相关探讨。最后针对两个工程案例,采用该计算模型对地表冻胀位移进行分析,得到地表冻胀位移随时间的变化规律,并与现场实测结果相比较,验证了模型的可靠性。该模型应用于隧道水平冻结施工前、冻结期内任意时刻的地表冻胀位移预测,可为工程冻结实施方案的合理确定提供有效依据。     

海堤地基固结系数反演与工后沉降分形模型预测

海堤工后沉降分析预测对工程维护和灾害防治具有重要意义。以海堤工后初期监测资料为基础,运用分层总和法和太沙基一维固结理论反演了土层固结系数,并通过对固结系数变化规律的分析,提出采用改进变维分形预测模型对具有较强趋势性的固结系数进行预测,进而对海堤工后可能发生的长期沉降进行预测。以浙江省洞头县北岙后二期东围堤工程的海堤监测数据作为研究对象进行了验证计算,结果表明,该方法反映了实际土体的固结沉降机制,能够较好地预测海堤工后沉降趋势,具有一定的工程实用意义。     

空沟对列车运行引起的地基振动隔振效果研究

为研究空沟对高速列车引起地基振动的隔振效果,在已有研究基础上改进出一种新的解析研究模型。模型中首次运用了饱和半空间模型来研究地基上隔振沟对高速列车的隔振效果;隔振沟则通过在饱和半空间土体上设置3个合适宽度、截面为矩形的弹性层来模拟;中间矩形弹性层为路堤,路堤上方放置了枕木与轨道。枕木与轨道分别通过纵向异性Kirchhoff薄板与Euler梁来模拟,饱和土地基采用Biot多孔饱和介质理论来描述。控制方程通过傅立叶变换与傅立叶级数展开,在变换域中进行求解。研究表明,随着列车运行速度的提高,空沟的隔振效果明显提高;饱和土体固-液相的耦合作用对隔振沟的隔振效果的影响明显,尤其当列车运行速度超过土体表面Rayleigh波速时,随着土体渗透系数的增加,空沟隔振效果显著降低。此外,列车运行速度超过土体表面波速时,饱和土地基上空沟的隔振效果明显优于相应单相弹性地基上空沟的隔振效果。     

五台山新太古代绿岩带中高镁安山岩（HMA）的化学特征

采用高镁安山岩（HMA）新的分类方法,描述了五台地区新太古代绿岩带中HMA的化学特征,并与实验的HMA以及自然界中典型的HMA（日本的Chichijima地区HMA和Setouchi地区HMA）进行了对比,显示该区HMA具有与上述HMA类似的化学特征。初步认为五台新太古代HMA产于洋俯冲带的环境,进而表明板块构造理论的适用时间范围可推演到新太古代。     

Stop making sense: the trials and tribulations of qualitative data analysis

An account is given of a research project with particular attention to the place of data analysis and interpretation. It shines some light into what often remains the black box of qualitative research practice, the mysterious process of generating a research narrative based on the data. The paper provides other apprentice researchers with an insight into what the data analytical experience can be like, warts and all. It raises issues to do with qualitative data analysis, the use of CAQDAS, grounded theory, the relationship between the narratives of the researcher and the researched, and about the training of qualitative researchers.     

CMT data inversion using a Bayesian information criterion to estimate seismogenic stress fields

We developed an inversion method to estimate the stress fields related to earthquake generation (seismogenic stress fields) from the centroid moment tensors (CMT) of seismic events by using Akaike's Bayesian information criterion (ABIC). On the idea that the occurrence of an earthquake releases some part of the seismogenic stress field around its hypocentre, we define the CMT of a seismic event by a weighted volume integral of the true but unknown seismogenic stress field. Representing each component of the seismogenic stress field by the superposition of a finite number of 3-D basis functions (tri-cubic B-splines), we obtain a set of linear observation equations to be solved for the expansion coefficients (model parameters). We introduce prior constraint on the roughness of the seismogenic stress field and combine it with observed data to construct a Bayesian model with hierarchic, highly flexible structure controlled by hyper-parameters. The optimum values of the hyper-parameters are objectively determined form observed data by using ABIC. Given the optimum values of the hyper-parameters, we can obtain the best estimates of model parameters by using a maximum likelihood algorithm. We tested the validity of the inversion method through numerical experiments on two synthetic CMT data sets, assuming the distribution of fault orientations to be aligned with the maximum shear stress plane in one case and to be random in the other case. Then we applied the inversion method to actual CMT data in northeast Japan, and obtained the pattern of the seismogenic stress field consistent with geophysical and geological observations.