高温作用下混凝土热-水-力耦合损伤分析模型

以混合物理论为基础建立了高温作用下混凝土的热-水-力耦合损伤分析模型．将混凝土视为由固体骨架、液态水、水蒸气、干燥气体和溶解气体共5种组分构成的混合物,模型的宏观平衡方程包括各组分的质量守恒方程、整体的能量守恒方程及动量守恒方程,模型所需的状态方程及本构关系全部给出,最后给出基于4个主要参数（固体骨架位移、气压力、毛细压力和温度）的控制方程．模型考虑了混凝土在高温作用下,水分的蒸发与冷凝、胶结材料的水化及脱水、溶解气的溶解与挥发等相变过程;从材料变形破坏过程中能量耗散特征入手,基于Lemaitre应变等价性假说和能量守恒原理得到力学损伤演化方程,并考虑了高温引起的热损伤对材料力学性能及力学损伤演化规律的影响,建立了热-力耦合损伤本构模型．     

基于∨-·运算的Fuzzy矩阵方程的摄动问题

定义基于∨ - .运算的 Fuzzy矩阵 A的同解简化矩阵 ,利用 A的同解简化矩阵 A(1) 给出基于∨ - .运算的 Fuzzy矩阵方程的矩阵解法 ,并研究这类 Fuzzy矩阵方程的摄动问题。     

随机双曲折现下的消费与最优投资策略探究

在本文中,我们主要讨论在随机双曲折现中,当投资者的消费行为是一个布朗运动时,其用于投资风险资产的最优投资额.基于汉密尔顿-雅克比-贝尔曼方程计算常数绝对风险厌恶型投资者的效用函数下的最优投资组合,并给出这组方程的近似解.在此基础上,我们分析当消费服从Wiener过程时风险资产投资的一些重要性质,研究消费行为与风险资产投资行为之间的关系.     

通过求解输运方程计算壁面距离

壁面距离在当代湍流模化中仍然扮演着关键角色,然而苦于遍历计算壁面距离的高昂代价,该文考虑了求解偏微分方程的途径．基于Eikonal方程构造出类Euler形式的输运方程,这样,可以直接利用求解Euler和Navier-Stokes方程的CFD程序使用的高效数值格式和部分代码．基于北航的MI-CFD(CFD for missles)数值平台,详尽地介绍了该输运方程在直角坐标下的求解过程．使用隐式LUSGS时间推进和迎风空间离散,发现该方程具有鲁棒快速的收敛特性．为了保证精度,网格度量系数必须也迎风插值计算．讨论了初始条件和边界条件的特殊处理．成功应用该壁面距离求解方法计算了几个含1-1对应网格和重叠网格的复杂外形．     

力学系统运动方程的矩阵表示法

本文阐述以矩阵形式建立力学系统运动方程的方法.采用矩阵形式可使动力学方程写得紧凑,以便于用电子计算机研究多自由度的力学系统.基于动力学基本规律、达朗倍尔-拉格朗日原理、哈密顿-奥斯塔拉格勒特斯基原理、以及高斯原理,列出了建立运动方程的方法.     

具有初应力的1-3型压电复合材料的横向共振分析

1-3型压电复合材料在极化处理过程中会出现极化残余应力,基于线性压电理论,采用解析方法研究了极化残余应力对1-3型压电复合材料的横向共振模态的影响．首先建立含有初始应力的平面波波动方程,基于Bloch波函数理论构建了位移和电势函数,从而最终得到方程的解．数值结果表明初始残余应力的出现降低了复合材料的横向振动频率,而材料的压电性则使得横向振动频率提高．初始剪应力对对横向振动的频率影响较小,可以忽略．第一阶及第二阶振动频率将随着初始应力的大小呈线性变化．得出的结论对超声换能器的制作和研究将提供有意义的指导作用．     

堆浸工艺中流动-反应-变形-传质耦合过程数值模拟及应用

基于堆浸过程中的孔隙变形,发展了一孔隙发育模型,并建立了堆浸工艺中流动-反应-变形-传质全耦合模型控制方程组,在喷淋强度不变和常水头的的情况下,给出了方程的数值解.数值结果描述了浸出过程矿堆中溶浸液饱和度的分布、溶浸剂浓度和浸出矿石浓度的分布,拟合结果表明:有用矿物浸出百分比与浸出时间之间呈3次关系,这一关系可用来预测堆浸中的有用金属的浸出回收率.     

分式方程中的字母参数

<正>初三总复习第二部分方程(组)与不等式(组),遇到这样一道题目:若关于x的方程(ax+1)/(x-2)=-1的解是正数,则a的取值范围是____.1题目的解答方法与过程1.1多数同学解题的方法与过程方法描述先将方程化简,因为题目中所给关于x的方程(ax+1)/(x-2)=-1的解是正数,所以想到先求出方程的解,再建立不等式,即这     

滞时Volterra积分方程数值方法的数值稳定性分析

本文给出数值方法解Volterra积分方程的稳定性分析,我们判定可约积分方法的数值稳定性基于如下试验方程其中τ是正常数,p和q是复值的。在上述试验方程的情况下,我们研究θ-方法及可约积分方法的稳定性。     

稳态Poisson-Nernst-Planck方程的残量型后验误差估计

针对稳态的Poisson-Nernst-Planck方程研究了一种残量型的后验误差估计子, 对方程的两个解-浓度和电势, 都分别给出了上界和下界估计. 数值实验表明, 基于这种后验误差估计子构造的自适应有限元算法对于稳态的Poisson-Nernst-Planck方程是有效的.     

电池系统建模中Butler-Volmer方程的同伦分析求解

