Beam188汇总.docx
《Beam188汇总.docx》由会员分享,可在线阅读,更多相关《Beam188汇总.docx(15页珍藏版)》请在冰豆网上搜索。
Beam188汇总
Beam188
3维线性有限应变梁单元
Beam188单元描述
Beam188单元适合于分析从细长到中等粗短的梁结构,该单元基于铁木辛哥梁结构理论,并考虑了剪切变形的影响。
Beam188是三维线性(2节点)或者二次梁单元。
每个节点有六个或者七个自由度,自由度的个数取决于KEYOPT
(1)的值。
当KEYOPT
(1)=0(缺省)时,每个节点有六个自由度;节点坐标系的x、y、z方向的平动和绕x、y、z轴的转动。
当KEYOPT
(1)=1时,每个节点有七个自由度,这时引入了第七个自由度(横截面的翘曲)。
这个单元非常适合线性、大角度转动和/并非线性大应变问题。
当NLGEOM打开的时候,beam188的应力刚化,在任何分析中都是缺省项。
应力强化选项使本单元能分析弯曲、横向及扭转稳定问题(用弧长法)分析特征值屈曲和塌陷)。
Beam188/beam189可以采用sectype、secdata、secoffset、secwrite及secread定义横截面。
本单元支持弹性、蠕变及素性模型(不考虑横截面子模型)。
这种单元类型的截面可以是不同材料组成的组和截面。
Beam188从6.0版本开始忽略任何实参数,参考seccontrols命令来定义横向剪切刚度和附加质量。
单元坐标系统(/psymb,esys)与beam188单元无关。
下图是单元几何示意图:
BEAM188输入数据
该单元的几何形状、节点位置、坐标体系如图“BEAMGeometry”所示,beam188由整体坐标系的节点i和j定义。
节点K是定义单元方向的所选方式,有关方向节点和梁的网格划分的信息可以参见ANSYSModelingandMeshingGuide中的GeneratingaBeamMeshWithOrientationNodes。
参考lmesh和latt命令描述可以得到k节点自动生成的详细资料。
Beam188可以在没有方向节点的情况下被定义。
在这种情况下,单元的x轴方向为i节点指向j节点。
对于两节点的情况,默认的y轴方向按平行x-y平面自动计算。
对于单元平行与z轴的情况(或者斜度在0.01%以内),单元的y轴的方向平行与整体坐标的y轴(如图)。
用第三个节点的选项,用户可以定义单元的x轴方向。
如果两者都定义了,那么第三节点的选项优先考虑。
第三个节点(K),如果采用的话,将和i、j节点一起定义包含单元x轴和z轴的平面(如图)。
如果该单元采用大变形分析,需要注意这个第三号节点紧紧在定义初始单元方向的时候有效。
梁单元是一维空间线单元。
横截面资料用sectype和secdata命令独立的提供,参见ANSYSStructuralAnalysisGuide的BeamAnalysisandCrossSections看详细资料。
截面与单元用截面ID号(SECNUM)来关联,截面号是独立的单元属性。
除了等截面,还可以用sectype命令中的锥形选项来定义锥形截面(参考DefiningaTaperedBeam)。
单元基于铁木辛哥梁理论,这个理论是一阶剪切变形理论;横向剪切应力在横截面是不变的,也就是说变形后横截面保持平面不发生扭曲。
Beam188是一阶铁木辛哥梁单元,沿着长度用了一个积分点,用默认的KEYOPT(3)设置。
因此,在i和j节点要求SMISC数值的时候,中间数值在两端节点均输出。
当KEYOPT
(1)设置为2,两个积分点作为延长的线性变量被运用。
Beam188/beam189单元可以用在细长或者短粗的梁。
由于一阶剪切变形的限制,只有适度的“粗”梁可以分析。
梁的长细比(GAL2/(EI))可以用来判定单元的适用性,这里:
G
剪切模量
A
截面积
L
长度
EI
抗弯刚度
需要注意的是这个比例的计算需要用一些全局距离尺寸,不是基于独立的单元尺度。
下面这个图提供了受端部集中荷载的悬臂梁的横向剪切变形的评估,这个例子可以作为一个很好的大致的指导。
我们推荐长细比要大于30。
Figure188.2TransverseShearDeformationEstimation
长细比(GAL2/(EI))
铁木辛哥/欧拉-伯努力( δTimoshenko/δEuler-Bernoulli)
25
1.120
50
1.060
100
1.030
1000
1.003
这些单元支持横向剪切力和横向剪切变应力的弹性关系。
你可以用seccontrols命令忽略默认的横向剪切刚度值。
无形变的状态决定了扭转作用引起的St.Venant翘曲变形,甚至可以用来定义屈服后的剪应力。
Ansys没有提供选项来不成型的结构重新计算,这种结构是由分析过程中的扭转剪切对横截面的作用和可能的截面的部分塑性屈服引起的。
正因为此,由扭转作用引起的非弹性大变形需要小心的来处理和较合。
在这样的情况下,要求用solid或者shell单元来替换。
Beam188/beam189单元支持“约束扭转”分析,通过定义梁节点的第七个自由度来实现。
Beam188单元默认的假设是截面的扭转是足够小的以至于可以忽略(KEYOPT
(1)=0)。
你可以激活它的扭转自由度通过定义KEYOPT
(1)=1。
当激活节点的扭转自由度的时候,每个节点有七个自由度:
UX,UY,UZ,ROTX,ROTY,ROTZ,andWARP。
当KEYOPT
(1)=1,双力矩合双弧线将被输出。
实际上,当两个“约束扭转”的单元以一个尖锐的角度组合在一起的时候,你需要耦合他们的唯一合转角,但是它们平面外的自由度解藕。
通过用两个节点在物理位置和运用合适的约束可以一般地实现。
这个过程很容易的(自动的)实现,通过ENDRELEASE命令,这个命令将两个临近横截面相交角度大于20度的单元的平面外扭转解耦。
Beam188允许改变横截面惯性属性来实现轴向伸长的功能。
默认的,截面面积改变来使得单元的体积变形后不变化。
这种默认的值对于弹塑性应用是适用的。
通过运用KEYOPT
(2),你可以选择使得横截面是恒定的或者刚性的。
单元的输出在单元积分位置和截面的积分点是可以用到的。
沿着梁长度的积分点(高丝点)如图Figure188.3:
"BEAM188ElementIntegrationStations"所示:
Figure 188.3 BEAM188ElementIntegrationStations
截面的应变和力(包括弯距)可以在这些积分点上得到。
单元支持输出选项来外推这些数值到单元的节点。
Beam188/beam189的在一些截面点的截面相关量(积分面积、位置、泊松比、函数倒数等)通过用sectype和secdata命令自动计算得到。
每个截面假定是由一系列预先决定的9节点元组合而成。
下面的图形描述了模型运用矩形截面亚类和槽型截面亚类。
每个截面单元有4个积分点,每个可能与独立的材料种类相关联。
Figure 188.4 Cross-SectionCells
Beam188/beam189提供在积分点和界面节点输出的选项。
你可以要求紧紧在截面的外表面输出。
(PRSSOL打印截面节点和截面积分点结果。
应力和应变在截面的截面打印,塑性应变,塑性作用,蠕变应力在截面的积分点输出。
当与单元相关的材料有非弹性的行为或者当截面的温度在截面中有变化,基本计算在截面的积分点上运行。
对于更多的普通的弹性的运用,单元运用预先计算好的单元积分点上的截面属性。
无论如何,应力和应变通过截面的积分点输出来计算。
如果截面指定为ASEC亚类,仅仅广义的应力和应变(轴力、弯距、横向剪切、弯曲、剪应力)能够输出。
3-D轮廓线和变形形状不能输出。
ASEC亚类紧紧可以作为细矩形来显示来定义梁的方向。
Beam188/beam189能够对组合梁进行分析,(例如,那些由两种或者两个以上材料复合而成的简单的实体梁)。
这些组件被假设为完全固接在一起的。
因此,该梁表现为一单一的元件。
多材料截面能力仅仅在梁的行为假定(铁木辛哥或者伯努力欧拉梁理论)成立的时候能运用。
用其他的话说,支持简单的传统铁木辛哥梁理论的扩展。
在这些地方可能应用到:
•双层金属带
•带金属加固的梁
•位于不同材料组成的层上的传感器
Beam188/beam189计算在截面刚度水平上的弯距和扭距的耦合。
横向的剪切也作为一个独立的量来计算。
这对于分层的组合物和夹层量可能会有很大的影响,如果街头处不平衡。
Beam188/189没有用高阶理论来计算剪切应力的变更贡献,如果这些作用必须考虑的话,就需要运用ANSYS实体单元。
要使beam188/beam189用于特殊的应用,作试验或者其他的数值分析。
在合适验证后使用对于组合截面的约束扭曲的选项
对于质量矩阵和一致荷载向量的赋值,比刚度矩阵使用的规则更高阶积分规则被使用到。
单元支持一致质量矩阵和集中质量矩阵。
用LUMPM,ON命令来激活集中质量矩阵。
一致质量矩阵时默认使用的。
每单位长度的附加质量将用ADDMAS截面控制来输入,参见"BEAM188InputSummary"。
在节点(这些截面定义了单元的x轴)上施加力,如果重心轴和单元的x轴不是共线的,施加的轴力将产生弯距。
如果质心和剪切中心不是重合的,施加的剪切力将导致扭转应力和弯曲。
因而需要设置节点在那些你需要施加力的位置。
可以适当的使用secoffset命令中的offsety和offsetz自变量。
默认的,ansys会使用量单元的质心作为参考轴。
单元荷载在NodeandElementLoads被描述。
压力可能被作为单元表面力被输入,就像Figure188.1:
"BEAM188Geometry"中带圈的数字所示。
正的压力指向单元内部。
水平压力作为单元长度的力来输入。
端部的压力作为力输入。
当keyopt(3)=0的时候(默认),beam188基于线性多项式,和其他的基于厄密多项式的单元(例如beam44)不同,一般来说要求网格划分要细化。
当keyopt(3)=2,ansys增加了一个中间积分点在内插值图标,有效的使得单元成为基于二次型功能的铁木辛哥梁。
这个选项迫切被要求,除非这个单元作为刚体使用,而且你必须维持和一阶shell单元的兼容性。
线性变化的弯距被经且的表现。
二次选项和beam189相似,有如下的不同:
•不论是否使用二次选项,beam188单元最初始的几何总是直线。
•你不能读取中间节点,所以边界条件/荷载不能在那些节点描述。
均布荷载是不允许描述偏移的。
不支持非节点的集中力。
用二次选项(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和isfile命令来定义初始应力状态。
要获取更多的信息,可以参考ANSYSBasicAnalysisGuide的InitialStressLoading。
可以替换的,你可以设置keyopt(10)=1来从用户的子程序ustress来读取出初始应力。
关于用户子程序的详细资料,参见ANSYSUserProgrammableFeatures的指南。
应力刚化作用在单元中没有自动计算,如果对应力刚化作用需要非对称矩阵,使用nropt,unsym。
在"BEAM188InputSummary"给出单元的输入总结。
BEAM188InputSummary
节点
I,J,K(K,方向点,可选但被要求)
自由度
UX,UY,UZ,ROTX,ROTY,ROTZifKEYOPT
(1)=0UX,UY,UZ,ROTX,ROTY,ROTZ,WARPifKEYOPT
(1)=1SectionControls
截面控制
TXZ,TXY,ADDMAS(SeeSECCONTROLS)(TXZandTXYdefaulttoA*GXZandA*GXY,respectively,whereA=cross-sectionalarea)TXZ和TXY默认分别是A×GXZ和A×GXY,这里A是截面面积MaterialProperties
材料属性
EX,(PRXYorNUXY),ALPX,DENS,GXY,GYZ,GXZ,DAMP
表面力
压力
face1(I-J)(-znormaldirection),
face2(I-J)(-ynormaldirection),
face3(I-J)(+xtangentialdirection),
face4(J)(+xaxialdirection),
face5(I)(-xdirection).
(用负数表示作用方向相反)
I和j是端节点
体力
温度
T(0,0),T(1,0),T(0,1)ateachendnode
特殊特征
Plasticity塑性
Viscoelasticity粘弹性
Viscoplasticity粘弹性
Creep蠕变
Stressstiffening应力刚化
Largedeflection大挠曲
Largestrain大应变
Initialstressimport初始应力引入
Birthanddeath(requiresKEYOPT(11)=1)单元的生死(要求keyopt(11)=1)Automaticselectionofelementtechnology自动选择单元技术。
支持下列用TB命令相关的数据表种类:
BISO,MISO,NLISO,BKIN,MKIN,KINH,CHABOCHE,HILL,RATE,CREEP,PRONY,SHIFT,CAST,andUSER.
Note
对于材料模型细节可以参见ANSYS,Inc.TheoryReference对于更多的关许单元技术选择的信息可以参见AutomaticSelectionofElementTechnologies和ETCONTROL
KEYOPT
(1)
扭转自由度
0--
默认;六个自由度,不限制扭转
1--
7个自由度(包括扭转),双力矩和双曲线被输出
KEYOPT
(2)
截面缩放比例
0--
默认;截面因为轴线拉伸效应被缩放;当大变形开关打开的时候被调用。
1--
截面被认为是刚性的(经典梁理论)
KEYOPT(3)
插值数据
0--
默认;线性多项式。
要求划分细致。
2--
二次型(对于铁木辛哥梁单元有效)运用中间节点(中点点用户无法修改)来提高单元的精度,能够精确的表示线性变化的弯距。
KEYOPT(4)
剪应力输出
0--
默认;仅仅输出扭转相关的剪应力
1--
仅仅输出弯曲相关的横向剪应力。
2--
紧紧输出前两种方式的组合状态。
KEYOPT(6)
在单元积分点输出控制
0--
默认;输出截面力、应变、和弯距
1--
和keyopt(6)=0相同,加上当前的截面单元
2--
和keyopt(6)=1相同加上单元基本方向(X、Y、Z)
3--
输出截面力、弯距和应力、曲率,外推到单元节点。
Note
仅仅当outpr,esol是激活状态的时候,Keyopt(6)通过keyopt(9)来激活。
当keyopt(6)、(7)、(8)和(9)都激活的时候,在单元输出中的应变是总应变。
这个“总”包括温度应变。
当单元材料是有塑性的时候,能够提供塑性应变和塑性作业。
在/post1,可替换的运用prssol命令。
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来读取初始应力。
Note
参考GuidetoANSYSUserProgrammableFeatures帮助用户书写子程序。
KEYOPT(11)
设置截面属性
0--
自动计算是否能够提前积分截面属性。
(默认)
1--
用户单元数值积分(在生/死功能的时候要求)
KEYOPT(12)
契型截面处理
0--
线性变化的契型截面分析;截面属性在每个积分点计算(默认),这种方法更加精确,但是计算量大。
1--
平均截面分析;对于契型截面单元,截面属性仅仅在中点计算。
这是划分网格的阶数的估计,但是,速度快。
Beam188的输出数据
这种单元用两种方式计算输出
•节点唯一和反应包括全部节点的计算。
•附加的单元输出在Table188.1:
"BEAM188ElementOutputDefinitions"描述。
在需要的地点,ansys要求keyopt(8)=2和keyopt(9)=2,参考ANSYSBasicAnalysisGuide来找到查看结果的方法。
要看beam188的3-D变形形状,运用OUTRES,MISC或者OUTRES命令,所有的静态和瞬态分析的命令。
要观察模态分析和特征值屈曲分析的3-D模态形状,必须用激活单元结果扩展模态(MXPAND命令Elcalc=YES的选项)
对于梁设计很常规的是使用轴力成分,轴力由轴向荷载和在各个端点的弯曲独立提供。
因此,beam188提供线性的应力输出作为它的SMISC输出命令的一部分,由下面的定义来指示:
SDIR是轴力引起的应力分量。
SDIR=FX/A,这里FX是轴力(SMISC的数值为1和14),A表示截面面积。
SBYT和SBYB是弯曲应力分量。
SBYT=-MZ*ymax/Izz
SBYB=-MZ*ymin/Izz
SBZT=MY*zmax/Iyy
SBZB=MY*zmin/Iyy
这里MY、MZ是弯距(SMISC数值是2、15、3、16)。
坐标ymax,ymin,zmax,和zmin是y和z坐标的最大和最小值。
数值Iyy和Izz是截面惯性距。
对于ASEC梁截面,ANSYS用最大和最小截面尺度,对于ASEC种类的截面,最大最小的Y和Z方向直接分别假定在+0.5到-0.5。
单元应力的相应定义:
EPELDIR=EXEPELBYT=-KZ*ymax
EPELBYB=-KZ*ymin
EPELBZT=KY*zmax
EPELBZB=KY*zmin
这里EX、KY和KZ是总应力和曲率(SMISC数值是7,8,9,20,21和22)
输出的应力仅仅对于单元的弹性行为严格有效。
Beam188总是组合应力来支持非线性材料的行为。
当单元和非线性材料相关的时候,组合应力最好作为线性近似来对待,应该谨慎的说明。
单元运用以下符号输出定义表格:
在name列的冒号表示该项目可以通过构成名字的方法来获得[ETABLE,ESOL]。
第0列表示该项有效的说明在文件Jobname.OUT中。
R列表示该项的结果显示在results文件中。
无论在0还是R列中,Y表示该项一直是可用的。
数值表示描述哪里该项是选择性提供的脚注,-表示该项不提供。
Table 188.1 BEAM188ElementOutputDefinitions
Name
Definition
O
R
EL
Elementnumber
Y
Y
NODES
Elementconnectivity
Y
Y
MAT
Materialnumber
Y
Y
C.G.:
X,Y,Z
Centerofgravity
Y
Y
AREA
Areaofcross-section
1
Y
SF:
Y,Z
Sectionshearforces
1
Y
SE:
Y,Z
Sectionshearstrains
1
Y
S:
XX,XZ,XY
Sectionpointstresses
2
Y
E:
XX,XZ,XY
Sectionpointstrains
2
Y
MX
Torsionalmoment
Y
Y
KX
Torsionalstrain
Y
Y
KY,KZ
Curvature
Y
Y
EX
Axialstrain
Y
Y
FX
Axialforce
Y
Y
MY,MZ
Bendingmoments
Y
Y
BM
Bimoment
3
3
BK
Bicurvature
3
3
Note
MoreoutputisdescribedonthePRSSOLcommandin/POST1
SeeKEYOPT(6)description
SeeKEYOPT(7),KEYOPT(8),KEYOPT(9)descriptions
SeeKEYOPT
(1)description
Table188.2:
"BEAM188ItemandSequenceNumbers"列出了通过etable命令用序列号方法提供的输出。
参见ANSYSBasicAnalysisGuide中的CreatinganElementTable和TheItemandSequenceNumberTable来获取更多的信息。
Table188.2:
"BEAM188ItemandSequenceNumbers"用到了下列符号:
Name
在Table188.1:
"BEAM188