结构稳定分析域外奇点法

利用等效横向荷载的概念，基于静力问题基本解，导出了平面杆系结构线性失稳分 析的域外奇点法. 该法避开使用形式复杂的稳定问题控制微分方程的基本解，使之 成为静力方法的一种直观推广，具有物理意义明确，便于实现等优点. 给出一数值算例，验证了方法的正确性与实用性.     

概率约束评估的功能度量法的混沌控制

功能度量法是基于可靠度的结构优化设计中评估概率约束的一种方法,其改进均值(AMV)迭代格式具有简洁、高效的优点,但对一些非线性功能函数搜索最小功能目标点时可能陷入周期振荡或混沌解,本文利用混沌反馈控制的稳定转换法对功能度量法的AMV迭代格式实施收敛控制.首先展示一些功能函数应用功能度量法AMV格式迭代计算产生了周期解和混沌解现象,并对迭代算法进行了混沌动力学分析.然后利用稳定转换法对功能度量法迭代失败的参数区间进行混沌控制,使嵌入周期和混沌轨道的不稳定不动点稳定化,获得了稳定收敛解,实现了迭代解的周期振荡、分岔和混沌控制.     

固定接触界面法向静弹性刚度

基于Hertz接触理论推导了两个微凸体之间互相作用的法向接触静弹性刚度.根据修正后的一个微接触点的平截面积尺寸分布,给出了界面的总法向接触静弹性条件刚度、总条件载荷的解析解.将法向静弹性刚度的解析解嵌入到有限元软件中,获得整机的理论模态.通过实验对解析解进行了定量验证.以机床结合部为研究对象,在理论振型与实验振型一致的...     

???????????????????????????

主要基于三维弹性力学和状态空间法给出了四边简支各向同性矩形层合板自由振动和 强迫振动问题的精确解. 首先基于三维弹性力学建立了层合板的基本方程,利用状态空间法 解决了层合板的自由振动问题,然后根据Lagrange动力学方程求解了层合板受到横向冲击时 的强迫振动响应.     

基于最优组合的长基线网络RTK三频载波模糊度快速解算

针对三频载波模糊度解算TCAR法因电离层、观测噪声等误差影响而导致长基线模糊度解算可靠性较低的问题,提出了一种改进的适用于长基线网络RTK的TCAR方法.该方法继承了传统TCAR法逐级固定整周模糊度的思想,充分根据模糊度解算中每一步骤不同特点及不同目的选择最优组合观测值.首先,利用组合限制条件确定超宽巷模糊度解算最优组合;然后,引入平滑思想提高电离层残差估计精度,同时结合求解窄巷模糊度的WL/NL最优组合快速获得窄巷模糊度;最后,通过线性无关宽巷与窄巷组合模糊度实现基础载波模糊度准确解算.算例分析表明,基于最优组合的改进TCAR法有效提高了模糊度解算的成功率,实现了长基线网络RTK三频载波模糊度的快速解算.     

蒸汽压力作用下板条的湿热粘弹性分析

本文分别采用Fick第二定律和Fourier热传导定律模拟了板条的湿气场和温度场,并利用分离变量法给出了板条湿气场和温度场的解析解;采用Wong模型的Kitano和Galloway形式,模拟了回流焊接过程中产生的蒸汽压力;基于Weitsman湿热粘弹性本构关系,利用半逆解法和Laplace变换法,给出了在蒸汽压力作用下湿热粘弹性板条位移场和应力场的解析解.借助解析解,研究了温度和湿度对蒸汽压力作用下粘弹性板条位移场和应力场时空分布的影响.     

基于无网格Galerkin方法的整体式柔性机构的多准则拓扑优化设计

基于一种新型的数值计算技术—无网格伽辽金法,提出了一种整体式柔性机构拓扑优化设计的新方法.利用移动最小二乘形函数和伽辽金的弱变分形式建立弹性问题的控制方程,用Lagrange乘子法增强本质边界条件.在优化问题中,同时综合考虑机构的柔性和结构的刚度要求,用折衷规划法建立了柔性机构拓扑优化的多准则优化模型,这样对于非凸的优化问题也能保证搜索到Pareto解集所有的有效解.基于SIMP密度函数惩罚模型和优化准则法,建立了一种设计变量的显式迭代格式.运用经典算例证明了文中方法的正确性和有效性.     

基于Erdogan基本解边界元法计算应力强度因子

