共查询到18条相似文献,搜索用时 125 毫秒
1.
形状设计灵敏度分析的改进的再生核质点法 总被引:3,自引:0,他引:3
基于物质导数概念和直接微分法,将再生核质点法应用于形状设计灵敏度分析(DSA)中。导出了基于无网格近似的灵敏度方程,特别强调了在考虑形状函数关于设计变量的物质导数时无网格方法与有限元法的不同。通过对RKPM形状函数及其物质导数进行矩式显式表述,提高了无网格方法的计算效率。对两个二维线弹性问题进行了位移灵敏度和应力灵敏度分析,计算结果与解析解吻合的很好;同时通过对通常的RKPM和改进的RKPM计算耗时的比较,显示了该方法不仅有效,而且可以显著地提高计算效率。 相似文献
2.
断裂力学是工程材料和结构的疲劳与断裂分析、损伤容限设计和结构完整性评定的理论基础. 应力强度因子作为线弹性裂纹尖端奇异场的单一表征参量和裂纹扩展驱动力, 在裂纹体的断裂力学分析中发挥着关键作用. 权函数法为复杂受载裂纹体的应力强度因子求解计算提供了强有力的解析工具, 不但具有远高于各类数值解法的计算效率, 而且精度可靠, 使用方便. 本文结合笔者团队在权函数法方面的长期研究工作, 对该方法自20世纪70年代初提出至今半个世纪以来, 国际断裂界在二维和三维权函数理论与应用方面的主要研究进展作了回顾和评述, 并对其未来发展提出了展望. 主要内容涵盖: 当前国际断裂界广泛应用的3种二维裂纹解析权函数法简介和以格林函数为基准的验证评价; 三维裂纹问题的片条合成权函数法和点载荷权函数法; 权函数法在复杂受载裂纹体的应力强度因子和裂纹张开位移等关键力学参量计算、内聚力/桥连等裂纹模型分析、共线多裂纹权函数理论及其在剩余强度预测等方面的应用, 以及复杂裂纹几何的工程化权函数分析和权函数法的反向应用问题. 相似文献
3.
基于ALE方法的3D充填流动模拟 总被引:1,自引:1,他引:0
基于任意拉格朗日-欧拉方法发展了三维充填流动的数值模拟方案.该方案采用ALE方法准确地追踪移动自由面的位置并避免了网格扭曲;基于移动最小二乘曲面拟合方法提出了移动自由面上网格节点重定位方法,将充填流动的网格更新过程简化为自由面附近的局部网格重划分过程,并通过分级多面体三角剖分实现,减小了网格划分的计算量,实现了实时网格生成.给出的数值算例结果表明了该数值模型对三维充填流动模拟的有效性. 相似文献
4.
三维计算断裂力学 总被引:2,自引:2,他引:0
断裂力学理论从1921年Griffith研究玻璃的脆性断裂问题开始,经历了从线弹性体系到弹塑性以及蠕变理论体系、从单参数到多参数体系和从理想的二维平面理论到实际三维含裂纹结构的三维断裂理论的发展历程。针对应力强度因子K和J积分以及C(t)积分的计算方法从理想化模型的理论计算发展到实际复杂工程结构裂纹体计算的各种商业软件平台以及专业的断裂理论分析平台。尤其是随着计算机技术的发展,对三维含裂结构的静态和扩展裂纹的计算模拟已经能够融入计算机辅助设计。结合本研究组近30年来在三维疲劳断裂理论和应用研究方面的体会,简述了三维计算断裂力学从裂纹体应力应变分析和断裂参数计算到三维蠕变断裂和疲劳裂纹扩展模拟的国内外进展,并对涉及的计算方法,包括原子尺度和跨尺度的计算模拟,以及目前面临的挑战性问题作了简要介绍和分析。 相似文献
5.
6.
高超声速气动热环境的数值计算对算法和网格的敏感度极高. 随着高超声速飞行器外形日益复杂, 生成高质量的结构网格时间成本呈指数增加, 难以满足工程应用的需求. 非结构/混合网格因具有很强的复杂外形适应能力, 为了缩短任务周期, 有必要在非结构/混合网格上开展高精度的气动热环境数值计算方法研究. 梯度重构方法是影响非结构/混合网格热流计算精度的重要因素之一. 本文通过引入多维梯度重构方法, 发展了基于常规的非结构/混合网格的高精度热流计算方法, 对典型的高超声速Benchmark算例(二维圆柱)进行了模拟, 并与气动力计算广泛采用的Green-Gauss类方法和最小二乘类方法进行了对比. 计算结果表明, 多维梯度重构方法能有效提高非结构/混合网格热流预测精度, 其鲁棒性和收敛性更好. 最后将多维梯度重构方法应用于常规混合网格的三维圆柱和三维双椭球绕流问题, 得到了与实验值吻合较好的热流计算结果, 展现了良好的应用前景. 相似文献
7.
本文所提算法适用于二维和三维多介质流体力学两步欧拉数值方法中输运计算的混合网格(包括自由面网格)界面处理。在一个混合网格中,界面被近似地看作直线(二维)或平面(三维)。整个方法分为三步:(1)第一步,用混合网格周围的八个网格的介质面积份额(二维)或二十六个网格的介质体积份额(三维)确定界面的法线方向;第二步,用混合网格的本身的介质面积份额(二维)或体积份额(三维)确定界面的方程(位置);第三步,用此直线方程求出通过网格边界的流以及下一时刻网格的面积份额(二维)或体积份额(三维)。最后给出了用此方法所做的一些数值计算及与SLIC算法的比较。 相似文献
8.
根据三维裂纹前沿局部坐标系中主法平面内的周向应力计算公式,可以确定主法平面内的最大周向应力σθ,并由最大周向拉应力裂纹破裂准则确定裂隙岩体的破裂载荷和裂纹的破裂方向以及破裂步长。在三维无网格方法的节点的影响域中,选择合适的节点密度,以确保计算精度,并对单向轴压条件下内置三维裂纹扩展进行了数值模拟研究。根据数值计算结果,分析了破裂载荷、裂纹扩展角和扩展步长,揭示了内置三维裂纹的破裂过程,结合插值法清楚地展示了三维裂纹的破裂曲面形态。研究结果表明,针对三维裂纹破裂过程和破裂曲面形态无网格方法数值计算结果与FLAC-3D计算结果和试验结果吻合较好。 相似文献
9.
10.
11.
12.
13.
《Journal of Fluids and Structures》2007,23(6):839-857
In this paper, we present a detailed derivation of the numerical method, Immersed Finite Element Method (IFEM), for the solution of fluid-structure interaction problems. This method is developed based on the Immersed Boundary (IB) method that was initiated by Peskin, with additional capabilities in handling nonuniform and independent meshes and applying arbitrary boundary conditions on both fluid and solid domains. A higher order interpolation function is adopted from one of the mesh-free methods, the Reproducing Kernel Particle Method (RKPM), which relieves the uniformity constraint of fluid meshes. Two 2-D example problems are presented to illustrate the capabilities of the algorithm. The accuracy in the numerical analysis demonstrates that the IFEM algorithm is a reliable and robust numerical approach to solve fluid and deformable solid interactions. 相似文献
14.
15.
The work is concerned with the modeling and simulation of large scale ductile fracture in plate and shell structures. A meshfree method – the reproducing kernel particle method (RKPM) – is used in numerical computations in order to enact dynamic crack propagation without remeshing. There are several novelties in the present approach. First, we have developed a crack surface approximation and particle split algorithm for three-dimensional through-thickness cracks. Second, to represent evolving crack surface in 3D shell structures, a 3D parametric visibility condition algorithm is proposed, which re-constructs the local connectivity map for particles near the crack tip or crack surfaces, so that the meshfree interpolation field can represent physical material separation in the computational domain. Third, the constitutive update formulas in explicit time integration by different versions of Gurson models and the rate-dependent Johnson–Cook model are implemented for 3D computations. Finally, the performance of different Gurson-type models are investigated and compared with the experimental data of large scale in-plate tear process. Numerical simulations of crack propagation in stiffened plates and shells demonstrate that the proposed method provides an effective means to simulate ductile fracture in large scale plate/shell structures with engineering accuracy. 相似文献
16.
基于Delaunay三角化的无网格法计算结果后处理 总被引:3,自引:0,他引:3
最近新发展起来的无网格方法由于不需要显式网格,节省了网格生成所需的大量时间,并且避免了网格畸变问题,所以在处理一些特殊问题如移动边界、大变形、高梯度等方面显示出特殊的优越性。但另一方面也使得计算结果的全域后处理遇到困难。提出了一种基于Delaunay三角形背景网格的实用无网格计算结果后处理方法,以无网格离散节点为顶点生成Delaunay三角形,将无网格法计算得到的节点应力值映射插值得到三角形内的应力场云图颜色。给出的二维线弹性应力分析算例表明方法可靠实用。 相似文献
17.
由于直接配点法在求解边值问题时边界上的求解精度较低,本文提出了Hermite梯度重构核近似配点法(HGCM)来改进边界求解精度。重构核近似是无网格法中一种常用的近似函数,但是其在求解高阶导数时格式复杂且非常耗时。HGCM采用梯度重构核近似构建形函数的任意高阶导数,提高了计算效率;通过Hermite配点法构建离散方程,提高了边界求解精度。这种方法在求解对应变系数四阶偏微分方程的功能梯度材料板的静力问题时精度高,计算效率高,并可进一步推广应用于高阶偏微分方程描述的边值问题。 相似文献
18.
《International Journal of Solids and Structures》1999,36(24):3709-3729
This paper is concerned with an estimation of the effective elastic properties of an anisotropic body permeated by ellipsoidal cracks following the Eshelbys method. The classical integral expression of the P-tensor, the symmetrised derivative of the Greens tensor, is given for a 3D defined crack embedded in an anisotropic medium. The numerical evaluation of the P-tensor is validated with several limiting cases of simplified geometry cracks. The interest of a 3D representation of the cracks is shown with several applications: influence of the ellipticity and crack thickness aspect ratios, growing cracks. To conclude, a comparison between theoretical results and experimental data is done for the load-induced change of the compliance tensor for a damaged composite material. 相似文献