峨山褐煤的分子结构和分子模拟

作为重要的化石能源，褐煤资源潜力巨大、分布广泛但综合利用率低。研究褐煤的分子结构模型，有助于预测褐煤在热解、液化和气化过程中的化学反应机理及反应路径，进而提高褐煤的综合应用水平。以云南峨山褐煤为研究对象，利用傅里叶变换红外光谱、13C核磁共振波谱及X射线光电子能谱等分析测试方法，获取了峨山褐煤的含碳、含氧及含氮结构参数。在此基础上，借助Gaussian 09计算平台，采用量子化学建模的方法构建并优化了峨山褐煤的分子结构模型。研究结果表明：峨山褐煤的芳碳率为39.20%，芳香碳结构主要为苯和萘，且芳香桥头碳与周边碳的比值χ_b为0.07；脂碳率为49.51%，脂肪碳结构主要为亚甲基，季碳和氧接脂碳；氧原子主要存在于羟基、醚氧、羰基和羧基结构中；含氮结构则以吡啶为主。基于元素分析、13C 核磁共振波谱分析，又经过热重实验消除褐煤中残余水分的影响后，计算出峨山褐煤的分子式为C₁₅₃H₁₃₇O₃₅N₂。依据分子式及分析结果计算出峨山褐煤的结构单元含量并构建出其初始结构模型，采用半经验法PM 3基组及密度泛函理论M06-2X/3-21G基组对初始分子构型进行优化。优化后的分子模型具有明显的三维立体特征，芳香环之间较为分散且在空间中排列不规则，芳香簇主要通过亚甲基、醚氧基、羰基、酯基和脂肪环连接，含氧官能团主要分布在分子边缘，脂肪族侧链较多。对优化后的分子模型进行振动频率计算进而获得了分子模型的模拟红外光谱，其与实验红外谱图吻合度良好，证明了峨山褐煤分子结构模型的准确性、合理性。分子结构模型的构建有利于直观地了解峨山褐煤的分子结构特征，从而有助于从微观分子角度研究峨山褐煤的宏观性质。同时，峨山褐煤分子结构模型可为其在热解、液化和气化等领域研究中提供理论指导。     

配合参数对单晶炉驱动系统的影响分析

光伏产业的发展使得对硅材料的需求日益增加,同时硅单晶生产行业竞争也日趋激烈。作为生产硅单晶的重要装备,单晶炉的稳定性和可靠性关系到硅单晶生产效率的提升和成本的下降,因此其驱动系统的设计和优化成为装备制造的关键环节。本文以NVT-HG2000-V1型硅单晶生长炉的驱动系统为研究对象,用SolidWorks三维建模实现虚拟装配,采用ADAMS建立其动力学仿真模型,并对驱动系统的运动过程进行仿真模拟。采用控制变量法定量分析了铜套与升降轴的配合间隙及丝杠参数对驱动力和驱动力矩的影响规律,进而在提高硅单晶生长炉装备稳定性和可靠性方面给出合理的技术建议。结果表明,铜套与升降轴的配合间隙达到0.071 mm后能有效降低驱动系统运行所需驱动力矩,丝杠倾斜度、螺纹螺距与螺纹间摩擦系数的增大均会导致驱动系统运行所需力矩大幅增加。     

一种排气式膛压模拟装置

为了在实验室条件下简洁、高效地获得与实际相符的膛压曲线，进而开展典型结构和材料膛压载荷响应特性研究，提出了压力舱内发射药燃烧同时发射药气体由排气件排出的膛压模拟装置。结合发射药燃烧理论和等熵流动模型，建立了排气式膛压模拟过程的数学模型。基于理想气体假设，利用Fluent软件模拟泄压过程质量流量规律，并与理论结果对比，确定了流量系数。分别根据76和155 mm火炮膛压曲线特点及小型化设计原则，对模拟装置性能参数进行了优化设计。优化结果表明，获得的压力曲线的增压速率和降压速率基本满足要求，峰值压力达到300 MPa，压力大于30 MPa历时10 ms以上。验证实验结果表明:压力曲线有良好的重复性，且与理论结果一致，装置工作可靠性高；以排放发射药气体方式模拟膛压曲线是可行的。     

