软粘土卸载特性计算模型研究

本文从三个方面对软粘土的卸载特性进行了理论研究。采用太沙基一维固结理论，编制了内嵌固结计算的修正剑桥模型程序．对软粘土卸载强度进行了分析。推导了基于修正剑桥模型的孔压表达式并编制了相应的程序对软粘土的卸载孔压特性进行了研究，结果显示，在伸长条件下．孔压先由零逐渐发展到负的最大值．然后又逐渐减小至零。突破零位后，发展成为正孔压并最后达到正的最大值。围压越大，产生的负孔压的极值越大，最终的正孔压越小。将传统边界面模型的线弹性卸载过程改进为弹塑性．建立了软粘土的边界面广义弹塑性模型，从而使边界面模型可用于分析软粘土卸载过程中的塑性变形问题。理论分析结果与试验进行了比较验证，证明理论研究方法是可行的。     

铝材轧制过程中润滑添加剂吸附特性研究

依据铝材轧制过程中润滑油添加剂浓度对摩擦因数的影响，同时以Ｔｅｍｋｉｎ吸附理论为基础，设定添加剂分子在摩擦表面的覆盖率与摩擦因数的联系，提出了铝材轧制过程中润滑添加剂吸附有的求解方法，还定义了铝材轧制时润滑添加剂饱和吸附最低浓度的概念，通过二辊轧机测定了在不同添加剂浓度时纯铝材轧制过程中的摩擦因数，并且计算了相应的润滑添加剂的吸附自由能及其饱和吸附最低浓度，表明理论分析结果与实测结果具有良好的一致     

GENERALIZED FINITE SPECTRAL METHOD FOR 1D BURGERS AND KDV EQUATIONS

A generalized finite spectral method is proposed. The method is of high-order accuracy. To attain high accuracy in time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme was used. To avoid numerical oscillations caused by the dispersion term in the KdV equation, two numerical techniques were introduced to improve the numerical stability. The Legendre, Chebyshev and Her-mite polynomials were used as the basis functions. The proposed numerical scheme is validated by applications to the Burgers equation (nonlinear convection- diffusion problem) and KdV equation (single solitary and 2-solitary wave problems), where analytical solutions are available for comparison. Numerical results agree very well with the corresponding analytical solutions in all cases.     

基于CUDA的有限元矩阵并行装配算法研究

构建航天飞行器的结构有限元模型是准确模拟飞行仿真、完成飞行器在轨飞行阶段结构故障监测和诊断的基础。采用细长体飞行器简化梁模型，提出新的基于CUDA（Compute Unified Device Architecture）的有限元单元刚度矩阵生成和总刚度矩阵组装算法。依据梁单元矩阵的对称性，结合GPU硬件架构提出并行生成算法并进行改进。为有效减少装配时间，在装配过程中采用着色算法，提出了基于GPU（Graphics Processing Unit）共享内存的非零项组装策略，通过在不同计算平台下算例对比，验证了新算法的快速性。数值算例表明，本文算法的求解效率较高，针对一定计算规模内的模型可满足快速计算与诊断的实时性要求。     

阶梯形变截面圆拱和圆环稳定计算的有效方法

利用等截面圆弧曲杆的稳定微分方程和本文导出的弯矩一般表达式，建立了静水压力作用下阶梯变截面圆拱和圆环稳定分析的代数方程组，得到了解析计算临界载荷的特征方程．算例表明本文方法简单有效     

剪切弯曲下短深梁位移数值计算精度的研究

分析了短深梁位移的基本解，评述了有限元解法中关于高次分布力简化为节点力的问题。使用虚功原理和位移线性插值推导了满足二次分布力作用的精确节点力分配公式。通过矩形截面悬臂梁自由终端受集中载荷的算例，阐明了提高短深梁位移计算精度的途径。     

用于先进半导体制程的光刻反向计算技术(ILT)

此篇论文将介绍一个用于半导体光罩上图样设计以及可用于实际生产的光刻反向计算技术（ILT）。在论文中将讨论有关ILT的最新发展，包括在超成像极限协助图样（SRAF）的生成，可增加制程宽容度的ILT，以及如何生成满足光罩生产标准的图样等方面。从内部的研究结果和客户的使用结果可以看出，ILT已经不再只是一种用于研究的工具，而是已经可以用于先进半导体制程的大规模生产。在对各个环节优化之后，ILT可以增加制程的宽容度，同时将光罩的成本控制在可以接受的水平。     

单向悬挂屋盖结构的风致气弹耦合效应数值模拟

综合运用计算流体力学和计算结构力学技术,建立了适用于索膜结构流固耦合风振分析的CFD数值模拟方法,并编制了相应的有限元计算程序。应用该程序对具有不同参数的大跨度单向屋盖结构进行了风振响应分析,探讨了来流风速、屋面质量和初始预张力等参数对结构流固耦合特性的影响。研究表明,由于屋盖前缘周期脱落的大尺度旋涡是诱导结构振动的主要原因,随着大尺度旋涡沿屋面向下游移动,其动能逐渐转化为结构变形能,使结构振动形态呈现明显的涡激振动特征。屋面质量、来流风速和初始预张力是影响结构流固耦合性能的主要参数。对于索膜结构,仅通过刚性模型来确定结构风压并进行风振响应分析,其结果是不可靠的。     

建筑风压场数值模拟中压力波动的平抑

在以同位网格为基础的简单流场压力计算中,通常采用动量插值方法来平抑流场中的压力波动现象;但是对于建筑风场等复杂的钝体绕流问题,由该平抑方法得到的收敛风压场仍可能存在小幅波动。为彻底解决同位网格格式下的压力波动,除采用动量插值方法外,本文提出了在压力校正方程的界面流速中添加压力梯度差值项的方法。算例分析表明,该方法计算得到的建筑风压场完全避免了压力波动现象,风压解与试验结果吻合良好。     

斜交连续梁桥数值分析研究

在分析了矩阵位移法和传递矩阵法的理论基础上，将二者组合而成矩阵混合法。矩阵混合法的实质就是将结构物的一端未知向量通过矩阵的简单乘法传递到另一端，再将已知条件代入求解的过程。在矩阵传递的过程中引入位移法中的前进代入的概念，即将未知向量不断地以新的未知向量替换，可以防止误差的传递和积累。将矩阵混合法运用于斜交连续梁桥内力计算，完成了矩阵混合法在斜交连续梁桥内力计算中的理论推导和程序的设计，给出了实际算例并与有限元法进行了比较。