2006年10月22日星期日

Diana的学习和应用

Diana的学习和应用

公司官方介绍:
Displacement analyser (Diana)
DIANA是由荷兰的TNO DIANA公司开发的通用有限元结构分析软件包。
  该软件以土木建筑工程为主,在世界各地拥有很高的声誉。 特别是在混凝土的裂缝开展分析,钢结构的局部稳定分析, 复合结构的极限承载力分析,地基的分段施工分析等非线性分析方面,得到了广大用户的高度评价。
2003年1月,建筑研究所负责开发DIANA软件的部门从荷兰应用科学研究所独立出来, 与英国的FEMSYS公司合并,成立了TNO DIANA公司. 该公司是荷兰应用科学研究所下属的百分之百的子公司。
分析类型:
线弹性静力分析(Linear Static Analysis)
特征值分析(Eigenvalue Analysis)
稳定分析(Euler Buckling Analysis)
非线性分析(Nonlinear Analysis)
静力问题 动力问题
材料非线性 几何非线性
动力分析(Dynamic Analysis)
频谱响应分析 时程响应分析
分段施工分析(Phased Analysis)
热传导分析、渗透分析:【线性、非线性问题、稳态分析、瞬态分析、与应力分析的连动】
流体-应力分析:[热应力分析、压密分析]
Lattice分析
管道线分析(Pipeline Analysis)
逆向分析(Parameter Estimation Analysis)
地基基础的液化分析(3种本 模型) (Liquefaction Analysis)
混凝土弹塑性断裂模型(Maekawa Model)


单元类型:
杆单元
梁单元
平面应力单元
平面应变单元
轴对称单元
弯曲板单元
平面壳单元
曲面壳单元
固体单元
界面单元
弹簧单元
结点质量单元
埋入式钢筋单元
流体单元
边界单元
埋入式钢筋单元

材料特性:
各向同性/各向异性材料
非线性弹性分析
弹塑性分析
粘塑性分析
粘弹性徐变分析
干燥收缩分析
裂缝分析(分散型裂缝模型,
离散型裂缝模型)
超弾性問題
界面单元的非线性材料特性
用户定义子程序功能
混凝土弹塑性断裂模型
钢材的移动硬化准则

载荷:
载荷
集中载荷
分布载荷
重力
离心力
位移载荷
钢筋单元的预应力
风载荷
水载荷
加速度
各种单元的初始应力

初始条件:
变形
自重应力(K0状态)
初始应力
初始速度
失稳后的変形

我的学习心得:
零零散散地学了一段时间的Diana,也不知道入门了没有
下面是我对Diana的一些肤浅的认识
优点:
1.对于钢筋的处理有其独到之处,可以比较简单的加入分离式的钢筋,并施加预应力。
2.对于混凝土的本构模型,有比较多的选择,从受压软化到受拉软化,模型较多,但也还是不够全面,比如W-W五参数模型就没有。
3.有lattice分析,不过它的建模比较麻烦。
4.有施工阶段分析,这点sap2000也有,但是对于复杂的实体非线性,sap2000还差一些。
5.命令流也挺过瘾,而且似乎还挺简洁的。
6.基本都是键盘操作,有时候可以很快,和早期的algor有点象。
7.分析的filo文件可以保存,挺不错的。
7.有些后处理操作做的不错。


不足
1.虽然键盘操作不错,但是鼠标功能还是差了些,整个界面用的是femsys的界面,太dos了一些。
2.命令流里面没有办法设置复杂一些的操作,比如循环,高级的参数变化,比起其他用python,lisp,apdl的要逊色一些。
3.混凝土有再多一些本构就好了。
4.帮助太过含糊,有许多不明之处,希望能象abaqus。

下面是我学习的第一个例子,按照学习光盘的录像来的,是前处理部分,计算之后可以得到如图的显示效果
不过和他们演示光盘中的结果略微不同。

有限元软件Diana入门例子:
[code]
!analysis of prestressed T beam
FEMGEN D9S
PROPERTY FE-PROG DIANA STRUCT_3D
yes
UTILITY SETUP UNITS LENGTH MILLIMETER
UTILITY SETUP UNITS MASS KILOGRAM
UTILITY SETUP UNITS TIME SECOND
GEOMETRY POINT 0
GEOMETRY SWEEP P1 TRANSLATE 0 100
EYE FRAME
GEOMETRY SWEEP L1 TRANSLATE 0 0 300
DRAWING DISPLAY
EYE FRAME
DRAWING DISPLAY
EYE FRAME
EYE LOCATE
GEOMETRY SWEEP L2 TRANSLATE 0 0 100
GEOMETRY SWEEP L6 TRANSLATE 0 0 -100
UTILITY UNDO
GEOMETRY SWEEP L6 TRANSLATE 0 -100
GEOMETRY SWEEP L7 TRANSLATE 0 100
MESHING DIVISION ELSIZE ALL 50
MESHING DIVISION ELSIZE ALL 50
LABEL GEOMETRY SURFACES S3
MESHING DIVISION ELSIZE ALL 50
EYE FRAME
LABEL GEOMETRY DIVISIONS
MESHING DIVISION LINE L1
MESHING DIVISION LINE L2
MESHING DIVISION LINE L5
EYE FRAME
CONSTRUCT SET OPEN
CONSTRUCT SET OPEN FACE1
CONSTRUCT SET APPEND ALL
CONSTRUCT SET CLOSE
MESHING TYPES ALL QU8 CQ40S
MESHING GENERATE
VIEW MESH CURRENT
VIEW GEOMETRY CURRENT

