BEAM188

BEAM188 — 3-D 线性有限应变梁

（基于Ansys 5.61的help）

MP ME ST PR PP ED

元素描述

BEAM188 适用于分析细长的梁。 元素是基于Timoshenko 梁理论的。 具有扭切变形效果。

BEAM188 是一个二节点的三维线性梁。 BEAM188 在每个节点上有6或7个自由度，（自由度）数目的变化是由KEYOPT(1)来控制的。当 KEYOPT(1) = 0时 (默认), 每节点有6个自由度。 分别是沿x,y,z的位移及绕其的转动。当 KEYOPT(1) = 1时,会添加第七个自由度 (翘曲量) 。

此元素能很好的应用于线性（分析），大偏转，大应力的非线性（分析）。

BEAM188包含应力刚度，在默认情况下，在某些分析中由 NLGEOM来打开。 在进行弯曲（ flexural）,侧向弯曲（ lateral）, 和扭转稳定性（ torsional stability）分析时，应力刚度应该是被打开的。

BEAM188 能够采用SECTYPE, SECDATA, SECOFFSET, SECWRITE,和 SECREAD来定义任何截面（形状）。. 弹性（elasticity）,蠕变（ creep）,和塑性（ plasticity） 模型都是允许的 (不考虑次截面形状)。

图1. BEAM188 3-D 线性有限应变梁

输入数据

（元素的）几何形状，节点为止，即元素坐标系图示于 BEAM188。BEAM188在模型坐标系中是由节点 I 和节点 J 来定义的。节点 K 是必需的元素方向点定义。有关方向点的相关信息详见 Generating a Beam Mesh With Orientation Nodes 在 ANSYS Modeling and Meshing Guide中。于 LMESH 和 LATT命令说明中可见节点 K 的自动定义的详细说明。

在空间中这是一个没有量纲的元素。截面形状是用SECTYPE 和 SECDATA 命令 (详见 ANSYS Commands Reference )来独立定宓摹Ｃ恳桓鼋孛嫘巫淳?囟ㄒ桓?ID 号(SECNUM)。 截面号是特定的元素属性。

梁元素是基于 Timoshenko 梁理论的，这是一个一阶切应变理论：横向切应变在截面中是常量；也就是说截面在变形后仍是平面。 BEAM188是一阶 Timoshenko 梁元素，它用一个点在长度上来（代替截面）。应此当在节点 I  和 J 上使用SMISC参数的话会显示每个端点节点的形心。  BEAM188 能被用于细长（slender）或粗壮（ stout？？？）的梁。因为一阶切应变理论的限制，自有适当厚度的梁能被分析。梁结构上的细长比 (GAL2/(EI)) 能够用来判断是否采用此元素：

G		切变模数
A		截面面积
L		构件长度
EI		弯曲刚度

在整体（偏移）距离而不是单个元素的情况下记录这个比值是重要的。 悬臂梁受向下的负载 提供了悬臂梁在受向下的负载的情况下横向切应变的一个估评。虽然这个结果不能外推到所有的情况， 但可以作为一个指导。 我们推荐细长比应大于30 。

图 2. 悬臂梁受向下的负载

细长比 (GAL2/(EI)>30)	Timoshenko/ Euler-Bernoulli
25	1.120
50	1.060
100	1.030
1000	1.003

元素能提供一个横向剪切力与横向切应变的弹性关系。你可以用实常量来定义横向剪切刚度。

 扭转变形的St. Venant 翘曲决定了一个综合状态，它可以使（材料）在屈服后的切应力变得平均。 ANSYS 不提供对横截面或可能出现塑性屈服的横截面上的扭切分布情况的换算。应此因扭转负载而引起的大的非弹性的变形应当进行讨论，（ansys）也会检查并给处警告。在这种情况下推荐用实体或壳模型来代替。

在默认情况下BEAM188 元素假设横截面上的弯曲很小可以被忽略(KEYOPT(1) = 0)。 你可以使用KEYOPT(1) = 1来打开弯曲度的自由度。 如果此自由度被打开那每个节点会有7个自由度： UX, UY, UZ, ROTX, ROTY, ROTZ, 和 WARP。

