三维弹塑性自然单元法算法实现

自然单元法是一种新兴的无网格数值计算方法,其实质是基于自然相邻插值(C∞)的伽辽金法。该方法计算精度与四边形或六面体单元有限元法相当,自然相邻插值函数比其他无网格法插值函数的计算速度快。由于自然相邻插值在凸域的边界上的相邻点之间是严格线性的,所以自然单元法在边界面的处理也相当简单。本文研究了在自然单元法中采用Von.Mises,Mohr-Coulomb和Drucker-Prager屈服准则解决三维弹塑性问题,并编制了相应计算程序,最后通过算例验证算法的正确性。     

无单元法求解正交各向异性板自由振动问题

无单元法是一种新出现的数值方法。本文对无单元法的数学基础—滑动最小二乘法进行了详细的研究,推导了无单元法的形函数,并对一些关键问题,如权函数的选取,正交基函数,边界条件,数值实现方法等得出了研究结论。用无单元法研究了正交各向异性板的自由振动问题,由动力学变分原理和滑动最小二乘法导出了正交各向异性板的无单元法质量矩阵和刚度矩阵,编制了相应的计算程序,通过计算实例验证了该方法的有效性。     

线性强化材料弹塑性分析的自然单元法

自然单元法(NEM)是一种求解偏微分方程的无网格数值方法,其形函数兼具无网格法的特点和传统有限元法的优点.本文基于塑性增量理论,将自然单元法应用于弹塑性问题的分析计算中.为实现近似函数在非凸边界上的线性变化,采用约束的自然单元法(C-NEM)进行形函数计算.给出了增量切线刚度法求解非线性控制方程的相关公式,并对加载状态的确定和过渡状态下比例因子的计算方法等问题进行了深入的研究.编制了Von-Mises屈服准则下线性强化材料模型的二维弹塑性分析计算程序.算例分析表明,用自然单元法分析弹塑性力学问题是可行的,具有前处理过程简单、可以方便地准确施加本质边界条件等优点.     

求解弹性力学问题的有理单元法

传统的位移有限元法采用多项式形式的位移试函数,对于边数大于4的多边形单元,构造满足单元间协调性要求的多项式形式位移插值函数是一件困难的工作。本文利用逆距离权插值的思想并考虑到单元节点的分布,建立了边数大于4多边形单元上的有理函数形式的形函数。利用有理试函数,采用Galerkin法推导出求解平面弹性力学问题的有理单元法。采用有理单元法求解弹性力学问题,求解区域根据需要可以划分为任意多边形单元,极大地提高了网格划分的灵活性。有理单元法不依赖等参变换,不同单元的形函数表达形式统一,方便计算程序的编写。     

面向对象的有限元程序设计

本文按照面向对象的程序设计方法,遵循有限元分析的本质,采用Ｃ＋＋语言,建立了有关描述有限元模型的类,给出了类的描述和它的实现方法。相关的类和方法包括处理矩阵的类、节点类、单元类、材料类、形函数类等。据此编制了有限元分析的数值计算程序,并给出了一个实例。     

多边形有限元研究进展

有限元法是数值求解偏微分方程边值问题的重要方法,采用不规则多边形单元网格, 可以方便有效地模拟材料的力学性能, 又使得区域网格剖分变得灵活方便. 特别是对于复杂的几何形状, 多边形单元网格具有更大的优势. 本文对国内外有关多边形有限元法的最新进展作了初步的总结和评述, 主要以基于位移法的多边形有限元为主.论述了多边形有限元的发展历史, 给出了多边形单元上的Wachspress插值、Laplace插值和重心坐标的一些最新研究成果. 与经典有限元法形函数为多项式形式不同, 多边形单元的形函数为有理函数或者无理函数形式. 多边形单元插值形函数满足线性完备性, 可以再现线性位移场, 像经典有限元法一样直接施加本质边界条件; 插值函数在多边形的边界上是线性的,确保不同单元间的自动协调. 不同单元的插值形函数表达公式形式统一, 方便混合单元网格计算的程序编写. 提出了多边形有限元法今后需要研究的问题.      

