ANSYS热结构耦合分析实例文档格式.docx

1、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

2、,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,

3、2,60,11,20,1,D_col NGEN,2,10,71,80,1,tf_colNGEN,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 !沿翼缘长度拷贝9个单元E,15,16,22,21,115,116,122,121 !定义腹板上第一个单元，编号为10EGEN,10,6,10 !拷贝单元10形成腹板E,71,72,82,81,171,172,18

4、2,181 !定义另一翼缘上一个单元，编号20EGEN,9,1,20 !拷贝单元20形成另一翼缘将第一层单元沿y方向拷贝60层，ANSYS自动生成所有的节点。节点编号每层加100。EGEN,60,100,1,28,1,0,（Dis_ver-H_beam/2）/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,

5、6190 !将顶层节点向上偏移梁的腹板净高度，NGEN,11,100,ALL,D_beam/10 !并生成10层单元 EGEN,11,100,1681,1708,1 NSEL,S,NODE,7101,7190 !将顶层节点向上偏移梁的翼缘厚度EGEN,2,100,1961,1988,1 将实体模型的柱向上延伸H_beam的高度，避免梁单元和实体单元在梁柱节点处切换NSEL,S,NODE,7201,7290 !生成6层单元，每层高H_beam/6NGEN,7,100,ALL,H_beam/6 EGEN,7,100,1989,2016,1 NSEL,ALL 第一根实体模型柱完成。共计：节点79层，

6、每层编号1-90， 逐层加100，顶层编号7801-7890单元78层，自动编号。每层28个，共28*78=2184个拷贝第一根柱生成第二根柱NGEN,2,10000,ALL,Dis_hor !拷贝所有的节点，节点号加10000EGEN,2,10000,1,2184,1 !拷贝所有的单元生成梁梁被夹在两根柱之间，实际长度为Dis_hor-H_col生成左侧形成一个梁截面的所有节点，梁的节点编号从20001开始N,20001,H_col/2,Dis_ver-H_beam/2,-W_beam/2 NGEN,4,1,20001,B_beam/3 N,20005,H_col/2,Dis_ver-H_b

7、eam/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 沿x方向偏移（Dis_hor-H_col）/100拷贝一层节点，节点编号加100NGEN,2,100,20001,20090,（Dis_hor-H_col）/100 生成梁的第一层截面单元两根柱单元总数为4368，故梁的单元编号从4369开始。E,20001,2

8、0002,20012,20011,20101,20102,20112,20111 !定义单元4369EGEN,7,1,4369 !拷贝单元4369形成一个翼缘E,20014,20015,20021,20020,20114,20115,20121,20120 !定义单元4376EGEN,10,6,4376 !拷贝单元4376形成腹板E,20071,20072,20082,20081,20171,20172,20182,20181 !定义单元4386EGEN,7,1,4386 !拷贝单元4386形成另一翼缘沿x方向拷贝100层生成整根梁，每层高度（Dis_hor-H_col）/100EGEN,10

9、0,100,4369,4392,1,（Dis_hor-H_col）/100 梁的实体模型完成总计：梁的节点为101层，每层编号1-88。从20001开始，逐层加100左端截面的节点为20001-20088； 右端截面的节点为30001-30088每层单元数为24个，总计24*100=2400个。单元编号为4369-6768-建立梁和柱连接处的耦合关系。CPINTF,all,0.002 FINISH /SOLU ANTYPE,TRANS !定义分析类型TUNIF,20 !定义初始温度-定义边界条件，并求解选择梁的上翼缘上表面，定义为HTUPNSEL,S,NODE,20081,30081,100N

10、SEL,A,NODE,20082,30082,100NSEL,A,NODE,20083,30083,100NSEL,A,NODE,20084,30084,100NSEL,A,NODE,20085,30085,100NSEL,A,NODE,20086,30086,100NSEL,A,NODE,20087,30087,100 CM,HTUP,NODE !定义以上所选节点为HTUP选择梁的下翼缘下表面，定义为HTDOWNNSEL,S,NODE,20001,30001,100NSEL,A,NODE,20002,30002,100NSEL,A,NODE,20003,30003,100NSEL,A,NOD