BEAM188 允许用一个轴向延伸率的函数来改变横截面的转动惯量。 默认情况下元素横截面的面积可以改变，但元素的体积在变形前后是相同的。此默认同样适用于elasto-plastic 情况。 使用 KEYOPT(2), 你能使横截面面积为一个常量或保持不变。

元素的输出在元素的积分位置和横截面的积分点上都是有效的。

梁在长度方向的积分点（高斯点）如（图） 积分位置  所示。

Figure 3. 3-D 线性有限应变梁元素的积分位置

截面的应力与力（包含弯矩）都是在积分点上获得的。 元素基本点的输出会外推到元素的节点。

BEAM188 的一些剖面关联量(面积的积分，位置，泊松比函数，函数的导数等等) ，在使用 SECTYPE 和 SECDATA命令定义截面时会自动分配到一个序列号。 每一个截面区域预定为由9个节点组成。 相交区域模型 举例说明矩形块和槽形块的模型情况。 每个单元有4个积分点。

Figure 4. BEAM188 相交区域模型

BEAM188 提供剖面积分点和节点的结果输出。但你仅能查看边界上的输出。 (PRSSOL 打印 BEAM188 剖面节点和积分点的解。应力和应变是在节点上的，塑性应力，塑性功，潜变应变则是在积分点上。)

当元素的材料具有非线性状态或有通过剖面的温度时，计算是在积分点上进行的。在大量通用弹性应用中]，元素采用剖面积分点的pre-calculated 特性。应此，应力与应变的输出均是经过了积分点的计算的。

如果截面分配了次截面 ASEC, 那么只有一般性的应力与应变 (轴向力, 弯矩, 切向应变, （弯曲）曲率, 和切应力)能够输出。 3-D 的轮廓图和变形显示图是不可用的。 ASEC 次截面只能被作为一个薄矩形块来显示验证梁的方向。

质量矩阵与负载向量的相容性的评估，相对于使用的刚度矩阵来说是一个高阶积分。元素提供包含相容性与集中的质量初矩阵。 使 LUMPM,ON 可以让质量矩阵（质量）集中。 （系统）默认使用相容性矩阵。 单位长度的质量可以用 ADDMAS 作为实常量来输入。详见 输入概述。

力是相加在节点上的(定义在元素主方向)。 如果形心轴不与元素主方向重合，那么附加的轴向力会引起弯曲。（同样）如果形心和扭转中心不重合的话，扭转力也会引起扭转变形和扭矩。 应次节点的定位应当与力的中心向重合。使用 SECOFFSET 命令可以适当的改变OFFSETY和 OFFSETZ  的幅角。默认情况下 ANSYS 用形心来定义元素的主轴。

在 节点和元素负载中有元素负载的描述。压力是作为一种面负载来作用在元素表面上的，（元素的面可见）图BEAM188.中带圈文字的显示。正向压力一压力（常规形式）输入。侧向压力以单位长度上的力来输入。 尾端压力以力（的形式）输入。

BEAM188 与ansys中的其他基于埃尔米特多项式（Hermitian polynomial）的元素(f比如说BEAM4).不同，它是基于线性多项式（linear polynomials）的。因此分布式（周延式）负载的偏移在说明中是不允许的。此外不支持非节点上的集中力。（必须加的话）推荐用加细元素的方法。 BEAM188 计算的准确性与收敛性与元素的细化程度相关。

温度作为一种体负载可以加在每个端点节点的三个方向上。 在端点上，加在元素主方向（x-axis）上的温度是 (T(0,0))， y 方向上为(T(1,0)),  z 方向上为 (T(0,1)).。第一个温度坐标T(0,0) 默认为 TUNIF。绻?龆ㄒ辶说谝桓鑫露龋?敲雌渌?木??衔?谝桓觥?如果仅在节点 I 上输入温度，那节点 J 默认对应于节点 I 。其他的输入形式如果未定以均默认为 TUNIF。

KEYOPT(10) = 1 用于从用户子程序中读入初始应力数据。用户子程序的详细叙述请见 ANSYS Guide to User Programmable Features 。

