有限元网格节点优化排序方法研究

本文采用的有限元网格节点优化排序方法是一种归纳,演绎方法,通过对在网格半带宽最小时节点排列的一些特征进行了分析,从其中归纳的出三条规律,通过大量的算例来验证一和可靠性,结果表明,对于任意划分有限单元网格,该方法通过对节点的重新排序,能够优化出最佳的半带宽值,大大减少了计算机内存有空间,提高了运算效率。     

一种有限元模型修正中的参数选择方法

有限元模型修正面临的首要问题是待修正参数的选择.目前主要是基于灵敏度分析的结果进行参数选择.本文从结构固有频率分析的能量法出发,论证了某阶固有频率主要受到同阶振型中发生弹性变形部位刚度的影响,据此提出一种待修正参数的选择方法,以振型中弹性变形为依据,选择弹性变形较大部位的参数作为待修正参数.用一个梁模型仿真验证了该方法的正确性,并将其应用于系杆拱桥的有限元模型修正中,修正后有限元计算结果与试验结果的最大误差缩小至3.3%.     

有限元计算中疏密网格过渡方法研究

工程计算中出于节省计算量的目的,往往需要在一个有限元模型中布置粗细不同的网格。为保证计算结果的准确性,必须保证网格突变情况下的位移协调问题。本文工作之一是在强天驰界面过渡单元的基础上,引入虚拟节点和子单元,在子单元中应用节理元思想,提出了基于最小势能原理的弹簧节理单元法。简化了积分运算,避免了精度要求极高的坐标转换,从而提高了方法的精度和实用性;二是提出了基于位移约束的主从自由度法,简便实用,只需简单的矩阵运算即可实现。两种方法均实现了不同尺寸网格间位移的协调性和刚度的匹配,从而使之满足有限元收敛准则,且生成的刚度阵具有对称性及带状性。算例证明两种方法精度良好,并可方便地应用于求解大规模工程问题。     

基于多层模拟方法的双层梁断裂有限元分析

本文提出一种用于含分层的双层梁线弹性断裂分析的有限元方法．将上下子梁均模拟为多个子层,采用只有平动位移自由度的新型梁单元,假设单元内的位移沿纵向和横向均线性变化,推导了该单元的单元刚度矩阵．将开裂部分和未开裂部分的子梁进行单元刚度矩阵组装,施加相应的等效结点力,得到整体平衡方程,并结合边界条件进行求解．为验证该方法的有效性和精度,开展非对称双悬臂梁(Asymmetric Double Cantilever Beam, ADCB)和单臂弯曲梁(Single Leg Bending, SLB)试样的断裂分析,利用虚拟裂纹闭合技术(Virtual Crack Closure Technique, VCCT)得到了试样的能量释放率及其分量,并将求得的结果与解析解和二维有限元解进行对比．计算结果表明,相对于传统双层模拟方法,该多层模拟方法能够精确、高效地计算各类梁试样的能量释放率及其分量,并且无需引入界面连续条件．     

并行有限元网格生成方法及其应用

大型工程数值仿真中,在前处理阶段需要生成千万甚至亿量级的网格,传统的串行网格生成方法由于内存和时间的限制,难以处理如此规模的网格。针对此问题,本文提出了一种大规模网格并行生成方法。首先基于推进波前法对几何模型进行初始体网格划分,接着利用图论理论进行区域分解,并通过表面单元恢复保持其几何精度,然后通过分裂法进行网格的并行生成。将所述方法应用到实际大型工程数值仿真前处理阶段,结果表明所述方法可以获得较好的并行效率,同时所产生的网格质量可以满足后续计算需要。     

保持特征的三角形有限元网格简化方法

基于二次误差的网格简化算法,提出了一种三角形有限元网格模型简化算法。在网格质量评估中考虑折叠代价,避免过于狭长的三角形出现,保证新生成单元的质量。同时,通过改进特征检测机制,加入用户控制,对模型上的重要细节特征进行精细定位,使算法适用范围更广泛,简化结果更能满足网格简化过程中细节特征的保持要求。     

一类混凝土力学参数的迭代多尺度有限元预测方法

根据多级配骨料混凝土特点,给出了一种迭代多尺度有限元方法预测其等效力学参数。本文首先介绍了多级配骨料混凝土的材料特点并给出了计算混凝土多尺度模型,然后基于多尺度方法介绍了计算混凝土力学参数的计算程序,最后针对小湾大坝中混凝土给出了算例说明了此种方法在预测混凝土等效力学参数的有效性。     

多尺度六面体单元网格连接方法研究

提出多结点六面体单元的结点形函数构建方法,解决了数字化分析时两个多尺度有限元网格在接触界面由于单元结点不匹配而导致的结点属性不能连续传递问题.首先将全局坐标下的六面体单元及其表面上多结点通过等参逆变换转成局部坐标下的规则六面体单元及其表面上多结点;在规则单元中,以每个结点为基点,分别沿三个正交的局部坐标方向在单元内寻找...     

复杂三维组合曲面的有限元网格生成方法

提出一种基于映射法的复杂三维组合曲面的有限元网格全自动生成方法。通过引入虚边界解决了闭合曲面在参数域中边界不完整的问题；通过调节虚边界提高了复杂组合曲面网格生成的质量。改进二维多边形区域的裁减算法，解决了闭合曲面在参数空间中的边界环形成问题。对曲面片公共边界进行统一离散化处理，以满足有限元网格的相容性要求。以边界表示(B—Rep)数据结构为基础，实现了组合曲面全自动网格剖分的总成算法．改进了曲面网格剖分布点算法，并结合局部连接、诊断交换等技术，优化了网格的整体质量。     

