动力学有限元
- 格式:docx
- 大小:155.21 KB
- 文档页数:10
有限元第9讲动力学问题有限单元法动力学问题是指研究物体在运动中的受力和受力作用下的运动状态,常见的应用是结构工程学中的振动分析。
有限单元法是解决结构工程学中动力学问题的常用方法之一。
本文将介绍动力学问题和有限单元法的基本概念,并介绍其应用。
动力学问题的定义动力学是研究质点或刚体运动情况的分支学科,在结构工程学中是指结构在做振动时所受的力和运动状态。
动力学问题可以分为两种类型:稳态动力学问题和非稳态动力学问题。
稳态动力学问题是指结构在振动状态下所受的恒定力,而非稳态动力学问题则是指结构所受的变化的力,例如冲击力或地震力。
动力学问题的求解包括两个方面:一是确定受力情况;二是求解结构的运动状态。
确定受力情况通常需要通过实验或计算确定,求解结构运动状态则可以通过有限单元法来解决。
在结构工程学中,动力学问题的应用非常广泛。
例如,建筑物抗震设计需要对建筑物在地震作用下的反应进行分析,桥梁工程需要对桥梁在行车作用或风力作用下的振动响应进行分析。
有限单元法的基本概念有限单元法是一种将结构离散成若干小单元的数值分析方法,将结构分割成细小的单元,每个单元内部假设为均匀且连续的,通过对单元本身的运动状态进行求解,进而推知整个结构的运动状态。
有限元法用于解决的问题包括静力学问题、动力学问题、热力学问题和流体问题等。
有限单元法求解动力学问题的步骤主要包括如下几个步骤:1.离散化:将连续结构离散化成有限的小单元,每个单元内部运动状态通过定义一定数量的节点来确定。
2.建立单元的动力学方程:根据单元的形状和材料性质,建立单元的动力学方程,并计算单元的振动特性,例如频率和模态。
3.组装单元的方程:将单个单元的方程组装成整个结构的方程。
4.边界条件的处理:利用结构的边界条件(例如支撑、铰支等),将结构自由度减少到实际问题所需要的自由度。
5.求解结构的运动状态:通过求解整个结构的方程,得到结构的运动状态。
6.后处理:根据求解结果,进行结果的可视化和分析。
多体动力学和有限元关系多体动力学和有限元关系多体动力学和有限元关系是两个在工程领域中被广泛应用的概念。
多体动力学主要描述了多个物体之间相互作用的力学行为,而有限元是一种数值分析方法,用于近似求解连续物体中的力学问题。
在本文中,将探讨多体动力学与有限元的关系以及它们在工程设计中的应用。
1. 多体动力学基本原理多体动力学是研究多个物体之间相互作用的力学学科。
在多体动力学中,物体被视为刚体或弹性体,它们之间通过力或力矩进行相互作用。
多体动力学的研究对象包括机械系统、流体系统和电路系统等。
通过分析物体之间的相互作用,可以得到系统的运动学和动力学方程,从而预测系统的运动和响应。
2. 有限元方法概述有限元方法是一种近似求解连续物体中力学问题的数值分析方法。
它将连续物体离散为有限数量的子区域,称为有限元。
每个有限元代表一个局部区域,在该区域内的物理行为被近似为一组简单的函数。
通过在每个有限元内应用力学原理,可以建立有限元方程组,并通过求解该方程组得到连续物体的近似解。
有限元方法的优势在于可以处理复杂几何形状和边界条件,并且可以灵活地模拟材料的非线性行为。
3. 多体动力学与有限元的关系多体动力学与有限元方法在某种程度上可以看作是相互补充的。
多体动力学主要关注物体之间的相互作用和运动规律,而有限元方法则更注重求解连续物体内部的力学问题。
在一些对物体之间的相互作用和约束较为复杂的情况下,可以将多体动力学与有限元方法相结合,以获得更准确的结果。
4. 多体动力学与有限元的应用多体动力学和有限元方法在工程设计中具有广泛的应用。
在机械系统设计中,可以使用多体动力学分析来评估机械系统的动态性能和稳定性,而有限元分析则可以用于优化机械结构的刚性和耐久性。
在车辆工程中,多体动力学可以用于模拟车辆的悬挂系统和转向系统的运动特性,而有限元分析可以用于优化车身结构的强度和刚度。
在建筑工程中,有限元方法可以用于评估结构的抗震性能,而多体动力学可以用于研究大楼在地震中的动态行为。
分子动力学的有限元长时间计算研究分子动力学的有限元长时间计算研究一、引言分子动力学(MD)是一种模拟和研究原子和分子在特定条件下运动规律的方法,广泛应用于材料科学、生物化学、药物设计等领域。
有限元长时间计算是指利用有限元方法对分子系统进行长时间的模拟和计算,以研究复杂的分子动力学行为和性能。
本文将对分子动力学的有限元长时间计算研究进行全面评估,并通过逐步深入的探讨,帮助读者更全面、深入地理解这一主题。
二、分子动力学的基本原理1. 分子动力学的基本公式与算法分子动力学模拟的基本公式是牛顿运动方程,根据原子间的相互作用力和势能函数,利用数值算法进行时间演化。
有限元方法是一种用数学方法将连续体划分成离散单元的方法,结合分子动力学,可以更精确地模拟原子和分子的行为。
2. 分子动力学的模拟条件和约束在进行分子动力学的模拟时,需要考虑温度、压力、边界条件等影响因素,并通过施加约束条件来模拟不同环境下的分子系统。
三、有限元长时间计算的原理和方法1. 有限元方法在分子系统中的应用有限元方法是一种数值计算方法,通过离散化和逼近的方式,可以有效地模拟复杂的分子系统,并计算其长时间行为。
2. 长时间计算的数值稳定性和精度在进行有限元长时间计算时,需要考虑数值稳定性和精度的问题,以确保模拟结果的准确性和可靠性。
四、分子系统长时间行为的研究1. 原子和分子的动力学行为利用有限元长时间计算方法,可以研究原子和分子在不同条件下的动力学行为,如振动、扭转、碰撞等。
2. 分子系统的热力学性质通过长时间计算,可以研究分子系统的热力学性质,如热容、热传导等,为材料科学和化学工程领域提供重要参考。
五、总结与展望本文通过对分子动力学的有限元长时间计算研究进行深入探讨,全面评估了其在原子和分子行为研究中的重要性和应用前景。
有限元方法的应用为分子系统的模拟和计算提供了更精确和可靠的手段,长时间计算的研究将为材料科学、生物化学和药物设计等领域的发展提供重要支持。
分子动力学方法与有限元
分子动力学方法和有限元方法是两种在物理和工程领域广泛使用
的数值模拟方法。
分子动力学方法是一种计算分子尺度下物质动力学行为的方法,
它基于牛顿力学和统计力学理论,通过求解运动方程模拟分子在时空
上的演化。
通过分子间相互作用力的计算,可以获得物质的力学、热学、电学、光学等性质。
分子动力学方法可以用于研究多种物质行为,如材料变形、热传导、气体扩散,以及生物分子的构象和动力学。
有限元方法是一种计算连续介质力学和热学问题的方法,其基本
思想是将连续介质离散化为有限个小的单元,对每个单元进行有限次
的近似,然后通过数值求解算法得到整体系统的数学模型。
有限元方
法广泛应用于工程领域,如结构力学、流体力学、热传导等问题的数
值模拟。
其优点是能够较为准确地计算复杂几何形状的结构的应力、
变形、热传导等问题,而且计算结果有较高的可信性。
尽管分子动力学方法和有限元方法的基本理念不同,但两种方法
都是通过建立数学模型和求解模型方程来描述物理过程的,因此它们
在模拟物质性质和行为方面都具有很大的优势。
在实际应用中,一些
复杂的问题需要同时运用两种方法进行研究,例如材料的多尺度模拟、纳米材料的性能研究等。
刚体动力学是研究刚体运动的力学学科。
刚体是指形状和大小在运动过程中保持不变的物体,刚体动力学研究刚体在受力作用下的运动规律和动力学特性。
刚体动力学主要包括以下几个方面:
运动学:研究刚体的位移、速度和加速度等与时间的关系,描述刚体的运动状态。
动力学方程:根据牛顿第二定律,建立刚体的动力学方程,描述刚体受到的力和加速度之间的关系。
转动运动:研究刚体绕固定轴进行转动的规律,包括转动惯量、角速度、角加速度等的计算和分析。
能量与动量守恒:研究刚体运动过程中的能量守恒和动量守恒定律,用于分析刚体的碰撞、旋转和平移等情况。
有限元方法(Finite Element Method,简称FEM)是一种数值计算方法,广泛应用于工程和科学领域,包括力学、结构分析、流体力学等。
有限元方法将连续的物体或结构分割成有限数量的小单元,通过求解这些小单元的力学方程,得到整个物体或结构的力学行为。
在刚体动力学中,有限元方法可以用于建立刚体的数学模型,通过将刚体分割成有限数量的单元,利用数值计算方法求解刚体的运动和力学响应。
这种方法可以有效地模拟复杂的刚体运动和受力情况,帮助分析和优化刚体系统的设计和性能。
有限元方法在刚体动力学中的应用包括刚体结构的动力学分析、碰撞和撞击的模拟、机械系统的优化等。
它提供了一种灵活、高效的数值计算工具,用于解决刚体动力学问题和工程实践中的设计和分析任务。
结构动力学中基于有限元方法的动力响应分析结构动力学是研究结构在外部载荷作用下的振动特性和动态响应的学科。
大型工程结构系统的复杂性和非线性特性给结构动力学分析提出了挑战,而有限元方法则成为求解这种非线性响应的一种重要手段。
在本文中,我们将探讨结构动力学中基于有限元方法的动力响应分析。
1. 有限元方法有限元法是一种现代数值计算方法。
它是把连续物体分割成多个单元,通过单元间的相互作用关系求解结构的内部应力、变形和各种响应的数值方法。
有限元法的基本思想是把复杂的整体结构分解成有限数量的小单元,并对每个小单元进行数学模型分析。
通过求解这些模型,可以推导出整个结构的力学特性和响应情况。
2. 结构动力学中的有限元方法在结构动力学中,有限元方法也是一种重要的分析方法。
一般来说,结构动力学的有限元模型应包括结构的物理性质、载荷和边界条件等。
在构建有限元模型之前,需要对结构几何形状进行测量和描述,然后将结构分割成有限数量的单元,每个单元都有一组节点和自由度,节点之间的相互作用关系是通过构建单元刚度矩阵来实现的。
在建立了完整的有限元模型后,可以采用不同的求解算法,如静力求解和动力求解进行解析求解。
3. 动力响应分析在有限元法中,一般需要对结构进行动力响应分析。
动力响应分析的主要目标是确定在特定载荷下结构的动态响应情况。
动态响应包括结构的位移、速度、加速度、应力和应变等。
这些响应都对结构的安全性、稳定性和寿命等方面产生影响,因此需要进行充分的动态响应分析。
在动力响应分析中,一般采用有限元模型接触外部载荷模拟结构振动情况。
通过分析结构的固有振动模态和相应的频率响应,可以计算出特定载荷下结构的动态响应。
在实际分析中,通常需要考虑多种载荷并结合计算机模拟技术实现更为准确的动态响应分析。
4. 结论本文简要介绍了结构动力学中基于有限元方法的动力响应分析。
有限元法是一种现代数值计算方法,它可以将结构分割成多个小单元,进行数值模拟,计算结构内部应力、变形和各种响应。
多体动力学和有限元关系多体动力学和有限元分析的关系多体动力学和有限元分析是两个重要的工程领域,它们在实际问题的求解和设计过程中起着至关重要的作用。
本文将探讨多体动力学和有限元分析之间的关系,并介绍它们在实践中的应用。
多体动力学是研究多体系统运动的学科,多体系统是由多个物体组成的力学系统。
多体动力学研究物体之间的相互关系和运动规律,包括质点的运动、刚体的运动和变形等。
多体动力学可以应用于机械工程、航空航天工程、土木工程等领域,用于分析和设计各种机械装置、结构和系统。
有限元分析是一种数值计算方法,用于求解复杂结构的力学问题。
有限元分析将连续体划分为有限数量的单元,然后在每个单元内对物体的运动和应力进行数值计算。
通过求解单元之间的相互作用,可以得到整个结构的运动和应力分布。
有限元分析适用于求解各种结构的静力学和动力学问题,包括弹性变形、稳定性分析、疲劳寿命预测等。
多体动力学和有限元分析之间有着密切的联系和相互依赖。
首先,多体动力学可以为有限元分析提供准确的质点和刚体的运动学参数。
在有限元分析中,需要将结构划分为单元,并对每个单元的运动状态进行计算。
多体动力学可以提供这些单元的位移和速度等数据,为有限元分析提供边界条件和初始条件。
其次,有限元分析可以为多体动力学提供准确的材料性能和载荷边界条件。
在多体动力学中,物体之间的相互作用和外部载荷会影响系统的运动状态。
有限元分析可以通过数值计算得出这些相互作用力和载荷,为多体动力学提供准确的力学参数。
最后,多体动力学和有限元分析的结果可以相互验证和校准。
在设计和分析过程中,可以通过对比两种方法的结果来评估其准确性和可靠性。
这种相互验证和校准可以提高工程设计的质量和可靠性,减少开发和生产过程中的风险和成本。
综上所述,多体动力学和有限元分析是两个相互依赖且相互补充的工程领域。
它们在工程设计和分析中发挥着重要的作用,可以提供准确和可靠的结果。
多体动力学和有限元分析的结合应用可以为各种工程问题的求解和设计提供有效的工具和方法。
近场动力学与有限元方法的区别近场动力学和有限元方法是两种常用的数值模拟技术,它们分别采用不同的算法和数学方法来对物理问题进行求解。
本文将从不同的角度介绍这两种方法的区别。
1. 原理近场动力学方法是一种基于分子动力学的数值模拟技术,通过对物质分子间的相互作用力的计算来模拟物体在微观层面上的运动。
其基本原理是根据分子级别上的物理成分,包括分子的质量、电荷、速度、能量等,建立分子的运动学方程和相互作用力的力学方程,并将其离散化求解。
该方法可以有效地模拟固体、液体、气体等基本物质的力学和热力学性质,并能够在微观层面上捕捉物质的瞬态状态和动态行为。
有限元方法是一种基于数学理论的数值模拟技术,主要用于求解物体在宏观尺度上的力学问题。
它基于连续介质力学的基本原理,将物体离散化成许多小部分,同时将其与周围环境的作用力相互作用,然后求解其静力学或动力学方程。
该方法适用于多种物理场问题,例如结构分析、流体力学、电磁场问题等。
2. 应用领域近场动力学方法主要应用于电子结构、材料科学、地球科学等领域,特别是在模拟材料缺陷、断裂、疲劳、热膨胀等行为方面具有非常好的效果。
它可以模拟区域尺度在几纳米到几微米之间的物质和介质,并可以考虑物质的非线性和非均匀性。
有限元方法在机械、土木、航空航天、汽车等工程领域的应用非常广泛,可以对如钢筋混凝土梁、飞机机翼、汽车车身等结构进行力学计算和模拟,并可以预测其强度、刚度、振动等性质。
3. 离散化方法在近场动力学方法中,物体被离散成小的粒子,这些粒子通过相互作用力发生相互作用,并在微观层面上发生运动。
而且,物体的运动受到长程协作以及短程反应力的影响。
因此,近场动力学方法注重离散效果和相互作用力的准确计算,以获得更精确的结果。
在有限元方法中,物体被分割成一系列小单元或网格,将力学问题转化为求解一个连续的偏微分方程,其中物体的行为由节点或区域内的变量来描述,通常是位移分量、应力分量或变形分量等。
6.2结构动力有限元法理论与模型一、基本原理在实际问题的求解中,应用最广的是基于位移的有限元素法。
此法的基本思想是把本来为连续的工程结构分割成在结点上相联的单元组合体。
取这些结点的位移为基本未知量,并假定每个单元中的位移用单元位移函数来描述,这实质上是假定了单元的模态。
在此基础上,利用能量变分原理进行单元分析的全结构分析,得到全结构的振动平衡方程,从而把连续体的动力学问题化为多自由度系统的振动问题。
有限元动力分析的基本过程是首先将工程结构离散化,通过选择合理的单元确定出分析模型,在此基础上选择位移函数,进行单元分析,确定单元的刚度、质量、阻尼、载荷矩阵,再经过坐标变换,通过能量变分原理,进行全结构分析,建立系统的振动平衡方程。
最后运用有限元数值方法进行方程的求解。
结构动力有限元法采用的单元位移函数与静力分析相同,基本原理和求解过程也与静力分析相同,不同之处仅在分析模型的确定与运动方程的建立方面。
二、动态分析模型的确定由于结构动态分析中除考虑弹性力外,还要考虑惯性力和阻尼力,其运动方程是常微分方程组,所以动态分析的复杂程度高,计算工作量大,有限元分析模型要尽量精炼、简单。
1.模型确定的基本原则•分析模型应与分析的目的相适应。
动力分析的目的各不相同,有的是为了提供固有特性计算动态响应或供控制系统用;有的是为了舱内提供振动环境。
不同的目的,通常要求不同的模态数与计算精度。
显然,用于估算基本固有频率的模型应当比计算冲击响应的模型简单。
用于设计计算的模型应当比用于校核计算的模型简单。
•分析模型要与选用的计算工具与计算条件相适应。
计算机软件种类日益丰富,选择分析模型要与所用程序、所用计算机容量相适应。
如对于容量大的计算机,可选用较为复杂的有限元模型,而对于容量小的计算机则在能反映结构动态性能的前提下尽量简化模型,使求解规模尽量小。
对于大模型,可选用子结构模型,采用模态综合方法求解。
应注意, 不一定模型愈精细精度就愈高。
模型愈复杂,往往带来了更繁杂的运算。
虽然模型误差小了,但计算误差加大。
不恰当的精细模型反而得到不佳的结果。
•模型应正确反映结构的实际特性。
一个具体结构的动态特性,主要取决于质量、刚度的大小与分布,取决于结构边界条件与阻尼特性。
因此,模型应尽量保持整个飞行器、甚至各部件的质量、质心位置不变,保证构造的刚度特性和传力路线基本不变。
真实反映分析对象的边界条件与阻尼特性。
因此,选取的单元应保持其几何形状、受力特性、变形特性与实际结构的几何形状特点、受力传力特点、变形特性相一致。
例如,对于薄壁结构翼面或舱段,不应选用板、壳单元来模拟蒙皮。
2.影响分析模型的主要因素(1)刚度分布。
飞行器属于三维物体,但根据动态行为特征对于具有互相正交的两个对称面的飞行器,有时可按二维处理;对于轴对称性强的导弹或发射器,当长细比很大时,往往可按一维模型处理。
但是当具有明显的局部非对称性或当因内部装载而引起严重耦合现象时,则需考虑采用三维模型。
刚度的具体模型可分为:⑻ 刚体一一刚度视为无穷大,如舱体内部设备在进行全弹固有特性分析时可视为刚体;(b)集中刚度模型一一例如运载器级间连接接头或发动机连接支座处、设备支承处等简化为集中线弹簧或集中扭转弹簧(如图6-1所示);图6-1 刚度集中分布模型示意图1 —板;2 —集中质量;3—线弹簧;4—梁。
(c)分段连续模型(有限元模型)——对于刚度变化复杂,又不连续变化的复杂结构,常采用分段连续体有限元模型(如图6-2所示)图6-2 刚度分段连续模型示意图1 —板;2—线弹簧。
(d)连续分布模型——例如有限元素法的单元可采用杆、梁、板、壳等连续体。
有限元模型实际上是上述各种典型刚度模型的复合。
如多级运载器可简化为分段连续、各段间用集中刚度连接的复合模型(如图6-3所示)。
图6-3 动态分析模型示意图 (a )集中参数模型;(b )有限元模型。
1 —板;2—梁;3 —集中质量。
此外,舱体内部压力、轴向载荷、定轴旋转等,将在飞行器数学模型中引起附加修正刚度。
质量分布。
将质量人为集中到选定的结点上。
质量矩阵是一个对角矩阵,其形式为当质量均匀分布时,可将质量平均分配给各相关结点上。
如果质量分布不均匀(如图6-4所示), 较简便的方法是酌情规定各结点所分担的区域,然后再把各区域质量分配给各结点。
集中质量模型图6-4 集中质量分配模型因为此种质量矩阵是一个对角线矩阵,便于一维存贮和进行对称分解。
对于结构本身集中质量占相当优势的系统,或网格划分较细的系统可采用此模型。
(b) —致质量模型。
此类单元质量矩阵为a式中「------- 质量密度;匸——单元位移函数矩阵。
它之所以被称为一致质量矩阵是因为建立刚度矩阵和质量矩阵所用的位移插值函数是一致的。
把所有单元质量矩阵叠加起来得到整个结构的总体质量矩阵。
一致质量矩阵总是正定满秩的。
飞行器往往存在非结构质量。
例如贮箱中液体晃动质量,对它无论是选用当量摆模型,还是采用弹簧质量系统,都要确定有效质量,一般以集中质量方式处理。
许多复杂结构常采用一致质量与集中质量的复合模型。
例如,图6-5所示为“土星V'所采用的复合模型。
质量模型的选择通常要结合模型网格的疏密、单元的种类及结构特点,酌情决定。
一般说来,单元划分较细、结构包含的集中质量(例如内部设备或配重块)多,则采用集中质量为宜。
如果能用位移函数矩阵丁组成动态时单元的真实变形,则采用一致质量矩阵得到的固有频率与振型比较可靠,而且接近真实值的上界。
与一致质量矩阵相比,集中质量矩阵的形成以及以后的运算都比较方便,存储量少,但给出的固有频率往往偏低,振型计算的误差较一致质量矩阵的大。
总之,两种质量矩阵在计算精度上的优劣很难得出明确的结论。
(3) 边界条件。
它包括全系统的边界条件和多级发射器级间的边界条件。
由于振型、频率对于边界条件很敏感,因而必须正确模拟边界条件,才能准确求得实际飞行器的动力分析结果。
自由、固支边界条件较易解决,难处理的是弹性支撑条件。
发射器或导弹与发射架之间、其各级之间连接边界均 属于弹性连接边界。
这些边界的支撑刚度对固有特性及动态响应的影响较大,一般均需通过地面静力 实验来解决。
总之,一个好的分析模型的确定,依赖于综合考虑分析对象、采用的解法以及精度要求,综合平 衡有效性与经济性,经过实验验证逐步修正,才能完善、合理地解决。
另外,对同一个飞行器,其飞 行阶段不同,任务要求不同,一般要采用不同的分析模型。
如图6-5所示,其中30个自由度的集中质 量一维模型,是在初步设计阶段和确定实验要求时采用的模型;而300个自由度的三维一维复合、集中 质量一致质量复合模型,则是为确立控制系统敏感元件的位置和计算固有特性而确定的分析模型。
今 后,我们将连续模型以外的其他模型称为离散模型,因为它们都是用离散参数来描述动态性能的。
图6-5 土星-V 在不同设计阶段的动力学模型(a ) 土星-V 火箭;(b ) 梁-杆模型;(c ) 梁-杆-1/4 壳模型;(d ) 1/4 壳模型;(e )三维模型。
(3) 合理地进行结点布置与单元划分在具体实施模型化工作时,必须正确布置结点或分割线、分割剖面的位置,确定有限元网格,实 现工程结构的离散化。
通常情况下,结点、分割线或分割剖面应置于结构几何形状、载荷以及材料特 性突变处。
因此,结点位置宜取在以下位置:•结构几何形状的拐点、结构开口处、厚度突变处; •载荷作用点; •应力集中点;•结构的约束点或支承处;•结构单元(部件或零构件)之间连接点; •结构中主要受力元件之间的相交点(连接点);riihiuAAA|心I 曲"由曲felhll Illi•要求输出位移或应力的点。
在进行网格划分时,首先应选择合理的单元。
可供选择的单元很多。
按形状分可分为一维、二维和三维单元,也可分为直线边单元和曲线边单元。
直线边单元计算简单,曲线边单元容易较好地拟合复杂几何形状,但计算工作量大。
另外,单元也有高低精度之分。
高精度单元虽然计算精度高,但自由度多,相应的计算规模大。
确定模型时,应该在满足精度要求的情况下,尽量选用简单的单元。
而在具体划分网格时应考虑以下内容:•结构特性(几何形状、材料性能、边界条件等)和外力变化剧烈部位网格应相对较密。
•二维和三维单元在各方向的长度比应尽量接近于1,避免畸形单元。
例如三角形单元的边长比不应超过1 : 3,多边形单元的相邻两边夹角不要接近180。
等。
•相邻单元的尺寸不宜相差太大。
•选取不同元素相组合时,应考虑同一结点上不同元素间自由度的协调。
•单元的疏密要视需要而定,除应力梯度大的区域外,不应片面追求密网格,以减小计算规模。
对于大型复杂结构,可考虑先进行整体粗网格计算,再根据粗算结果,对应力梯度大的区域进一步加密网格进行“细算”。
•在条件合适时,应充分利用结构的对称性确定网格与模型,以减小计算规模。
•网格划分的结果应保证各单元的尺寸与单元的力学特性相符合。
例如,对于普通的板单元,其边长应远大于板厚;对于梁单元,通常梁的长度应远大于其截面尺寸。
如果出现了粗短梁,则应在单元分析中考虑剪切效应。
三、运动方程的建立采用瞬时最小势能原理建立运动方程。
1.单元分析取单元内任意一点的位移向量丿、为基本未知量,它是位置」和时间:的函数。
由有限元素法知曲,f"(如f)而且cr= D E 式中J ——单元位移插值函数矩阵,即位移函数矩阵;%")—单元结点位移向量;「[丁——分别为应变、应力变量;(6-1 ) (6-2 ) (6-3 )根据瞬时最小势能原理,系统在外力作用下,在满足位移边界条件和协调条件的所有位移场中, 真实位移场必满足平衡条件并使总势能取极值,即有变分式(511=0而单元总势能-二为(6-5)式中 ---------- 弹性体应变能;——外载荷势能。
由弹性力学的基本理论,有S 1 fff 小刖Q = 4/ (川 s j DBS )免(6-6)3如皿(6-7)式中的第一项为体力项,第二项为外激励力项。
式中」一 ------ 体力向量;人 ----- 外激力向量。
而且丄丄心「丁( 6-8)式中L ---- 单元内力向量;〔:——阻尼系数;L ,------ 体积密度;-1' -1' —— 依次为速度向量与加速度向量。
将(6-1)式代入,考虑到%仅仅是时间一的函数,可得3t D分别为几何矩阵和弹性矩阵。
(6-4)匕=-衬(出矿加亠jff N T CN ^-n n川矿吨3)-训M 晟(6-6 )式代入式(6-5)得单元总势能的表达式:C e = [jj CiV T NdQ17M , = \\I PN TN ^Q并依次称食」「•上丁约为单元的刚度矩阵、质量矩阵、阻尼矩阵、结点力向量。