首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 109 毫秒
1.
形状设计灵敏度分析的改进的再生核质点法   总被引:3,自引:0,他引:3  
基于物质导数概念和直接微分法,将再生核质点法应用于形状设计灵敏度分析(DSA)中。导出了基于无网格近似的灵敏度方程,特别强调了在考虑形状函数关于设计变量的物质导数时无网格方法与有限元法的不同。通过对RKPM形状函数及其物质导数进行矩式显式表述,提高了无网格方法的计算效率。对两个二维线弹性问题进行了位移灵敏度和应力灵敏度分析,计算结果与解析解吻合的很好;同时通过对通常的RKPM和改进的RKPM计算耗时的比较,显示了该方法不仅有效,而且可以显著地提高计算效率。  相似文献   

2.
应用再造核质点法(RKPM)进行了结构裂纹计算问题的研究。将不连续处理技术一可视准则和衍射方法应用于RKPM来模拟裂纹附近场函数,实现二维和三维裂纹体的分析。避免了有限元方法中裂纹附近复杂的网格剖分工作。应用面向对象技术在计算程序中实现了二维和三维裂纹体的应力场计算分析,并且将结构离散工作和裂缝网格构造工作分离,从而提高计算效率。简单的算例表明应用RKPM方法在二维和三维含裂纹结构计算是有效的。  相似文献   

3.
基于计算力学中的结构优化思想,应用一种新型的显式几何更新算法,自行编制C++程序,实现地下管道形状设计的自动优化。管道内的流体假设为牛顿不可压缩流,并考虑惯性项。优化区域主要为管道竖直方向和水平方向的过渡段。形状优化的设计变量是几何边界的有限元节点坐标,优化目标是实现流体黏性能耗散的最小化。优化过程基于形状梯度,即通过形状敏感度分析来求解目标函数相对于设计变量的偏导数。所使用的显式几何更新算法既可以通过网格清晰描述形状,也可以大范围地自动更新网格。详细介绍了地下管道自动形状优化过程的关键步骤。通过数值算例探讨了不同注入速度、密度和黏度对其最优形状的影响。  相似文献   

4.
基于径向基函数的无单元法求解力学问题误差分析   总被引:1,自引:0,他引:1  
径向基函数形状参数的选择在无单元法数值计算中一直是一个热门的问题,现在已总结出许多确定形状参数的经验公式. 但还没有相关研究表明这些形状参数是如何随着影响域尺寸而变化的. 本文研究了MQ(multi-quadrics) 径向基函数中形状参数对无单元法计算误差的影响. 首先,从理论上分析了形函数导数随着形状参数值的变化趋势,和以计算点为中心节点对称布置与不对称布置的形函数导数的变化规律;然后分析了影响域尺寸对误差的影响,得到了在不同影响域尺寸下,误差随形状参数值变化的规律;在此基础上,给出了影响域范围值.  相似文献   

5.
配点类无网格法需要计算近似函数的二阶导数,因而在移动最小二乘(MLS)近似中至少要采用二次基函数。本文利用Voronoi图对双重点移动最小二乘近似法进行了改进,建立了基于Voronoi图的双重点移动最小二乘近似(VDG),并利用加权最小二乘法离散微分方程,导出了双重点最小二乘配点无网格法(MD GLS)。该方法将求解域用节点离散,并以节点为生成点建立Voronoi图,取Voronoi多边形的顶点为辅助点。近似函数及其二阶导数的计算过程可分解为两个步骤:首先用场函数节点值拟合辅助点处近似函数的一阶导数,再以辅助点处近似函数的一阶导数值拟合节点处近似函数的二阶导数。由于在每一步中只需计算MLS形函数及其一阶导数,这种近似方法需要较少的影响点和较小的影响域。同时借助于Voronoi结构的优良几何性质,可以快速地搜索影响点。研究表明,与基于MLS的加权最小二乘无网格法(MWLS)相比,这种方法可以显著提高计算效率,并且在精度和收敛性方面也有所改善。  相似文献   

6.
有理多项式技术有工程结构可靠度分析中的应用   总被引:3,自引:1,他引:2  
给出了在相关空间中改进JC方法计算结构可靠指标的迭代格式,并采用有理多项式技术计算功能函数的偏导数,以模拟实际工程中常见的功能函数不能明确的可靠度计算问题。因此,功能函数无论是线性或非线性,显式或隐式,该方法都简单且适用。数值结果表明本文方法具有较好的效率和精度。  相似文献   