Butler-Volmer方程是电化学系统中描述电极动力学过程的本构方程,具有强非线性．为了对这一方程(耦合两个Ohm方程)进行解析求解,在同伦分析方法的框架下,发展了满足简单条件的广义非线性算子的算法,以取代原同伦分析中的非线性算子．该广义非线性算子的构造保证了高阶形变方程的线性特征．这一方法的有效性通过一些算例得到了验证．最后通过同伦分析方法对Butler-Volmer方程进行了求解,结果显示过电位和电流密度的级数解析解与数值解吻合很好,并有很好的收敛效率．     

一种求解变系数渗压固结微分方程的数值方法

对变系数渗压固结微分方程的求解过程进行了深入研究,提出一种精细积分半解析数值方法,首先对渗压固结微分方程在空间离散,建立起对于时间的常微分方程组,然后对时程积分,利用矩阵指数函数可以在计算机字长范围内精确计算的特点,得到精密解答.当指数函数用Taylor展开式的一阶近似替代时,精细积分转化为差分方程.用matlab语言编写程序进行求解,得到孔隙比在固结过程中的分布规律,并通过模型试验进行了验证.     

两股票情形下欧式期权定价新方法

假定股票价格服从布朗运动驱动的随机微分方程,从随机动力学的角度出发考虑欧式期权定价问题.由Fokker-Planck-Kolmogrov得到了股票价格过程的概率转移密度函数,基于此,可以求得两股票情形下各种欧式类型未定权益的定价公式.为欧式期权定价提供了一个新方法.     

求解热传导反问题的一种正则化Newton型迭代法

讨论热传导方程求解系数的一个反问题．把问题归结为一个非线性不适定的算子方程后,考虑该方程的Newton型迭代方法．对线性化后的Newton方程用隐式迭代法求解,关键的一步是引入了一种新的更合理的确定(内)迭代步数的后验准则．对新方法及对照的Tikhonov方法和Bakushiskii方法进行了数值实验,结果显示了新方法具有明显的优越性．     

Numerical Analysis of Time Fractional-Order Unsaturated Flow and Its Application北大核心CSCD

非饱和渗流过程的数值模拟对土质边坡稳定性分析、地下污染物迁移模拟等众多领域有着重要的意义。Richards方程由于其普遍适用性被广泛地应用,然而Richards方程所描述的渗流过程并未考虑在自然环境和实验中存在的反常扩散现象。针对这一问题,该文结合Caputo导数得到了具有更广泛渗流意义的时间分数阶Richards方程,采用有限差分法得到其离散格式并采用Picard法迭代求解,以及对分数阶参数和土水特征曲线进行了敏感性分析。最后,结合土柱入渗实验数据,比较了不同土水特征曲线下时间分数阶Richards方程得到的数值解。结果表明,VGM模型的时间分数阶Richards方程与实测数据具有更好的拟合效果,能够更好地描述地下水在非饱和土中的渗流过程。     

基于跳扩散过程的欧式交换期权定价

考虑了股票价格服从带时滞泊松跳的跳扩散模型的欧式交换期权定价问题,运用无套利理论推导出期权价值微分方程,利用变换计价单位的方法,得到交换期权的显示定价公式.     

Numerical Treatment of Homogeneous Semi-Markov Processes in Transient Case–a Straightforward Approach

This paper presents the numerical solution of the process evolution equation of a homogeneous semi-Markov process (HSMP) with a general quadrature method. Furthermore, results that justify this approach proving that the numerical solution tends to the evolution equation of the continuous time HSMP are given. The results obtained generalize classical results on integral equation numerical solutions applying them to particular kinds of integral equation systems. A method for obtaining the discrete time HSMP is shown by applying a very particular quadrature formula for the discretization. Following that, the problem of obtaining the continuous time HSMP from the discrete one is considered. In addition, the discrete time HSMP in matrix form is presented and the fact that the solution of the evolution equation of this process always exists is proved. Afterwards, an algorithm for solving the discrete time HSMP is given. Finally, a simple application of the HSMP is given for a real data social security example.     

Dynamic Modeling of Spatial Cooperation Between Dual-Arm Mobile Manipulators北大核心CSCD

移动机械臂进行空间协作时会产生复杂的非线性耦合,使得采用Lagrange方程或Newton-Euler法直接进行建模极为繁琐。针对双移动机械臂空间协作问题,提出了一种结合Udwadia-Kalaba (U-K)方法与Lagrange方程建立动力学模型的方法。在建模过程中,将负载简化为连杆,选择负载中心断开的方式对系统进行分解,从而避免了机械臂末端关节断开导致的末端关节转角与连杆转角的约束信息缺失问题;将分割形成的两个子系统通过Lagrange方程进行建模,得到了子系统的动力学模型;再将协作系统的固有几何关系通过约束形式引入,应用U-K方法得到了协作系统动力学模型,减少了建立动力学模型所需要的计算量;最后通过数值仿真验证了该方法所得到的动力学模型的准确性。     

随机利率下期权定价

本文在随机利率的基础上,考虑股票价格过程和利率过程分别为扩散过程和Ito过程,并且在相关的假设下,运用鞅方法推导出欧式期权价值过程所满足的微分方程;以及利率满足一种特殊方程时,运用最优停止的鞅方法,得到了随机利率下美式期权的价格和最优停时.     

悬臂梁大挠度问题的双参数摄动解

利用拟线性分析方法中的一阶导数代换对基本方程进行了简化处理,将基本的微分积分方程转化为非线性代数方程组,并结合双参数摄动研究了悬臂梁的大挠度问题．与已有的研究结果比较表明: 将拟线性方法用于研究悬臂梁的大挠度问题,计算较为简便, 同时又具有良好的精度．