粘弹性力学的对应原理及其数值反演方法

积分变换是处理粘弹性混合边值问题的重要数学工具,积分变换的应用使粘弹性混合边值问题在象空间与相应弹性混合边值问题对应起来,从而使粘弹性混合边值问题的求解可以继承和借鉴弹性问题的求解方法,再利用积分反演方法就可求得时间域粘弹性边值问题的解．本文结合国内外的研究成果,就粘弹性力学中存在的各种对应原理及数值反演方法进行了归类和总结．结合在求解粘弹性边值问题中的应用,对各类方法的特点进行了评述,并指出存在的问题及发展新的数值方法的研究重点．      

平面粘弹性问题的辛求解方法

将弹性力学辛对偶求解方法与Laplace变换相结合,提出了一个求解粘弹性平面问题的新方法。首先利用Laplace变换,将粘弹性平面问题转化为一个准弹性问题,在辛弹性力学的框架下,利用分离变量和辛本征展开法对其进行求解,然后由逆变换得到原问题的解。为证明方法的有效性,求解分析了矩形域平面粘弹性圣维南问题,得到了令人满意的结果。     

考虑损伤的粘弹性梁的纯弯曲

根据粘弹性损伤理论,分析了带损伤粘弹性矩形梁在受纯弯曲时损伤对应力的影响,得到了在Laplace变换域内损伤场和应力场的分布.利用Laplace数值逆变换,分别得到了损伤弹性梁和损伤粘弹性梁的最大应力和最大损伤值,分析了材料的粘性对梁内应力和损伤的影响.     

弹性力学的复变量无网格流形方法

为克服无网格流形方法配点过多、计算速度慢、容易形成病态方程组等缺点,将复变量移动最小二乘法与无网格流形方法相结合,提出了弹性力学的复变量无网格流形方法。分别采用线性基本与二次基进行计算,并与无网格流形方法相比。研究表明该方法计算量小、精度高。     

平面裂纹问题的h, p, hp型自适应无网格方法的研究

无网格方法以其独特的优点：不需“网格”（即节点间的连接信息）划分,特别适合自适应的分析,在分析中只需要高梯度域简单地插入离散点（ｈ型）或保持模型节点数、分布、覆盖大小均不变,中增加高误差覆盖上的函数的多项式阶次（ｐ型）,便可以得到更高精度的数值模型。针对平面弹性问题发展和推导一种显式后验误差指示公式,对平面裂纹实例进行了ｈ型,ｐ型,ｈｐ型三种不同类型的无网格自适应分析,数值分析结果表明了这种自适应     

无网格方法求解稳定渗流问题

使用无单元伽辽金(EFG)方法求解圆形油藏中心井和矩形油藏裂缝井两种稳定渗流模型。在中心井模型计算过程中，观察到采用对数等分布置节点是最有效的；计算结果与理论解和有限元解相比较，表明无网格方法是一种比有限元更为精确的方法。裂缝井模型通过在初始节点基础之上加密节点，获得了比较好的结果，并且给出了等压力曲线图。     

SHPB实验中粘弹性试件内部应力波的传播

建立描述SHPB实验中线性粘弹性试件内部应力波传播的控制方程组.根据试件两端与入射杆及透射杆接触的应力波特征关系给出耦合边界条件.对方程组和定解条件进行Laplace变换,求得试件内部应力在变换域像函数的表达式.采用数值反变换技术进行反Laplace变换,获得试件两端的应力时程曲线.对现有的固定Talbot反变换算法进行改进:将入射波像函数分解为基本部分和延迟部分,利用固定Talbot算法对基本部分入射波作用下的波动问题求解,其他部分的解通过延迟定理得到,最终解为两部分的叠加.采用这种改进算法得到的不同入射波下粘弹性试件的内部应力解与传统的基于特征线数值模拟方法的结果吻合.在此基础上探讨了粘弹性试件的几何参数和材料本构参数对透射波波形的影响.     

无网格局部Petrov-Galerkin方法在弹塑性断裂力学问题中的应用

