《ANSYS-Beam188单元应用》由会员分享,可在线阅读,更多相关《ANSYS-Beam188单元应用(10页珍藏版)》请在金锄头文库上搜索。
1、Beam188/189单元基于Timoshenko梁理论(一阶剪切变形理论:横向剪切应变在横截面上是常数,也就是说,变形后的横截面保持平面不发生扭曲)而开发的,并考虑了剪切变形的影响,适合于分析从细长到中等粗细的梁结构。该单元提供了无约束和有约束的横截面的翘曲选项。Beam188是一种3D线性、二次或三次的2节点梁单元。Beam189是一种3D二次3节点梁单元。每个节点有六个或者七个自由度,包括x、y、z 方向的平动自由度和绕x、y、z 轴的转动自由度,还有一个可选择的翘曲自由度。该单元非常适合线性、大角度转动或大应变非线性问题。beam188的应力刚化选项在任何大挠度分析中都是缺省打开的,从
2、而可以分析弯曲、横向及扭转稳定问题(进行特征值屈曲分析或(采用弧长法或非线性稳定法)破坏研究)。Beam188/beam189单元支持弹性、塑性,蠕变及其他非线性材料模型。这种单元还可以采用多种材料组成的截面。该单元还支持横向剪力和横向剪应变的弹性关系,但不能使用高阶理论证明剪应力的分布变化。下图是单元几何示意图:该单元的几何形状、节点位置、坐标体系和压力方向如图所示,beam188 由整体坐标系的节点i 和j 定义。 对于Beam188梁单元,当采用默认的KEYOPT(3)=0,则采用线性的形函数,沿着长度用了一个积分点,因此,单元求解量沿长度保持不变;当KEYOPT(3)=2,该单元就生成
3、一个内插节点,并采用二次形函数,沿长度用了两个积分点,单元求解量沿长度线性变化;当KEYOPT(3)=3,该单元就生成两个内节点,并采用三次形函数,沿长度用了三个积分点,单元求解量沿长度二次变化;当在下面情况下需要考虑高阶单元内插时,推荐二次和三次选项:1) 变截面的单元;2) 单元内存在非均布荷载(包含梯形荷载)时,三次形函数选项比二次选项提供更好的结果。(对于局部的分布荷载和非节点集中荷载情况,只有三次选项有效);3) 单元可能承受高度不均匀变形时。(比如土木工程结构中的个别框架构件用单个单元模拟时)Beam188单元的二次和三次选项有两个限制:1) 虽然单元采用高阶内插,但是beam18
4、8的初始几何按直线处理;2) 因为内节点是不可影响的,所以在这些节点上不允许有边界(或荷载或初始)条件。由于这些限制,所以如果beam189模型的中间节点作用有边界(或荷载或初始)条件或者中间节点不在单元中点时,需要注意beam188的二次选项和beam189的差异。同样,beam188的三次选项不同于传统三次(Hermitian)梁单元。未变形的状态决定了计算扭转作用的St.Venant 翘曲函数,该翘曲函数用来定义屈服后的剪应变。Ansys在 没有提供选项来重新计算在分析过程中变形状态的横截面扭转剪力分布和横截面可能的部分塑性屈服。因此,由扭转作用引起的大的非弹性变形需要谨慎处理和验证。在
5、这样的情况下,推荐采用solid 或者shell 单元来模拟。Beam188 可以在没有方向节点的情况下被定义。在这种情况下,单元的x 轴方向为i 节点指向j 节点。对于两节点的情况,默认的y 轴方向按平行x-y 平面自动计算。对于单元平行与z 轴的情况(或者斜度在0.01%以内),单元的y 轴的方向平行与整体坐标的y 轴(如图)。用第三个节点的选项,用户可以定义单元的x 轴方向。如果两者都定义了,那么第三节点的选项优先考虑。第三个节点(K),如果采用的话,将和i、j 节点一起定义包含单元x 轴和z 轴的平面(如图)。如果该单元采用大变形分析,需要注意这个第三号节点紧紧在定义初始单元方向的时候
6、有效。Beam188/beam189 提供在积分点和界面节点输出的选项。你可以要求紧紧在截面的外表面输出。(PRSSOL 打印截面节点和截面积分点结果。应力和应变在截面的截面打印,塑性应变,塑性作用,蠕变应力在截面的积分点输出。当与单元相关的材料有非弹性的行为或者当截面的温度在截面中有变化,基本计算在截面的积分点上运行。对于更多的普通的弹性的运用,单元运用预先计算好的单元积分点上的截面属性。无论如何,应力和应变通过截面的积分点输出来计算。如果截面指定为ASEC 亚类,仅仅广义的应力和应变(轴力、弯距、横向剪切、弯曲、剪应力)能够输出。3-D 轮廓线和变形形状不能输出。ASEC 亚类紧紧可以作为
7、细矩形来显示来定义梁的方向。Beam188/beam189 能够对组合梁进行分析,(例如,那些由两种或者两个以上材料复合而成的简单的实体梁)。这些组件被假设为完全固接在一起的。因此,该梁表现为一单一的元件。多材料截面能力仅仅在梁的行为假定(铁木辛哥或者伯努力欧拉梁理论)成立的时候能运用。用其他的话说,支持简单的传统铁木辛哥梁理论的扩展。在这些地方可能应用到:& #8226; 双层金属带& #8226; 带金属加固的梁& #8226; 位于不同材料组成的层上的传感器Beam188/beam189 计算在截面刚度水平上的弯距和扭距的耦合。横向的剪切也作为一个独立的量来计算。这对于分层的组合物和夹层
8、量可能会有很大的影响,如果街头处不平衡。Beam188/189 没有用高阶理论来计算剪切应力的变更贡献,如果这些作用必须考虑的话,就需要运用ANSYS 实体单元。要使beam188/beam189 用于特殊的应用,作试验或者其他的数值分析。在合适验证后使用对于组合截面的约束扭曲的选项对于质量矩阵和一致荷载向量的赋值,比刚度矩阵使用的规则更高阶积分规则被使用到。单元支持一致质量矩阵和集中质量矩阵。用LUMPM,ON 命令来激活集中质量矩阵。一致质量矩阵时默认使用的。每单位长度的附加质量将用ADDMAS 截面控制来输入,参见BEAM188 Input Summary。在节点(这些截面定义了单元的x
9、 轴)上施加力,如果重心轴和单元的x 轴不是共线的,施加的轴力将产生弯距。如果质心和剪切中心不是重合的,施加的剪切力将导致扭转应力和弯曲。因而需要设置节点在那些你需要施加力的位置。可以适当的使用secoffset 命令中的offsety 和offsetz 自变量。默认的,ansys 会使用量单元的质心作为参考轴。单元荷载在Node and Element Loads 被描述。压力可能被作为单元表面力被输入,就像Figure 188.1: BEAM188 Geometry中带圈的数字所示。正的压力指向单元内部。水平压力作为单元长度的力来输入。端部的压力作为力输入。当keyopt(3)=0 的时候
10、(默认),beam188 基于线性多项式,和其他的基于厄密多项式的单元(例如beam44)不同,一般来说要求网格划分要细化。当keyopt(3)=2,ansys 增加了一个中间积分点在内插值图标,有效的使得单元成为基于二次型功能的铁木辛哥梁。这个选项迫切被要求,除非这个单元作为刚体使用,而且你必须维持和一阶shell 单元的兼容性。线性变化的弯距被经且的表现。二次选项和beam189 相似,有如下的不同:& #8226; 不论是否使用二次选项,beam188 单元最初始的几何总是直线。& #8226; 你不能读取中间节点,所以边界条件/荷载不能在那些节点描述。均布荷载是不允许描述偏移的。不支持
11、非节点的集中力。用二次选项(keyopt (3)=2 当单元大和契型截面相关。温度可以作为单元的体力在梁的每个端部节点的三个位置输入,单元的温度在单元的x 轴被输入(T(0,0),和在离开x 轴一个单元长度的y 轴(T(1,0), 和在离开x 轴一个单元长度的z 方向(T(0,1)。第一坐标温度T(0,0) 默认是TUNIF。如果所有的温度在第一次以后是没有指明的,那么它们默认的就为第一次输入的温度。如果所有i 节点的温度均输入了,j 节点的都没有指明,那么j 节点的温度默认的是等于i 节点的温度。对于其他的输入模式,没有指明的温度默认的是TUNIF。你可以对该单元通过istress 和isf
12、ile 命令来定义初始应力状态。要获取更多的信息,可以参考ANSYS Basic Analysis Guide的Initial Stress Loading。可以替换的,你可以设置keyopt(10)=1 来从用户的子程序ustress 来读取出初始应力。关于用户子程序的详细资料,参见ANSYS User Programmable Features 的指南。应力刚化作用在单元中没有自动计算,如果对应力刚化作用需要非对称矩阵,使用nropt,unsym。在BEAM188 Input Summary给出单元的输入总结。BEAM188 Input Summary节点I, J, K (K, 方向点,可
13、选但被要求)自由度UX, UY, UZ, ROTX, ROTY, ROTZ if KEYOPT(1) = 0 UX, UY, UZ, ROTX, ROTY, ROTZ, WARP if KEYOPT(1) = 1 Section Controls截面控制TXZ, TXY, ADDMAS (See SECCONTROLS) (TXZ and TXY default to A*GXZ and A*GXY, respectively, where A = cross-sectional area) TXZ 和TXY 默认分别是AGXZ 和AGXY,这里A 是截面面积Material Properti
14、es材料属性EX, (PRXY or NUXY), ALPX, DENS, GXY, GYZ, GXZ, DAMP表面力压力face 1 (I-J) (-z normal direction),face 2 (I-J) (-y normal direction),face 3 (I-J) (+x tangential direction),face 4 (J) (+x axial direction),face 5 (I) (-x direction).(用负数表示作用方向相反)I 和j 是端节点体力温度T(0,0), T(1,0), T(0,1) at each end node特殊特征Pl
15、asticity 塑性Viscoelasticity 粘弹性Viscoplasticity 粘弹性Creep 蠕变Stress stiffening 应力刚化Large deflection 大挠曲Large strain 大应变Initial stress import 初始应力引入Birth and death (requires KEYOPT(11) = 1) 单元的生死(要求keyopt(11)=1) Automatic selection of element technology 自动选择单元技术。支持下列用TB 命令相关的数据表种类: BISO,MISO, NLISO, BKIN, MKIN, KINH, CHABOCHE, HILL, RATE, CREEP, PRONY,SHIFT, CAST, and USER.Note对于材料模型细节可以参见ANSYS, Inc. Theory Reference 对于更多的关许单元技术选择的信息可以参见Automatic Selection of Element Technologies 和ETCONTROLKEYOPT(1)扭转自由度0 -