弹性矩形板方程在非线性边界条件下整体解的存在唯一性

考虑了在非线性边界条件下的弹性矩形板方程.利用Galerkin方法,首先证明了该方程在非线性边界(a)及初值w⁰∈W,w¹∈W的条件下初边值问题存在唯一整体弱解w(t).其次证明了该方程在非线性边界(b)及初值w⁰∈W₁,w¹∈W₁的条件 关键词： 弹性矩形板方程 非线性边界条件 初边值问题 整体解     

平面撞击流偏斜振荡的实验研究与大涡模拟

采用实验和大涡模拟对喷嘴出口雷诺数(Re= U₀ hρ/μ, 其中 U₀为出口平均速度, h为平面喷嘴出口狭缝高度, ρ和 μ分别为流体的密度与动力黏度)为25–10000, 喷嘴间距 L为4h–40h范围内的平面撞击流偏斜振荡特性进行了研究. 通过对平面撞击流模拟和实验的结果进行比较, 验证了数值模拟的可靠性, 并对平面撞击流发生偏斜振荡的无因次参数(喷嘴间距 L/h与出口雷诺数 Re)范围进行划分, 重点考察了湍流平面撞击流的偏斜振荡周期及速度-压力变化特征. 研究结果表明大涡模拟能对平面撞击流的偏斜振荡进行有效预报; 当平面撞击流发生周期性偏斜振荡时, 特定点的压力与速度也发生周期性变化, 且变化周期与偏斜振荡周期一致, 偏斜振荡本质上是由速度-压力的周期性变化和转换引起的. 关键词： 平面撞击流 偏斜振荡 大涡模拟     

高压下CaPo弹性性质和热力学性质的第一性原理研究

利用密度泛函理论的平面波赝势方法预测研究了CaPo从岩盐结构(B1结构)到氯化铯结构(B2结构)的相变以及B1结构CaPo高压下的弹性性质以及热力学性质等.通过等焓原理发现B1→B2的相变压力为22.8GPa. 同时计算了B1结构CaPo高压下的弹性常数以及剪切模量、杨氏模量等相关弹性参数,结果发现当压力超过20GPa时,B1结构CaPo开始不稳定了,这和等焓原理所得结果相符合. 最后通过Debye模型成功获取了B1结构C 关键词： 相变 弹性性质 热力学性质 CaPo     

氘含量对DKDP晶体横向受激拉曼散射增益系数的影响

本文测量了不同氘含量K(H_1-xD_x)₂PO₄晶体(DKDP晶体)在Z(XX)Y散射配置下的自发拉曼散射光谱, 并详细分析了氘含量对与横向受激拉曼散射(TSRS)增益系数有关的拉曼频移、半峰宽和散射强度的影响. 然后通过与去离子水拉曼散射对比得出了不同氘含量DKDP晶体的TSRS增益系数. 结果表明随着氘含量的增加DKDP晶体的TSRS增益系数先减小至KDP晶体的40.1%, 后增大至68.9%; 本文认为掺氘后拉曼半峰宽的变化是引起TSRS增益系数随氘含量变化的主要原因.     

描述高自旋弱耦合系统的广义积算符

本文提出了Z算符与角动量积算符相结合的广义积算符,用于描述高自旋弱耦合系统NMR实验谱。给出了高自旋系统在自旋耦合作用下的演变公式,Z算符在脉冲作用下的变化公式。以氘(D)为例,计算了重聚INEPT,DEPT和同核COSY实验谱,并对计算方法作了讨论。     

自锻破片形成的二维数值模拟

自锻破片装药(简称SFF)是从聚能破甲装药派生出来的一种新的反坦克装药结构。描述SFF的形成过程,我们采用二维轴对称弹性-理想塑性-流体力学方程组。基本数值方法采用M.L.Wilkins的迴路积分法和Legrange四边形网格。数值模拟得到的SFF完整的形成过程和试验中x射线高速摄影照片基本相符。数值模拟中的若干方法简化了计算,节省了机时。     

二氧化铀电子结构和弹性性质的第一性原理研究

本文采用第一性原理的方法系统研究了UO₂的晶体结构、电子结构和弹性性质. 在计算中采用广义梯度近似结合Hubbard U项描述电子的局域强关联效应. 首先通过计算能带带隙大小并与理论值比较的方法, 得到了合理的有效库仑相关作用能(U_eff)的取值, 同时通过态密度的计算, 进一步验证了U_eff取值的合理性. 计算得到UO₂中U原子的U_eff值为3.30 eV (U_eff=U-J, U=3.70 eV, J=0.40 eV). 应用此参数计算得到的UO₂晶格常数为5.54 Å, 带隙宽度为2.17 eV. 该结果优于目前现有的研究结果, 同时在同样的U_eff值条件下计算所得到的弹性常数与实验值也符合得较好. 相较于之前的基于实验测量并分析得到的U_eff值, 我们所采用的方法在对UO₂性质描述上更为准确. 不同的有效库仑相关作用能取值下的态密度结果表明, 有效库仑相关作用能的大小可以影响铀原子5f电子轨道的分布.     

