基于光滑有限元法的热-弹相场断裂研究

基于相场理论，论文构建了热力耦合作用下热-弹性材料的复杂断裂控制方程，采用光滑有限元法对该模型进行了数值离散；设计编写了相应的Matlab计算程序；开展了热-力载荷下的二维拉伸断裂模拟，将结果和有限元计算结果对比，吻合一致.同时采用该模型开展热冲击下陶瓷板的多裂纹扩展过程模拟，计算结果和实验结果吻合良好，验证了论文方法的有效性.     

气流式喷嘴雾化过程的分形特征及液滴分裂模型

研究了在雾化等过程中液滴粒度分形特征，提出了—个液滴分裂服从均匀分布概率的数学模型，通过计算机模拟，计算出同代液滴的结构分维和特征分维，并用对异代液滴计算得到的结果，验证粒度分布中的结构分维，结果吻合很好．实验测定了液滴分裂的分形维数，与模拟结果相符。     

气流式喷雾化过程的分形特征及液滴分裂模型

研究了在雾化等过程中液滴粒度分形特征，提出了一个液滴分裂服从均匀分布概率的数学模型，通过计算机模拟，计算出同代液滴的结构分维和特征分维，并用对异代液滴计算得到的结果，验证粒度分布中的结构分维，结果吻合很好，实验测定了液滴分裂的分形维数，与模拟结果相符。     

高层建筑测压风洞试验校准与不确定性的量化

高层建筑风洞试验可靠性常基于多家机构的CAARC标准高层建筑风洞试验结果进行对比检验，然而其风洞试验结果的不确定性难以量化。本研究对D类地貌下缩尺比为1∶400的CAARC标准高层建筑进行刚性模型测压风洞试验，首先研究了风洞试验参考位置对试验结果的影响，结果表明测压试验的参考位置宜选择受湍流影响较小的模型上方风洞顶部；其次，研究了采样次数对平均和极值风压系数的影响规律，为获得较稳定的极值风压系数，建议最小试验采样次数15次；最后，利用试验结果与国内外多家研究机构的结果进行对比分析，验证了本研究风洞试验的可靠性，量化了试验结果的不确定性，结果表明本研究高层建筑风洞试验的风压系数建模不确定性偏差为0.98,变异系数为0.14。     

二维翼型亚临界定常流动的有限元计算

本文使用了曲四边形Ｈｅｒｍｉｔｅ等参元，通过变分有限元法，计算了二维翼型亚临界定常位势流动中的压力分布，数值计算结果与试验结果符合较好。文中还使用了有限元法中的局部线性化理论，插值出元素中心点密度去处理迭代过程中元素系数阵的计算，其结果与通常的有限元法计算相一致，但计算时间却大大减少了。     

中柱失效工况下单层组合框架受力理论分析

基于已有试验，本文提出了轴向弹簧和转动弹簧约束的组合框架简化计算模型，推导了在集中荷载作用下刚性、半刚性连接平面组合梁双跨钢框架移除中柱后的荷载-跨中挠度曲线计算公式，考虑了压拱效应和连接刚度的影响，将公式计算结果与试验结果进行对比证明了推导过程和计算公式的可靠性．在对比结果的基础上，根据已有试验和模拟结果对已有悬索阶段的计算公式进行修正，修正后的公式与试验结果吻合较好．     

二维高亚音速Laval喷管流场的有限元计算

本文使用了八节点曲四边形等参元，通过变分有限元法，对二维高亚音速Ｌａｖａｌ喷管位势流场进行了计算，结果是满意的，并和一维结果进行了对比。文中还使用了有限元法特有的局部线化理论，取单元中心点密度去处理迭代过程中单元系数阵的计算，其结果与通常的有限元法计算相一致，但计算时间却大大减少了。     

气体—颗粒激波松弛流动特性研究

提出了一种解气体－颗粒粘性流动方程组的方法。在此方法中，计算颗粒阻力时采用Ｈｅｎｄｅｒｓｏｎ阻力系数公式，得到两相激波松弛流动数值结果，并与实验结果作了比较，详细分析和讨论了计算结果，揭示了两相激波松弛流动特性随颗粒容积比分的变化规律     

采用NND方法计算三维喷管气流场

本文运用NND显式差分格式，计算了三维喷管气流场。气流场计算的基本方程为一般贴体坐标系下三维守恒型的欧拉方程。采用了时间分裂法和Steger-Warming矢通量分裂技术。在喷管内沿周向的每个由轴线和壁面构成的子午面上根据泊松方程生成贴体网格。本文运用三维程序计算了轴对称JPL喷管，同时与实验结果和前人采用轴对称二维程序所计算的结果做了对比。最后，本文还计算了三维矢量喷管，计算结果与现有的实验结果一致。通过轴对称JPL喷管和三维矢量喷管的计算考核，表明建立的算法和编写的计算程序是正确的。文中提出了采用子午面形式的贴体网格时奇性轴的处理方法。计算结果表明在喷管壁面处，马赫数与压强的计算结果与实验值吻合较好，而在喷管轴线处，只有当网格较密时，才能得出与实验结果接近的计算结果。     

整装式液体随行装药燃烧的实验与数值模拟研究

报道了在整装式液体随行装药技术研究方面取得的最新实验结果，观测了随行药在不同点火延迟时间下燃烧对膛压曲线的影响，并建立了随行装药的内弹道经典模型，计算结果与实测结果有较好的一致性。     

入水后爆炸问题的数值模拟技术

