简单命令流的纠错

西风独自凉... 发表于 2012-3-13 09:38:20 | 显示全部楼层 | 阅读模式
悬赏0堡币
本帖最后由 西风独自凉 于 2012-3-13 09:42 编辑

当我输入egen,60,100,1,28,1,,,,,,,(Dis_ver-H_beam/2)/60 时,ansys崩溃,出现下面提示
1.JPG
请各位大侠指点,纠正命令流中错误,可适当加堡币
fini
/cle
/PREP7
/TITLE,Part 1:heat transfer analysis
ET,1,SOLID70 !定义单元类型
!-----------------------------------------------------------------
!定义参数
!-----------------------------------------------------------------
W_col=0.2039 !柱截面宽度
H_col=0.2062 !柱截面高度
tf_col=0.0125 !柱翼缘厚度
tw_col=0.0083 !柱腹板厚度
B_col=(W_col-tw_col)/2 !柱翼缘伸出长度
D_col=H_col-2*tf_col !柱腹板净高
W_beam=0.165 !梁截面宽度
H_beam=0.3038 !梁截面高度
tf_beam=0.0102 !梁翼缘厚度
tw_beam=0.0061 !梁腹板厚度
B_beam=(W_beam-tw_beam)/2 !梁翼缘伸出长度
D_beam=H_beam-2*tf_beam !梁腹板净高
Dis_hor=5.5 !框架水平跨间距
Dis_ver=3.0 !框架竖向层高
pp=(W_col-W_beam)/2
!----------------------------------------------------------------------------------
!定义热分析材料特性
!----------------------------------------------------------------------------------
!热分析需要定义的材料特性包括导热性,比热,密度等
MPTEMP,,20,800,900,1000 !定义随温度变化的钢材的导热性
MPDATA,KXX,1,,53.334,27.36,27.36,27.36
MPTEMP !清除温度场
MPTEMP,,20,100,180,260,380 !定义随温度变化的钢材的比热
MPDATA,C,1,,439.8,487.62,522.33,550.75,596.52
MPTEMP,,500,600,640,720,735
MPDATA,C,1,,666.5,759.92,798.67,1388,5000
MPTEMP,,750,830,900,1000
MPDATA,C,1,,1483,725,650,650
MP,DENS,1,7850 !定义钢材的密度
!---------------------------------------------------------------------------------
!建立分析模型
!---------------------------------------------------------------------------------
!采用直接生成节点和单元的办法建立实体模型。框架除防火间以外的部分不参与传热
!反应。因此,仅建立防火间的分析模型
!生成第一根柱
N,1,-H_col/2,,-W_col/2 !产生构成柱截面的节点
N,2,-H_col/2,,-W_col/2+pp
NGEN,4,1,2,,,,,(B_col-pp)/3
N,6,-H_col/2,,tw_col/2
NGEN,4,1,6,,,,,(B_col-pp)/3
N,10,-H_col/2,,W_col/2
NGEN,2,10,1,10,1,tf_col
NGEN,10,6,15,16,1,D_col/10
NGEN,2,60,11,20,1,D_col
NGEN,2,10,71,80,1,tf_col
NGEN,2,100,all,,,,(Dis_ver-H_beam/2)/60 !将已生成的节点沿y方向偏移
!(Dis_ver-H_beam/2)/60拷贝一层
!节点号加100
E,1,2,12,11,101,102,112,111 !生成第一层单元
EGEN,9,1,1
E,15,16,22,21,115,116,122,121
EGEN,10,6,10
E,71,72,82,81,171,172,182,181
EGEN,9,1,20
egen,60,100,1,28,1,,,,,,,(dis_ver-h_beam/2)/60!将第一层单元沿y方向拷贝60层
!生成第一根柱
!将柱的模型继续向上延伸梁的截面高度,生成梁柱节点。
!节点处单元尺寸尽量与梁的截面单元尺寸一致。
NSEL,S,NODE,,6001,6090,1 !沿梁的翼缘厚度生成一层单元
NGEN,2,100,ALL,,,,tf_beam
EGEN,2,100,1653,1680,1
NSEL,ALL
NSEL,S,NODE,,6101,6190 !沿梁的腹板高度生成10层单元
NGEN,11,100,ALL,,,,D_beam/10
EGEN,11,100,1681,1708,1
NSEL,ALL
NSEL,S,NODE,,7101,7190 !继续沿梁的翼缘厚度生成一层单元
NGEN,2,100,ALL,,,,tf_beam
EGEN,2,100,1961,1988,1
NSEL,ALL
!第一根实体模型柱完成。
!共计:节点79层,每层编号1-90, 逐层加100,顶层编号7801-7890
!单元78层,自动编号。每层28个,共28*78=2184个
!生成第二根柱及梁柱节点
NGEN,2,10000,ALL,,,Dis_hor !从第一根柱拷贝所有的节点
!节点号加10000
EGEN,2,10000,1,2184,1 !从第一根柱拷贝所有的单元
!生成梁
!梁被夹在两根柱之间,实际长度为Dis_hor-H_col
N,20001,H_col/2,Dis_ver-H_beam/2,-W_beam/2
!生成梁的截面的所有节点
!梁的节点编号从20001开始
NGEN,4,1,20001,,,,,B_beam/3
N,20005,H_col/2,Dis_ver-H_beam/2,tw_beam/2
NGEN,4,1,20005,,,,,B_beam/3
NGEN,2,10,20001,20008,1,,tf_beam
NGEN,10,6,20014,20015,1,,D_beam/10
NGEN,2,60,20011,20018,1,,D_beam
NGEN,2,10,20071,20078,1,,tf_beam
NGEN,2,100,20001,20090,,(Dis_hor-H_col)/100 !沿x方向偏移(Dis_hor-H_col)/100
!拷贝一层节点
E,20001,20002,20012,20011,20101,20102,20112,20111 !生成梁的第一层截面单元
!两根柱单元总数为4032
!故梁的单元编号从4033开始
EGEN,7,1,4033
E,20014,20015,20021,20020,20114,20115,20121,20120
EGEN,10,6,4040
E,20071,20072,20082,20081,20171,20172,20182,20181
EGEN,7,1,4050
EGEN,100,100,4033,4056,1,,,,,,(Dis_hor-H_col)/100 !沿x方向拷贝100层生成整根梁
!梁的实体模型完成
!总计:梁的节点为101层,每层编号1-88。从20001开始,逐层加100
!左端截面的节点为20001-20088; 右端截面的节点为30001-30088
!每层单元数为24个,总计24*100=2400个。单元编号为4369-6768
!------------------------------------------------------------------------------
!建立梁和柱连接处的耦合关系
!------------------------------------------------------------------------------
!自动耦合所有节点坐标重合的节点。梁翼缘的节点和柱的侧面完全重合,可以自动耦合。
!梁腹板的节点距离柱侧面相应节点的距离为(tw_col-tw_beam)/2=0.0011
!因此,设置耦合误差为0.002时,也能自动耦合。
CPINTF,all,0.002
FINISH

您需要登录后才可以回帖 登录 | 注册

本版积分规则

 
 
  • QQ:56984982
  • 点击这里给我发消息
    电话:13527553862
    站务咨询群桥头堡站务咨询桥梁专业交流群:
    中国桥梁专业领袖群
    工作时间
    8:00-18:00