梁单元公式
abaqus第四讲：应用梁单元解析
（a）梁截面无偏置
（b）梁截面有偏置
北京怡格明思工程技术有限公司
Innovating through simulation
如图所示的工字型梁附着在一个1.2单位厚的壳上。通过定义梁的节点从I截面的底部的偏移量，梁截面的定位可以如图所示。在这种情况下，偏移 量为0.6，亦即壳厚度的一半。
壳截面厚度1.2
梁单元曲率
梁单元的曲率是基于梁的n2方向相对于梁轴的取向。如果n2方向不与 梁轴正交（即，梁轴的方向不与切向t一致），则认为梁单元有初始弯曲。 要模拟曲梁结构，可能需要使用两种方法直接定义n2方向，它允 许你更好地控制对曲率进行模拟： 一种是给出n2矢量的分量作为节点坐标的第4、第5和第6个数据 值； 另一种是使用*NORMAL选项直接地指定法线方向（添加该选项可 以通过ABAQUS/CAE中的Keywords Editor（关键词编辑器））
图 工字型梁用作壳单元的加强件
北京怡格明思工程技术有限公司
Innovating through simulation
你也可以指定形心和剪切中心的位置；这些位置也可以从梁的节点偏置， 从而使你很容易地模拟加强件。 另外也可以分别定义梁节点和壳节点，并在两个节点之间采用一个刚 性梁的约束连接梁和壳。
北京怡格明思工程技术有限公司
Innovating through simulation
二、计算公式和积分
在ABAQUS中的所有梁单元都是梁柱类单元，这意味着它们可以产 生轴向、弯曲和扭转变形。Timoshenko梁单元还考虑了横向剪切变形 的影响。
北京怡格明思工程技术有限公司
Innovating through simulation
剪切变形 ：
线性单元（B21和B31）和二次单元（B22和B32）是考虑剪切变形的 Timoshenko梁单元；因此，它们既适用于模拟剪切变形起重要作用的深梁 又适用于模拟剪切变形不太重要的细长梁。 ABAQUS假设这些梁单元的横 向剪切刚度为线弹性和常数 。 三次单元，称为Euler-Bernoulli梁单元（B23和B33），它们不能模拟 剪切变形。这些单元的横截面在变形过程中与梁的轴线保持垂直 ，因此， 应用三次梁单元模拟相对细长构件的结构更为有效。 对于静态分析，常常可用一个三次单元模拟一个结构构件，而对于动态 分析，也只采用很少数量的单元。这些单元假设剪切变形是可以忽略的。
abaqus第四讲：应用梁单元共20页文档
如图所示的工字型梁附着在一个1.2单位厚的壳上。通过定义梁的节点从I截面的底部的偏移量，梁截面的定位可以如图所示。在这种情况下，偏移 量为0.6，亦即壳厚度的一半。
壳截面厚度1.2 图 工字型梁用作壳单元的加强件
北京怡格明思工程技术有限公司
Innovating through simulation
北京怡格明思工程技术有限公司
Innovating through simulation
扭转响应翘曲
结构构件经常承受扭矩，几乎所有的三维框架结构都会发生这种情况。 在一个构件中引起弯曲的载荷，可能在另一个构件中引起了扭转，如图所 示
北京怡格明思工程技术有限公司
Innovating through simulation
剪切变形 ：
线性单元（B21和B31）和二次单元（B22和B32）是考虑剪切变形的 Timoshenko梁单元；因此，它们既适用于模拟剪切变形起重要作用的深梁 又适用于模拟剪切变形不太重要的细长梁。 ABAQUS假设这些梁单元的横 向剪切刚度为线弹性和常数 。
北京怡格明思工程技术有限公司
Innovating through simulation
ABAQUS中常用的横截面形状：
当梁的轮廓与梁的截面特性相关时，可以指定或者是在分析过程中计算 截面的工程性质，或者是让ABAQUS预先计算它们（在分析开始时）。
选择前者可以应用于线性或者非线性的材料行为（例如，截面刚度因非
另一种是使用*NORMAL选项直接地指定法线方向（添加该选项可 以通过ABAQUS/CAE中的Keywords Editor（关键词编辑器））
北京怡格明思工程技术有限公司
Innovating through simulation
梁单元的几何刚度
梁单元的几何刚度全文共四篇示例，供读者参考第一篇示例：梁单元是有限元分析中常用的一种元素，用于模拟结构中的梁元件。
在有限元分析中，每个梁单元由两个节点、一个横截面和一系列物理性质组成，如材料的弹性模量、截面的面积和惯性矩等。
梁单元的几何刚度是评估结构在受力情况下的扭曲和弯曲变形能力的重要参数之一。
梁单元的几何刚度反映了梁元件在受力情况下的抗弯能力，具有重要的物理意义。
在实际的工程应用中，梁元件的几何刚度可以通过梁单元的有限元模拟来评估，帮助工程师更好地了解结构的受力性能，制定合理的结构设计方案。
在计算梁单元的几何刚度时，需要考虑横截面的形状、尺寸和材料的物理性质等因素。
一般来说，梁单元的几何刚度与截面的几何形状密切相关，例如矩形梁和圆形梁的几何刚度相差较大。
材料的弹性模量、截面的高度和宽度等参数也会影响梁单元的几何刚度。
第二篇示例：梁单元是有限元分析中常用的一个元素，用于模拟实际物体中的横向力和弯曲力。
在有限元分析中，主要包括四个基本力学元素：杆单元、梁单元、壳单元和体单元。
梁单元是用来模拟梁的弯曲变形、传递弯曲载荷和抗弯刚度。
梁单元的几何刚度指的是梁在其几何形状和尺寸的影响下对弯曲应变的抵抗能力，也可以理解为梁在受到外力作用时对弯曲变形的抵抗程度。
梁单元的几何刚度与梁的材料性质、截面形状和尺寸等因素密切相关。
一般来说，梁的几何刚度随着横截面积的增大而增加，随着长度的增大而减小。
这是因为较大的横截面积可以承受更大的弯曲力，而较长的长度则会导致梁在弯曲过程中发生更明显的变形，从而减小梁的抵抗能力。
在设计梁单元时，需要综合考虑这些因素，以确保梁具有足够的几何刚度来承受外部载荷。
在有限元分析中，梁单元的几何刚度通常通过弯曲刚度矩阵来描述。
弯曲刚度矩阵包括四个弯曲刚度分量，分别表示梁在x、y和z方向上的弯曲刚度以及横截面的剪切刚度。
这些弯曲刚度分量可以通过梁单元的几何形状和尺寸来计算，从而得到梁单元的整体几何刚度矩阵。
汽车结构有限元分析03单元类型及单元分析
1.一维单元分析
主要有：杆单元、梁单元、管单元等 。
1.1杆单元---最简单的两节点一维单元， 用于杆件承受轴向力分析。
设杆单元横截面积为Ａ，长度为l，轴 向分布载荷q为(x) 。单元2个节点的位移 向量为： e ui u j T
由空间弹性力学几何方程，得应变表达式：
{} [B]{ }e [[B1 ][B2 ][B20 ]]{ }e
由空间弹性力学物理方程，单元内的应力可 以表示成：
[ ] [D][ ] [D][B]{ }e [S]{ }e
单元刚度矩阵为 ：
[k]e
[B]T [D][B]dV
[k1e1
[k
e 21
这其中设定单元位移模式，利用虚功 原理建立单元节点力与节点位移关系并组建 单元刚度矩阵的过程，我们将其称为单元分 析。
为了使有限元法的解在单元尺寸逐步趋 小时能够收敛于精确解，所构造的单元位移 函数必须满足以下三方面的条件：
1）位移模式中必须包括反映刚体位移的项；
2）位移模式中必须包括反映常应变的线性位 移项；
这样空间梁单元就由3个节点组成i,，j,k 点必
须在一个平面内，但不能共线。i节点到j节
点为单元坐标系的x轴，y轴(或z轴)在节点i、
j和k构成的平面上且与x轴垂直，应用右手定
则可以确定另一坐标iz, 轴j, k(或y轴)。
三点
确定后，单元坐标系即确定，梁单元的截面
方位也就完全确定下来。所增加的一个用于
] ]
[k1e2 ]
[k
e 22
]
[k1e20
[k
第5章 杆单元和梁单元
1 u2 E (2) A(2) (2) 2 u3 l
1 1 u2 1 1 1 u 2 R2 3
u1 在这里，把表达成整体位移矢量 u 2 的函数，如下： u 3
5.1 杆件系统的有限元分析方法
(1) (1) (1)
F3 10N
，进行相应的单元应力计算。得到的结果如下：
0 u1 4 u2 2.5 10 m u 7.5 10 4 m 3
(2) ( x) 5 103 (1) 0.05MPa (2) = 0.1MPa
第五章 杆单元和梁单元
第5章 杆单元和梁单元
本章主要介绍利用杆单元及梁单元进行结构静力学的有限 元分析原理。首先介绍了杆单元的分析方法，详细给出了采用 杆单元进行有限元分析的整个过程；紧接着介绍了平面梁单元 ，以一个平面悬臂梁力学模型为分析实例，分别采用材料力学 、弹性力学解析计算以及有限元法进行了分析与求解，以加深 读者对有限元法的理解。
E (2) A(2) (2) u2 1 u2 l 0 F3 (2) (2) E A u3 2 u3 l (2)
5.1.1 一维杆单元
u2 由最小势能原理，势能函数对未知位移 求变分，满足 u3 的条件是 ，得如下方程式 0, 0
P 1 , u1
E e , Ae , l e
1
图 5-2 杆单元
P2 , u2
2
对于两个节点的杆单元，存在如下节点力和节点位移的关 系式 u P 1 e 1 (5.1) k
P2
u2
其中， k e 称为单元刚度矩阵
5.1.1 一维杆单元
第十七章LS-DYNA的隐式求解.
第十七章 LS-DYNA 的隐式求解LS-DYNA作为著名的显式求解程序只能求解瞬态动力问题，对于时间历程较长的静、动力问题， LS-DYNA的显式中心差分法有它的局限性，而一些与瞬态动力分析紧密相关的问题要求LS-DYNA也能够求解，如：冲压成型后的回弹计算应力初始化冲击后常时间低频动力响应静力分析特征值分析实际上从950版本开始，LS-DYNA已增加了隐式求解功能。
刚开始的应用主要在冲压成型后的回弹计算，经过960版到970版的发展，LS-DYNA的隐式求解功能已大增强，已经能满足以上的求解需要。
17.1显式与隐式的区别：17.1.1 LS-DYNA显式求解：采用中心差分方法进行显式时间积分Man fn方程非耦合，可以直接求解(显式)但需要常小的时间步保持稳定状态不需要求解刚度矩阵 ext fnint适合冲击、穿透等高频非线性动力响应问题17.1.2 LS-DYNA隐式求解：采用Newmark隐式时间积分M an K un1 f对于线性问题，无条件稳定可以采用大的时间步 extn1fintn Man对于非线性问题，需要一系列线性逼近（Newton-Raphson）叠代求解需要求解刚度矩阵适合静力问题、低频动力问题及特征值分析。
17.2 LS-DYNA中隐式分析的激活及相关关键字在LS-DYNA中，缺省的求解是显式求解，可以通过下面的关键字来激活隐式求解：*CONTROL_IMPLICIT_GENERAL*CONTROL_IMPLICIT_GENERAL$ imflag dt0 iefs nstepsb igso1 0.01 0 0 0其中参数imflag=1激活全隐式求解imflag=0为缺省的显式求解。
imflag=2为显式求解后无缝进行隐式回弹求解。
该关键字对于所有隐式求解分析来说都是必需的。
与隐式求解相关的其它关键字：*CONTROL_IMPLICIT_LINEAR(v960版本改为*CONTROL_IMPLICIT_SOLVER)*CONTROL_IMPLICIT_NONLINEAR(在v960版后改为*CONTROL_IMPLICIT_SOLUTION)*CONTROL_IMPLICIT_AUTO*CONTROL_IMPLICIT_STABILIZATION*CONTROL_IMPLICIT_DYNAMICS*CONTROL_IMPLICIT_EIGENVALUE*CONTROL_IMPLICIT_BUCKLE（v970）*CONTROL_IMPLICIT_MODES（v970）使用*CONTROL_IMPLICIT_GENERAL激活隐式求解后，还需要与上面的相关关键字进行联合，从而对不同的问题设置不同的关键字进行求解。
梁单元分布弯矩
梁单元分布弯矩介绍梁单元分布弯矩是结构力学和有限元分析中一个重要的概念。
在建筑和工程领域中，梁承受着各种静力和动力荷载，弯矩是描述梁在受力过程中产生弯曲变形的一种指标。
梁单元分布弯矩的计算和分析对于结构的设计、安全评估和优化具有重要意义。
梁单元的基本概念梁单元是有限元分析中的一种元素，用于建模和分析实际梁结构。
梁单元通常由两个节点和一个连接节点的单元组成。
梁单元的长度可以根据实际情况进行调整，以满足模型的需求。
梁单元通常具有两个自由度，即横向位移和转角，可以通过梁单元力反力关系来求解。
梁单元分布弯矩的定义梁单元分布弯矩是指在梁单元的每一个点上，所受到的弯矩大小和方向。
梁单元分布弯矩可以根据材料力学和结构力学的基本原理进行计算。
在有限元分析中，往往通过施加不同类型的荷载来计算梁单元分布弯矩，例如均布荷载、集中荷载和扭转荷载等。
梁单元分布弯矩的计算方法梁单元分布弯矩的计算方法主要有两种：解析解和数值解。
解析解解析解是指通过数学公式和理论推导得到的精确解。
对于简单的梁结构，可以通过解析解的方法来计算梁单元分布弯矩。
解析解的优点是精确度高，可以直接得到梁单元各个点的分布弯矩值。
但是，解析解的适用范围有限，对于复杂结构和非线性问题，解析解往往难以得到。
数值解数值解是通过有限元分析等数值方法计算得到的近似解。
在数值解的计算中，梁结构被离散为许多小单元，然后通过数值方法进行求解。
数值解的优点是适用范围广，可以处理各种形状和边界条件的梁结构。
在计算梁单元分布弯矩时，数值解可以通过迭代计算每个节点上的弯矩值，得到整个梁单元的分布弯矩图。
梁单元分布弯矩的影响因素梁单元的分布弯矩受到多种因素的影响，主要包括以下几点：荷载类型和大小梁单元的分布弯矩与荷载类型和大小密切相关。
不同类型的荷载会对梁单元产生不同的弯矩分布。
例如，在均布荷载下，梁单元的弯矩分布呈现抛物线形状；而在集中荷载下，梁单元的弯矩分布呈现尖峰状。
梁单元的几何形状和长度梁单元的几何形状和长度对分布弯矩具有重要影响。
梁的有限元分析原理 - 考虑剪切变形影响的梁单元
代人
比较：弯曲梁 单元中的单刚
得到：
等截面梁单元有限元分析
8
长沙理工大学
小结
剪切变形的影响通过系数b反映在刚度矩阵中，使刚度减弱。 对矩形截面：
，当l >>h，b趋于0，可以忽略剪力变形的影响。
等截面梁单元有限元分析
9
长沙理工大学
Timoshenko梁单元
铁木辛柯梁单元——采用两个独立变量 挠度 w
几何关系，曲率
对比
等截面梁单元有限元分析
3
最小势能原理
长沙理工大学
k为截面剪切校正因子
1.经典梁单元 2.铁木辛柯梁单元
——C1型单元 ——C0型单元
等截面梁单元有限元分析
4
长沙理工大学
在经典梁单元基础上引入剪切变形的影响. 挠度叠加
结点位移
其中
采用不考虑剪切变形梁单元的w相同的Hermite插值； 采用2结点的Lagrange插值，即线性插值。
解决方法
假设剪切应变
代替插值函数
计算泛函的剪切应变能时，θ采用低一 阶，和dw/dx同阶插值函数代替原插值 函数
18
等截面梁单元有限元分析
长沙理工大学
等截面梁单元有限元分析
——考虑剪切变形的梁单元
2014.4.13
1
长沙理工大学
介绍.
轴力构件 axial elements 杆单元
受弯构件 flexural elements 梁单元
考虑剪切变形的梁单元
等截面梁单元有限元分析
2
长沙理工大学 假设：梁内的横向剪切力Q所产生的剪切变形将引起梁的附加挠度， 并使原来垂直于中面的截面变形后不再和中面垂直，而且发生翘曲。 考虑剪切变形的梁单元 但在这里，假设原来垂直于中面的截面变形后仍保持为平面。 几何描述