地震作用下结构弹塑性动力分析算例Word格式.docx

资源描述

地震作用下结构弹塑性动力分析算例Word格式.docx

《地震作用下结构弹塑性动力分析算例Word格式.docx》由会员分享，可在线阅读，更多相关《地震作用下结构弹塑性动力分析算例Word格式.docx（13页珍藏版）》请在冰豆网上搜索。

地震作用下结构弹塑性动力分析算例Word格式.docx

步长0.02s,

我时程分析的时候用0.001s

！

在森特罗波里（常用地震波文件）

当在结构中安装摩擦耗能装置后，安装的摩擦耗能装置情况如下：

从一层到6层，耗能支撑的刚度：

（KN/MM）

2671.22185.32261.11245.11800.3589.6

耗能装置的起滑位移，从一层到6层：

5.66.98.17.42.92.1

周福霖的结果如下（地震加速度峰值为8度大震400gal,与上面相同）

各层位移反应（从一层到6层）（mm）

15.48.616.078.448.448.05

我的结果如下：

26.69.646.289.8216.78.03

欢迎有兴趣的朋友帮我验证一下，

谢谢,我的qq:

47932480,基本上在线，

你塑性分析的恢复力模型是取的哪种模型？

三折线模型吗？

另外你未给出耗能装置的耗能起滑力，是否直接依据耗能装置的刚度与起滑位移来确定起滑力！

我取的是退化四线型,为了比较,简化为理想弹塑性后与周福霖比较,

起滑力是根据刚度和起滑位移相乘得到,

没有加耗能装置的时候,

ansys的命令流文件如下:

/BATCH

/COM,ANSYSRELEASE5.7UP2000120810:

20:

4310/06/2002

/input,menust,tmp,'

,,,,,,,,,,,,,,,1

/GRA,POWER

/GST,ON

/PLO,INFO,3

/COL,PBAK,ON,1,BLUE

/NOPR

/PMETH,OFF,0

KEYW,PR_SET,1

KEYW,PR_STRUC,1

KEYW,PR_THERM,0

KEYW,PR_FLUID,0

KEYW,PR_MULTI,0

/GO

/COM,

/COM,PreferencesforGUIfilteringhavebeensettodisplay:

/COM,Structural

/PREP7

ET,1,COMBIN40

KEYOPT,1,1,1

KEYOPT,1,3,0

KEYOPT,1,4,0

KEYOPT,1,6,2

R,1,583970000,0,1057142.857,0,8350771,0,

R,2,583264000,0,952040.8163,0,7115820.8,0,

R,3,583264000,0,952040.8163,0,7290800,0,

R,4,474254000,0,952040.8163,0,4789965.4,0,

R,5,474254000,0,952040.8163,0,4789965.4,0,

R,6,474254000,0,952040.8163,0,4315711.4,0,

R,6,454756000,0,625510.2041,0,3819950.4,0,

R,5,474254000,0,952040.8163,0,4315711.4,0,

N,0,3.6,0,0,,,,

NDELE,1

N,1,3.6,0,0,,,,

N,1,4.0,0,0,,,,

N,2,7.6,0,0,,,,

N,3,11.2,0,0,,,,

N,4,14.8,0,0,,,,

N,5,18.4,0,0,,,,

N,6,22,0,0,,,,

N,7,25.6,0,0,,,,

TYPE,1

MAT,

REAL,1

ESYS,0

SECNUM,

TSHAP,LINE

FLST,2,2,1

FITEM,2,1

FITEM,2,2

E,P51X

REAL,2

FITEM,2,3

REAL,3

FITEM,2,4

REAL,4

FITEM,2,5

REAL,5

FITEM,2,6

REAL,6

FITEM,2,7

FINISH

/SOLU

ANTYPE,4

TRNOPT,FULL

LUMPM,0

FLST,2,1,1,ORDE,1

D,P51X,,,,,,ALL,,,,,

*SET,NT,751

*SET,DT,0.02

*dim,ac,,NT

/input,aaaa,txt

41:

4910/06/2002

RESUME

5910/06/2002

42:

0110/06/2002

/EXIT,NOSAV

43:

0010/06/2002

KEYOPT,1,4,1

N,1,0,0,0,,,,

N,2,3.6,0,0,,,,