应用各种数值计算方法及计算处理技术,编制程序实现了对弹体入水后爆炸问题全过程的数值模拟.其中弹体与水,爆轰产物气体与弹壳之间的相互作用通过流固耦合技术来描述;水面与空气、爆轰产物与水、空气之间的相互作用采用VOF方法(Volume-of-Fluid)来描述;采用了刚体-柔体转换、单元失效删除等计算技术以更高效、更好地模...     

带约束多体系统动力学方程的隐式算法

研究了带约束多体系统隐式算法,用子矩阵的形式推导出了多体系统正则方程的Ｊａｃｏｂｉ矩阵,它适用于多种隐式算法并给出了隐式Ｒｕｎｇｅ－Ｋｕｔｔａ算法,最后用一算例表明了隐式算法的计算效率和精度明显优于算法。     

具有奇异位置的多体系统动力学方程的隐式算法

本文研究了在运动过程中具有奇异位置的多体系统动力学方程的隐式算法，给出了隐式算法所用的Ｊａｃｏｂｉ矩阵，并建立了该矩阵中各子矩阵间的计算关系，提高了计算效率，计算结果表明隐式算法的计算速度和精度明显优于显式算法。     

空化、超空化流动的数值模拟方法研究

基于结构化网格，运用可压缩流N—S方程及k-ε湍流模型对流场进行求解，在低压区域引入一种基于混合密度函数的空化模型对轴对称体的空化、超空化流动进行了数值模拟．通过将半球圆柱的计算结果与实验数据和前人的计算结果进行对比，验证了所发展的数值方法的可靠性．最后，采用非定常的数值方法，研究了钝头体射弹的空化、超空化流动特性，并模拟了其超空泡的发展过程．     

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.     

Timoshenko梁谐响应求解新方法探索

在空间域上采用只与结点有关的无网格方法离散,在时间域上采用精细积分方法求 解. 无网格离散过程中,利用伽辽金积分等效弱形式代替微分形式的控制方程,并 用修正变分原理满足位移边界条件,采用移动最小二乘法求解离散的形函数,把形 函数代入等效积分弱形式得到离散的二阶方程;精细积分过程中非齐次项采 用Romberg积分. 同时给出了两种不同边界条件的谐响 应求解的两个数值算例,得到了精确的数值结果.     

利用多小波自适应格式求解流体力学方程

高阶计算格式的高精度、高分辨率对提高复杂流场的计算水平有重要的意义,为了提高AUSMPW格式对流场计算中激波等间断的分辨率,减小数值振荡,在原有AUSMPW格式的基础之上,利用多小波对函数进行多尺度分解,并采取阈值的方法生成自适应网格,提出了一种新的基于多小波自适应算法的AUSMPW格式,理论上可以达到任意阶精度. 将所得的压强、密度与原格式、TVD格式及WENO格式的计算结果进行了比较分析. 结果表明改进后的AUSMPW格式较原格式具有更高的分辨率、更强的捕捉间断的能力及更低的数值耗散.      

混凝土中爆炸模拟的数值方法比较

应用AUTODYN商业软件对装药在混凝土靶介质中的爆炸过程进行数值模拟,对比了混凝土中爆炸模拟的不同方法。将各数值方法模拟得到的混凝土毁伤特征尺寸与相关实验数据进行比较,研究了各数值方法的准确性及描述混凝土毁伤特征的可靠性参数。结果发现：Euler算法计算的压力最大,SPH和Lagrange耦合算法计算的压力最小;Euler算法和SPH算法计算的压力衰减过程会发生比较显著的扰动。以不同的平均网格尺寸对数值模拟结果的影响分析了各种数值方法的稳定性。比较靶体上同一测点的压力时间历程表明网格变化对测点处压力影响不大。通过应用各种数值方法模拟再现混凝土中的爆炸过程,揭示了各种方法的优缺点,明确了现有商业软件模拟混凝土中爆炸的适用性。     

The study of quasi wavelets based numerical method applied to Burgers' equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｏｎｃｅｐｔｏｆｗａｖｅｌｅｔｗａｓｆｉｒｓｔｉｎｔｒｏｄｕｃｅｄｂｙｇｅｏｐｈｙｓｉｃｉｓｔＭｏｒｌｅｔ[1].Ａｗａｖｅｌｅｔｆｕｎｃｔｉｏｎｔｈａｔｉｓｆｉｎｉｔｅ＿ｅｎｅｒｇｙｆｕｎｃｔｉｏｎｓｗｉｔｈｗｅｌｌｌｏｃａｌｉｚａｔｉｏｎｐｒｏｐｅｒｔｉｅｓｃａｎｂｅｄｅｓｉｇｎａｔｅｄｂｙｄｉｌａｔｉｏｎａｎｄｔｒａｎｓｌａｔｉｏｎ .Ｔｈｅａｂｉｌｉｔｙｏｆｗａｖｅｌｅｔｂａｓｅｄｍｅｔｈｏｄｔｏａｎａｌｙｓｅｆｕｎｃｔｉｏｎｓｉｎｔｅｒｍｏｆｔｈｅｉｒｌｏｃａ…     

贯流风机流场的实验和数值研究

贯流风机流场的分析研究是对其进行结构调整、性能优化的前提。本文运用烟迹和闪光照相的方法对贯流风机的内部流动进行了流场显示研究，得到了瞬态流场照片，并通过多点总压和静压测量得到了速度分布。本文还对贯流风机的流场进行了数值模拟。结果表明：数值计算与实验结果符合得较好，这两者的结合是研究贯流风机流动和改进其性能的有效手段。本文最后就贯流风机的噪声进行了初步研究，得出了一些有价值的结果。