ansys中的Beam188单元中文说明

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. 悬臂梁受向下的负载

Timoshenko/

细长比 (GAL2/(EI)>30)

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 section1Y

SF Section forces1Y

SE Section strains1Y

S Section point stresses2Y

E Section point strains2Y

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获取节点解的顺序号

IL n获取中间位置n 的解的顺序号

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

Name Item I J

Axial Force SMISC114

Bending Moment My SMISC215

Bending Moment Mz SMISC316

Torque Mx SMISC417

Shear Force in XZ Plane SMISC518

Shear Force in XY Plane SMISC619

Axial Strain SMISC720

Curvature Kyy SMISC821

Curvature Kzz SMISC922

Torsion curvature Kxx SMISC1023

Transverse Shear Stain (XZ)SMISC1124

Transverse Shear Strain (XY)SMISC1225

Area of Cross Section SMISC1326

Bimoment SMISC2729

Bicurvature SMISC2830

假定与限值

梁长度不能为零。截面划分不能多于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, BEAM4and BEAM44are 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.