11、E,20004,30004,100NSEL,A,NODE,20005,30005,100NSEL,A,NODE,20006,30006,100NSEL,A,NODE,20007,30007,100CM,HTDOWN,NODE !定义以上所选节点为HTDOWN选择梁的左侧，定义为HTLEFTNSEL,S,NODE,20011,30011,100NSEL,A,NODE,20012,30012,100NSEL,A,NODE,20013,30013,100NSEL,A,NODE,20014,30014,100NSEL,A,NODE,20020,30020,100NSEL,A,NODE,20026,30

12、026,100NSEL,A,NODE,20032,30032,100NSEL,A,NODE,20038,30038,100NSEL,A,NODE,20044,30044,100NSEL,A,NODE,20050,30050,100NSEL,A,NODE,20056,30056,100NSEL,A,NODE,20062,30062,100NSEL,A,NODE,20068,30068,100NSEL,A,NODE,20074,30074,100NSEL,A,NODE,20071,30071,100NSEL,A,NODE,20072,30072,100NSEL,A,NODE,20073,30073

13、,100CM,HTLEFT,NODE !定义以上所选节点为HTLEFT选择梁的右侧，定义为HTRIGHTNSEL,S,NODE,20015,30015,100NSEL,A,NODE,20016,30016,100NSEL,A,NODE,20017,30017,100NSEL,A,NODE,20018,30018,100NSEL,A,NODE,20021,30021,100NSEL,A,NODE,20027,30027,100NSEL,A,NODE,20033,30033,100NSEL,A,NODE,20039,30039,100NSEL,A,NODE,20045,30045,100NSEL,

14、A,NODE,20051,30051,100NSEL,A,NODE,20057,30057,100NSEL,A,NODE,20063,30063,100NSEL,A,NODE,20069,30069,100NSEL,A,NODE,20075,30075,100NSEL,A,NODE,20076,30076,100NSEL,A,NODE,20077,30077,100NSEL,A,NODE,20078,30078,100CM,HTRIGHT,NODE !定义以上所选节点为HTRIGHT施加热边界条件并求解*DO,tm,60,180,60 !定义时间参数tm从60到600（秒）Time,tm !当

15、前时间为tmDELTIM,20 !定义初始时间步长AUTOTS,ON !打开自动步长控制Temp=20+345*LOG10（8*tm/60+1 !计算环境空气温度SF,HTUP,CONV,25,Temp !对受热边界施加对流作用SF,HTDOWN,CONV,25,TempSF,HTLEFT,CONV,25,TempSF,HTRIGHT,CONV,25,TempSF,HTUP,RDSF,0.7,1 !定义热辐射场SF,HTDOWN,RDSF,0.7,2SF,HTLEFT,RDSF,0.7,3SF,HTRIGHT,RDSF,0.7,4STEF,5.6696E-8 !定义Stefan-Boltzma

16、nn 常数TOFFST,273 !定义绝对温度偏差SPCTEMP,1,Temp !定义热辐射场的温度为火的温度SPCTEMP,2,TempSPCTEMP,3,TempSPCTEMP,4,TempSOLVE*ENDDOFINISHANSYS 热-结构耦合分析实例（2）结构分析/PREP7/TITLE,Part 2: structural analysisET,1,SOLID45,1,1 !对应于SOLID70的结构单元为SOLID45ET,2,BEAM188单元类型2为BEAM188-定义结构分析材料特性fy=275E+6常温下屈服应力exx=2.1E+11常温下杨氏模量MPTEMP清除旧的温度

17、场MPTEMP,20,100,200,300,400定义随温度变化的杨氏模量MPDATA,EX,1,exx,exx,0.9*exx,0.8*exx,0.7*exxMPTEMP,500,600,700,800,900MPDATA,EX,1,0.6*exx,0.31*exx,0.13*exx,0.09*exx,0.0675*exxMP,NUXY,1,0.3定义泊松比MP,ALPX,1,1.4E-5定义热膨胀系数定义随温度变化的应力-应变关系TB,MISO,1,10,3共10个温度，每个温度时的应力-应变由3个点描述TBTEMP,2020度时的应力-应变关系TBPT,fy/exx,fyTBPT,0.