输入概述中给出了元素输入的一个概括说明。

BEAM188 输入概述

Element Name		BEAM188
Nodes		I, J, K
Degrees of Freedom		UX, UY, UZ, ROTX, ROTY, ROTZ if KEYOPT(1) = 0 UX, UY, UZ, ROTX, ROTY, ROTZ, WARP if KEYOPT(1) = 1
Real Constants		(Blank), TYZ1, TXZ1, ADDMAS
Material Properties		EX, EY, EZ, (PRXY, PRYZ, PRXZ, or NUXY, NUYZ, NUXZ), ALPX, ALPY, ALPZ, DENS, GXY, GYZ, GXZ, DAMP
Surface Loads		Pressure -  　 face 1 (I-J) (-z normal direction), face 2 (I-J) (-y normal direction), face 3 (I-J) (+x tangential direction), face 4 (I) (+x axial direction), face 5 (J) (-x direction). 反方向请用负值。
Body Loads		Temperatures -  　 T(0,0), T(1,0), T(0,1) at each end node
Special Features		塑性，潜变，应力强化，大变形，大应力，初始应力输入（Plasticity, Creep, Stress stiffening, Large deflection, Large strain, Initial stress import）。 提供如下的TB 命令项：BISO, MISO, NLISO, BKIN, MKIN, KINH, CHABOCHE, and CREEP. 详见 ANSYS Theory Reference  。
KEYOPT(1)		0 -  　 默认; 六 DOF, 无弯曲 1 -  　 七 DOF (包括弯曲)
KEYOPT(2)		0 -  　 默认; 截面随轴向变长的函数而变化，但须 NLGEOM,ON 。 1 -  　 截面尺寸假定为常量 (经典梁理论)
KEYOPT(5)		0 -  　 默认; 对称的压载荷刚度 1 -  　 非对称的压载荷刚度 2 -  　 忽略压载荷刚度
Note -		只有在 OUTPR,ESOL 激活时KEYOPT(6) 才能通过KEYOPT(9)激活。当 KEYOPTs 6, 7, 8, 和 9 都被激活时， 元素输出中的应力是一个总的应变。 "总" 意味着同时包含了热应变。如果元素的材料定义是定义的塑性材料，那么可以包含塑性应变和塑性功。 可以在 /POST1  中用PRSSOL显示。
KEYOPT(6)		元素积分点输出控制 0 -  　 默认; 输出截面力，截面应变，弯矩 1 -  　 与 KEYOPT(6) = 0 相同，增加截面面积 2 -  　 与 KEYOPT(6) = 1 相同，增加元素基本方向（x,y,z） 3 -  　 输出截面向单元节点外推的 力/力矩，应变/曲率
KEYOPT(7)		截面积分点输出控制 (在截面压尺寸= ASEC时无效) 0 -  　 默认；没有输出 1 -  　 （输出）最大最小应力/应变 2 -  　 与KEYOPT(7) = 1 相同，增加每个截面节点的应力应变输出
KEYOPT(8)		截面节点输出控制 (在截面压尺寸= ASEC时无效) 0 -  　 默认；没有输出 1 -  　 （输出）最大最小应力/应变 2 -  　 与KEYOPT(8) = 1 相同，增加沿截面外边界的应力应变输出 3 -  　 与KEYOPT(8) = 1 相同，增加每个截面节点的应力应变
KEYOPT(9)		元素节点和截面节点的外推值的输出控制 (在截面压尺寸= ASEC时无效) 0 -  　 默认；没有输出 1 -  　 （输出）最大最小应力/应变 2 -  　 与 KEYOPT(9) = 1 相同，增加沿截面外边界的应力应变输出 3 -  　 与KEYOPT(9) = 1  相同，增加每个截面节点的应力应变
KEYOPT(10)		0 -  　 没有用户子程序来提供初始应力（默认） 1 -  　 用USTRESS 从用户子程序中读取初始应力(详见 ANSYS Guide to User Programmable Features 关于用户子程序的章节).

1. 切向应变刚度

输出数据

单元解的输出有两种：

• 节点位移解包含于节点解中。
• 元素附加解请见 元素输出定义