基于Caputo导数的分数阶非线性振动系统响应计算

研究了含分数阶Caputo导数的非线性振动系统响应的数值计算方法。首先，由Caputo分数阶导数算子的叠加关系，得到含分数阶导数项非线性振动系统状态方程的标准形式。其次，基于Caputo导数与Riemann-Liouville导数和Grunwald-Letnikov导数间的关系，推导计算了Caputo导数的一般数值迭代格式。本文方法不要求状态方程中各分数阶导数阶数相等，弱化了已有算法中对分数阶导数阶数的限制，并可推广到多自由度的情形。随后，选择若干有解析解的算例验证了本文方法的正确性。最后，以多吸引子共存的分数阶Duffing振子系统为例，比较Caputo和GL两种算法所得结果，说明了用GL算法求解存在的问题。     

菱形颗粒冲蚀磨损特性试验及仿真研究

本文中针对单个硬质角形颗粒冲击金属材料表面的过程,设计了弹射试验装置,研究菱形颗粒冲击行为及冲蚀机理.采用高速摄像机,捕捉不同冲击速度v_i、冲击角度α_i和方位角度θ_i下颗粒的运动轨迹.建立了基于拉格朗日法的FEM-SPH耦合数值计算模型,借助于模型进一步分析了角形颗粒的运动学行为和变形凹坑形态.结果表明:冲击角α和方位角θ是决定颗粒旋转的关键因素,在某一固定冲击角αi下存在一个临界方位角θcr_i,当θiθ_(cri)时颗粒冲击后发生前旋旋转,当θ_iθ_(cri)时颗粒冲击后发生后旋旋转;冲击诱导的颗粒旋转对冲蚀机理的影响较大,颗粒前旋旋转对金属材料产生"耕犁"作用,后旋旋转对金属材料产生"撬起剔除"作用.颗粒的动能损失受到冲击角α_i和方位角θ_i的影响较大,临界方位角θ_(cri)下颗粒的动能损失最大,凹坑变形最严重.     

浅水区海底管道周围海床渗透压力计算研究

我国海上油田开采起步较晚，大部分油田处于浅水区，因此，在设计管道时，应充分考虑由浅水区波浪引起的管道周围海床渗流力。根据浅水波相关假设，考虑自由水面非线性影响，推导出椭圆余弦波的波面方程，在此基础上进一步得到一个关于速度势的表达式，并根据该表达式得出作用于海床表面的波压公式。考虑海床土的压缩性，推导出一阶近似椭圆余弦波作用下浅水区埋置管道周围海床的渗流压力解析解，最后将计算结果与大型水槽试验及以往研究成果作对比。结果表明，在椭圆余弦波的作用下，由一阶椭圆余弦波理论得到的计算结果与试验结果规律基本一致，与相似工况下的现有理论成果数值基本相同，具有一定的可行性和工程价值。     

双侧下颌支矢状劈开截骨术后患者可施加合适咬合力的数值研究

双侧下颌支矢状劈开截骨术在临床手术中广泛用于治疗缩颌、凸颌和颌面不对称患者的颌面畸形症状。本文用数值方法研究双侧下颌支矢状劈开截骨术并前移下颌骨后患者可施加的合适咬合力。首先，用缩颌患者的CT断层扫面数据建立下颌骨、关节盘、关节窝和相关软组织的高度几何相似性的几何模型；然后，用有限元模型模拟双侧下颌支矢状劈开截骨术后患者的咬合过程。通过分析发现，术后关节盘上的位移处于正常的范围，而应力却远远超过正常范围，因此，将应力因素作为确定术后患者可施加的合适咬合力的标准。当施加正常咬合力的25%时，关节盘上的应力分布和术前没有显著性差异，同时下颌骨上的应力水平也处于正常范围。     

求解矩阵方程AXB+CXD=F参数迭代法的最优参数分析

本文针对求矩阵方程AXB+CXD=F唯一解的参数迭代法，分析当矩阵A，B，C，D均是Hermite正（负）定矩阵时，迭代矩阵的特征值表达式，给出了最优参数的确定方法，并提出了相应的加速算法.