18、02,fyTBPT,0.15,fyTBTEMP,100100度时的应力-应变关系TBTEMP,200200度时的应力-应变关系TBPT,0.807*fy/（0.9*exx）,0.807*fyTBTEMP,300300度时的应力-应变关系TBPT,0.613*fy/（0.8*exx）,0.613*fyTBTEMP,400400度时的应力-应变关系TBPT,0.420*fy/（0.7*exx）,0.420*fyTBTEMP,500500度时的应力-应变关系TBPT,0.360*fy/（0.6*exx）,0.360*fyTBPT,0.02,0.780*fyTBPT,0.15,0.780*fyTBTE

19、MP,600600度时的应力-应变关系TBPT,0.180*fy/（0.310*exx）,0.180*fyTBPT,0.02,0.470*fyTBPT,0.15,0.470*fyTBTEMP,700700度时的应力-应变关系TBPT,0.075*fy/（0.130*exx）,0.075*fyTBPT,0.02,0.230*fyTBPT,0.15,0.230*fyTBTEMP,800800度时的应力-应变关系TBPT,0.050*fy/（0.090*exx）,0.050*fyTBPT,0.02,0.110*fyTBPT,0.15,0.110*fyTBTEMP,900900度时的应力-应变关系TB

20、PT,0.0375*fy/（0.0675*exx）,0.0375*fyTBPT,0.02,0.060*fyTBPT,0.15,0.060*fy-定义梁和柱的截面特性SECTYPE,1,beam,I,column定义柱截面为截面类型1SECDATA,W_col,W_col,H_col,tf_col,tf_col,tw_col定义柱的截面尺寸SECTYPE,2,beam,I,beam定义梁截面为截面类型2SECDATA,W_beam,W_beam,H_beam,tf_beam,tf_beam,tw_beam定义梁截面尺寸用梁单元建立框架的剩余部分的模型-K,1,Dis_ver+H_beam*1.5

21、定义生成框架的关键点K,2,2*Dis_verK,3,3*Dis_verK,4,Dis_hor,Dis_ver+H_beam*1.5K,5,Dis_hor,2*Dis_verK,6,Dis_hor,3*Dis_verK,7,Dis_hor+H_col/2,Dis_verK,8,2*Dis_horK,9,2*Dis_hor,Dis_verK,10,2*Dis_hor,2*Dis_verK,11,2*Dis_hor,3*Dis_verK,12,3*Dis_horK,13,3*Dis_hor,Dis_verK,14,3*Dis_hor,2*Dis_verK,15,3*Dis_hor,3*Dis_ve

22、rK,100,-3,3定义用于确定梁的主轴方向的关键点K,200,5,20生成线L,1,2线1-10用来生成柱单元L,2,3L,4,5L,5,6L,8,9L,9,10L,10,11L,12,13L,13,14L,14,15L,2,5线11-18用来生成梁单元L,3,6L,7,9L,5,10L,6,11L,9,13L,10,14L,14,15定义线的属性LSEL,S,LINE,1,10,1定义线1-10 （柱）的属性LATT,1,2,100,1LSEL,ALLLSEL,S,LINE,11,18,1定义线11-18（梁）的属性LATT,1,2,200,2LSEL,ALL划分单元LESIZE,ALL,0.3定义单元尺寸LMESH,ALL划分单元-建立耦合与约束关系CPINTF,ALL,0.002自动耦合实体模型部分实体模型和线模型之间有三个接口：两个柱端的连接，以及底层中跨的梁左端连接到第二根实体柱的侧面建立关键点1和第一根柱柱端的连接N1=NODE（0,Dis_ver+H_beam*1.5,0）找到对应于关键点1的节点号，标记为N1num=0