N,2,4.0,0,0,,,,

N,3,7.6,0,0,,,,

N,4,11.2,0,0,,,,

N,5,14.8,0,0,,,,

N,6,18.4,0,0,,,,

N,7,22,0,0,,,,

*SET,NT,1501

NSUBST,20,0,0

ALPHAD,00.452709

BETAD,0.004295

TIME,0.02

NCNV,0,0,0,0,0

RESCONTRL,DEFINE,ALL,1,1

*do,i,1,751

acel,ac（i）,0,0

time,i*0.02

solve

*enddo

/POST26

GPLOT

NSOL,2,2,U,X,

PRVAR,2,,,,,,

TIME,15.02

AUTOTS,-1

DELTIM,0.001,,,1

KBC,0

TSRES,ERASE

*dim,acc,NT

*dim,acc,,NT

acel,acc（i）,0,0

NSOL,3,3,U,X,

PRVAR,3,,,,,,

*STAT

*STATUS,ACC,1,1501,1,,1,,

*dim,a,,1501

*STATUS,A,1,1501,1,,1,,

*dim,acc,1501

*dim,acc,,1501

*do,i,1,751,1

NSOL,4,4,U,X,

PRVAR,4,,,,,,

NSOL,5,5,U,X,

PRVAR,5,,,,,,

NSOL,6,6,U,X,

PRVAR,6,,,,,,

NSOL,7,7,U,X,

PRVAR,7,,,,,,

PLVAR,6,,,,,,,,,,

PLVAR,7,,,,,,,,,,

PLVAR,5,,,,,,,,,,

PLVAR,4,,,,,,,,,,

PLVAR,3,,,,,,,,,,

PLVAR,2,,,,,,,,,,

PLVAR,1,,,,,,,,,,

QUOT,5,1,2,,,,,1,1,

ADD,5,-1,-1,,,,,2,3,0,

ADD,7,1,-1,,,,,2,3,0,

NSOL,3,3,U,X,UX

NSOL,2,2,U,X,UX

ADD,7,1,-1,,,,,3,2,0,

ADD,7,3,2,,,,,1,-1,0,

运行结果:

各层位移（相对地面的位移）

3.716.918.2714.916.316.6

我的程序结果:

3.716.908.2814.9816.3716.60

我发觉ansys的阻尼矩阵是时时更新,

第一次,我的程序结果是阻尼矩阵的计算都采用第一刚度系数,

这一次,我把阻尼矩阵时时更新,结果和ansys差不多,,

我用sap2000的结果如下,（没有加耗能装置）

各层位移:

3.596.498.7710.812.0812.60

我把sap2000的相关文件放在论坛上（sap2000）专栏,

由于sap2000在这方面的一些帮助文件我没有看明白,

因此,这个也许是错的哦

大家一起看一下,

当安装了耗能装置后:

各层安装情况如下:

支撑刚度

KN/mm）

5571.05829.85717.32141.22242.94487.4

起滑位移（mm）

10.84.43.15.33.11.7

ansys的结果:

各层位移（mm）

4.868.8013.3421.6526.9527.93

我的程序（mm）

4.969.0013.3721.6626.9727.93

ansys命令流文件与上基本相同,只是改一下combin40单元的实参数.

加了这样的耗能装置后

一层到6层（KN/mm）

这样的设置与上贴的区别是,此时原结构也进入了弹塑性状态:

用ansys分析的时候,在每个节点之间用两个combin40单元连接.

ansys的结果如下:

2.533.403.9647.5254.5859.25

我的程序结果（mm）

2.533.403.9747.5554.6859.24

回答如下：

1。

我用的是wilson-theta法，

虽然时程分析的方法很多，但是都是求微分方程的解，因此，结果应该差不多。

2、理想弹塑性模型－应该算是双线性模型的特例吧。

第二刚度为0

3、阻尼矩阵的alpha和beta由初始刚度确定然后保持不变。

阻尼矩阵根据当时的K生成。

4。

滞回模型各控制点的数值已经假定已知。

已给出。

以您现在给出的条件是可以计算的了。

因为我的程序是针对砌体结构编制的，需要改造，所以没能为您计算这道题，有时间我会用“高阶单步法”（一个比威尔逊精度高点的计算方法）给您算一下的。