线性状态方程的精确计算

以解析法为基础,采用新的精确计算方法求解线性系统的状态方程．该方法属于一种２Ｎ类算法,它是把每一时间步长再进一步细分为２Ｎ个精细步长,然后在精细步长内利用泰勒展开式求解矩阵函数．与传统的Ｒｕｎｇｅ－Ｋｕｔａ方法比较,这种方法兼有解析法与数值法的优点,稳定性好,可以采用较大的时间步长以提高计算效率．     

非线性动力学方程的自适应精细积分

将定常结构动力方程的精细积分算法推广应用于非线性动力学问题的求解．对非线性项的线性化处理使该方法的计算精度对时间步长非常敏感,为此将龙贝格积分法引入该方法,提出了由此而产生的指数矩阵的快速精细算法,从而使时间步长的选择具有了自适应性,计算精度和效率均得到提高。     

空调动态负荷变步长计算方法

多层墙体传热反应系数／Ｚ传递系数的异步长关系，为发展变步长的反应系数法／Ｚ传递函数法开辟了一条可行的途径．本文进一步提出了多层墙体传热反应系数／Ｚ传递系数异步长关系的一种修正形式，将异步长关系由单向拓展为双向，并在此基础上，结合线性系统叠加原理，提出空调动态负荷计算的反应系数法和变步长Ｚ传递函数法，可以有效地降低计算量．实例计算表明效果良好     

用度量矩阵加速修改AHP中的判断矩阵

通过分析判断矩阵、一致性矩阵、导出矩阵及度量矩阵的关系，提出一种修改判断矩阵的预测加速修正法．当判断矩阵的一致性较差时，基于度量矩阵中偏离大的元素对判断矩阵一致性的影响较大，通过度量矩阵得出加速修正的步长．每次修改判断矩阵的一对元素即可进行判断矩阵的修正．实例分析表明，预测加速修正法是可行的，且可根据问题的性质，灵活确定修正的步长．     

龙格—库塔法和传输矩阵法求解非均匀长光纤光栅问题的比较

通过对几个具体的例子，分析研究龙格-库塔法和传输矩阵法求解非均匀光纤光栅的异同。发现在计算过程中矩阵法的运算速度快，精确度高，特别适合求解非均匀长光栅问题。而龙格-库塔法的计算结果则受步长的影响，当步长取得足够小时，两种方法计算结果基本吻合。     

直流输电系统的暂态仿真

应用结构矩阵法建立了双桥换流站的数学模型．为提高计算精度，提出以余弦外推公式预报换相电压，以抛物线拟合阀电流以判断阀的截止时刻，采取变步长与定步长相结合的方法来简化处理因阀通断造成的间断点．还讨论了交直流电网的接口问题．算例表明本方法是可行的．     

翅片管式蒸发器结霜问题的数值分析与实验研究

翅片管蒸发器结霜会增加传热热阻和流动阻力．为了合理确定化霜周期,必须对结霜过程作深入了解．文中应用动量、能量、质量输运方程模拟翅片管式蒸发器的结霜过程,并用数值求解方法预测了霜沉降量、蒸发器换热量随时间的变化过程．翅片管结霜是一个瞬变问题,为了计算方便,将其简化为准稳态问题,即在每个时间步长内,认为该过程是稳定的．然后把所得的霜层厚度以及霜层表面温度作为下一时间步长内传热传质的边界条件．在研究中考虑了霜层阻力引起风量下降这一因素,并通过测定空气除湿量来确定霜沉降量,以验证计算结果,两者吻合得很好     

结构非线性动力学方程的自适应精细积分算法

将定常结构动力方程的精细时程积分算法推广应用于非线性动力学问题时,对非线性项的线性化处理使该方法的计算精度对时间步长非常敏感.为此本研究中将龙贝格积分法引入该方法,提出了由此而产生的指数矩阵的快速精细算法,从而使时间步长的选择具有了自适应性,且计算精度和效率均得到提高.文中的数值算例给出了该方法的计算精度和效率.     

结构动力响应的无条件稳定计算格式

以三次Ｂ样函数作为速度时域函数，应用配点法，在任一时间步长中列出两个时刻的结构物运动方程式残值为零的条件，并使结果满足了稳定条件，得到一个计算结构物动力员应无条件稳定的计算格式。     

求解紊流流动方程组的预估-校正法

SIMPLE类方法是一种有效的求解Navier-Stokes方程组的原始变量法．然而，对非恒定流，尽管SIMPLE类方法采用隐式差分离散，允许较大时间步长，但对每一时间步都需要反复的迭代计算，求解相当费时，从分步法的概念出发，本文讨论了用预估-校正法求解隐式离散的Navier-Stokes方程组和k-ε紊流方程的过程，并与SIMPLE方法进行了比较．     

能力谱法在腋撑转换剪力墙结构抗震性能设计中的应用

