基于流固耦合作用的抽水引起地面沉降计算模型

过量开采地下承压水会引起地基土体的变形，进一步导致地面沉降和地裂缝等地质灾害，抽水引起地层的变形过程是土层中孔隙水和土颗粒相互作用的复杂过程．两步法和流固耦合法是计算地面沉降常用的两种思路，它们的适用性各有不同．本文基于理想的地层模型，分别建立两步法计算模型和流固耦合计算模型，利用FLAC3D6.0有限差分软件对模型进行求解，对比分析两类方法计算结果的差异．结果表明流固耦合法符合力学机制，能够更好地反映土体变形的发展过程，但计算时间更长．而两步法适用于稳定流，在渗流达到稳定流时两步法计算的孔压分布和地面沉降能够满足工程精度要求．本文的研究成果为合理选取地面沉降计算方法提供依据．     

复振型导数计算的移位多次模态加速法

本文研究了含粘阻尼结构的复振型导数计算问题，将导数计算问题成是一个简谐激谐的振的响应计算问题，采用多次模态加速法和移位法，导出了复振型导数计算的移位多次模态加速法，该方法具有确的数学和物理意义，可导出了已有的各种计算方法。算例表明本方法计算复振型导数只需用很少几个模态即可保证精度，计算量大大减少。     

复振型导数计算的移位多次模态加速法

本文研究了含粘性阻尼结构的复振型导数计算问题，将导数计算问题看成是一个简谐激振的响应计算问题，采用多次模态加速法和移位法，导出了复振型导数计算的移位多次模态加速法。该方法具有明确的数学和物理意义，可导出已有的各种计算方法。算例表明本方法计算复振型导数只需用很少几个模态即可保证精度，计算量大大减少。     

层状土中群桩沉降分析

本文提出一个半解析的方法分析单桩和群桩基础的沉降，地基土可以是任意层状的，并且可以是有限深的，单桩沉降采用荷载传递法计算。群桩中桩－桩相互作用采用近解析法计算。通过与其它数值分析方法结果的对比，表明本文提出的方法与更严密的计算方法有很好的一致性，能大大减少计算量，可分析桩数较多的群桩基础。     

特征向量导数计算各种模态法的比较和发展

在结构动力学修改以及结构动态设计中，灵敏度分析是十分重要的一环，往往也是主要工作量所在。如何在保证精度前提下减少灵敏度分析工作量和计算时间，具有重要意义。本文对特征向量灵敏度分析中计算特征向量导数的模态法作了综合评述，除了经典模态法和修正模态法之外，还对新近发展的迭代模态法和移位模态法进行了介绍。最后通过典型结构实例的计算仿真，对各种模态法作了比较。     

非等桩长或非等桩径的桩基沉降计算

本文提出了一种非等桩长和非等桩径的桩基沉降计算方法，计算实例表明，该法简便，实用，可用于桩基优化设计计算。     

基于虚拟激励法的结构随机振动声辐射分析

本文基于有限元法、边界元法和虚拟激励法，对随机激励下结构振动声辐射问题进行研究。提出了一种计算随机激励下结构振动声辐射问题的新方法，其中，有限元法用于计算结构谐振响应，边界元法用于计算结构振动声辐射，虚拟激励法结合有限元和边界元计算随机激励下结构振动声辐射问题。 数值算例表明，本文方法在计算精度上与传统方法等价，且更具高效性。     

热结构稳态响应的耦合灵敏度分析方法

研究结构稳态热变形和热应力的灵敏度分析方法，给出了直接法和伴随法两种算法。考虑了温度场的耦合作用，在直接法中需要计算温度场对设计变量的导数，在伴随法中需要计算热载荷对温度场的导数。对尺寸和形状两类设计变量的灵敏度分析算例，验证了本文方法的精度。伴随法在应用程序中的实现，为大型结构优化提供了高效率的灵敏度计算方法。     

三维板料成形过程的显式有限元分析

本文采用基于随动坐标系的假定应变域壳单元及显式有限元格式求解三维板料成形问题。板材料服从Ｈｉｌ各向异性弹塑性准则，板料与模具之间的接触界面由主仆面接触搜寻法处理，接触力由罚参数法计算。文中给出了几个三维成形过程的计算实例。数值算例表明，本文方法具有较高的计算精度和计算效率，可在微机上分析中等复杂程度的成形过程     

三维高频声波的矩阵压缩边界节点法模拟

本文采用边界节点法(Boundary Knot Method, BKM)求解三维高频声场．由于高频赫姆霍兹方程的解是振荡的，极大影响了数值求解的精确度，需要在计算区域增加离散点，这会增加计算量．同时对于大规模声学问题，依靠边界节点法形成的插值矩阵为满秩，导致计算量过高和存储量过大．所以本文采用矩阵压缩技术(Matrix Compression, MC)，在有效继承边界节点法高精确度的基础上减少计算内存需求和时间，从而提高计算效率．数值实验表明，MC-BKM 求解精度高、收敛速度快、计算时间少，在高频大规模声波问题中应用前景广泛．     

