基于FDM-DEM耦合的冲击损伤大理岩静态断裂力学特征研究

为探究循环冲击损伤后大理岩的静态断裂力学特征,基于有限差分（finite difference method,FDM）-离散元（discrete element method,DEM）耦合的建模技术构建了三维分离式霍普金森压杆（split Hopkinson pressure bar,SHPB）数值模型,其中杆件系统和岩石试件分别采用FLAC3D和PFC3D程序建模。利用该模型对中心直切槽半圆盘（NSCB）试样进行了恒定子弹速度下的循环冲击,随后对受损试样进行静态三点弯曲断裂实验。通过编写Fish程序,提取试样断裂面数据,对断裂面进行重构并定量计算表面粗糙度。通过与相关室内实验结果的对比分析,验证了本文数值分析的合理性与可靠性。模拟结果表明,随着循环冲击次数的增加,试样内部微裂纹、破碎颗粒均增加。连接力场分布混乱,部分力链发生断裂。力链的变化是试样力学性能劣化的根本原因。在静态三点弯曲断裂实验中,冲击5次后试样的静态断裂韧度较天然试样产生一定程度的降低。试样在静载过程中产生的微裂纹和碎块的数量随循环冲击次数的增加而增加,断裂面粗糙度随循环冲击次数的增加而增加。     

有红利美式看跌期权定价的Crank-Nicolson有限差分法

本文基于B-S微分方程,采用Crank-Nicolson差分格式(简称C-N差分格式)求解支付固定红利的美式看跌期权价值,给出实证分析,并对C-N差分格式和隐含的差分格式进行了比较.结果表明,用C-N差分格式可以得到更加精确、有效的数值解.     

二维变系数反应扩散方程的紧交替方向差分格式

1 引言 在研究热传导过程、气体扩散现象和电磁场的传播等问题时,常常遇到抛物型偏微分方程。用有限差分方法研究这类问题的数值解法目前已经有了许多工作。对于二维、三维抛物方程的数值求解比较理想的方法是交替方向法。     

考虑热流固耦合作用的多孔介质孔隙尺度两相流动模拟

为了研究深层油气资源在岩石多孔介质内的运移过程, 使用一种基于Darcy-Brinkman-Biot的流固耦合数值方法, 结合传热模型, 完成了Duhamel-Neumann热弹性应力的计算, 实现了在孔隙模拟多孔介质内的考虑热流固耦合作用的两相流动过程. 模型通过求解Navier-Stokes方程完成对孔隙空间内多相流体的计算, 通过求解Darcy方程完成流体在岩石固体颗粒内的计算, 二者通过以动能方式耦合的形式, 计算出岩石固体颗粒质点的位移, 从而实现了流固耦合计算. 在此基础上, 加入传热模型考虑温度场对两相渗流过程的影响. 温度场通过以产生热弹性应力的形式作用于岩石固体颗粒, 总体上实现热流固耦合过程. 基于数值模型, 模拟油水两相流体在二维多孔介质模型内受热流固耦合作用的流动过程. 研究结果表明: 热应力与流固耦合作用产生的应力方向相反, 使得总应力比单独考虑流固耦合作用下的应力小; 温度的增加使得模型孔隙度增加, 但当注入温差达到150 K后, 孔隙度不再有明显增加; 温度的增加使得水相的相对渗流能力增加, 等渗点左移.      

有限图上高阶Yamabe型方程的非平凡解

该文研究了以下高阶Yamabe型方程Lm,pu - g|u|p-2 u = λ f|u|α-2 u在有限图上的非平凡正解的存在性,其中Lm,p是一个2m阶差分算子,它是一种p次(-Δ)m算子更一般化,α≥p≥2,g＞0和f ＞0是定义在G的所有顶点上的实函数,m≥1是一个整数.     

二维光子晶体异质结构的带隙分析

提出异质结构的二维光子晶体模型,用时域有限差分法计算其带隙,数值计算结果表明,其带隙宽度变大.     

变温场中具损伤粘弹性矩形板的非线性动力响应分析

基于热粘弹性理论、Von Karman板理论和连续损伤力学,导出了二维状态下各向同性材料的变温粘弹性本构方程,建立了含损伤效应的各向同性粘弹性矩形板在变温场中的非线性运动控制方程,且应用有限差分法对问题进行求解.算例中,讨论了损伤演化及温度场等因素对粘弹性矩形板非线性动力学行为的影响,得出一些有意义的结论.     

数字图像频域相关分析法及其应用

提出了一种以灰度为载体的数字图像频域相关(DFC)计量的方法,通过傅里叶变换和滤波,直接提取物体位移信息。与传统的激光散斑法比较,无需干涉光源,可用白光照明,因此对实验环境的要求很低;与数字图像相关法比较,计算速度更快,对照明要求更低,不受刚体转动影响。新方法不仅设备简单、容易实现,还具有灵敏度可调,测量过程无需人工干预的良好特征,由于充分利用了计算机数字图像处理的功能,抗干扰性强,是一种适合工业现场测量的有发展前途的位移全场测试新技术。     

有限深水域中三维频域自由面格林函数及其导数的计算方法

三维频域自由面格林函数及其偏导数的数值计算是海洋结构物水动力分析的难点.本文对有限水深格林函数及其偏导数的数值算法作出两点改进:①对原函数的级数表达式进行变形,提出适用于数值计算的形式,解决高频率情况下级数解的计算失真问题;②对Gauss-Laguerre积分法作出改进,解决高频率大水深情况下积分解的数值溢出问题.计算...