欧洲计算流体力学——个人观点

介绍了70年代初到20世纪末欧洲科学家对CFD发展的主要贡献,并重点关注了对国际流体力学数值方法会议(ICNMFD)和国际计算力学学术会议(ISCFD)的成功举办具有重大贡献的科学家.ICNMFD始于1969年(新西伯利亚),由参与了美国和苏联太空竞赛的两国科学家联合创办,而ISCFD则作为ICNMFD在亚太地区的同类会议初创于1995年(东京).由Oshima教授和本人提议,得益于上述两个会议全体科学委员会成员的共同努力,ICNMFD和ISCFD自2000年开始,合并成为每两年举办一次的国际计算流体力学会议(ICCFD).绝大部分本文讨论所涉及的学者通常都在受邀进行大会报告后,成为了前述三大国际会议中的一个或多个科学委员会的成员.作者对文中讨论所涉及的研究人员的圈定负有全权责任,上述人员的选择更多地基于应用而非理论计算流体力学,因此尚存在不完善之处.此外由于文中所讨论的学者大都是作者在40多年计算流体力学研究生涯中的私交好友,因此并没有将年轻一代的科学家考虑进来,但这并非由于他们的贡献不重要.另外值得注意的是欧洲的计算流体力学同美国的计算流体力学有非常紧密的联系,本文讨论所涉及的研究学者有些最早从欧洲开始涉足计算流体力学而现在则生活和工作于美国(例如Glowinski,Roe,Van Leer等人),而有些则是在美国取得了博士学位(例如Launder,Napolitano),他们都同美国的科学家进行过合作,因此很难清晰地对欧洲和美国的计算流体力学进行界定.对那些现已移居美国的欧洲科学家们,本文优先考虑了他们移居之前的贡献.最后但并非是最不重要的,作者通过电子邮件一一联系了文中所涉及的科学家,希望他们本人能够给出其自身贡献的评价草稿,从而帮助作者以最权威的视角在文中相应章节对其工作进行评述.本文向提供自述评价的诸位致以谢意,同时也希望其他同仁能够对本人不恰当的擅自评述给予谅解.     

复合材料结构非线性分析软件功用简介

简介了复合材料结构非线性有限元分析的功能，用途以单元处理技术。本软件可对复合材料结构进行应力应变分析，轴对称分析，空间体分析和空间板壳分析，其非线性分析功能包括；材料非线性分析，小位移小应变分析，大位移小应变分析和大位移大应变分析，同时，还具有较好的前处理和后处理功能。     

固体火箭发动机燃烧室喷管统一流场计算

采用ＳＩＭＰＬＥ算法的可压缩形式求解固体火箭发动机燃烧室喷管内的全速度统一流场;在非交错的配置网格基础上用有限体积法离散Ｎ－Ｓ方程,采用强隐式算法（ＳＩＰ）求解线性代数方程组;采用边界标志法实现复杂几何形状下的多区域统一计算;采用喷管外型逼近法保障计算稳定性;通过与同类计算结果的对比表明：上述算法的组合能够成功地获取固体火箭发动机内流场的整体结构。     

固体火箭发动机药柱裂纹腔内三维流场瞬态特性分析

以固体发动机药柱内存在的楔形裂纹为研究对象,采用三维流场控制方程,应用有限体积法计算了发动机点火启动阶段裂纹腔内的对流燃烧过程。在裂纹腔侧壁被点燃前,裂纹腔内的燃气压力基本呈均匀分布,且约等于燃烧室燃气压力;在裂纹腔侧壁被点燃后,燃气压力逐渐呈现出上部低、下部高的分布,且腔内平均压力远高于燃烧室内燃气压力;裂纹腔侧壁开口边缘处的推进剂首先达到点火温度开始燃烧,燃面迅速向内推进,燃气以非常高的速度向外流出裂纹腔。     

计算辐射冷却推力室不稳定温度场的隐式格式

导出了计算液体火箭发动机辐射冷却推力室不稳定温度场的隐式差分格式,并作了线性化处理。据此编制通用程序,可以计算任意辐射冷却推力室的不稳定温度场。     

SRM壳体/绝热层脱粘的XFEM

将扩展有限元法(XFEM)用于研究固体火箭发动机(SRM)壳体/绝热层的脱粘。将层隙型和紧贴型脱粘等效为双材料界面裂纹,建立了XFEM的双材料界面裂纹分析模型,并给出了其中核心算法——应变矩阵的求解,通过相互作用积分求出应力强度因子和能量释放率。算例结果表明,XFEM的计算结果与理论值的差异较小,可用于有效求解脱粘问题。     

花瓣铺层C/C喷管扩张段外缠绝热层工艺固化过程应力-应变分析

用有限元法对花瓣铺层C／C喷管扩张段外缠绝热层工艺固化过程的应力－应变进行了分析；并分析了C／C材料的层间缺陷对热应力的影响。通过分析可知，由于花瓣铺层C／C材料与绝热层材料线膨胀系数不一致，导致它们在固化过程后产生热应力；而C／C扩张段中少数缺陷的存在，使得缺陷附近的层间剪切应力大幅增加，从而可能导致其结构的完整性被破坏。文中的分析可为花瓣铺层C／C扩张段的设计和工艺提供参数。     

复合材料机翼翼梁与蒙皮间应力分析

本文采用三维有限元模式，用“ＡＢＡＱＵＳ”分析程序，进行了复合材料机翼翼梁和蒙皮结构的应力分析，对所获得的缘条和蒙皮之间六个应力分量的分布曲线，与采用二维有限元模式的结果做了对比分析，得出了一些有益的结论。     

Analytical solution methods for eigenbuckling of symmetric cross-ply composite laminates

Based on the first-order shear deformation theory, this paper explores the analytical methods for eigenbuckling of symmetric cross-ply rectangular composite laminates with a pair of parallel edges simply supported and the remaining two edges arbitrarily constrained. The main con-tribution of present paper lies in two aspects: one is to present a simple and effective analytical method, namely, the separation-of-variables method, which can generate the closed-form buckling solutions without any computational difficulty;the other is to incorporate the accurate computation method of exponential matrix into the state space technique to avoid the inevitable numerically ill-conditioned problems reported in several literatures. The results obtained via both analytical meth-ods are identical, and a good agreement with their counterparts in literature is observed. The separation-of-variables method can generate exact solutions within 1 s, which is impossible if the state space method is employed. Besides, the combination of the accurate computation method of exponential matrix and the state space method greatly improves the computational efficiency and gives correct results compared with the straightforward use of state space method.     

基于双量程推力传感器的精确测力技术

针对双量程推力传感器在使用中测试精度低的问题,结合传感器固有特性,开展了传感器动态特性分析及误差补偿方法的研究。通过有限元分析和力锤冲击试验,结合结构分析获得了双量程力传感器等效二自由度和三自由度动力学模型。利用近似求积数值法对冲击力试验输出信号进行补偿,获得了接近标准信号的测试结果,对该补偿方法进行了有效性验证。结果显示:本方法对发动机试验的推力信号进行了补偿,推动了单室双推力发动机推力的精确测量,为双推力发动机的推力测试提供了有效途径。