多数情况下，我们推荐 KEYOPT(8) = 2 和 KEYOPT(9) = 2。 详见 ANSYS Basic Analysis Guide f 。

要在结构静态或瞬态分析中观察 BEAM188 的3-D不变形的形状,可以用 OUTRES,MISC 或 OUTRES,ALL 。 要在屈曲分析中观察3D模型，必须进行模态扩展(Elcalc = YES on MXPAND) 。

元素输出表格中所用的符号：

冒号（：）表示可以由[ETABLE, ESOL]的形式获取。 O 列表是存在于 Jobname.OUT中。  R 列表示存在于结果文件中。

在 O 或 R 列中, Y 表示该项肯定有, 数字则表示注释条件下获得， " -- " 表示不存在此项。

表 1. BEAM188 元素输出描述

Name	Definition	O	R
EL	Element number	Y	Y
NODES	Element connectivity	Y	Y
MAT	Material number	Y	Y
VOLU	Volume	Y	Y
XC, YC, ZC	Location where results are reported	Y
AREA	Area of cross section	1	Y
SF	Section forces	1	Y
SE	Section strains	1	Y
S	Section point stresses	2	Y
E	Section point strains	2	Y

1. 见 KEYOPT(6)的描述
2. 见 KEYOPT(7), KEYOPT(8), KEYOPT(9)的描述
3. 自由当形心用 *GET 获取

使用 ETABLE 和 ESOL 命令的项目和顺序号(KEYOPT(9) = 0) 列出了输出通过 ETABLE 命令采用顺序号获取的方法。在 The General Postprocessor (POST1) 于 ANSYS Basic Analysis Guide 和 The "Item and Sequence Number" Table 中有详细介绍。下面是 使用 ETABLE 和 ESOL 命令的项目和顺序号(KEYOPT(9) = 0)的说明：

Name		定义于 元素输出定义中
Item		ETABLE 命令中的项目名称
E		获取单值的元素解的顺序号
I,J		获取节点解的顺序号
ILn		获取中间位置 n 的解的顺序号

表 2. BEAM188 使用 ETABLE和 ESOL 命令的项目和顺序号

Name	Item	I	J
Axial Force	SMISC	1	14
Bending Moment My	SMISC	2	15
Bending Moment Mz	SMISC	3	16
Torque Mx	SMISC	4	17
Shear Force in XZ Plane	SMISC	5	18
Shear Force in XY Plane	SMISC	6	19
Axial Strain	SMISC	7	20
Curvature Kyy	SMISC	8	21
Curvature Kzz	SMISC	9	22
Torsion curvature Kxx	SMISC	10	23
Transverse Shear Stain (XZ)	SMISC	11	24
Transverse Shear Strain (XY)	SMISC	12	25
Area of Cross Section	SMISC	13	26
Bimoment	SMISC	27	29
Bicurvature	SMISC	28	30

假定与限值

梁长度不能为零。截面划分不能多于250个。 This restriction is applicable to the standard library sections that allow user specification of section model size, and for the MESH section subtype. By default (KEYOPT(1) = 0), the effect of warping restraint is assumed to be negligible. Cross section failure or folding is not accounted for.

It is a common practice in civil engineering to model the frame members of a typical multi-storied structure using a single element for each member. Because of cubic interpolation of lateral displacement, BEAM4 and BEAM44 are well-suited for such an approach. However, if BEAM188 is used in that type of application, be sure to use several elements for each frame member. BEAM188 includes the effects of transverse shear.

This element works best with the full Newton-Raphson solution scheme (that is, the default choice in solution control). For nonlinear problems that are dominated by large rotations, we recommend that you do not use PRED,ON.

Note that only moderately "thick" beams may be analyzed. See the Input Data section for more information.

The treatment of material nonlinearities at the section points is based on uncoupled behavior between the stress components. Axial and the torsional shear stress components behave independently. As a result of this approximation, the equivalent stress output by PRSSOL may at times exceed yield stress and yet no active yielding is reported. Under such circumstance, use OUTPR and set KEYOPT(7) = 2 to see actual section point shear stress values.

Product Restrictions

There are no product-specific restrictions for this element.