solid65.docx

solid65.docx

《solid65.docx》由会员分享，可在线阅读，更多相关《solid65.docx（17页珍藏版）》请在冰豆网上搜索。

solid65

solid65

SOLID65单元性质

SOLID65单元描述：

SOLID65单元用于含钢筋或不含钢筋的三维实体模型。

该实体模型可具有拉裂与压碎的性能。

在混凝土的应用方面，如用单元的实体性能来模拟混凝土，而用加筋性能来模拟钢筋的作用。

当然该单元也可用于其它方面，如加筋复合材料（如玻璃纤维）及地质材料（如岩石）。

该单元具有八个节点，每个节点有三个自由度，即x,y,z三个方向的线位移；还可对三个方向的含筋情况进行定义。

本单元与SOLID45单元（三维结构实体单元）的相似，只是增加了描述开裂与压碎的性能。

本单元最重要的方面在于其对材料非线性的处理。

其可模拟混凝土的开裂（三个正交方向）、压碎、塑性变形及徐变，还可模拟钢筋的拉伸、压缩、塑性变形及蠕变，但不能模拟钢筋的剪切性能。

有关SOLID65单元的更细节的描述请参见《ANSYS理论手册》。

SOLID65的几何模型图

SOLID65输入数据：

关于单元几何图形、节点位置、单元坐标系请见上图。

单元性质为八节点各向同性材料，单元包括一种实体材料和三种钢筋材料，用命令MAT输入对混凝土材料的定义，而有关钢筋的细则需在实常数中定义，包括材料号、体积率、方向角（THETA,PHI），钢筋的方向角可通过命令/Eshape以图示方式校验。

体积率是指钢筋的体积与整个单元体积的比，钢筋的方向通过单元坐标系中的两个角度（度制）来定义。

当钢筋的材料号为0或等于单元的材料号时则不考虑它的作用。

另外，有关混凝土的材料定义，如剪切传递系数，拉应力，压应力都应在数据表中给出，详细描述见表“SOLID65混凝土材料数据表”。

通常剪力传递系数为0～1.0，0表示平滑的裂缝（完全丧失剪力传递作用），1表示粗糙的裂缝（几乎没有失去剪力传递作用）。

这就有利于对裂缝开裂与闭合进行描述。

有关单元荷载的描述见“节点单元荷载”（ANSYS帮助中专有一节）。

压力作为面荷载作用在单元表面如“SOLID65的几何模型图”中带圈数字所示。

主动力作用在单元内。

温度和热流量可在节点上作为单元体荷载输入。

节点I的温度T（I）默认为TUNIF，如其它节点温度没有被指定，则它们默认为T（I）。

对于其它的输入模型未指定温度时默认值都为TUNIF。

对热流量的设定除用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,PRXYorNUXY,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）

热流量-- 

FL（I）,FL（J）,FL（K）,FL（L）,FL（M）,FL（N）,FL（O）,FL（P）

特性：

Plasticity（塑性）

Creep（蠕变\徐变）

Cracking（开裂）

Crushing（压碎）

Largedeflection（大变形）

Largestrain（大应变）

Stressstiffening（应力强化）

Birthanddeath　（单元生死）

Adaptivedescent

KEYOPT

（1）大变形控制:

0 -- 

考虑大变形

1 -- 

不考虑

KEYOPT（5）混凝土线性解的输出控制:

0 -- 

只打印质心的线性解

1 -- 

给出每个积分点的解

2 -- 

输出节点应力

KEYOPT（6）混凝土非线性解输出控制:

0 -- 

只打印质心的线性解

3 -- 

同时还给出积分点的解

KEYOPT（7）开裂后应力松驰考虑选项:

0 -- 

不考虑拉伸应力松驰

1 -- 

考虑应力松驰，有助于计算收敛

SOLID65单元信息：

列在“SOLID65混凝土材料数据表”中的数据项是通过TB命令来建立的。

表中未重新赋值的项默认为0或其特定默认值。

数据表首先用TB（withLAB=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输出数据：

与单元相关的解答输出项有以下两方面：

●所有节点的节点位移；

●其它输出项见“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

热流量

作用在各节点上的热流量FL（I）,FL（J）,FL（K）,FL（L）,FL（M）,FL（N）,FL（O）,FL（P）

Y

Y

XC,YC,ZC

Locationwhereresultsarereported

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

三轴应力与屈服面应力之比

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单元状态表）对应项含意：

●Crushed：

混凝土被压碎

●Opend：

混凝土开裂，裂缝张开

●Colesed：

混凝土开裂，但裂缝闭合

●Neither：

混凝土即未开裂也未压碎

3、钢筋求解项给出每种钢筋情况