采用无网格局部Petroy-Galerkin方法来分析弹塑性断裂力学问题.这种无网格方法采用移动最小二乘法(MLS)来构造近似试函数和采用Heaviside函数作为加权残值法中的权函数,由于近似函数不满足KroneckerDelta条件,因此采用直接插值法来施加本质边界条件.如果不考虑体力,所形成的整体刚度矩阵只包含局部边界积分,而不包含局部域积分和奇异积分.采用增量Newton-Raphson迭代法来求解弹塑性增量形式的局部Petrov-Galerkin方程.数值算例结果表明,该文方法对于弹塑性断裂力学问题的求解是可行的和有效的,并且所得到的结果具有较好的精度.     

Laplace变换法及其在粘弹性波中传播的应用

针对波传播分析理论的发展历程进行了简要的综述,详细介绍了几种处理粘弹性波传播问题的分析方法,重点讲解Laplace变换法以及Laplace变换在粘弹性波中的应用,对比分析几种方法在各自应用上的优劣,由于Laplace变换法能准确地描述应力波在任意时刻、任意点的波动情况,在处理大尺寸混凝土类构件中应力波传播问题时具有其独特的优势。     

固体力学中的无网格方法

简要地介绍了目前在无网格方法中主要使用的近似方案:移动最小二乘法、核函数法和单位分解法,并对这些方案的联系及各自的一致性条件进行了讨论;此外,对于无网格方法在数值计算中的离散方案、积分方案、边界条件的引入以及如何处理场函数或其导数的不连续性加以论述．      

Galerkin meshless methods based on partition of unity quadrature

Introduction Meshlessmethodsarenewmethodsofnumericalcomputationwhichhavebeendeveloped rapidlyinrecentyears.Inthesemethods,onlynodesareneeded,meshinformationistotally unnecessary.Thiscanavoidorpartlyavoidthedifficultyofmeshgeneration.Duetohigh accuracyandstability,Galerkinmeshlessmethodsareappliedbroadly,butitisunavoidable tocomputetheintegrationoverthewholephysicaldomaininGalerkinweakform,whichisa greatchallengeforGalerkinmeshlessmethodsbecauseoftheabsenceofmesh.TocarryouttheintegrationinGal…     

混合变换法在无网格伽辽金方法中的应用

移动最小二乘近似的非插值特征给无网格伽辽金方法的应用带来了一定的的困难，本文将再生核质点法中的混合变换法与无网格伽辽金方法相结合，通过对移动最小二乘近似进行局部修正，实现了无网格伽辽金方法中本质边界条件的精确施加。对权函数、影响半径、积分阶等对计算精度的影响进行了有益的探讨。数值计算结果表明了方法的可行性。     

Correct Solutions of Plane Elastic Problems for a Half-Plane

Continuous and integrable solutions and one-to-one relationships between boundary forces and displacements are found through the direct integration of the differential equations of the plane elastic problem for a half-plane with boundary conditions for either forces or displacements or with mixed boundary conditions. The necessary equilibrium conditions for forces and the compatibility conditions for displacements that ensure the correctness of the solutions are formulated     

自然单元法在Winkler地基薄板计算中的应用

自然单元法是一种基于Voronoi图及Delaunay三角形剖分图,以自然邻接点插值为试函数的一种无网格数值方法.本文以目前该方法中自然邻接点的Laplace插值形函数为基础,求出了其一阶及二阶导函数,建立了Winkler地基上薄板弯曲挠度的自然单元法求解控制方程,并编制了相应的计算程序,通过算例分析表明了本文方法的可行性和有效性.     

基于顺序雨流法和临界平面法的多轴变幅低周疲劳寿命评价

在多轴变幅疲劳寿命预测过程中,合适的雨流计数法对复杂加载历程分析非常重要,但是大多数雨流计数过程往往无法保持原始的加载顺序特性,进而会导致非保守的疲劳损伤和寿命预测.本文提出一种考虑加载顺序效应并基于临界面概念的多轴实时顺序雨流法,该方法既具有实时顺序计数特点,同时与Bannantine-Socie多轴雨流法结合,可以实现对主要通道内的载荷历程实时的雨流计数.基于Morrow模型,提出一种新的考虑加载顺序的线性损伤累积方法.相对于传统雨流计数法需要得到完整的载荷数据后才能进行分析的特点,新方法计算效率更高,实用性更强.通过对316L不锈钢的多轴疲劳试验数据的分析,验证了该方法在多轴疲劳寿命预测过程中的有效性.     