GEOMETRY SWEEP FACE1 FACE2 40 TRANSLATE 4000
EYE FRAME
GEOMETRY SWEEP FACE2 FACE3 12 TRANSLATE 1200
EYE FRAME
CONSTRUCT SET OPEN BEAM
CONSTRUCT SET APPEND ALL
CONSTRUCT SET CLOSE
!EYE ZOOM .778 .384 .7819999 .378
EYE ZOOM .828 .385 .957 .167
LABEL GEOMETRY LINES ALL
GEOMETRY MOVE L37 TRANSLATE 0 0 80
yes



EYE DIRECTION 0 -1 0
EYE FRAME
CONSTRUCT SET OPEN TMP
CONSTRUCT SET APPEND CURSOR POLYGON .899 .555 .9 .499 .969 .504 .965 .562 .887 .544
CONSTRUCT SET CLOSE
GEOMETRY MOVE TMP TRANSLATE 0 0 -40
yes

EYE FRAME
MESHING TYPES ALL NONE
MESHING TYPES ALL HE20 CHX60
MESHING GENERATE
VIEW MESH CURRENT
VIEW HIDDEN SHADE
VIEW HIDDEN SHADE
VIEW GEOMETRY CURRENT
VIEW HIDDEN FILL
VIEW HIDDEN OFF
EYE FRAME

EYE FRAME
CONSTRUCT SET OPEN
CONSTRUCT SET OPEN TOP
CONSTRUCT SET APPEND CURSOR POLYGON .576E-1 .5470001 .588E-1 .523 .714 .524 .9359999 .517 .938 .542 .711 .554 .42E-1 .542
CONSTRUCT SET CLOSE
VIEW GEOMETRY +TOP RED
EYE FRAME
EYE ZOOM .113 .692 .241 .538
PROPERTY BOUNDARY CONSTRAINT P3 PINNED
PROPERTY BOUNDARY CONSTRAINT P4 PINNED
EYE FRAME
EYE ZOOM .6580001 .531 .786 .315
LABEL GEOMETRY LINES L14
PROPERTY BOUNDARY CONSTRAINT L14 Y Z
LABEL MESH CONSTRNT

EYE FRAME
GEOMETRY POINT 0 50 350
GEOMETRY POINT 2000 50 100
GEOMETRY POINT 4000 50 350
GEOMETRY POINT 5200 50 200
REINFORCE BAR SECTION P31 P32 P33
REINFORCE BAR SECTION P33 P34
REINFORCE BAR BAR1 RE1 RE2

PROPERTY MATERIAL BE1 ELASTIC ISOTROP 35000 0.2
PROPERTY MATERIAL BE1 MASS DENSITY 2400E-9 0
PROPERTY MATERIAL BE1 STATNONL CONCBRIT CRACK CONSTA TENSIO1 ULTIMATE TAUCRI1 NONE 2.5 0.0003 0.2

PROPERTY MATERIAL RP1 ELASTIC REINFORC NOBOND 2.1E5 ""
PROPERTY MATERIAL RP1 STATNONL REINFORC VMISES NONE 500
PROPERTY PHYSICAL PH1 GEOMETRY REINFORC BAR 1000

PROPERTY ATTACH BEAM BE1
PROPERTY ATTACH BAR1 RP1
PROPERTY ATTACH BAR1 PH1
LABEL GEOMETRY MATERIALS


PROPERTY LOADS PRESTRESS BAR1 100
PROPERTY LOADS GRAVITY ALL -9.8 Z
PROPERTY LOADS PRESSURE LO3 2 TOP -0.05 Z

save
yes

UTILITY WRITE DIANA D9S.DAT
yes
ANALYSE D9S
[/code]

前段时间曾经把所有的混凝土本构都学习了一遍,也做了些学习笔记,有空再整理一下。

图1是我自己做的textpad的高亮文件,查错还是挺不错的:)

图2是例子的结果抓图