7.
无网格RKPM法模拟平板轧制过程   总被引:3,自引:0,他引:3  
将无网格再生核质点法(RKPM)用于刚塑性可压缩材料轧制过程的模拟,采用罚函数满足本质边界条件,选用矩形影响域的张量积核函数,利用有限元网格作为积分单元,对求解域内和边界上采用不同的高斯积分方案。并将计算结果与刚塑性有限元计算结果和文献中的实验数据相比较,说明RKPM方法用于刚塑性可压缩材料轧制过程的可行性和正确性。  相似文献   

8.
黄娟  姚林泉 《力学季刊》2007,28(3):461-470
无网格法是求解微分方程定解问题的一种新数值方法.移动最小二乘近似只要求近似函数在各节点处的误差的平方和最小,对近似函数导数的误差没有任何约束.而广义移动最小二乘近似要求近似函数及其导数在所有节点处的误差的平方和最小.为了降低计算工作量,本文构造了要求近似函数在全部节点处和任意阶导数在部分节点处误差的平方和最小的改进广义移动最小二乘近似.数值计算显示本文提供的方法关于函数值和各阶导数值都具有很高的精度.  相似文献   

9.
采用核近似光滑粒子流体力学SPH(Smooth Particle Hydrodynamics)方法计算一元函数二阶导数和多元函数二阶偏导数,对避免核函数导数算法进行了改进,提出了对核函数光滑长度处处具有o(h)精度的算法,并分别推导出一元和多元函数的修正公式。用不同的粒子间距和不同的光滑长度进行了二阶导数的计算和热传导问题的模拟,进行了修正方法与原方法的误差分析。结果表明,本文提出的修正措施在提高精度、减少误差及加快收敛速度等方面起到很大的作用。  相似文献   

10.
本文提出了任意四边形薄板弯曲单元以解决非正交边界薄板弯曲的问题。文中给出了二维坐标变换的二阶导数的雅可比矩阵及其逆矩阵显式、形函数、荷载列阵及内力矩阵。通过四个计算实例说明精度是相当高的,并与三角形、矩形单元及函数解进行了比较。  相似文献   

11.
将重构核粒子法和势问题的边界积分方程方法结合,提出了势问题的重构核粒子边界无单元 法. 推导了势问题的重构核粒子边界无单元法的公式,研究其数值积分方案,建立了重构核 粒子边界无单元法的离散化边界积分方程,并推导了重构核粒子边界无单元法的内点位势的 积分公式. 重构核粒子法形成的形函数具有重构核函数的光滑性,且能再现多项式在插值点 的精确值,所以该方法具有更高的精度. 最后给出了数值算例,验证了所提方法的有效性 和正确性. }  相似文献   

12.
A novel ellipsoidal acoustic infinite element is proposed. It is based a new pressure representation, which can describe and solve the eUipsoidal acoustic field more exactly. The shape functions of this novel acoustic infinite element are similar .to the Burnett‘s method, while the weight functions are defined as the product of the complex conjugates of the shaped functions and an additional weighting factor. The code of this method is cheap to generate as for 1-D element because only 1-D integral needs to be numerical. Coupling with the standard finite element, this method provides a capability for very efficiently modeling acoustic fields surrounding structures of virtually any practical shape. This novel method was deduced in brief and the conclusion was kept in detail. To test the feasibility of this novel method efficiently, in the examples the infinite elements were considered, excluding the finite elements relative. This novel eUipsoidal acoustic infinite element can deduce the analytic solution of an oscillating sphere. The example of a prolate spheroid shows that the novel infinite element is superior to the boundary element and other acoustic infinite elements. Analytical and numerical results of these examples show that this novel method is feasible.  相似文献   

13.
This study proposes an elastodynamic time-reversal imaging method for the shape reconstruction of flaws with the traction-free boundary. The proposed method is a generalization of the synthetic aperture focusing technique (SAFT) implemented with the aid of time-reversal back-propagation computation. The reconstruction formula is derived from the boundary condition applied on the flaw boundary without using a diffraction theory yielding an inherently normalized objective function. No restriction is imposed on the wave source or the material properties as long as the wave medium is linearly elastic and lossless. There are two major advantages with the proposed method. First, it is free from ray tracing. The method can thus work with an arbitrary source and the wave modes including inhomogeneous waves. Second, signal conditioning is not required to focus a scattered field to the correct location in the reconstructed image. Numerical examples are presented to demonstrate the performance and benefits of the proposed method. In the numerical examples, the shape of cavities and cracks in an elastic layer is reconstructed from synthetic data. The results show that the method works well with the reverberating incident fields in which inhomogeneous and body waves coexist. Also demonstrated is that the curvature of the cavities is resolved well at the scale of the incident wave length. Through the numerical examples, the method is shown to be a versatile imaging method potentially useful for ultrasonic nondestructive characterization of small flaws.  相似文献   

