道路曲线放样方法的实验

对曲线放样的常见方法进行了比较分析,提出了改进方法,结合实例对自由设站极坐标法、DGPS-RTK法曲线测设方法进行了深入讨论,给出了相应的数据处理方法和实验方法,实验结果表明DGPS-RTK技术使曲线放样精度更高,实现更便捷。     

基于差分法的半参数模型粗差估计

运用差分法,得到半参数模型中参数向量和非参数向量的估计,并且当粗差的位置已经确定时,可以利用基于最小二乘原理的粗差探测方法,估计出半参数模型中的观测量所含粗差。通过算例分析,证明了算法的有效性,最后提出了粗差定位与定值时应注意的问题。     

国家海洋经济发展试点政策效应及异质性评价——基于双重差分法的实证检验

国家海洋经济发展试点政策是我国为建设“海洋强国”的重要举措,2010—2013年期间,我国已设立5个试点省市。基于2005—2016年我国11个沿海省市的省级面板数据,运用双重差分法(DID)对国家海洋经济发展试点政策的海洋经济增长效应进行实证研究,并进一步探讨海洋科技创新水平在试点政策对海洋经济增长影响中的调节作用。研究发现：(1)国家海洋经济发展试点政策能够显著促进海洋经济增长;(2)试点省市海洋经济实力越强,对区域的拉动作用越大,国家海洋经济发展试点政策对环渤海经济区试点和泛珠三角经济区试点的海洋经济促进作用更大;(3)海洋科技创新水平具有显著的正向调节作用,海洋科技创新水平提高会加强国家海洋经济发展试点政策对海洋经济增长的正向作用。     

多腔钢-混凝土组合壁式柱弱轴方向抗震性能试验研究

用于异形柱框架结构的各种截面型式的多腔钢-混凝土组合柱具有良好的使用功能和可装配建造性能,与异形柱面外有侧肢不同,保证矩形截面壁式柱弱轴方向的抗震性能更为重要。为了探究多腔钢-混凝土组合壁式柱弱轴方向的抗震性能与受力机制,运用低周往复荷载试验与有限元分析相结合的方法开展研究。结果表明：压弯荷载作用下,多腔钢-混凝土组合壁式柱弱轴方向的滞回曲线饱满,其全过程受力可划分为稳定承载、转化过渡和峰后二阶效应显现三个阶段,试件对应的两个特征分界点的层间位移角分别为1/57和1/27;试件屈服点和失效破坏点对应的层间位移角分别为1/85和1/23,均大于规范限值,延性系数为3.73;两端钢管、中间钢板和腔内混凝土可协同受力,混凝土未见明显压碎和开裂,钢管和钢板未见明显损伤;所得有限元模拟结果与试验结果吻合度高,骨架曲线平均相对误差为4.6%。上述成果深化了对多腔钢-混凝土组合壁式柱弱轴方向抗震性能的认知,为工程应用提供了参考。     

探地雷达GPR正演模拟的时域有限差分实现

近年来,随着数字处理技术和电子技术的飞速发展,探地雷达（GPR）的实际应用范围迅速扩大,现已覆盖考古、矿产资源勘探、水文地质调查、岩土勘查、无损检测、工程建筑物结构调查、军事等众多领域,解决了很多工程实际问题,成为浅层勘探的有力工具.而探地雷达的理论研究与实际的应用相比,具有明显的滞后性.但是解释人员要达到精确地对探地雷达实际资料的进行解析,必须事先了解地质体的雷达反射剖面的特征,所以作为反演与解释基础的复杂地电模型的探地雷达正演模拟技术,就成了探地雷达理论研究的主要内容之一.本文以麦克斯韦两个旋度方程为基本出发点,运用K.S.Yee的空间网格模型理论和时域有限差分法的基本原理,推导出二维空间的探地雷达正演方程组,并详细地分析了差分格式中半空间步长与半时间步长的实现方法,及其雷达波电场与磁场分量在计算机上相互关系的C程序实现.然后讨论了数值频散关系及其产生原因,通过同时考虑时域有限差分法及Yee氏网格的特点,推导出了符合探地雷达实际传播规律的理想频散关系,作者自制了探地雷达正演程序,并分别计算了Mur超吸收边界条件及无边界条件下的雷达地电模型,通过对比可知,超吸收边界条件可利用,大大地减少截断边界处的干扰波,达到用有限区域达到在无限空间传播的效果.最后作者利用自制程序,对“V”字形和同一斜面上的五个圆的两个典型的探地雷达地电模型进行了正演模拟,得到了正演剖面图,消除了边界反射后的雷达剖面能很好地指导工作人员对雷达实测剖面的地质解释,同时使正演研究更符合实际的地质情况.     

频率域电磁剖面有限差分法2.5维正演数值模拟

