基础滑移隔震体系的地震响应计算

利用随机振动理论求出了基础滑移隔震房屋在无频谱过程的白噪声过程作用下的响应特性，探讨了基础滑移隔震的振动时最佳阻尼系数取值，提出了基础最大滑移量在白噪声过程作用下的计算公式，并进行了基础滑移隔震位移量的试验测定。计算结果与试验结果表明，本文提出的基础滑移隔震体系的计算公式是可行的，为滑移隔震房屋设计证结构安全、合理提供了理论依据。     

无剪力分配新法解剪力墙框架混合结构

本文将无剪力分配原理(Moment distribution without shearing force)推广到多层多跨剪力墙框架结构(Frame shear wall structures),使转角弯矩分配同时剪力平衡,以加速收敛.此法便于模似手算操作进行微机计算(Calculation with pocket computer);也便于用集体分配原理(Moment integration distribution)迅速得出手算成果。     

HIGH PERFORMANCE SPARSE SOLVER FOR UNSYMMETRICAL LINEAR EQUATIONS WITH OUT-OF-CORE STRATEGIES AND ITS APPLICATION ON MESHLESS METHODS

A new direct method for solving unsymmetrical sparse linear systems(USLS) arising from meshless methods was introduced. Computation of certain meshless methods such as meshless local Petrov-Galerkin (MLPG) method need to solve large USLS. The proposed solution method for unsymmetrical case performs factorization processes symmetrically on the upper and lower triangular portion of matrix, which differs from previous work based on general unsymmetrical process, and attains higher performance. It is shown that the solution algorithm for USLS can be simply derived from the existing approaches for the symmetrical case. The new matrix factorization algorithm in our method can be implemented easily by modifying a standard JKI symmetrical matrix factorization code. Multi-blocked out-of-core strategies were also developed to expand the solution scale. The approach convincingly increases the speed of the solution process, which is demonstrated with the numerical tests.     

阻抗式和简单调压室甩荷时水位波动的显式计算方法

目前对调压室水位波动的研究主要有图解法和数值计算。尚未看到考虑阻力项后调压室水位波动随时间变化的显函数计算方法。本文从调压室的基本微分方程出发，采用泰勒级数展开式推导了阻抗式和简单式调压室甩荷时水位波动的显式计算方法。运用这种方法不仅可以计算调压室水位波动的过程，而且可以方便的计算调压室的最高涌浪和第二振幅。     

强流长距离加速过程中的束包络数值计算

直线感应加速器束包络的计算编码（ABC）能给出束流从注入器经加速器直到最后聚焦的包络图，可以方便地研究各种磁场位形及隙压参数对束流输运的影响。作为算例，本文给出了束流为3kA，在长达40余米直线感应加速器的多种磁场位形下的束包络。束包络的计算对调试加速器磁场位形保证束流的正常输运及抑制束流不稳定性可提供帮助。     

球形封闭容器内一个简单的煤粉燃烧爆炸模型

在分析了大量球形封闭容器内煤粉燃烧爆炸实验数据基础上，考虑了煤粉燃烧爆炸机理所涉及的湍流燃烧、相变、各种化学反应动力学过程等复杂因素，并且对球形封闭容器内由于煤粉混合不均匀造成的燃烧不充分给予了考虑，得到了球形封闭容器内煤粉燃烧爆炸特征的数值计算结果，计算的压力－时间曲线与实验结果符合较好。     

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

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

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

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

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

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

用同伦方法反演非饱和土中溶质迁移参数

非饱和土中溶质迁移参数反演问题可以归结为非线性算子方程的求解问题. 将同伦方法 引入该问题的求解，通过构造线性同伦将原问题转化为求解同伦函数最小值的无约束优化问 题. 同时在分析了同伦参数正则化效应的基础上，提出一种两段同伦参数修正方法. 即在求 解的初始阶段，根据拟Sigmoid函数调整同伦参数，以追踪同伦路径，保证计算稳定地进行； 在迭代的后期，采用与残差相关的同伦参数修正方法，以抵抗观测噪声对求解的影响. 数值 算例为求解带有平衡及非平衡吸附效应的一维非饱和土中溶质迁移模型参数反演问题，计算 结果表明了该方法的大范围收敛性及较强的抵抗观测噪声的能力.