引入含裂纹问题基本解(Erdogan基本解)，提出了基于Erdogan基本解的样条虚边界 元法，并阐述了该法在实施过程中的特点与具体做法. 采用该方法详细分析了若干 典型裂纹问题，全面考察了方法的计算精度和收敛情况，以及在求解复杂裂纹问题方面 的能力. 结果显示，该方法具有精度高、收敛快、计算能力强等优点，是裂纹问题分析中 一种具有竞争力的通用计算方法.     

Matlab编程实现求解基于振型叠加法的多自由度动力学系统的时域和频域响应

给出了求解多自由度动力学系统响应的M atlab程序,这些程序基于振型叠加法可用于求解由质量矩阵M和刚度矩阵K以及常见阻尼矩阵描述的线性离散系统的时域和频域解.对于无阻尼系统,用户可以选择数值解或符号解析解(以时间或频率表示),并利用复模态叠加法计算了阻尼系统的数值解.总结了模态叠加方法下动力学响应的求解,并在简短的M atlab程序中实现.以三自由度系统和悬臂梁模型为例说明了程序的应用.这些程序也可用于工程应用中,通过对商用有限元软件包产生的质量和刚度进行后处理,产生感兴趣的时域和频域响应.     

四边简支各向同性矩形厚板的插值矩阵法解

针对强厚度矩形板四边简支情况,论文根据状态变量法思想,基于三维弹性理论基本方程,以3个位移分量及3个应力分量按双三角级数展开,将三维弹性力学控制方程转化为常微分方程边值问题.尽管一些各向异性弹性矩形厚板早已由状态空间法获得分析解,可是各向同性厚板的分析解至今难以获得,因为状态空间解法中特征方程有重根问题而不易于收敛.论文提出采用插值矩阵法直接对常微分方程进行求解,获得各向同性矩形厚板在四边简支边界条件下三维理论的位移和应力解,并与有限元精细结果进行比较,证明了本文解的准确性.     

Elliptic averaging methods using the sum of Jacobian elliptic delta and zeta functions as the generating solution

The present paper describes an improved version of the elliptic averaging method that provides a highly accurate periodic solution of a non-linear system based on the single-degree-of-freedom Duffing oscillator with a snap-through spring. In the proposed method, the sum of the Jacobian elliptic delta and zeta functions is used as the generating solution of the averaging method. The proposed method can be used to obtain the non-odd-order solution, which includes both even- and odd-order harmonic components. The stability analysis for the approximate solution obtained by the present method is also discussed. The stability of the solution is determined from the characteristic multiplier based on Floquet’s theorem. The proposed method is applied to a fundamental oscillator in a non-linear system. The numerical results demonstrate that the proposed method is very effective for analyzing the periodic solution of half-swing mode for systems based on Duffing oscillators with a snap-through spring.     

基于GPS/DR的车载导航设备的改进及算法

将 GPS 定位航位推算(DR)技术结合起来,可以用于车辆的定位和导航中。但是现有的航位推算设备都需要用角度传感器测量车辆的行驶方向,在实际应用中存在很多不便而且性价比不高。为解决这一问题,通过研究设计了只使用速度传感器就可以进行航位推算的设备,并对设备的工作原理、体系结构和基本算法进行了介绍。通过对模拟实验数据进行分析,证明了这种航位推算方法可以很好地弥补 GPS 导航定位的不足。     

A hybrid analytical–numerical solution for 3D notched plates

This study used a hybrid analytical and numerical method to analyze three-dimensional (3D) elastic bodies with sharp-V notches. The proposed method separates the 3D equilibrium equation into primary and shadow parts, where the solution of the primary part is the analytical solution under the generalized plane-strain theory, and the shadow part is solved numerically using a weak form based on the finite element theory. A least-squares method is then used to find the multiplication factors of these primary and shadow modes using 3D finite element results. Numerical simulations indicate that the proposed method can accurately simulate the singularities near a sharp V-notch. The major advantage of this method is that a 3D whole displacement field with the singular effect based on the theoretical solution near the notch can be obtained for anisotropic materials under arbitrary loads.     

Time-dependent motion of a two-dimensional floating elastic plate

