ansys中的Beam188单元中文说明.doc
如有侵权,请联系网站删除,仅供学习与交流ansys中的Beam188单元中文说明【精品文档】第 8 页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包含应力刚度,在默认情况下,在某些分析中由 来打开。 在进行弯曲( flexural),侧向弯曲( lateral), 和扭转稳定性( torsional stability)分析时,应力刚度应该是被打开的。 BEAM188 能够采用, , , ,和 来定义任何截面(形状)。. 弹性(elasticity),蠕变( creep),和塑性( plasticity) 模型都是允许的 (不考虑次截面形状)。图1. BEAM188 3-D 线性有限应变梁输入数据(元素的)几何形状,节点为止,即元素坐标系图示于 。BEAM188在模型坐标系中是由节点 I 和节点 J 来定义的。节点 K 是必需的元素方向点定义。 有关方向点的相关信息详见 在 中。于 和 命令说明中可见节点 K 的自动定义的详细说明。在空间中这是一个没有量纲的元素。截面形状是用 和 命令 (详见 )来独立定宓摹恳桓鼋孛嫘巫淳囟桓?ID 号()。 截面号是特定的元素属性。梁元素是基于 Timoshenko 梁理论的,这是一个一阶切应变理论:横向切应变在截面中是常量;也就是说截面在变形后仍是平面。 BEAM188是一阶 Timoshenko 梁元素,它用一个点在长度上来(代替截面)。 应此当在节点 I 和 J 上使用SMISC参数的话会显示每个端点节点的形心。 BEAM188 能被用于细长(slender)或粗壮( stout?)的梁。因为一阶切应变理论的限制,自有适当厚度的梁能被分析。 梁结构上的细长比 (GAL2/(EI) 能够用来判断是否采用此元素:G切变模数A截面面积L构件长度EI弯曲刚度在整体(偏移)距离而不是单个元素的情况下记录这个比值是重要的。 提供了悬臂梁在受向下的负载的情况下横向切应变的一个估评。 虽然这个结果不能外推到所有的情况, 但可以作为一个指导。 我们推荐细长比应大于30 。图 2. 悬臂梁受向下的负载细长比 (GAL2/(EI)>30)Timoshenko/ Euler-Bernoulli251.120501.0601001.03010001.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 的一些剖面关联量(面积的积分,位置,泊松比函数,函数的导数等等) ,在使用 和 命令定义截面时会自动分配到一个序列号。 每一个截面区域预定为由9个节点组成。 举例说明矩形块和槽形块的模型情况。 每个单元有4个积分点。Figure 4. BEAM188 相交区域模型BEAM188 提供剖面积分点和节点的结果输出。但你仅能查看边界上的输出。 ( 打印 BEAM188 剖面节点和积分点的解。 应力和应变是在节点上的,塑性应力,塑性功,潜变应变则是在积分点上。)当元素的材料具有非线性状态或有通过剖面的温度时,计算是在积分点上进行的。在大量通用弹性应用中,元素采用剖面积分点的pre-calculated 特性。应此,应力与应变的输出均是经过了积分点的计算的。如果截面分配了次截面 ASEC, 那么只有一般性的应力与应变 (轴向力, 弯矩, 切向应变, (弯曲)曲率, 和切应力)能够输出。 3-D 的轮廓图和变形显示图是不可用的。 ASEC 次截面只能被作为一个薄矩形块来显示验证梁的方向。质量矩阵与负载向量的相容性的评估,相对于使用的刚度矩阵来说是一个高阶积分。元素提供包含相容性与集中的质量初矩阵。 使 ,ON 可以让质量矩阵(质量)集中。 (系统)默认使用相容性矩阵。 单位长度的质量可以用 ADDMAS 作为实常量来输入。详见 。力是相加在节点上的(定义在元素主方向)。 如果形心轴不与元素主方向重合,那么附加的轴向力会引起弯曲。 (同样)如果形心和扭转中心不重合的话,扭转力也会引起扭转变形和扭矩。 应次节点的定位应当与力的中心向重合。使用 命令可以适当的改变OFFSETY和 OFFSETZ 的幅角。默认情况下 ANSYS 用形心来定义元素的主轴。在 中有元素负载的描述。压力是作为一种面负载来作用在元素表面上的,(元素的面可见)图.中带圈文字的显示。正向压力一压力(常规形式)输入。侧向压力以单位长度上的力来输入。 尾端压力以力(的形式)输入。BEAM188 与ansys中的其他基于埃尔米特多项式(Hermitian polynomial)的元素(f比如说).不同,它是基于线性多项式(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 NameBEAM188NodesI, J, KDegrees of FreedomUX, UY, UZ, ROTX, ROTY, ROTZ if KEYOPT(1) = 0UX, UY, UZ, ROTX, ROTY, ROTZ, WARP if KEYOPT(1) = 1Real Constants(Blank), TYZ1, TXZ1, ADDMASMaterial PropertiesEX, EY, EZ, (PRXY, PRYZ, PRXZ, or NUXY, NUYZ, NUXZ), ALPX, ALPY, ALPZ, DENS, GXY, GYZ, GXZ, DAMPSurface LoadsPressure - 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 LoadsTemperatures - T(0,0), T(1,0), T(0,1) at each end nodeSpecial Features塑性,潜变,应力强化,大变形,大应力,初始应力输入(Plasticity, Creep, Stress stiffening, Large deflection, Large strain, Initial stress import)。 提供如下的 命令项:BISO, MISO, NLISO, BKIN, MKIN, KINH, CHABOCHE, and CREEP. 详见 ANSYS Theory Reference 。KEYOPT(1)0 - 默认; 六 DOF, 无弯曲1 - 七 DOF (包括弯曲)KEYOPT(2)0 - 默认; 截面随轴向变长的函数而变化,但须 ,ON 。1 - 截面尺寸假定为常量 (经典梁理论)KEYOPT(5)0 - 默认; 对称的压载荷刚度1 - 非对称的压载荷刚度2 - 忽略压载荷刚度Note - 只有在 ,ESOL 激活时KEYOPT(6) 才能通过KEYOPT(9)激活。当 KEYOPTs 6, 7, 8, 和 9 都被激活时, 元素输出中的应力是一个总的应变。 "总" 意味着同时包含了热应变。如果元素的材料定义是定义的塑性材料,那么可以包含塑性应变和塑性功。 可以在 /POST1 中用显示。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。 详见 f 。要在结构静态或瞬态分析中观察 BEAM188 的3-D不变形的形状,可以用 ,MISC 或 ,ALL 。 要在屈曲分析中观察3D模型,必须进行模态扩展(Elcalc = YES on ) 。元素输出表格中所用的符号:冒号(:)表示可以由, 的形式获取。 O 列表是存在于 Jobname.OUT中。 R 列表示存在于结果文件中。在 O 或 R 列中, Y 表示该项肯定有, 数字则表示注释条件下获得, " - " 表示不存在此项。表 1. BEAM188 元素输出描述NameDefinitionORELElement numberYYNODESElement connectivityYYMATMaterial numberYYVOLUVolumeYYXC, YC, ZCLocation where results are reportedY AREAArea of cross sectionYSFSection forcesYSESection strainsYSSection point stressesYESection point strainsY1. 见 KEYOPT(6)的描述 2. 见 KEYOPT(7), KEYOPT(8), KEYOPT(9)的描述 3. 自由当形心用 获取 列出了输出通过 命令采用顺序号获取的方法。在 于 和 中有详细介绍。下面是 的说明:Name定义于 中Item 命令中的项目名称E获取单值的元素解的顺序号I,J获取节点解的顺序号ILn获取中间位置 n 的解的顺序号 表 2. BEAM188 使用 和 命令的项目和顺序号NameItemIJAxial ForceSMISC114Bending Moment MySMISC215Bending Moment MzSMISC316Torque MxSMISC417Shear Force in XZ PlaneSMISC518Shear Force in XY PlaneSMISC619Axial StrainSMISC720Curvature KyySMISC821Curvature KzzSMISC922Torsion curvature KxxSMISC1023Transverse Shear Stain (XZ)SMISC1124Transverse Shear Strain (XY)SMISC1225Area of Cross SectionSMISC1326BimomentSMISC2729BicurvatureSMISC2830假定与限值梁长度不能为零。截面划分不能多于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, and 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 ,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 may at times exceed yield stress and yet no active yielding is reported. Under such circumstance, use and set KEYOPT(7) = 2 to see actual section point shear stress values.Product RestrictionsThere are no product-specific restrictions for this element.