SOLID65单元描述:SOLID65单元用于含钢筋或不含钢筋的三维实体模型。该实体模型可具有拉裂与压碎的性能。在混凝土的应用方面,如用单元的实体性能来模拟混凝土,而用加筋性能来模拟钢筋的作用。当然该单元也可用于其它方面,如加筋复合材料(如
SOLID65单元描述:
SOLID65单元用于含钢筋或不含钢筋的三维实体模型。该实体模型可具有拉裂与压碎的性能。在混凝土的应用方面,如用单元的实体性能来模拟混凝土,而用加筋性能来模拟钢筋的作用。当然该单元也可用于其它方面,如加筋复合材料(如玻璃纤维)及地质材料(如岩石)。该单元具有八个节点,CAE,有限元分析每个节点有三个自由度,即x,y,z三个方向的线位移;还可对三个方向的含筋情况进行定义。
本单元与SOLID45单元(三维结构实体单元)的相似,只是增加了描述开裂与压碎的性能。本单元最重要的方面在于其对材料非线性的处理。其可模拟混凝土的开裂(三个正交方向)、压碎、CAE,有限元分析塑性变形及徐变,还可模拟钢筋的拉伸、压缩、塑性变形及蠕变,但不能模拟钢筋的剪切性能。有关SOLID65单元的更细节的描述请参见《ANSYS理论手册》。
SOLID65的几何模型图
SOLID65输入数据:
关于单元几何图形、节点位置、单元坐标系请见上图。单元性质为八节点各向同性材料,单元包括一种实体材料和三种钢筋材料,用命令MAT输入对混凝土材料的定义,CAE,有限元分析而有关钢筋的细则需在实常数中定义,包括材料号、体积率、方向角(THETA, PHI),钢筋的方向角可通过命令/Eshape以图示方式校验。
体积率是指钢筋的体积与整个单元体积的比,钢筋的方向通过单元坐标系中的两个角度(度制)来定义。当钢筋的材料号为0或等于单元的材料号时则不考虑它的作用。
另外,有关混凝土的材料定义,如剪切传递系数,拉应力,压应力都应在数据表中给出,详细描述见表“SOLID65混凝土材料数据表”。通常剪力传递系数为0~1.0,0表示平滑的裂缝(完全丧失剪力传递作用),1表示粗糙的裂缝(几乎没有失去剪力传递作用)。这就有利于对裂缝开裂与闭合进行描述。
有关单元荷载的描述见“节点单元荷载”(ANSYS帮助中专有一节)。压力作为面荷载作用在单元表面如“SOLID65的几何模型图”中带圈数字所示。主动力作用在单元内。温度和影响(fluences一词不知如何译好)可在节点上作为单元体荷载输入。节点I的温度T(I)默认为TUNIF,如其它节点温度没有被指定,则它们默认为T(I)。对于其它的输入模型未指定温度时默认值都为TUNIF。对影响(fluence)的设定除用0取代TUNIF外与温度的设定是相同的。
用命令TREF和BETAD分别用来设定整体的基准温度和阻尼值。用MAT命令指定与单元相关的基准温度值(MP,REFT)或阻尼值(MP,DAMP),但不能对钢筋的材料号进行以上定义。
KEYOPT(1)用于设定是否考虑大变形,KEYOPT(5)和KEYOPT(6)则提供是多种单元输出选项(详见单元解答)。
KEYOPT(7)是与是否考虑应力松弛相关的项,当KEYOPT=1时表示考虑,目的是加速裂缝即将开裂时计算的收敛(在混凝土材料数据表的第9个系数中的输入值即为拉伸应力松弛的折减系数)。应力松弛并不能反应因为次生裂缝的产生而引起的应力应变关系的变化。松弛系数在裂缝处为零,因此,相应的开裂面上的刚度也是零。
在几何非线性分析时可用“SOLCONTROL,,,INCP”命令设定考虑抗压刚度的影响。抗压刚度的影响在线性屈曲分析中会被自动考虑。
SOLD65单元输入总结:
节点:
I,J,K,L,M,N,O,P
自由度:
UX,UY,UZ
实常数:
MAT1, VR1, THETA1, PHI1, MAT2, VR2,
THETA2, PHI2, MAT3, VR3, THETA3, PHI3
(这里的 MATn 是材料号, VRn 是体积率, 而THETAn 和 PHIn 方向角,代表了三种钢筋材料。)
材料性质:
EX, ALPX, PRXY or NUXY, DENS (用于混凝土)
EX, ALPX, DENS (用于每种钢筋)
还可通过MAT命令设定阻尼与基准温度。
表面荷载:
压力―――
面1(J-I-L-K),面2(I-J-N-M),面3(J-K-O-N),
面4(K-L-P-O),面5(L-I-M-P),面6(M-N-O-P)
体荷载:
温度―――
T(I), T(J), T(K), T(L), T(M), T(N), T(O), T(P)
影响(Fluences)--
FL(I), FL(J), FL(K), FL(L), FL(M), FL(N), FL(O), FL(P)
特性:
Plasticity(塑性)
Creep(蠕变\徐变)
Cracking(开裂)
Crushing(压碎)
Large deflection(大变形)
Large strain(大应变)
Stress stiffening(应力强化)
Birth and death (生死)(单元生死?)
Adaptive descent
KEYOPT(1) 大变形控制:
0 --
考虑大变形
1 --
不考虑
KEYOPT(5) 混凝土线性解的输出控制:
0 --
只打印质心的线性解
1 --
给出每个积分点的解
2 --
输出节点应力
KEYOPT(6) 混凝土非线性解输出控制:
0 --
只打印质心的线性解
3 --
同时还给出积分点的解
KEYOPT(7) 开裂后应力松驰考虑选项:
0 --
不考虑拉伸应力松驰
1 --
考虑应力松驰,有助于计算收敛
SOLID65单元信息:
列在“SOLID65混凝土材料数据表”中的数据项是通过TB命令来建立的。表中未重新赋值的项默认为0或其特定默认值。数据表首先用TB(with LAB=CONC)命令生成,再用命令TBTEMP定义温度,用命令TBDATA定义其余表项。最多可定义六种温度,相应与每种温度都可用TBDATA建立C1~C9各项的值。
SOLID65混凝土材料数据表
常数
含义
1
裂缝张开剪力传递系数.
2
裂缝闭合剪力传递系数
3
单轴抗拉强度
4
单轴抗压强度
5
双轴抗压强度
6
围压大小
7
围压下双轴抗压强度
8
围压下单轴抗压强度
9
拉应力折减系数 当 KEYOPT(7) = 1时有效 (默认为 0.6).
注:当变量3(4)被设为-1时表示混凝土无开裂(无压碎),如只输入前1-4项,则5-8项取默认值,如5-8项中给定了其中一项的值,则其余3项也必需给出。
SOLID65 输出数据:
与单元相关的解答输出项有以下两方面:
l 所有节点的节点位移;
l 其它输出项见“SOLID65单元输出数据说明表”
一些细则的说明可见“SOLID65应力输出图”。单元应力的方向平行于单元坐标系,只有当非线性特性被考虑时才有相应的输出,当然也只有对钢筋参数进行过定义,才有关于它的输出。如可能发生开裂或压碎,那具体情况也会在积分点上输出,因为开裂和压碎可能发生在任一积分点上。在POST1中用命令PLCRACK可显示各积分点的状态。对输出结果更普便的说明可参见“结果输出说明”,对各种结果数据(图形)的查看方法请见《ANSYS入门指南》。
SOLID65应力输出图
单元输出说明表的有关事项解释如下:
表中第一列给出了各输出项的名称,用命令ETABLE(POST1)及ESOL(POST26)可定义这些变量用于查询。第三列表示某一变量值是否在输出文件中给出,第四列某一变量值是否在结果文件中给出。
无论是第三还是第四列,“Y”表示可以输出,列中的具体数值则表示在满足特定条件时才输出,而“-”则表示不输出。
SOLID65单元输出数据说明表
Name
Definition
O
R
EL
单元号
Y
Y
NODES
节点: I, J, K, L, M, N, O, P
Y
Y
MAT
材料号
Y
Y
NREINF
钢筋数量
Y
-
VOLU:
体积
Y
Y
PRES
压力作用面P1由节点 J, I, L, K围成; P2由I, J, N, M围成; P3由 J, K, O, N围成; P4由K, L, P, O围成; P5由L, I, M, P围成; P6由M, N, O, P围成。
Y
Y
TEMP
作用在各节点上的温度 T(I), T(J), T(K), T(L), T(M), T(N), T(O), T(P)
Y
Y
FLUEN
作用在各节点上的影响 FL(I), FL(J), FL(K), FL(L), FL(M), FL(N), FL(O), FL(P)
Y
Y
XC, YC, ZC
Location where results are reported
Y
6
S:X, Y, Z, XY, YZ, XZ
应力
1
1
S:1, 2, 3
主应力
1
1
S:INT
应力强度
1
1
S:EQV
等效应力
1
1
EPEL:X, Y, Z, XY, YZ, XZ
弹性应变
1
1
EPEL:1, 2, 3
弹性主应变
1
-
EPEL:EQV
等效弹性应变 [7]
1
1
EPTH:X, Y, Z, XY, YZ, XZ
平均热应变
1
1
EPTH:EQV
等效热应变 [7]
1
1
EPPL:X, Y, Z, XY, YZ, XZ
塑性应变均值
4
4
EPPL:EQV
等效塑性应变 [7]
4
4
EPCR:X, Y, Z, XY, YZ, XZ
蠕变均值
4
4
EPCR:EQV
等效蠕变 [7]
4
4
NL:EPEQ
等效塑性应变均值
4
4
NL:SRAT
Ratio of trial stress to stress on yield surface
4
4
NL:SEPL
应力应变曲线上的平均等效应力(混凝土的)
4
4
NL:HPRES
静水压力
-
4
THETCR, PHICR
裂缝的方向角THETA和 PHI
1
1
STATUS
单元状态
2
2
IRF
钢筋号
3
-
MAT
材料号
3
-
VR
体积率
3
-
THETA
X-Y平面内的方向角
3
-
PHI
X-Y平面外的方向角
3
-
EPEL
单轴弹性应变
3
-
S
单轴弹性应力
3
-
EPEL
单轴弹性应变均值
5
5
EPPL
单轴塑性应变均值
5
5
SEPL
应力应变曲线上的平均等效应力(钢筋的)
5
5
EPCR
单轴蠕变均值(钢筋)
5
5
1、 混凝土求解项目(输出每个积分点(如KEYOPT(5)=1)和质心)
2、 单元状态表(SOLID65单元状态表)对应项含意:
l Crushed:混凝土被压碎
l Opend:混凝土开裂,裂缝张开
l Colesed:混凝土开裂,但裂缝闭合
l Neither:混凝土即未开裂也未压碎
3、 钢筋求解项给出每种钢筋情况
4、 混凝土非线性积分点解(如KEYOPT(6)=3 and the element has a nonlinear material)
5、 钢筋非线性积分点解(如KEYOPT(6)=3 and the rebar has a nonlinear material)
6、 仅在质心利用*GET命令可得的选项。
7、 等效应变用一个有效的泊松比:对于弹性和热量问题这个值通过“MP,PRXY”命令设定,对于塑性和蠕变问题这个值被设为0.5。
SOLID65混合单元输出项
Description
Names of Items Output
O
R
Nodal Stress Solution
TEMP, S(X, Y, Z, XY, YZ, XZ), SINT, SEQV
1
-
1、如果KEYOPT(5)=2,则输出每个节点的值。
SOLID65单元状态对应表
状态
方向 1的状态
方向2的状态
方向 3的状态
1
Crushed
Crushed
Crushed
2
Open
Neither
Neither
3
Closed
Neither
Neither
4
Open
Open
Neither
5
Open
Open
Open
6
Closed
Open
Open
7
Closed
Open
Neither
8
Open
Closed
Open
9
Closed
Closed
Open
10
Open
Closed
Neither
11
Open
Open
Closed
12
Closed
Open
Closed
13
Closed
Closed
Neither
14
Open
Closed
Closed
15
Closed
Closed
Closed
16
Neither
Neither
Neither
“SOLID65项目和序号表”中列出了在后处理中可通过ETABLE命令加参数及数字序号的方法定义可列表察看的有关变量的细则。详细参见《ANSYS基本分析指南》中有关“The General Postprocessor (POST1)”和“The Item and Sequence Number Table”部分。下面是表格的一些使用说明:
Name
指在“SOLID65单元输出数据说明表”中的有关变量。
Item
命令ETABLE中使用的参数。
I,J,...,P
节点I,J,...,P所对应的数字序号。
IP
积分点对应的数字序号
SOLID65项目和序号表
Output Quantity Name(变量名)
ETABLE和 ESOL 命令输入项
Item(参数)
Rebar 1
Rebar 2
Rebar 3
EPEL单轴弹性应变均值
SMISC
SIG
SMISC
EPPL单轴塑性应变均值
NMISC
EPCR单轴蠕变均值
NMISC
SEPL应力应变曲线上的平均等效应力
NMISC
SRAT
NMISC
Output Quantity Name(变量名)
ETABLE 和 ESOL命令输入项
Item
(参数)
I
J
K
L
M
N
O
P
P1 (压力作用面1)
SMISC
P2 (压力作用面2)
SMISC
P3 (压力作用面3)
SMISC
P4 (压力作用面4)
SMISC
P5 (压力作用面5)
SMISC
P6 (压力作用面6)
SMISC
S:1 主应力1
NMISC
S:2 主应力2
NMISC
S:3 主应力3
NMISC
S:INT 应力强度
NMISC
S:EQV 等效应力
NMISC
FLUEN作用在各节点上的影响
NMISC
Output Quantity Name
(变量名)
ETABLE 和 ESOL命令输入项
Item
Integration Point
STATUS 单元状态
NMISC
Dir 1
THETCR 平面内的方向角
NMISC
PHICR 平面外的方向角
NMISC
Dir 2
THETCR 平面内的方向角
NMISC
PHICR 平面外的方向角
NMISC
Dir 3
THETCR 平面内的方向角
NMISC
PHICR 平面外的方向角
NMISC
SOLID65的假定和限制:
l 不允许使用零体积单元;
l 单元可以如“几何模型图”编号,也可将图中的IJKL面与MNOP面的编号交换。同时单元不能被扭转导致形成两个单独体,在对单元不恰当编号时这种情况最可能发生。
l 所有的单元应有八个节点
l 当K与L,O与P节点重合时单元形状便成为棱柱体形,当然也可退化为四面体。其它的形状会被四面体自动替换。
l 当考虑单元的钢筋作用时,钢筋被假定分散在整个单元中,且所有钢筋的总体积率不能起过1.0。
l 单元是非线性的故要求迭代求解。
l 当同时考虑混凝土的开裂与压碎时,应注意要缓慢加载,以免在实际可承受荷载通过闭合裂缝传递前出现混凝土的假压碎现象。这种现象通过泊松效应常常发生在与大量开裂应变垂直的未开裂的方向上。同样也会在压碎的积分点上出现,输出的塑性和蠕变应变值来自于先前子步的收敛。而且,当裂缝已经产生,则弹性应变的输出量就包含了开裂应变。单元开裂或压碎后失去的抗剪作用将不能被传递到钢筋上,因为钢筋没有抗剪刚度。
l 在考虑开裂或压碎的材料非线性问题中以下两项最好不要考虑:
n 应力强化效应。
n 大应变、大变形。
否则,结果可能不收敛或不正确,特别在有大转角情况下。 |