对于实际生产中遇到的2.5维问题,从麦克斯韦方程组出发,把地电参数变化小的走向方向转化到波数域,用一系列波数模拟三维源的特征,并在波数域中,得到2组关于Hx (k_x,y,z)和Ex (k_x,y,z)的偏微分方程。选取适当的kx值,用有限差分法在y-z平面的网格中求解,再通过反傅里叶变换得到空间域中的电磁场。在验证了算法的正确性之后,对不同埋深的直立异常体、倾斜异常体及断陷模型进行了数值模拟,其结果直观地显示了异常体磁异常的形态,同时研究了程序对低阻异常体的横向分辨率。对层状大地模型,用井间收发方式进行了模拟研究,并取得了较好的结果,这对今后实际勘探应用奠定了基础。      

变荷载下基于指数形式渗流的一维固结分析

建立了考虑指数形式渗流以及变荷载条件下的一维固结微分方程,采用相对稳定的Crank-Nicolson差分格式获得控制方程的差分解答并验证了计算程序的可靠性。结果表明,当指数大于1时,较小时间因子下固结速率比达西渗流快,较大时间因子下固结速率比达西渗流慢;而当指数小于1时,较小时间因子下固结速率比达西渗流慢,较大的时间因子下固结速率比达西渗流快。在土层厚度相同的情况下,指数大于1时作用于土层的荷载越小,固结速率越慢;基于指数形式渗流,传统一维固结理论中室内土样固结与实际地基土层固结之间的相似关系不再成立;加荷速率越快,则土层的固结速率越快。     

基于三次插值函数算法的时间积分方案与二阶时空余差数值模式——以原始大气运动方程与理想全球模拟个例为例

从大气运动原始方程和欧拉算符出发,用泰勒级数展开,给出二阶时空微商余项预报方程。进而讨论用三次插值函数——双三次曲面拟合求上游点的准拉格朗日时间积分方案与相应的二阶时空余差数值模式——“双三次模式”。则双三次模式是通过实现各个大气物理量场的二阶可导,从而可对预报方程做空间非线性(“三次”)时间离散积分,成为“双三次曲面拟合——时间步积分——双三次曲面拟合——……”一种新算法数值模式。讨论双三次数值模式的数学基础:三次插值函数及其数值分析极性定律用于数值模式。指出:双三次模式和谱模式都具有数学“收敛性”;而Coons双三次曲面具有对变量场拟合二阶可导“最优性”;和Hermite双三次曲面片具有对网格变量场二阶可导运算“等价性”。又指出:有限差分模式的中央差近似斜率和曲率,分别是三次样条斜率和曲率作“三点平滑”。双三次模式适合采用原始大气运动方程,适合采用准拉格朗日时间积分方案,并给出一个理想全球模拟个例。因大气运动本质上是非线性的,理论上可按变量场双三次曲面曲率判断,以采用符合物理诠释的局域或单点平滑,以保持模式时间积分稳定性。且未来容易实现全球多重/时变套网格双三次数值模式。      

基于二阶锥规划有限元增量加载法的条形浅基础极限承载力分析

地基极限承载力分析是土力学研究中的一个经典课题。基于Hellinger-Reissner混合变分原理和有限元方法,将岩土体弹塑性问题构造成基于有限元框架的二阶锥规划（second-order cone programming,SOCP）问题,进而提出一种基于二阶锥规划理论的增量有限元法,即FEM-SOCP法。将岩土体弹塑性问题构造成二阶锥规划的数学优化问题,可以避免采用传统弹塑性计算中复杂的应力点积分等算法和屈服面棱角的平滑处理。此外,对于二阶锥规划问题,可以采用具有原始?对偶内点求解法的标准数学规划求解器MOSEK进行求解。将增量加载FEM-SOCP法应用于经典的基底粗糙的条形浅基础地基极限承载力分析中,分别考虑了关联和非关联塑性条件下的Mohr-Coulomb屈服准则。数值结果表明：所提出的增量加载FEM-SOCP法获得的地基承载力系数及地基承载力与传统FEM计算结果基本一致,而与常规有限元计算结果相比,基于增量加载的FEM-SOCP法所获得的屈服区更加平滑。     

丽江-清镇剖面上地壳速度结构及其与鲁甸MS6.5级地震孕震环境的关系

2014年8月3日,云南省昭通-鲁甸地区发生M_S6.5级地震,造成了重大的人员伤亡和财产损失.鲁甸震区位于扬子块体的西缘,小江断裂带的东侧北东向的昭通-莲峰断裂带内.由于至今没有穿越该断裂带的人工源深地震测深剖面,而丽江-攀枝花-清镇650 km长深地震测深剖面距离鲁甸主震区不超过50 km,利用宽角地震资料的初至波震相,通过有限差分反演揭示该地区上地壳速度结构,可以为鲁甸震区的地震定位、地震孕育机制等提供深部速度模型.速度剖面显示：剖面结晶基底厚度平均为2 km左右;小江断裂带速度较低,东西两侧的速度较高;因此小江断裂带区域地壳强度比较低,加上断裂两侧的应变速率很高,所以小江断裂带和旁边的鲁甸-昭通断裂带,未来具有发生较大地震的可能,值得关注.