边界单元法在断裂力学中的应用

本文从弹性力基本方程出发,叙述了用边界单元法求解断裂力学问题的基本原理和数值计算方法.在微机上用常量单元计算了承拉双边裂纹板和中心裂纹板的应力强度因子K_1和裂纹张开位移,论证了边界单元法用于断裂问题的有效性.     

动态裂纹迅速扩展的数值分析

本文运用有限元法对双悬臂梁(DCB)试件动态裂纹迅速扩展进行了分析,指出当采用线性位移场的单元时,线性释放节点约束力是合理的。实例计算表明:动态裂纹尖端应力强度因子计算值与实验值能较好地符合,同时对文中提出的两参数——动态裂纹尖端张开角及材料止裂时的裂纹尖端张开角——的计算结果表明:它们与动态裂纹尖端应力强度因子 K_(1D)及材料止裂韧性 K_(1A)等价。建议在动态裂纹迅速扩展及止裂分析中也可采用此参数。     

交叉裂纹各裂尖应力强度因子同时快速求解的Williams单元

为了研究交叉裂纹各裂尖应力强度因子之间的相互影响,建立了交叉裂纹各裂尖应力强度因子同时快速求解的Williams单元,以十字交叉裂纹为例,分别对正交十字裂纹或斜十字裂纹各裂尖应力强度因子相关参数进行研究,分析了薄板尺寸、裂纹夹角与奇异区尺寸等相关参数对十字交叉裂纹各裂尖应力强度因子的影响,其中K_Ⅰ、K_(Ⅱ)为裂尖Ⅰ、Ⅱ型应力强度因子。算例分析表明,当板宽与裂纹长度满足W/a≥13时,能够忽略板的尺寸对正交十字裂纹裂尖应力强度因子的影响,可视为无限大板;对于斜十字裂纹情况,当水平裂纹长度不变时,随着斜裂纹长度参数d的增大,K_(Ⅰ,A)、K_(Ⅱ,C)和K_(Ⅰ,D)逐渐递减,而K_(Ⅱ,A)、K_(Ⅰ,B)、K_(Ⅱ,B)和K_(Ⅰ,C)逐渐递增,对于K_(Ⅱ,D)则出现先减后增的趋势;当斜裂纹长度参数d=0时,随着裂纹夹角γ的变化,K_(Ⅰ,A)、K_(Ⅰ,B)大小相等,符号相同,K_(II,A)、K_(Ⅱ,B)则大小相等,符号相反,且各裂尖应力强度因子对裂尖奇异区尺寸不敏感。证明了Williams单元能够同时快速求解交叉裂纹各裂尖应力强度因子,且具有较高的计算精度。     

高阶园奇异单元法计算混合型应力强度因子K_Ⅰ K_Ⅱ

在裂纹尖端单元和常规单元的连接面上,用矩阵方法使结点位移协调一致.给出了离心力的广义载荷列阵.通过广义结点位移解就可求得K_Ⅰ K_(Ⅱ).计算实例表明,精度和计算速度是满意的.     

斜载荷作用下带直裂纹无限大板问题

本文对于受斜载荷作用下带直裂纹无限大板问题采用复变函数的幂级数解法求得应力、应变、位移的分析解和应力强度因子K_Ⅰ、K_Ⅱ的表达式;并指出H·Tada、P·C·Paris等人在(2)、(3)文献中在这个问题上的错误。     

热交换管内壁半椭圆裂纹应力强度因子的数值研究

使用ANSYS软件建立了热交换管内壁横截面上含半椭圆裂纹的有限元分析模型,用来模拟裂纹尖端应力奇异性.通过对裂纹前缘进行离散,得到不同扩展的裂纹前缘.在不同热机械耦合载荷下,计算并比较了不同裂纹参数下有限元模型的应力强度因子K_Ⅰ值.讨论了K_Ⅰ值分布规律以及影响其分布情况的因素.     

裂纹表面线性分布力作用下I型应力强度因子的解析解

