基于Novozhilov理论带隔板连续弯箱位移参数的动态Bayes估计

推导了基于Novozhilov理论的简支弯箱有限曲条控制方程,并结合柔度法解决了带隔板连续弯箱的求解问题。首次建立了带隔板连续弯箱位移参数的动态Bayes误差函数,推导了相应的动态Bayes均值和方差表达式,提出步长的一维自动寻优方案后,结合共轭梯度法研究了带隔板连续弯箱位移参数的动态Bayes估计方法,同时给出了带隔板连续弯箱位移参数具体估计步骤。通过典型算例分析,总结了带隔板连续弯箱位移参数先验信息准确性判定方法及位移参数动态Bayes随机估计的其它重要结论。     

金属丝网——油液介质耦合减振器抗冲击动态力学特性研究

基于金属丝网－油液介质耦合减振器原理样机的抗冲击动态力学特性进行研究，样机的多参数匹配冲击试验研究表明该减振器具有强非线性动态特征，并能实现抗大冲击；在此基础上结合流体理论和库伦摩擦理论建立了原理样机非线性冲击动态特性模型，研究减振器抗冲击力学特性的设计可控性，该建模方法直接引入结构参数，为直接设计具体的器件建立了理论基础。     

基于共轭梯度法带隔板连续薄壁直箱位移参数的动态Bayes估计

推导了简支薄壁直箱有限条元控制方程,并结合柔度理论解决了带隔板连续薄壁直箱的超静定问题.基于Bayes定理,建立了带隔板连续薄壁直箱位移参数的动态Bayes误差函数,推导了相应的位移参数计算公式,并应用共轭梯度法研究了带隔板连续薄壁直箱位移参数的动态Bayes估计方法,同时给出了连续薄壁直箱位移参数的具体估计步骤,编制了相应的计算程序.通过典型算例分析可知,联合有限条元法和共轭梯度法进行带隔板连续薄壁直箱位移参数的动态Bayes估计是可行的,只需利用若干测点位移印能同时进行多参数估计,且先验信息合理时,迭代过程稳定地收敛于位移参数实际值.     

基于模态参数考虑边界条件变异的桥梁结构损伤识别

根据桥梁结构的实际工程特性,分析其边界条件变异、结构损伤及其参数变化,采用约束优化理论,建立以实测和理论模态参数误差平方和最小为目标函数的优化反演问题。基于矩阵摄动理论引入与结构动力方程对应的特征值和特征向量的一阶、二阶摄动量,将优化反演问题简化为非线性最小二乘法优化反演问题。针对桥梁结构边界条件对模态参数影响显著的实际情况,实施桥梁结构边界条件预识别,采用单元模态应变能方法预定位损伤,提出考虑边界条件变异的桥梁结构损伤识别具体流程。以一磁浮轨道梁方案为例,采用数值模拟进行边界条件变异及损伤的识别验证,结果表明:该方法能够有效识别边界条件的变异及构件损伤,识别参数的相对误差最大为12.48%,具有较高的识别精度。     

基于mbS模式的有限单元法

mbS模式及其有限元法是在固体和结构分析模型中引入薄膜、弯曲和剪切理论，且采用纯拉压、纯弯和纯剪单元进行分析的数值方法。在时空系中剖分物质单元和时间单元上构造以指数函数和贝塞尔函数为插入函数且按Lagrange插值条件的薄膜、弯曲和剪切等基本位移函数，由此得到更加完备和耦合的固体和结构实体单元的变形模式，根据能量泛函变分原理得到静动力有限元基本方程的一致格式。研究表明，mbS模式及其有限元法可用于梁柱和板壳等结构的静动力分析及屈曲分析。     

基于渐近传递函数解的截锥壳单元

普通截锥壳单元是分析旋转壳结构的常用单元,但应力计算的精度较差;而渐近传递函数解在圆锥壳的应力分析方面具有很高的计算精度。本文针对一般截锥壳单元应力计算精度不高的缺点,将传递函数法与有限元法进行结合,以圆锥壳的渐近传递函数解为插值函数,直接构造了一种高精度的截锥壳单元,该单元位移插值模式满足相容性和完备性要求,并具有力学概念清楚、计算精度高等特点。数值算例表明,采用该单元进行圆锥壳的内力和自由振动     

用于超弹性分析的高效杂交应变--EAS固体壳单元的研究

提出了一个用于橡胶材料分析的弱变分方程，基于一个可压缩Neo—Hookean模型，将EAS模式和杂交元法有机地结合起来，推导了一个高效、稳定的杂交应变——EAS固体壳单元，通过巧妙地选择应交插值函数和本文中提出的一个精化措施，克服了橡胶材料所表现出的大应变超弹性本构为杂交元正交化法的实施所带来的困难，保证了整个单元列式都仅采用低阶高斯积分，显著提高了计算效率，确保橡胶材料不可压缩性计算的顺利进行，并克服了固体壳元的厚度自锁问题。     

大转动曲边柱壳非线性3-D动力学建模及分析

