复杂载荷三维裂纹分析双重边界元法

提出可用于高温、高转速状态下的热动力机械三维含裂构件热弹性分析方法——双重边界元法.首先建立了考虑温度及离心载荷的双重边界积分方程组,并对边界积分方程组的选取及适用范围进行了讨论。然后提出角非快调元模型离散技术。接着提出超奇异积分方程分析去除奇异性方法及数值积分技术.数值算例表明计算结果与有关权函数解十分吻合,说明了用双重边界元法计算复杂载荷条件下三维应力强度因子的有效性.还讨论了有关热应力强度因子权函数解的适用范围.     

扁球壳的边界元几何非线性分析

本文从壳体位移的三个微分方程出发，采用付立叶积分变换的基本解，利用加权残值法推导了几何非线性边界积分方程。这种基本解的壳体边界元法类似于板的非线性边界元法，各种变量物理意义明确，能方便地处理各种复杂边界条件及有开口情况。文末算例说明本文方法的可行性、收敛性和精确性，并与二变量边界单元法或有限元结果相比较，吻合较好。     

各向同性和横观各向同性材料的统一点力解

本文根据横观各向同性弹性力学的通解获得了无限体的点力解，由它可以直接退化到各向同性情形的Ｋｅｌｖｉｎ解，利用这个点力解编制的边界元法程序，适用于横观各向同性材料也适用于各向同性材料，因此是真正的统一点力解。还用边界元法计算了两个数值例题。     

平板网架结构分析的超级有限元法

本文给出一种二维矩形超级单元，用以分析平板网架结构，其做法是将大型平板网架结构离散成一系列矩形超级单元，考虑弯曲、剪切、挤压、拉压、翘曲等多种非经典变形效应，采用构件自由度向超级元自由度的转换把大型多构件问题的求解变为二维问题的求解，可大大减少未知量，又能保证精度，并可方便地与一般有限元软件连接。     

Experimental and FEM study of windshield subjected to high speed bird impact

Both experimental and finite element model (FEM) results are presented for the dynamic strength behavior of windshield subjected to bird impact. The experimental data taken from a series of high speed photographs are compared with the numerical results predicted by using FEM in which the windshield was modeled entirely with solid elements and the bird body was approximately simulated by an elastic-plastic material with failure element behavior. Effective plastic strain and element pressure were adopted as the failure criteria and once the pressure or the effective plastic strain of an element reached the critical value, the element would lose the tensile resistance capability completely. The deflection and stress distribution in the windshield were obtained. It is shown that the result from the finite element analysis agrees with those from the full-scale bird impact test. The project supported by the National Natural Science Foundation of China (10272011)     

用均匀化理论分析蜂窝结构的等效弹性参数

在线弹性范围内,根据均匀化理论,并结合有限元方法推导出适用于二维周期性结构的均匀化的有限元格式（Homo FEM）,计算出不同相对密度下的规则蜂窝结构的等效弹性模量Ee和泊松比νe.同时,利用蜂窝结构的代表胞元模型,用常规的有限元方法（FEM）模拟计算出相应的等效弹性参数.最后将两种数值计算结果与己有的理论公式进行了比较和分析讨论.结果表明:在考察的相对密度全场范围内（0～0.4）,HOmO FEM得到的蜂窝结构的 Ee和νe 与 FEM使用平面实体单元模拟计算得到的结果一致吻合,反映出 Homo FEM数值方法的客观准确性和可行性.而 Gibson公式和 W-K得到的等效弹性模量值 Ee只是在较小相对密度的情况下（小于0.15）与数值计算结果吻合.当结构的相对密度较大时,必须考虑胞棱附近区域由应力集中导致的复杂的应力和应变分布的影响.     

动态J积分及Rice J积分一致性比较与分析

首先回顾了动态J积分的发展进程,然后给出了两种常用动态J积分和Rice J积分的数学计算式,并对它们的一致性加以讨论,由此提出一种广义J积分的定义。然后通过算例,基于ANSYS软件,建立平面应变情形下纯I型中心裂纹模型,并利用编写的APDL程序计算了各个J积分,通过结果比较,进一步验证了在动态加载条件下,动态J积分比Rice J积分具有更好的路径守恒性,但仍可以看出它们的取值具有一定的一致性。最后,通过J积分间接法和裂尖位移外推法分别计算裂纹动态应力强度因子,所得结果吻合很好,验证了编写程序的正确性。     

Application of ‘Operational Quadrature Methods’ in Time Domain Boundary Element Methods

The usual time domain Boundary Element Method (BEM) contains fundamentalsolutions which are convoluted with time-dependent boundary data andintegrated over the boundary surface. Here, a new approach for theevaluation of the convolution integrals, the so-called OperationalQuadrature Methods developed by Lubich, is presented. In thisformulation, the convolution integral is numerically approximated by aquadrature formula whose weights are determined using the Laplacetransform of the fundamental solution and a linear multisep method. Tostudy the behaviour of the method, the numerical convolution of afundamental solution with a unit step function is compared with theanalytical result. Then, a time domain Boundary Element formulationapplying the Operational Quadrature Methods is derived. For thisformulation only the fundamental solutions in Laplace domain arenecessary. The properties of the new formulation are studied with anumerical example.     

透平机械离心压缩机和离心泵叶轮动力分析的三维旋转循环对称CN群算法

以离心压缩机叶轮为例,建立三维有限元模型,采用循环对称CN群算法计算了它的动力特性,详细讨论离心力作用下叶轮的应力状态,分析与轴连接的不同方式对叶轮应力的影响。应用循环对称算法能大大减少建模与计算的工作量,有很好的实用价值。计算结果表明,离心压缩机叶轮在旋转中,叶轮外缘应力比较高,轮盖上的应力较轮盘上的应力大,叶片根部容易产生应力集中,应采用适当的园角避免应力集中,过盈配合的叶轮内孔处应力比较大,需要确定合理的过盈量,既要保证叶轮不松动也要保证叶轮内孔处的强度条件。     

正交各向异性位势问题边界元法中几乎奇异积分的解析算法

几乎奇异积分的计算困难阻碍了边界元法的工程应用。本文针对二维正交各向异性位势问题边界元法中近边界点的几乎奇异积分,采用分部积分法,导出一种直接的解析计算公式。该解析公式可以精确计算线性单元上的几乎奇异积分。对二次单元,可将其细分为几个线性元,采用该解析公式近似计算其边界积分。当内点离当前积分单元较远时,仍保持常规高斯数值积分模式;而当内点离其较近时,因常规高斯积分结果失效,则采用该解析积分取代高斯数值积分。数值算例证明了该算法的有效性和精确性。二次元计算结果比线性元计算结果更精确。