离子碰撞过程中电荷转移的理论研究

用经典轨道蒙托卡罗(CTMC)方法计算离子与氢原子碰撞过程电荷转移截面,并计算俘获电子的n(n为主量子数)壳层分布。对于多电荷离子的碰撞过程,利用约化变量得到了截面的q标度规律。计算结果与实验测量符合得很好。     

大规模宽频弹性动力学分析的快速定向压缩边界元法

发展一种大规模宽频弹性动力学分析的快速定向压缩边界元法.证明弹性动力学核函数具有定向低秩特性，为采用快速定向压缩算法提供理论基础.根据S波的波数，将节点之间的相互作用划分为低频相互作用和高频相互作用，并将后者进一步划分为与多个楔形区.在楔形区上，可以采用核函数的定向低秩特性进行快速计算.低频相互作用与核无关快速多极边界元法中计算方法相同，不同方向楔形区上的变换矩阵可以采用坐标系旋转的方法进行快速计算.可对任意频率进行快速谐响应分析.数值算例表明：该方法可以将宽频弹性动力学问题计算复杂度降低到O（N log^αN）.与卷积求积边界元法相结合，也可以应用于弹性动力学瞬态分析.     

阶矩阵及其在传统预处理方法中的应用

本文应用矩阵元素阶和阶矩阵概念,讨论了ICCG和MICCG这两种传统的预处理方法在实用中的一些问题。为什么ICCG(s,t)在s+t固定时取(s,t)=(1,1),(1,2),(1,3),(2,4),(3,5),…有较高的收敛速度?为什么MICCG(m)当m>3时迭代次数不变?ICCG和MICCG的填入方式如何系统化?MICCG是否总比ICCG收敛速度高?本文拟作一个初步的讨论。通过LU分解的阶矩阵,本文给出了按阶递增的填入原则,将ICCG和MICCG系统化为P阶ICCG和P阶MICCG,并讨论了MICCG原有填入方式存在的问题。应用误差阵的阶矩阵,本文讨沦了MICCG迭代参数选取中存在的问题,给出了合理的参数选取方法。通过不同算例,本文还比较了ICCG和MICCG的计算效率。     

L12型铝合金的结构、弹性和电子性质的第一性原理研究

运用第一性原理方法研究了L1₂型铝合金相Al₃Sc和Al₃Zr的晶体结构、电子结构和弹性.结合能和形成能的计算表明,两种合金具有较强的合金化能力,且Al₃Zr较Al₃Sc具有更强的结构稳定性.电子结构分析表明,费米能级以下较多的价电子数决定了Al₃Zr具有较强的结构稳定性.计算并分析比较了两种合金相的单晶弹性常数（C₁₁,C₁₂和C₄₄）以及多晶弹性模量（体弹性模量B、剪切模量G、杨氏模量Y、泊松比ν和各向异性因子A）.通过对比实验和其他理论计算结果,进一步分析和解释了两种合金相的力学性质. 关键词： 铝合金 第一性原理 结构和电子性质 弹性     

通孔尺寸对铜互连应力迁移失效的影响

在200℃温度下进行了700h双层铜互连(M1/M2)的应力迁移加速老化试验, 结合有限元分析和聚焦离子束(focused-ion-beam,简称FIB)技术研究了通孔直径分别为500和350nm的铜互连应力诱生空洞失效现象, 探讨了应力诱生空洞的形成机理, 并分析了通孔尺寸对铜互连应力迁移的影响. 结果表明,M1互连应力和应力梯度在通孔底部边缘处达到极大值. 应力梯度在应力诱生空洞成核过程中起主导作用, 由张应力产生的过剩空位在应力梯度作用下沿Cu M1/SiN界面作扩散运动并在应力梯度极大值处成核生长成空洞. 由于M1互连应力沿横向方向变化较快, 因此应力诱生空洞的横向生长速率较大. 当通孔直径增大时,互连应力和应力梯度值增大, 并导致应力诱生空洞的生长速率上升. 关键词： 铜互连 应力迁移 应力诱生空洞 失效     

超高速碎片云撞击下薄板变形与撕裂建模研究