多重子结构法在非线性控制系统中的应用

多重子结构法是计算结构力学程序系统采用的最广泛的方法之一，由于计算结构力学与最优控制存在模拟关系，本文同样基于此法对非线性控制系统进行计算，并编制了相应的程序系统。本文工作进一步说明了计算结构力学与最优控制理论存在模拟关系。     

组合网架分析的康托洛维奇－加权残数法

本文给出的组合网架的一种实用计算方法是用康托洛维奇法把二维平板问题转化成一维问题，而后用加权残数法得到问题的解，计算结果表明，本文方法计算简便且有相当的精度，亦可用于一般网架结构的设计．     

组合网架分析的康托洛维奇－加权残数法

本文给出的组合网架的一种实用计算方法是用康托洛维奇法把二维平板问题转化成一维问题，而后用加权残数法得到问题的解，计算结果表明，本文方法计算简便且有相当的精度，亦可用于一般网架结构的设计．     

三维热权函数法和多虚拟裂纹扩展技术

热权函数法直接利用温度场与热权函数的乘积的积分来求应力强度因子(SIF)的过渡过程，它可以免除对每一时刻进行热弹性力学有限元或边界元应力分析，计算效率大大提高。本文给出了三维热权函数法的基本方程，并提出了求解三维热权函数法基本方程的多虚拟裂纹扩展法(MVCE法)。在MVCE法中，可以引入无穷多个虚拟裂纹扩展模式；虚拟裂纹扩展模式与应力强度因子的插值直接相联系，所得到的方程组的系数矩阵是一个三对角矩阵，具有良好的计算性能。它对于裂纹前缘SIF分布急剧变化的情况，有良好的数值模拟能力。实例计算表明，MVCE法具有权高的计算效率，并具有很高的计算精度。     

任意支承条件下连续板系结构的有限板块解法

本文基于薄板小挠度弯曲理论，构造出板元内部解析，边界挠度和边界法向弯矩以带补充项的付氏级数逼近，同时考虑域内多点支承作用的板元位移函数，给出了一处适用于任意支承条件下连续板系结构的有限板块法求解格式。数值计算结果表明：本文的方法具有良好的计算精度和计算效率，适于工程应用。     

腐蚀管道的可靠度分析方法研究

对腐蚀管道的模糊可靠度计算方法进行了深入研究，提出了确定性可靠度和模糊可靠度两个基本概念；提出了基于断裂失效判据、基于FAD失效评估图和基于剩余强度3种方法计算模糊可靠度的方法；提出了3种新的方法计算确定性可靠度，即改进的JC法、GA—JC法和改进的Monte-Carlo法．以胜利油田某试验区已加入缓蚀剂的注水管道为例，采用本文建立的计算模糊可靠度的3种计算方法和确定性可靠度的3种方法，分别计算注水管道可靠度随着时间的变化趋势．在模糊可靠度的3种计算方法中，基于剩余强度失效方法比较适中，是计算模糊可靠度较好的方法；在确定性可靠度的3种计算方法中，改进的Monte-Carlo法较好．     

关于复杂管系气柱固有频率计算的若干问题

本文论述了复杂管道系统气柱固有频率计算中存在的若干问题，在模型简化、计算方法等方面提出一个基于频响函数改进措施，并增加了验根过程。通过和实验数据、以往传递矩阵法计算结果的比较和实际管系的应用计算，证明本文方法是令人满意和有效的。     

有限元线法在解FPK方程中的应用

本文首次将有限元线法（ＦＥＭＯＬ）引入到非一性随机振动领域，给出了相应的单元特性的计算格式，探索了求ＦＰＫ方程的方法，并以船舶侧摆运动和Ｄｕｆｆｉｎｇ为例进行了必要的讨论。计算结果表明，用有限元线法求解非线性随机振动问题，较有限 法有所需内存小、计算精度好的效率高的优点。     

非线性随机过程的有限元分析

本文根据实际复杂工程结构独立随机变量多，计算工作量大的特点，提出了用广义随机变量替代一般的独立的变量以减少基本随机变量的数目，用响应曲面法替代敏度法或摄动法以提高计算效率，算例证明，以上措施取得了良好的效果，有助于非线性随机有限元法在实际工程可靠性分析中的应用。     

非线性随机过程的有限元分析

本文根据实际复杂工程结构独立随机变量多，计算工作量大的特点，提出了用广义随机变量替代一般的独立随机变量以减少基本随机变量的数目，用响应曲面法替代敏度法或摄动法以提高计算效率。算例证明，以上措施取得了良好的效果，有助于非线性随机有限元法在实际工程结构可靠性分析中的应用。