4、混凝土非线性积分点解（如KEYOPT（6）＝3　andtheelementhasanonlinearmaterial）

5、钢筋非线性积分点解（如KEYOPT（6）＝3　andtherebarhasanonlinearmaterial）

6、仅在质心利用＊GET命令可得的选项。

7、等效应变用一个有效的泊松比：

对于弹性和热量问题这个值通过“MP，PRXY”命令设定，对于塑性和蠕变问题这个值被设为0.5。

SOLID65混合单元输出项

Description

NamesofItemsOutput

O

R

NodalStressSolution

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基本分析指南》中有关“TheGeneralPostprocessor（POST1）”和“TheItemandSequenceNumberTable”部分。

下面是表格的一些使用说明：

Name

指在“SOLID65单元输出数据说明表”中的有关变量。

Item

命令ETABLE中使用的参数。

I,J,...,P

节点I,J,...,P所对应的数字序号。

IP

积分点对应的数字序号

SOLID65项目和序号表

OutputQuantityName（变量名）

ETABLE和ESOL命令输入项

Item（参数）

Rebar1

Rebar2

Rebar3

EPEL单轴弹性应变均值

SMISC

1

3

5

SIG

SMISC

2

4

6

EPPL单轴塑性应变均值

NMISC

41

45

49

EPCR单轴蠕变均值

NMISC

42

46

50

SEPL应力应变曲线上的平均等效应力

NMISC

43

47

51

SRAT

NMISC

44

48

52

OutputQuantityName（变量名）

ETABLE和ESOL命令输入项

Item

（参数）

I

J

K

L

M

N

O

P

P1　（压力作用面1）

SMISC

8

7

10

9

-

-

-

-

P2　（压力作用面2）

SMISC

11

12

-

-

14

13

-

-

P3　（压力作用面3）

SMISC

-

15

16

-

-

18

17

-

P4　（压力作用面4）

SMISC

-

-

19

20

-

-

22

21

P5　（压力作用面5）

SMISC

24

-

-

23

25

-

-

26

P6　（压力作用面6）

SMISC

-

-

-

-

27

28

29

30

S:

1　主应力1

NMISC

1

6

11

16

21

26

31

36

S:

2　主应力2

NMISC

2

7

12

17

22

27

32

37

S:

3　主应力3

NMISC

3

8

13

18

23

28

33

38

S:

INT　应力强度

NMISC

4

9

14

19

24

29

34

39

S:

EQV　等效应力

NMISC

5

10

15

20

25

30

35

40

热流量作用在各节点上的影响

NMISC

109

110

111

112

113

114

115

116

OutputQuantityName

（变量名）

ETABLE和ESOL命令输入项

Item

IntegrationPoint

1

2

3

4

5

6

7

8

STATUS　单元状态

NMISC

53

60

67

74

81

88

95

102

Dir1

THETCR　平面内的方向角

NMISC

54

61

68

75

82

89

96

103

PHICR　平面外的方向角

NMISC

55

62

69

76

83

90

97

104

Dir2

THETCR　平面内的方向角

NMISC

56

63

70

77

84

91

98

105

PHICR　平面外的方向角

NMISC

57

64

71

78

85

92

99

106

Dir3

THETCR　平面内的方向角

NMISC

58

65

72

79

86

93

100

107

PHICR　平面外的方向角

NMISC

59

66

73

80

87

94

101

108

SOLID65的假定和限制：

●不允许使用零体积单元；

●单元可以如“几何模型图”编号，也可将图中的IJKL面与MNOP面的编号交换。

同时单元不能被扭转导致形成两个单独体，在对单元不恰当编号时这种情况最可能发生。

●所有的单元应有八个节点

●当K与L，O与P节点重合时单元形状便成为棱柱体形，当然也可退化为四面体。

其它的形状会被四面体自动替换。

●当考虑单元的钢筋作用时，钢筋被假定分散在整个单元中，且所有钢筋的总体积率不能起过1.0。

●单元是非线性的故要求迭代求解。

●当同时考虑混凝土的开裂与压碎时，应注意要缓慢加载，以免在实际可承受荷载通过闭合裂缝传递前出现混凝土的假压碎现象。

这种现象通过泊松效应常常发生在与大量开裂应变垂直的未开裂的方向上。

同样也会在压碎的积分点上出现，输出的塑性和蠕变应变值来自于先前子步的收敛。

而且，当裂缝已经产生，则弹性应变的输出量就包含了开裂应变。

单元开裂或压碎后失去的抗剪作用将不能被传递到钢筋上，因为钢筋没有抗剪刚度。

●在考虑开裂或压碎的材料非线性问题中以下两项最好不要考虑：

⏹应力强化效应。

⏹大应变、大变形。

否则，结果可能不收敛或不正确，特别在有大转角情况下。