基于有限元弹性变形运算的全六面体网格生成方法研究

以映射法为基础并结合网格划分经验,提出了原理简单的六面体网格生成新办法.该方法根据物体轮廓选择初始网格,设定表面结点强制变形到目的曲面,经由有限元弹性计算确定内部节点的位置.在检查全体单元质量以后,调整畸形单元从而生成目的网格.通过为一个复杂的马头门模型构建全六面体网格,最后证明了本文所述方法的可行性.     

有限元模型阻尼特性的复参数修正方法研究

阻尼对于结构动力学响应具有重要的影响,但有限元模型一般很难对阻尼特性进行精确建模.基于实测频响函数,研究了一种有限元模型阻尼特性的复参数修正方法.以待修正区域各单元质量、刚度矩阵的比例修正系数为复修正参数,建立了单元矩阵比例修正的灵敏度方程直接算法,并对比分析了复修正参数与不同阻尼特性之间的数学关系.以六自由度集中参数模型和25杆平面桁架模型为例,验证了复参数修正方法在阻尼特性修正中的有效性.     

结构电磁有限元分析的网格依赖性研究

分析了基于有限元方法的超材料电磁性能仿真分析结果对有限元网格的依赖性.通过对仿真结果的分析,发现基于S参数的材料性能参数反演结果对S参数的计算准确度的依赖性巨大,特别是处在谐振频域的S参数误差会导致等效电磁参数的性态误差.针对基于求解频率的网格自适应技术对远离求解频率的谐振频域的分析结果的准确度改善不够这个问题,需要确定合适的求解频率来执行网格自适应.提出了一种改进电磁性能仿真分析准确度的方法,其基本思想是首先给定一个初始的求解频率值,根据Kramers-Kronig关系反演等效电磁参数,进而获得谐振频率的近似值,并把此值作为求解频率来实现网格自适应以获得合适的网格剖分,最后经扫频获得S参数并进行等效电磁参数的反演.数值结果表明,该方法可有效提高仿真分析和参数反演的准确度.     

三维有限元六面体网格几何自适应再生成方法

基于几何自适应的六面体网再生成方法,提出了一种以栅格法为基本方法,对有限元网格再划分过程中的网格再划分标准、几何形状的获取以及新旧网格之间物理量的传递等关键问题进行了研究.重点介绍了基于几何自适应的六面体网格再生成算法,首先对旧网格实体模型进行识别,根据实体模型几何特征建立加密源点信息场;然后采用栅格法生成核心网格并对...     

多尺度材料模拟的自适应有限元方法研究

研究了材料模拟中一类新型耦合多尺度的自适应有限元方法. 采用微观分子动力学耦合宏观有限元的桥尺度方法来模拟材料破坏的前期行为,其中宏观有限元计算推广到了一般非结构三角形网格. 材料破坏形成后,停止微观尺度的计算,它的进一步发展和演化通过一个宏观模型来描述,采用自适应有限元方法来求解这一宏观模型. 其中,后验误差估计的基础是变分多尺度理论,即自适应网格加密是基于粗尺度上残差分布和细尺度上单元Green's函数. 计算中采用了破坏准则来模拟材料的断裂. 数值实验表明了方法的有效性.      

基于实验模态参数的结构多层结合部参数识别

本文从一般情况出发，提出了一种利用结构实验模态参数与子结构有限元模型，识别结构多层结合部刚度与阻尼参数的方法。当结合部模态向量分别为可测，部分可测或不可测时，给出了相应的识别算法，特别地，当结合部模态向量为可测时，其识别算法无任何矩阵求逆运算。算例表明，本文方法具有满意的识别精度。     

有限元网格全自动生成的四分法

本文论述了有限元网格全自动生成的四分法。这种方法适用于任意形状的平面结构,且形体内部可以带有若干个任何形状的孔洞;自动生成的单元以四点元为主,边界辅以三角形元;该法还允许有限元网格局部加密,加密区交接处引入五点元,网格加密部位可以有若干处。有限元网格全自动划分算法在CAD工作站上实现,用人机交互会话方式控制输入/输出操作,方便、灵活、高效,可做为前处理器与任何一个有限元分析程序系统相联接。本文所采用的方法和数据结构,可以推广到带孔洞的在意形状空间形体,实现对更广泛的三维结构的有限元网格全自动生成,     

一种微机械有限元仿真的多重网格预处理方法

在微机械结构设计中常用有限元方法进行仿真,然而,通用的三维有限元仿真算法在进行有限元分析时耗时较长,影响了仿真效率。针对常规MEMS工艺仅有二维成形能力的特征,提出了将二维多重网格方法融入条件共轭梯度法的预处理的方法。并以一种三明治式微加速度计的ANSYS模态分析为具体算例,提出了用多重网格思想对网点分布进行插值以求在不太提高条件数的情况下提高对这类低频非局部振型的模态分析问题的求解精度的二维预处理算法。所提出的多重网络的网络划分预处理算法将预处理前后的同精度分析计算时间由206 min减少到38 min,但相对精度从0.018%降低到了0.038%。这种算法可用于二维成形结构的静力分析和模态分析并提高设计效率。     

循环对称结构的模态参数计算方法

利用循环对称结构部件模态的等同性，推求了循环对称结构中的混合型动态子结构法的综合方程。     

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

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