按照刚度等效、内力等效的原则,对剪力墙结构进行模拟框架等效代换,从而能够通过分析杆件单元塑性铰出现的特点、次序来研究剪力墙结构的弹塑性性能。针对高层建筑转换结构使用中的要求,提出了一种腋撑式的新型梁式转换结构,并对该种结构进行了剪力墙结构算例的Push-over分析。采用能力谱法评估了该结构的整体抗震性能,表明该方法是对带转换层剪力墙结构进行弹塑性分析、评估结构抗震安全性的一种有效方法。     

多层厂房模型振动特性的理论分析与实验研究

运用有限元法与实验模态分析相结合的方法，以一个典型厂房的缩微模型为研究对象，研究多屡工业厂房的动态特性，给出了模型的低阶固有频率与振型。     

基于FMM的Krylov子空间IGMRES(m)新算法及其应用

研究了Krylov子空间GMRES(m)算法的基本理论,提出一种基于FMM的Krylov子空间截断型IGMRES(m)新算法.给出三物体弹性摩擦接触算例,计算结果表明,所提出算法在保证计算精度的前提下,可以大大减少迭代次数,显著提高计算效率.     

填充墙对隔震框架结构地震反应的影响分析

以一幢7层框架结构建筑为例,分别考虑常规抗震方案、以橡胶支座或铅芯橡胶支座为隔震元件的基底隔震方案。用非线性时程分析法探讨了填充墙对上述两种方案地震反应的影响。研究结果表明,对常规方案,考虑填充墙作用,诸层加速反应增加;而对于隔震方案,诸层加速度反应减小。     

输电网的中长期规划(Ⅰ)──支路追加法

输电网中长期规划的任务是：规划何时、何地、架设多少条何种电压等级的线路，经济、可靠地消除系统各阶段原有线路出现的潮流过载现象。它可分为单阶级段的静态决策和多阶段的动态决策两个子问题。本文利用追加连支路修改节点阻抗矩阵的思想，将静态决策问题用一个０—１整数规划模型表达，求解较为迅速。同时它具有解决含孤立电源接入系统的输电网规划的能力。计算实例表明了方法的功能和实用性。     

巷道锚杆各向异性作用效果的数值分析

将巷道围岩锚固区视为宏观上的各向异性连续介质,由锚杆支护参数确定其等效力学参数,利用快速拉格朗日元法 Ｆ Ｌ Ａ Ｃ 程序,分析了等效模型与实际模型的巷道变形特点。结果表明,等价效果主要依赖于锚杆的支护密度,密度越大,二者越接近;反之,误差增大。此外,锚杆长度对等价效果的影响很小。     

摩擦接触系统多极边界元法二次建模与数值优化

讨论了摩擦接触系统的建模与求解问题.将系统方程组的多极边界元法求解归结为变分问题并证明其等价性.在系统最优化数学模型的基础上,通过引入严格凸二次函数建立系统的二次规划模型,证明其最优解存在且唯一,并给出模型求解的最优性条件.     

剪力墙结构的静动力分析

以西北工业大学南村4#、5#小高层住宅为例,进行了剪力墙结构的壳元模型在不同受力状态的静力分析、模态分析、动力分析．分析剪力墙结构的基本振型、结构的侧移和层间侧移、结构的等效应力等,结果表明：该结构的最大层间侧移为二层、最大应力值在二层底部、底部的一到五层应力相对较大,一般在工程中需要加强．     

速度型减震技术在抗震水准提高的改造工程中的应用

福州琅岐岛某酒店公寓楼改造为中学宿舍,抗震设防类别由标准设防(丙类)提升为重点设防类(乙类),采用黏滞阻尼器对其进行抗震加固.对加设黏滞阻尼器模型和传统加固模型进行了罕遇地震下的抗震验算,对结构层间位移角、楼层位移、楼层剪力、地震输入能量分布和结构构件损伤等进行抗震性能分析,以及经济性对比.采用黏滞阻尼器进行抗震加固能有效将位移控制在2～3倍的弹性层间位移角限值,大幅提高结构抗震超强能力,替代主体结构构件耗散50%～60%地震能量,减轻结构构件的损伤.加设阻尼器加固后的结构抗震性能明显优于传统的抗震加固方式,且大幅降低结构在改造加固中的工程造价(降低约30%),可作为类似工程的参考.     

强震作用下钢-混凝土结构弹塑性损伤分析

通过LS-DYNA程序二次开发了钢材的弹塑性损伤本构模型,并分别建立了钢框架和混凝土核心筒的损伤准则.对强震作用下某20层的钢-混凝土混合结构的层间位移、框架与核心筒之间剪力分配等进行了数值分析,结果表明混凝土核心筒变形能力较钢框架差,强震作用下容易在薄弱层处产生变形集中破坏;对两种结构体系损伤发展过程分析表明,损伤指数能很好地跟踪两种结构体系的抗震能力退化过程,该模型和损伤准则可以用于强震作用下钢-混凝土混合结构的弹塑性损伤分析.