含裂纹体的单位分解扩展无网格方法

不连续体的数值模拟尤其是动态裂纹的追踪问题一直是工程界研究的热点和难点问题。无网格方法仅仅需要结点信息,非常适合于求解这类问题。基于单位分解思想,在移动最小二乘近似函数(MLS)中根据裂纹面的不连续位移增加一个Heaviside函数,在裂尖则增加四个扩展函数描述渐进裂纹位移场;应用Galerkin方法推导了平衡方程的离散线性方程,并给出了求解裂纹问题应力强度因子的计算公式。与其他类型的扩展无网格相比,在裂尖处近似函数不需要使用可视准则,很容易生成r1/2奇异;另一个优势是影响域并没有因为裂纹的存在而改变,不会降低方程的稀疏性,求解效率较高。数值算例表明,该方法能方便有效地模拟不连续问题,具有十分广阔的应用空间。     

Numerical-analytical method to determine the stress state of an elastic rectangular plate

A numerical-analytical method for solving the plane problem of elasticity is proposed. Systems of nonorthogonal functions are used. The method involves the minimization of a quadratic form that is equal to the integral of the sum of squared residuals of the solution and given forces. An explicit expression for stresses is derived. Bessel's inequality and the convergence of the solution are proved. The accuracy of the boundary conditions is estimated. The stress and strain distribution in the plate depending on the maximum magnitude of distributed forces and the size of their localization area is analyzed numerically. New quantitative and qualitative features of the stress distribution in the plate are established __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 1, pp. 90–98, January 2008.     

加权残数法通用程序包MWRAP的开发技术与1程应用

要想使加权残数法广泛应用于工程实际问题,需要解决两方面的问题,一是精度问题,二是必须研制和开发通用软件。本文阐述了加权残数法通用程序包MWRAP的结构设计和功能设置及其前后处理技术,详细阐述了其前后处理程序的程序设计流程;在该软件中计算模块,所用的计算方法为加权残数和有限元耦合法,为此推导了加权残数和有限元耦合法的公式,从计算结果来看该方法具有较高的精度,可以满足工程需要。最后,运用MWRAP计算了钢管混凝土柱肩梁结构的位移和应力,与实验结果相比吻合较好,说明了MWRAP具有较高的精度,加上其完善的前后处理功能,可以应用于工程实际问题。     

数值流形方法在流固耦合谐振分析中的应用

数值流形方法(流形法)是石根华博士利用现代数学中流形分析的有限覆盖技术建立起来的新的数值分析方法,统一解决了连续和非连续变形的力学问题,具有广阔的应用前景。本文将流形法应用于交界面耦合的流固振动分析,采用平面矩形数学网格,针对无粘、无旋、不可压缩流体和无阻尼的固体结构,提出分析流固耦合系统简谐振动的高阶流形法公式,其中,采用拉格朗日乘子法引入流场的已知边界条件。本文还初步研究了在特殊的无限远流场中采用解析解覆盖函数的实现技术。文中算例体现了流形法网格划分的方便性和计算的高精度,显示出流形法在数值解和解析解联合运用上的优势。     

涡量修正方法在切向炉内流场数值模拟中的应用

涡量是流体运动的一个基本物理量，旋涡是流体运动中常见的一种基本形态，又是湍流的一种基本结构，它在理论上和工程实践中占有相当重要的位置。本文在吴镇远（Ｊ．Ｃ．Ｗｕ）方法的基础上，提出了一种新的数值计算方法——涡量修正方法。对炉内速度场及涡量场进行了计算，其速度场的计算结果与实验基本符合，在网格划分基本相同的条件下与ＳＩＭＰＬＥ算法所得结果相比，涡量修正算法更接近于实验结果，为流场的数值模拟提供了一条新的途径。