为研究多层板结构中薄板在碎片云作用下的变形与破坏问题,开展了超高速撞击多层板实验。实验结果表明,薄板在高速碎片云冲击下的典型破坏特征为中央穿孔及环孔凹陷变形与花瓣型撕裂。在此基础上,考虑弯矩和膜力作用,建立了描述薄板在轴对称分布强冲击载荷作用下大变形的理想刚塑性环板模型,据此可以计算环板变形过程的横向与径向速度场,结合Grady破碎理论,可以计算花瓣型撕裂的花瓣数,理论计算值与实验比较吻合。研究结果可以为多层板结构在超高速弹丸撞击下的毁伤评估提供理论基础。     

K3[Cr(CN)6]中Cr3+电子跃迁几率(I)

本文首先计算了K₃[C_r(CN)₆]中C_r³⁺在外加直流磁场H绕晶体形态轴(a,b,c)情况下能级的分裂,然后计算相应的波函致,最后计算了电子极大受激相对跃迂几率。     

Al的微观组态与LaNi3.75Al1.25的结构和弹性第一性原理研究

在全电子水平上，采用广义梯度近似密度泛函理论和全势能线性缀加平面波方法并结合二维立方拟合方法，对LaNi_3.75Al_1.25合金的晶体结构与弹性性质进行了理论研究.计算结果给出合金的晶格常数a=b=0.5137　nm，c=0.4018　nm，Al原子在晶胞中的微观分布为同时占据部分3g和2c等价格位，弹性常数C₁₁+C₁₂=281.2，C₁₃=82.3，C₃₃=227.3，以及体弹性模量B＝124.5、切变模量G＝68.2　GPa.还对态密度、能带结构和电荷密度进行了计算分析，并给出材料LaNi_3.75Al_1.25的电子线性比热系数23.45　mJ/molK². 关键词： 3.75Al_1.25')" href="#">LaNi_3.75Al_1.25 储氢合金 全势能线性缀加平面波 弹性常数     

中子输运方程中的α本征值计算

用多群Sn方法(离散纵标法)独立编制了α本征值计算程序,在一定程度上解决了次临界情形下α本征值计算遇到的困难,并对含空腔系统给出了一个简化处理模型.对标准问题的计算表明,该程序算法可靠,计算结果可信.     

第一性原理计算α-ScDx(D=H,He)的弹性常数

本文采用第一性原理方法研究了α相钪氢化合物(α-ScH _x )和α相钪氦化合物(α-ScHe _x )的弹性常数,其中x=0,1/4,1/8,1/32表示H原子和He原子在α-ScD_x(D=H,He)中的浓度.研究结果表明,对于钪,氢和氦对晶体的弹性性质的影响截然不同.钪氢化合物晶体的弹性常数基本上随着晶体中氢的浓度增加而增加,这跟实验测量得出的结论一致;然而,钪氦化合物体系的弹性常数几乎随着氦浓度的 关键词： 第一性原理 弹性常数 x')" href="#">α-ScH _x x')" href="#">α-ScHe _x     

连续态数值波函数归一化的近似计算及TRK求和中余项的积分处理

采用改进的Numerov格式,对连续态波函数及其归一化问题作了细致的计算.指出了Cowan在ε→0时状态波函数的归一化中存在的问题,给出了更加准确、有效的归一化计算目的.为了检验这一算法,计算了氢原子各状态(l=0~3)下的Thomas-Reiche-Kuhn(TRK)和.利用最小二乘法,给出了TRK求和中无穷项的渐近函数形式,并由此采用积分目的得到了余项的近似计算结果.     

第一性原理计算MAlH4(M=Na,K)的结构和弹性性质

基于密度泛函理论研究并比较了两种储氢材料NaAlH₄和KAlH₄的晶格参数,弹性性质和电子结构.计算结果表明NaAlH₄和KAlH₄都是绝缘体,Al—H之间是共价键,M(M=Na,K)与AlH₄之间是离子键.NaAlH₄和KAlH₄的弹性常数都比较小而且NaAlH₄的弹性常数要大于KAlH₄的弹性常数,对此给予了解释. 关键词： 配位金属氢化物 电子结构 弹性性质     

空间碎片超高速撞击充气压力容器前壁准静态破坏分析

针对空间碎片超高速撞击充气压力容器前壁发生准静态破坏问题,将其简化为受双向拉应力的圆孔边双裂纹的线弹性断裂问题进行处理;并在数值模拟及理论分析的基础上建立了充气压力容器前壁发生准静态破坏的预报模型,得到了当球形弹丸撞击速度为7.0 km/s时、壁厚为1.0 mm的Al5754圆柱形压力容器前壁发生准静态破坏的临界应力曲...