对于I型裂纹,考虑到裂纹表面受线性分布约束应力作用,利用复变函数方法,从两个基本的解析函数出发,推导出应力强度因子(SIF)和裂纹张开位移(COD)的解析解.数值计算了3种形式的应力分布,即右梯形分布(情况I)、均匀分布(情况II)和左梯形分布(情况III).通过对3种形式的应力强度因子和裂纹张开位移进行对比,发现裂纹表面约束应力的分布形式和位置对应力强度因子和裂纹张开位移有很大影响;随着约束应力区远离裂纹中心,应力强度因子减小,而裂纹中心张开位移随之增大.     

平板表面裂纹应力强度因子和应力分布规律的确定

基于有限元法(FEM)和节点位移方法,研究了在中间裂纹与侧裂纹处裂纹尖端的应力强度因子和应力场.使用8节点四边形等参元和1/4节点退化单元,对网格密度和裂纹长度对计算精度的影响进行了研究.改善裂纹尖端周围子区域的节点密度,在保证结果精度的同时,可以节省计算时间.为了计算和分析裂纹尖端的应力强度因子、应力场与位移场,利用Matlab/Simulink编写了关于平板Ⅰ型和Ⅱ型应力强度因子的代码,并比较了有限元方法与精确方法的计算结果.结果表明,所提出的方法有助于提高计算精度与收敛速度,且算法合理,可以提高仿真精度并指导工程设计.     

Ⅰ型裂纹尖端圆弧对应力强度因子影响的数值研究

应力强度因子表征了裂纹尖端奇异应力场的强度,它是研究裂纹扩展规律和带裂纹构件强度的基础。采用有限元法,对受均布荷载作用存在边缘Ⅰ型裂纹的平面板进行了数值分析。研究了裂纹尖端圆弧对应力强度因子的影响,分别计算了具有不同裂尖圆弧r_0的Ⅰ型裂纹的应力强度因子。采用应力法计算不同半径处的表观应力强度因子,插值到裂尖圆弧而得。根据计算结果,绘制K_Ⅰ-r_0曲线,利用最小二乘法拟合至尖裂纹(r_0=0)即得理想尖裂纹的应力强度因子,与解析解相差仅0.7%。该曲线为带圆弧裂纹的应力强度因子测试和带圆弧裂纹构件强度计算提供了依据。     

高强度轴承钢滚动接触疲劳的Ⅱ型裂纹扩展驱动力

本文探讨了滚动接触产生的残余应力对Ⅱ型应力强度因子幅度(ΔK)的影响,此ΔK_Ⅱ是指承受重复滚动接触轴承圈的次表层裂纹扩展驱动力。对于高强度轴承钢,分别在短寿命(N=2—10~7),产生较小周向残余应力的情况下及在长寿命(N=10~8—10~(10))产生较大周向残余应力的情况下进行了计算。在计算时还把赫芝应力、裂纹相对深度、裂纹倾斜角度及裂纹表面间的摩擦力等作为变量。在下列前提下进行计算:二维的、重复的、无摩擦的纯滚动接触;符合线弹性断裂力学条件;存在的是平面的、小的次表面裂纹。最后将计算结果与实际观察的裂纹生长剥落情况对比。研究结果表明:除裂纹倾角极小的特殊情况外,一般说在较长寿命下所形成较大的残余压应力能明显地使次表层粗糙裂纹的扩展驱动力(ΔK_Ⅱ)减小。     

基于IMM-IME结构的藏文输入法研究

IMM-IME(输入法管理器-输入法编辑器)结构是Windows多语言支持的核心组件之一,为输入法的开发提供了强有力的支持.通过讨论中文Windows操作系统平台下输入法的基本工作原理,分析了中文Windows平台下基于IMM-IME结构的输入法的构成、接口、设计思想以及一些关键技术,并结合这些技术实现了基于IMM-IME结构的"央金藏文输入法".     

介绍一种改良的抗酸染色法

本文介绍并通过实验证明了一种建立的改良抗酸染色法与常规染色法具有同样的可靠性。     

求解指标1的微分代数方程组的一类新方法

给出求解指标1的微分代数方程组的一类新的计算方法.将微磁学仿真的方法推广到求解微分代数方程组,并给出方法的收敛性和相容性分析.利用与伴随法相复合的方法,提高方法的收敛阶.并将方法应用于晶体管放大器的模型中.数值实验表明方法是有效的.     