14.
AN IMPROVED HYBRID BOUNDARY NODE METHOD IN TWO-DIMENSIONAL SOLIDS   总被引:3,自引:2,他引:1  
The hybrid boundary node method (HBNM) is a promising method for solving boundary value problems with the hybrid displacement variational formulation and shape functions from the moving least squares(MLS) approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the latter. Following its application in solving potential problems, it is further developed and numerically implemented for 2D solids in this paper. The rigid movement method is employed to solve the hyper-singular integrations. Numerical examples for some 2D solids have been given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method are studied through numerical examples.  相似文献   

15.
采用径向基函数配点法分析考虑剪切效应的梁板弯曲问题,该方法利用径向基函数作为近似函数,基于配点法离散方程,通过最小二乘法求解。径向基函数配点法在离散和计算过程中不需要任何形式的网格划分,是一种真正的无网格法;径向基函数可以用一元函数来描述多元函数,存在明显的储存和运算简单的特点;而基于配点法求解不需要积分,提高了计算效率。分析考虑剪切效应的薄梁板问题时,传统的有限元法或无网格法求解均会存在剪切锁闭问题,而径向基函数在全域内存在无限连续性,能够准确地满足Kirchhoff约束条件,因此径向基函数配点法能够消除剪切锁闭现象,而且不会出现应力波动。该方法的优势在于,其不仅易于离散、精度高,而且具有指数收敛率,计算效率高。数值算例验证了上述结论和该方法的稳定性。  相似文献   

16.
The reproducing kernel particle method(RKPM) has been efficiently applied to problems with large deformations, high gradients and high modal density. In this paper, it is extended to solve a nonlocal problem modeled by a fractional advection–diffusion equation(FADE), which exhibits a boundary layer with low regularity. We formulate this method on a moving least-square approach. Via the enrichment of fractional-order power functions to the traditional integer-order basis for RKPM,leading terms of the solution to the FADE can be exactly reproduced, which guarantees a good approximation to the boundary layer. Numerical tests are performed to verify the proposed approach.  相似文献   

17.
提出了一种光滑积分伪弱形式,将光滑积分拓展至被积函数非偏导项求解。结合光滑应变技术和伪弱形式,可实现有限元系统方程统一光滑积分求解,即对刚度矩阵和质量矩阵中的应变矩阵和形函数矩阵均可进行光滑积分处理,并转化为光滑子域的边界积分。光滑积分伪弱形式与光滑应变技术比较,增加了形函数矩阵不定积分处理过程,且没有降低有限元求解对形函数连续性的要求。不过,伪弱形式改变了单元积分的求解形式,连续质量矩阵求解也无需坐标映射和雅可比矩阵计算。以轴对称二维问题为研究对象,结果表明极度不规则三角形和四边形单元光滑积分伪弱形式在静态和动态有限元方程求解中也具有很好的精度。  相似文献   

18.
<正>A three-dimensional(3-D)approach based on the state space method is proposed to study size-dependent mechanical properties of ultra-thin plate-like elastic structures considering surface effects.The structure is modeled as a laminate composed of a bulk bounded with upper and bottom surface layers,which are allowed to have different material properties from the bulk layer.State equations,including the surface properties of the structure,can be established on the basis of 3-D fundamental elasticity to analyze the size-dependent static characteristics of the thin plate-like structure.Compared with two-dimensional plate theories based size-dependent models for thin film structures in literature,the present 3-D approach is exact,which can provide benchmark results to assess the accuracy of 2-D plate theories and various numerical approaches. To show the feasibility of the proposed approach,a 3-D analytical solution for a simply supported plate-like thin structure including surface layers is derived.An algorithm is proposed for the calculation of the state equations obtained to ensure that the numerical results can reveal the surface effects clearly even for extremely thin surface layers.Numerical examples are carried out to exhibit the surface effects and some discussions are provided based on the results obtained.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号