边坡弹塑性稳定性分析的无网格法

将无单元伽辽金法(EFGM) 推广到求解弹-- 黏塑性问题, 自行编制了相应的计算程序, 并应用于土质边坡的弹塑性稳定分析. 通过无单元伽辽金法计算结果与有限元法的计算结果的对比分析, 可以发现用黏弹塑性问题的无单元伽辽金法计算程序去求解弹塑性问题是方便可行的, 本文编制的计算程序稳定性好, 收敛速度快.     

自然单元法研究进展

自然单元法是一种基于Voronoi图和Delaunay三角化几何结构,以自然邻点插值为试函数的一种新型数值方法.其既具有无网格方法和经典有限元方法的优点,又克服了两者的一些缺陷,是一种发展前景广阔的求解微分方程的数值方法.自然单元法的形函数满足插值性质,可以像有限元法一样直接施加本质边界条件,不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题.由于自然单元法是无网格方法,可以方便处理有限元方法较难处理的一些问题,例如移动边界和大变形等问题.自然单元法与其他数值方法的最根本区别于其插值格式的不同.将自然邻点插值用于Galerkin过程,就得到基于Voronoi结构的自然单元Galerkin法.自然邻点插值有自然邻点Sibson插值和Laplace插值(非Sibson插值)两种.Laplace插值比Sibson插值在计算上要简单的多,并且不论对凸的或非凸的区域都能精确施加本质边界条件.以Laplace插值为试函数的自然单元法在数值实施上比以Sibson插值为试函数的自然单元法简单.本文对基于Voronoi结构的自然邻点插值和自然单元法的基本思想作了介绍,综述了国内外关于自然单元法的研究成果,总结了自然单元法的优点和尚需解决的问题.      

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

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

面向对象的钢筋混凝土有限元非线性分析程序设计

采用面向对象的程序设计方法，结合钢筋混凝土有限元非线性分析模型，利用MFC(Microsoft Foundation Class Library)建立了有关描述钢筋混凝土有限元非线性分析的类，包括混凝土单元类、钢筋单元类、粘结单元类、非线性计算类，并给出了这些类的描述和它的实现方法，并为钢筋混凝土有限元非线性分析程序采用更合理的面向对象的方法提供了思路。     

面向对象有限元程序的类设计

传统有限元程序设计中的许多问题可以通过面向对象程序设计方法加以解决。本文介绍了我们尝试采用ＢｏｒｌａｎｄＣ＋＋２．０进行的一种面向对象有限元程序设计，着重进行有限元程序的类划分与定义，并将几个典型的类设计成类等级。面向对象语言的继承性、封装性、多态性等特性在我们的程序设计中均得到了运用。与传统的有限元程序（通常采用Ｆｏｒｔｒａｎ）相比，面向对象有限元程序更加结构化、更易于编写、更易于维护和扩充，程序代码的可重用成分更大，它为开发大型有限元分析软件提供了一条新途径     

面向对象有限元方法研究进展

综述了面向对象有限元方法的研究现状。基于国内外大量相关文献，首先评述了有限元方法的面向对象分析过程（OOA），包括对象的识别，对象属性，方法和关系的确定等等。其次涉及有限元方法的面向对象设计过程（OOD），具体讨论了人机交互，任务管理和数据管理等内容。再次讨论了有限元方法的面向对象编程过程（OOP），如编程语言的选择，面向对象有限元的应用和集成等问题。此外还延拓到面向对象方法增强有限元软件的表现，面向对象有限元的扩展等。最后展望该方法的发展动态。     

OBJECT-ORIENTED FINITE ELEMENT ANALYSIS AND PROGRAMMING IN VC++

ＩｎｔｒｏｄｕｃｔｉｏｎＳｉｎｃｅｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓｐｒｏｇｒａｍｓｄｅａｌｗｉｔｈｔｈｅｔｈｅｏｒｉｅｓａｎｄｍｅｔｈｏｄｓｏｆｔｈｒｅｅｆｉｅｌｄｓｏｆｍｅｃｈａｎｉｃｓ,ａｐｐｌｉｅｄｍａｔｈｅｍａｔｉｃｓａｎｄｃｏｍｐｕｔｅｒｓｃｉｅｎｃｅ,ｔｈｅｙａｒｅｃｏｍｐｌｉｃａｔｅｄａｎｄｅｒｒｏｒ＿ｐｒｏｎｅ .Ｄｕｒｉｎｇｄｅｖｅｌｏｐｍｅｎｔｐｅｏｐｌｅａｌｌｔｈｒｏｕｇｈｓｅｅｋｔｈｅｍｅｔｈｏｄｔｏｉｍｐｒｏｖｅｐｒｏｇｒａｍｓ.Ｔｈｅｏｂｊｅｃｔ＿ｏｒｉｅｎｔｅ…     

弹性接触问题的形状优化设计方法

设计模型复杂和结构的高度非线性等问题使得目前结构优化技术中很少涉及接触问题。面向对象的有限元法利用了面向对象编程的优点，有效地提高了编程的效率，并便于程序的扩充。在Visual C 软件平台上，开发编制了面向对象的有限元法程序，并对含接触问题的结构进行形状优化设计，重点推导并给出了含接触问题有限元法的灵敏度全解析分析的表达式，包括位移和应力的灵敏度分析、接触条件的灵敏度分析及相应的边界支配方程等。最后通过算例，验证了所提出的优化方法的正确性。     

连续胶体材料化学-力学耦合有限元理论及其在ABAQUS中UEL程序实现

基于哈密顿原理,得到水凝胶的化学-力学耦合控制方程的等效积分形式和有限元形式。在整体坐标系下推导出用形函数表示的化学-力学耦合应变矩阵和单元刚度矩阵,并且得到在局部坐标系下的离散化形式。结合ABAQUS软件,编制了用户单元子程序UEL,通过数值算例验证了所开发单元的正确性,为在ABAQUS软件中实现各种耦合问题的有限元UEL编程提供了参考依据。     

基于PVM的网络并行有限元面向对象的设计

子结构是有限元并行计算常用的一种方法,本文采用面向对象的方法,首先对子结构进行了面向对象的设计,得到了其类层次结构图;然后针对工作站网络有限元并行计算环境。提出了基于PVM消息传递平台上的Shadow—Mirror数据传输模型,该模型在有限元并行计算数据传输时,充分发挥数据面向对象的特性,采用设置数据缓冲区、短消息合并等方法以缩短数据通信时间,并据此编制了相应的程序。计算结果表明,使用文中提出的面向对象的Shadow—Mirror数据传输模型可以得到较为理想的并行加速比,而且随着问题规模增大,并行加速比增高。本文研究内容为进一步开展基于工作站网络的并行有限元研究提供了一个可参考的基础。     

含缺陷结构的塑性极限分析

结合极限分析中的数学规划理论和有限元技术,提出了三维含结构极限分析的数学规划方法,并采用罚函数法引入塑性不可压条件,对于考虑多组独立变化载荷作用的情况,提出了加载路径射线辐射求解方案,并基于这种射线射状的加载路径,推导了多组载荷联合作用下结构塑性极限上限分析的数学规划格式,编制了相应的有限元程序,文中的数值了该方法的正确性与有效性。     

空间网格结构几何非线性有限元分析方法的研究

在对 T.L/ U.L.法和 C.Oran梁柱单元有限元法进行系统研究比较的基础上 ,推导了结合以上两种理论的几何非线性有限元列式。C.Oran梁柱单元在分析轴力为主的结构时具有非常高的效率 ,但在分析纯弯曲问题时却存在很大困难。本文采用 VC 6 .0编制了面向对象的非线性有限元程序 ,对其进行了考题验算分析 ,得出了有用的结论。