Three methods are presented to determine the motion of a two-dimensional finite elastic plate floating on the water surface, which is released from rest and allowed to evolve freely. The first method is based on a generalized eigenfunction expansion and it is valid for all water depths. The second method is based on an integral equation derived from the Fourier transform, and it is valid for all water depths, although computations are made only for water of infinite depth. These two methods are both based on the frequency-domain solution—however no obvious connection exists between the two methods. The third method is valid only for shallow water, and it expresses the solution as the sum over decaying modes. We present a new derivation of the integral equation for a floating plate based on the Fourier transform of the equations of motion in the time domain. The solution obtained by each method is compared in the appropriate regime, and excellent agreement is found, thereby providing benchmark solutions. We also investigate the regime of validity of the infinite and shallow-depth solutions, and show that both give good results for a quite wide range of depths.     

Operational modal filter and its applications

In the paper, the procedure for the estimation of modal filter coefficients from output-only data is presented. The basic concept of the procedure consists in frequency response functions synthesis based on the knowledge of an operational modal model. A method of operational mode shapes scaling is described. The method is then compared with the classical modal filter and with modal filtration of responses spectra, which is sometimes used as a solution for modal filtration based on the output-only data. Each solution is applied to load identification and damage detection. The study shows the method verification on data obtained from laboratory experiment.     

结构电磁有限元分析的网格依赖性研究

分析了基于有限元方法的超材料电磁性能仿真分析结果对有限元网格的依赖性.通过对仿真结果的分析,发现基于S参数的材料性能参数反演结果对S参数的计算准确度的依赖性巨大,特别是处在谐振频域的S参数误差会导致等效电磁参数的性态误差.针对基于求解频率的网格自适应技术对远离求解频率的谐振频域的分析结果的准确度改善不够这个问题,需要确定合适的求解频率来执行网格自适应.提出了一种改进电磁性能仿真分析准确度的方法,其基本思想是首先给定一个初始的求解频率值,根据Kramers-Kronig关系反演等效电磁参数,进而获得谐振频率的近似值,并把此值作为求解频率来实现网格自适应以获得合适的网格剖分,最后经扫频获得S参数并进行等效电磁参数的反演.数值结果表明,该方法可有效提高仿真分析和参数反演的准确度.     

K-S函数集成局部性能约束的结构拓扑优化二阶逼近解法

应用K-S (Kreisselmeier–Steinhauser)函数,对结构拓扑优化问题中的局部性能如应力、疲劳寿命等进行集成然后求解。首先针对互逆规划的单目标多约束模型（称为s方模型）及多目标单约束模型（称为m方模型）,应用结构拓扑优化ICM方法,分别建立了基于K-S函数集成处理的优化模型,推导了集成化的约束（对s方模型）或目标（对m方模型）函数的一阶及二阶导数,采用序列二次规划模型对所建立的优化模型进行迭代求解,依据K-T条件给出了二次规划模型的迭代求解公式。然后基于K-S函数阐述了s方模型的集成迭代解法,亦即集成方法。最后,阐述了基于K-S函数的s方模型和m方模型交替融合的迭代解法,亦即集成-集成方法。结果表明集成-集成方法比单纯的集成方法收敛更快。     

结构杆件撤除后简化分析一般方法

本文提出结构杆件撤除后简化分析的一般方法,视杆件撤除后新结构的解为原结构的解与原结构上加以关联荷载解的叠加,利用原结构的解及总体刚度矩阵的逆矩陈求解新结构,无需进行结构重分析,使计算量大大减小,文中基于有限元法导出求关联荷载的一般算式。该法适用于空间网格结构和刚架结构一根杆件撤除问题,也适用于多根杆件同时撤除问题。     

求解非比例阻尼体系复模态的实模态摄动法

根据工程结构的实际情况,建立了非比例阻尼结构体系复模态特性的近似求解方 法------实模态摄动法. 这一方法以复Ritz向量展开原理为基础, 把非比例阻尼结构体系复模态特性的分析过程分解为两个基本步骤,首先以结构体系的实模 态向量构建复Ritz向量的求解子空间,然后通过非线性复代数方程组的求解代替扩阶后的复 特征值方程的求解,从而简化了计算过程.通过两个算例表明: 这一方法不仅计算简便，而且具有较高的计算精度和执行效率, 对于复杂的非比例阻尼系统是很适用的，具有一定的工程应用价值.     

边界元法中线性方程组的拟波阵分块求解

以边界单元法形成的线性方程组为例提出了一种分块求解具有非线称满系数矩阵方程组的拟波阵方法。这种在于高斯消去法基本原理的分块求解法,采用了拟波阵求解技术,使常驻计算机内存的仅为二行分块的系数子矩阵。