钢筋混凝土分析，各位高手帮助我成长！

coinage

各位高手：
      小弟以前一直用ansys做结构分析，但是越来越发现ansys分析混凝土的结果不太理想，在网上看到许多高手推荐abaqus，所以我也在学了。刚学了几天，有一点收获，对于我要分析的问题写了一段inp文件，发上来给各位看看，帮我查查错误，也希望能够多提意见。
      我要分析的问题是混凝土结构植筋的问题，简单来说就是在混凝土上钻孔，灌植筋胶，插入钢筋，使得钢筋和混凝土粘结工作。
      我的思路是：将混凝土和钢筋分开建模，混凝土C3D8单元，钢筋用Truss单元，然后再embed。
      请各位高手帮助小弟成长!谢谢!

:(:(:(
忙了几天，在原来的基础上改了一下
不过运行的时候有错误

Unknown keyword "fileformate". The keyword may be misspelled, obsolete, or invalid.
in keyword *STATIC, file "TRY1.inp", line 108: Unknown part instance name 0
Analysis Input File Processor exited with an error.

Boundary conditions are specified on inactive dof of 3 nodes. The nodes have been identified in node set WarnNodeBCInactiveDof.

修改后的inp文件
大家帮我看看吧

*HEADING
**JOB NAME:ANCHOR GROUP   REBAR DIAMETER=18 SAPCE=8D DEPTH=15D
*PREPRINT,ECHO=YES,HISTORY=YES,MODEL=YES
*RESTART,WRITE,FREQ=1
*FILE FORMATE,ZERO INCREMENT
**UNIT,KN,METER
**==========混凝土建模===========
*PART,NAME=CONCRETE
*NODE
1,0.,0.,0.
21,1.,0.,0.
421,0.,1.,0.
441,1.,1.,0.
10001,0.,0.,0.5
10021,1.,0.,0.5
10421,0.,1.,0.5
10441,1.,1.,0.5
*NGEN,NSET=NBT1
1,21,1
*NGEN,NSET=NBT2
421,441,1
*NFILL,NSET=NBUTTOM
NBT1,NBT2,20,21
*NGEN,NSET=NTP1
10001,10021,1
*NGEN,NSET=NTP2
10421,10441,1
*NFILL,NSET=NTOP
NTP1,NTP2,20,21
*NFILL,NSET=NALL
NBUTTOM,NTOP,10,1000
*ELEMENT,TYPE=C3D8
1,1,2,23,22,1001,1002,1023,1022
*ELGEN,ELSET=EALL
1,20,1,1,20,21,21,10,1000,1000
*SOLID SECTION,ELSET=EALL,MATERIAL=M-1
1.0
*END PART
**===========钢筋建模==============
*PART,NAME=REBAR
*NODE
1,0.,0.,0.
6,0.,0.,0.27
7,0.144,0.,0.
12,0.144,0.,0.27
13,0.288,0.,0.,
18,0.288,0.,0.27
*NGEN
1,6,1
7,12,1
13,18,1
*NSET,NSET=NL
6,12,18
*ELEMENT,TYPE=T3D2
1,1,2
2,2,3
3,3,4
4,4,5
5,5,6
7,7,8
8,8,9
9,9,10
10,10,11
11,11,12
13,13,14
14,14,15
15,15,16
16,16,17
17,17,18
*ELSET,ELSET=REBARALL,GEN
1,5,1
7,11,1
13,17,1
*SOLID SECTION,ELSET=REBARALL,MATERIAL=M-2
0.0002545
*END PART
**============ASSEMBLY=============
*ASSEMBLY,NAME=ASSEMBLY
*INSTANCE,NAME=CONCRETE-1,PART=CONCRETE
0.,0.,0.
*END INSTANCE
*INSTANCE,NAME=REBAR-1,PART=REBAR
0.356,0.5,0.23
*END INSTANCE
*EMBEDDED ELEMENT,HOST ELSET=CONCRETE-1.EALL
REBAR-1.REBARALL
*END ASSEMBLY
**===========材料性质==============
*MATERIAL,NAME=M-1
*ELASTIC
3.0E10,0.2
*CONCRETE DAMAGED PLASTICITY
38.0
*CONCRETE COMPRESSION HARDENING
2.0E7,0.0
3.0E7,0.0015
*CONCRETE TENSION STIFFENING,TYPE=DISPLACEMENT
3.33E6,0.0
3.33E4,7.447447447E-5
*MATERIAL,NAME=M-2
*ELASTIC
2.9E11
*PLASTIC
3.2E10
*BOUNDARY
CONCRETE-1.NBUTTOM,1,3
*STEP,INC=1000,NLGEOM=YES
*STATIC,RIKS
0.001,1,,0.02,0.1,0.02,3,
*BOUNDARY
REBAR-1.NL,3,3,5E-3
*OUTPUT,FIELD,VARIABLE=PRESELECT
*END STEP

[ 本帖最后由 coinage 于 2007-3-8 22:29 编辑 ]

coinage

ABAQUS的前处理好像比ansys差了一些，ansys写好的命令流贴进命令窗口后按回车就可以运行了，而abaqus好像不行，而且调用inp文件不能生成cae文件，这在模型检查时会遇到问题，起码不方便，不知道各位怎么看