对曲边柱壳受轴向非均匀内压作用下的大转动几何非线性3-D动力学行为进行了研究.基于Nayfeh and Pai[1]非线性壳体理论,给出了考虑几何非线性的3-D混合型（含内力与位移）动力学模型.为了克服该强非线性模型难以求解的问题,依据分析获得的结构静动态变形关系,采用Lagrange方程推导建立了基于结构静态解的曲边柱壳多自由度3-D动力学方程,并对其进行了线性化与降阶处理,结合差分法获得了一套高效的求解算法.与LS-DYNA有限元结果的吻合,验证了本文方法的正确性.最后分析了单元数和计算时间步分别对有限元模型和本文方法的影响,发现求解精度随着计算时间步的减小不断提高直至趋于稳定.同时对采用本文方法获得的曲边柱壳动态变形模式的分析表明:结构动态响应与其所受内压载荷沿轴向的分布形式关系紧密,可以通过改变或者设计内压轴向分布形式来影响以及控制结构的动态变形模式,从而应用于曲边柱壳结构设计及优化的工程实际中.     

复变量移动最小二乘法及其应用

提出了复变量移动最小二乘法，并详细讨论了基于正交基函数的复变量移动最小二乘 法. 然后，将复变量移动最小二乘法和弹性力学的边界无单元法结合，提出了弹性力学的复 变量边界无单元法，推导了相应的公式，并给出了数值算例. 基于正交基函数的复变量移动 最小二乘法的优点是不形成病态方程组、精度高，所形成的无网格方法计算量小. 复变量边 界无单元法是边界积分方程的无网格方法的直接列式法，容易引入边界条件，且具有更高的 精度.     

固体破坏理论的若干问题

本文列举和初步探讨了位于当代固体破坏理论的前沿上的十五个问题。它们包括：固体破坏理论的封闭、宏微观非线性力学、尺度效应、多层次计算、细微观实验力学技术、断裂过程区的描述、异质体的强度与韧性、材料的多层次抗破坏设计、高速破坏过程、非均匀介质的动态断裂与前兆、疲劳的起源、破坏元技术、内损伤破坏与微结构形貌演化失稳的识别、材料的微结构演化、电致失效力学。     

大学物理《力学》中的自由度分析

对《力学》中的物体自由度进行多方面分析,以深化教学、提高学生正 确分析物理问题的能力.使用实际教学分析的研究方法,在《力学》范围内讨论自由度与坐标、 自由与约束的关系并得以下结论: (1) 同一物体的自由度随其所在的``空间'不同而不同, 不因坐标系的选取不同而 异, 在同类参考系中不因参考系的动静而有别;(2)自由度遵循叠加原理. 讨论了质点系的总自由度及相关计算问题,并指出研究《力学》中自由度的意义.     

Development and design of new methods of measurement of state of strain of structures

The present paper deals with development and design of new methods utilizing Wiedemann's effect for determination of state of strain in building structures. Wiedemann's effect and some features of torsional strain of magnetic field are the basis of new experimental method. Especially the point electromagnetic strain gages using the effect of pure torsion of electromagnetic field to enable universal examination. For strain-gage measurements, almost all physical quantities are used which can be related to the variation in length of the structures. From the electric strain measurements, the most commonly used methods are the measurements by resonance-wire strain gages or by electric-resistance strain gages. In this paper, electromagnetic strain gages are discussed using the Wiedemann effect, and the author describes some new measuring equipment and his own suggestions and methods based on an application of this effect.     

Method of constructing estimates of the solution of equations of filtration of an aerated liquid

The exact solutions of the nonlinear equations of filtration of an aerated liquid have been obtained in [1–3]. In [4] the system of equations of an aerated liquid have been reduced to the heat-conduction equation under certain assumptions. An approximate method of computing the nonsteady flow of an aerated liquid is given in [5], where the real flow pattern is replaced by a computational scheme of successive change of stationary states. In [6] the same problem is solved by the method of averaging. In the present article estimates of the solution of the equations for nonstationary filtration of an aerated liquid in one-dimensional layer are constructed under certain conditions imposed on the desired functions. These estimates can be used as approximate solutions with known error or for the verification of the accuracy of different approximate methods. We note that the use of comparison theorem for the estimate of solutions of equations of nonlinear filtration is discussed in [7–9]. The methods of constructing estimates of solutions of various problems of heat conduction are given in [10, 11]     

Calculation of aerodynamic characteristics of arrays of profiles of arbitrary shape

It is well known that the problem on nonseparating potential flow of an incompressible fluid about an array of profiles reduces to an integral equation for a certain real function, determined on the contours of the profiles of the array. As such a function one can take, as was done, for instance, in [1–5], the relative velocity of the fluid on the profiles of the array. For arrays of profiles of arbitrary shape it is necessary to solve the corresponding integral equation numerically. In the particular examples of the calculation of aerodynamic arrays that are available [1–3] the numerical methods used were based on the approximate evaluation of contour integrals by rectangle formulas. As investigations showed, sizeable errors arose thereby in the approximate solution obtained, these being especially significant in the case of curved profiles of relatively small bulk. In the present paper a method for the numerical solution of the integral equation obtained in [5] is proposed. The method is based on the replacement of a profile of the array with an inscribed N polygon, the length of whose sides is of the order N^–1 and whose internal angles are close to . Convergence with increasing N of the numerical solution to an exact solution of the integral equations at the reference points is demonstrated. Examples of the calculation are given.Novosibirsk. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 105–112, March–April, 1972.