Sramek公式及修正式计算每搏量与导纳法的比较

本文应用阻抗—导纳式心输出量机对40名正常成年人(男、女各半)进行每搏量(SV)的测定,并把Sramek公式及Sramek修正式和导纳法3种不同方法测得的SV进行比较。结果表明,应用Sramek公式和Sramek修正式测得的SV均大于标准对照值,而导纳法的SV值在标准对照值范围内。提示导纳法测得的SV值可能更为可靠。     

三维快速多极虚边界元配点法

以三维弹性力学问题为研究背景,提出了一种三维快速多极虚边界元配点法的求解思想,即将三维快速多极展开的基本思想和广义极小残值法运用于求解传统虚边界元配点法方程.文中将三维弹性问题的基本解推导为适合于虚边界元快速多极算法的展开格式,经数值计算格式的演变,使求解方程的计算量和储存量与所求问题的计算自由度数成线性比例,以达到数值模拟大规模自由度问题的目的.算例说明了该方法的可行性、计算效率和计算精度.此外,该方法的思想具有一般性,应用上具有扩展性.     

应用FDM与AWE技术快速求解导体目标宽带RCS

渐近波形估计(AWE)技术是分析目标宽带有效数值方法之一,但该方法需要多次存储高阶频率导数阻抗矩阵,内存消耗大.本文将快速偶极子法与AWE技术相结合,只需存储近区场阻抗矩阵及其高阶频率导数阻抗矩阵,并且大大加速了在迭代求解过程中的矩阵矢量乘积运算.与传统AWE技术相比,计算时间和内存消耗都得到了有效缩减,数值结果证明了本方法的有效性和精确性.     

移动最小二乘配点的Petrov-Galerkin局部无网格方法

基于加权残值法和移动最小二乘(MLS)法并结合局部Petrov-Galerkin无网格方法(MLPG)的灵活性,将移动最小二乘配点法应用到无网格方法当中,建立了MLS配点无网格法的基本方程.在局部子域上利用Petrov-Galerkin原理给出了微分方程局部弱形式,通过惩罚因子引入本质边界条件;将局部弱对称形式进行离散化后,推导出移动最小二乘配点的Petrov-Galerkin局部无网格系统的刚度矩阵、载荷矩阵.通过数值算例证明该方法具有很高精确性、有效性和实用性.     

凹凸边界形状热弹性问题的自适应无网格法

为了解决具有凹凸边界形状的平面非定常拟静态热力耦合问题,采用无网格伽辽金算法(EFG)进行求解,使用移动最小二乘法构造形函数及拉格朗日乘子法处理本质边界条件(第一类边界条件),通过引入Voronoi邻接准则和后验误差式,对后续结果进行自适应优化;构建了一种新的适用于非定常拟静态热力耦合问题的EFG法自适应计算模型。为了验证计算模型的可行性,分别计算在二维混合边界条件下光滑与凹凸边界形状平面的温度场以及位移场的分布,并与有限单元法的计算结果进行对比,表征了有限单元法和无网格法计算结果的差异,验证了非定常拟静态热力耦合问题的EFG法计算模型的有效性和精确性。     

The notched filtering characteristics of stratified volume holographic grating

Utilizing the tool of beam propagation method (BPM) to calculate the zeroth order diffraction beam intensity, we find SVHG displays notched diffraction response as a function of the readout wavelength. Using the method of SA and considering the variance of refractive index as the readout wavelength changes, a practiced notch filter can be designed and the period of the filter is discussed. Foundation item: Supported by the National Natural Science Foundation of China (69607005) Biography: LI Xiang-yang (1974-), male, Master candidate: Research direction: message optics.     

微分求积方法及其在力学应用中的若干新进展

简单地介绍了微分求积法（DQM）的基本原理及发展、影响微分求积方法数值解精度的因素、与有限差分方法和有限元方法比较而言DQM的优点,同时重点介绍了近年来DQM在固体力学和流体力学相关领域应用中的若干新进展和展望.