ABAQUS在混凝土温度场和徐变温度应力场计算中的应用
基于ABAQUS的地下室基础底板温度场及温度应力分析_马向红

表1
材料参数
材料
密度 导热系数 比热 放热系数 kg/m3 KJ/m·h·℃ KJ/kg·℃ KJ/m2·h·℃
混凝土 土壤
2 400 1 750
2.710 0.586
0.963 1.005
34.73
江苏建筑 2011 年第 6 期(总第 145 期) 2.2 边界条件和初始条件
地下室底板与地基相连接,计算温度场时,土壤底部和
基土共划分 6 300 个单元。
图 1 有限元模型
图 2 浇 筑 后 30 d 三 维 温 度 场 云 图 (图 中 NT11 为 温 度,单位为℃)
如图 2 所 示 , 节点 3276、9828、16308 分 别 为从 上 至 下
的 3 个 点 , 坐 标 为 (10.7,400E -3,0)、 (10.7,0,0) 和 (10.7, -
混 凝 土 底 板 的 平 面 尺 寸 取 为 42.8 m×23.1 m, 厚 度 为 800 mm,浇 筑 在 土 壤 地 基 上 ,计 算 中 受 混 凝 土 影 响 的 土 壤 深度随时间增加而增加,但当其深度超过 1.5 m 后其影响已 很小, 故 取 土壤 的 计 算 深 度 为 1.5 m[5]。 混 凝 土 和 土 壤 的 热 力参数见表 1。
土壤尺寸为 21.4 m×11.55 m×1.5 m。 具体模型见图 1。 对结
构先进行热分析,得到其温度场,计算应力场时,读入热分
析结果,将其 作为 结 构 荷载 。 混 凝 土取 C40,混 凝 土弹 性 模
量随齡期的变化为:
E=(1-exp(-0.09×i))×3.25×104 N/mm2
(3)
参考文献
111
[1] 梁 绍 连 . 新 加 坡 的 住 房 保 障 制 度 [EB/OL].http://www. /p-57119048.html. [2] 常州市区公共租赁住房管理办法.(2009-07-16). [3] 张 瑞 雪.保 障 性 住 房 建 设 存 在 的 问 题 及 相 关 对 策 建 议 [J].河 南 建 材 ,2009,30(6):39-40. [4] 刘 玉 峰,陈 刚.完 善 公 租 房 制 度 的 政 策 建 议— — — 以 重 庆 市 为 例[J].建 筑 经 济,2010(07). [5] 李 宇 中 .浅 谈 利 用 REITs 发 展 公 租 房 建 设— — — 促 进 我 国 保 障 性 住 房 体 系 进 一 步 完 善 [J].现 代 商 业 ,2010(11). [6] 强 真 . 我 国 住 宅 供 应 模 式 优 化 及 价 格 形 成 机 制 的 研 究— —— 基于公共租赁房的启示[J]价格理论与实践.2009(07).
基于ABAQUS的预应力型钢混凝土结构实验温度场分析

基于ABAQUS的预应力型钢混凝土结构实验温度场分析摘要:在火灾和荷载共同作用的基础上,对预应力型钢混凝土简支梁进行耐火性能试验研究,分析了构件内部温度场分布规律。
通过试验研究,揭示了高温作用下预应力型钢混凝土简支梁力学性能劣化机理。
并利用ABAQUS大型非线性有限元分析软件,基于试验得到的升温曲线,对预应力型钢混凝土梁进行了传热分析,得到了构件内部不同时刻温度分布情况,并将测点处的温度时间关系曲线与试验实测值进行了比较。
关键词:预应力型钢混凝土,温度场,试验研究,非线性分析Abstract: On the basis of fire and load function, we have researched prestressed steel reinforced concrete and supported refractory performance test. This paper analyzes the internal temperature field distribution component. According to the test results, and reveals the effect of prestressed steel under high temperature simply-supported reinforced concrete mechanics performance degradation mechanism.By using ABAQUS large nonlinear finite element analysis software, based on the test to get heating curves of prestressed steel reinforced concrete beams are heat transfer analysis, get the internal components at different temperature distribution, and will point the temperature of the place time relationship curve and the experimental values are compared.Key Words: prestressed steel reinforced concrete, the temperature field, the experimental research, nonlinear analysis1 引言近年来,随着我国经济的迅速发展,对建筑功能的要求越来越高,建筑设计正朝着体型复杂、功能多样的综合性方向发展,相应的结构形式也变得复杂多样,大跨度、重载结构不断涌现,对结构体系提出了新的更高的要求。
ABAQUS在大体积混凝土徐变温度应力计算中的应用

ABAQUS在大体积混凝土徐变温度应力计算中的应用
王建;刘爱龙
【期刊名称】《河海大学学报(自然科学版)》
【年(卷),期】2008(036)004
【摘要】对ABAQUS进行了二次开发,建立了弹性徐变本构模型,提出了一种考虑表面散热对水管冷却影响的等效算法,并进行了数值测试和分析.结果表明,在ABAQUS中增加上述模型和算法后,可以实现大体积混凝土温度场和徐变应力的计算,拓宽ABAQUS软件在土木工程中的应用范围.
【总页数】6页(P532-537)
【作者】王建;刘爱龙
【作者单位】河海大学水利水电工程学院,江苏,南京210098;西北勘测设计研究院,陕西,西安710065
【正文语种】中文
【中图分类】TV315
【相关文献】
1.大体积混凝土温度徐变应力的ANSYS计算模块 [J], 李立峰;谢攀
2.基于Ansys平台增量法求解大体积混凝土早期弹性徐变温度应力 [J], 关战伟
3.基于Ansys平台增量法求解大体积混凝土早期弹性徐变温度应力 [J], 关战伟
4.大体积混凝土结构随机温度徐变应力计算方法研究 [J], 刘宁;刘光廷
5.普通工业玻璃退火与应力的几个问题Ⅲ在退火温度范围内连续冷却时玻璃中温差应力松弛的计算——变温下应力松弛的计算 [J], 林亢
因版权原因,仅展示原文概要,查看原文内容请购买。
基于ABAQUS的自由辊温度场及热应力场分析

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!(2)多种零件排样选择part1、part2、part3,其中,part1的长宽不变为100×70,part2的长宽为200×100,part3的长宽为150×90,板材大小仍然为设800×500,X 间距、Y 间距和零件与板材的边距都设为8,排样对话框设置和在AutoCAD 中的排样结果如图4、图5。
5结论通过排样程序可以看出,ObjectARX 编程环境提供的与MFC 相关的用户界面类AdUi 和AcUi 使开发的应用程序能与AutoCAD 风格界面很好地融合在一起,能在同一个界面对AutoCAD 和MFC 对话框进行操作,同时应用Visual.C ++设计较为复杂的对话框。
虽然用ObjectARX 开发应用程序功能强大,但掌握ObjectARX编程方法并不容易,尤其是运行到AutoCAD 界面出现的错误,因为没有具体的错误提示,很难从程序中找到错误,需要经过不停的调试才能找到原因,这给程序设计带来很大困难。
[参考文献][1]秦洪现,崔惠岚,孙剑,等.Autodesk 系列产品开发培训教程[M ].北京:化学工业出版社,2008.[2]江思敏,曹默,胡春江.AutoCAD2000开发工具———ObjectARX开发工具与应用实例[M ].北京:人民邮电出版社,1999.[3]刘蓉梅,姜秀萍,华徐勇,等.ObjectARX 二次开发及应用实例[J ].机械设计与制造,2002(3):27-29.(编辑昊天)作者简介:谢友宝(1968-),男,教授,硕士研究生导师,主要研究方向为机电一体化设备研制、数控技术、CAD/CAM 技术、计算机软硬件系统开发等。
收稿日期:2009-06-18图3排样结果图2设置对话框图5排样结果图4对话框设置基于ABAQUS 的自由辊温度场及热应力场分析杨桂芳1,罗会信1,林刚2,代宗岭2(1.武汉科技大学机械自动化学院,武汉430081;2.中冶京诚工程技术有限公司,北京100081)自由辊是连铸机中重要的零件,在结晶器、支撑导向段、扇形段中都有使用。
Abaqus焊接模拟分析程序(包括应力场和温度场)

【我的硕士论文的一部分】求解温度场!上表面上没有对流换热边界条件!单位制:米、秒、摄氏度!/CLEAR,START/FILNAME,temp,0/COM,ANSYS RELEASE 10.0 UP20050718 00:09:52 11/26/2007/CONFIG, NRES, 5000/PREP7/VIEW,1,1,2,3/ANG,1/REP,FAST!*!===============================================================================================!指定单元ET,1,SOLID70!*!*!===============================================================================================!材料属性!=====================================================================================!316LMPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,7850MPTEMP,,,,,,,,MPTEMP,1,20MPTEMP,2,300MPTEMP,3,900MPTEMP,4,1400MPTEMP,5,2000MPDATA,KXX,1,,18.6MPDATA,KXX,1,,21.4MPDATA,KXX,1,,28.4MPDATA,KXX,1,,33.9MPDATA,KXX,1,,48MPTEMP,,,,,,,,MPTEMP,1,20MPTEMP,2,600MPTEMP,3,800MPTEMP,4,1400MPTEMP,5,2000MPDATA,C,1,,502MPDATA,C,1,,612MPDATA,C,1,,635MPDATA,C,1,,659MPDATA,C,1,,670MPTEMP,,,,,,,,MPTEMP,1,20MPTEMP,2,1300MPTEMP,3,1410MPTEMP,4,1440MPTEMP,5,1550MPTEMP,6,2000MPDATA,ENTH,1,,7.88e7MPDATA,ENTH,1,,6.131e9MPDATA,ENTH,1,,7.347e9MPDATA,ENTH,1,,9.145e9MPDATA,ENTH,1,,1.03e10MPDATA,ENTH,1,,1.272e10!===============================================================================================!定义常量WidthBase=0.025 !宽度HeightBase=0.02 !基底高度Length=0.09 !长度WidthClad=0.0015 !宽度HeightDeposition=0.00375 !覆层高度Layer=15 !层数HeightClad=HeightDeposition/layerdt=0.0001 !小量Temp=20 !环境温度InitTemp=300 !初始温度CoffConv=30 !对流换热系数!===============================================================================================!定义常量Velocity=0.003 !扫描速度StepDis=0.0015 !每个载荷步位移LaserPower=700 !激光功率Radius=0.0015 !激光光斑半径Area=3.14159265*(Radius**2) !激光光斑面积Factor=0.3 !吸收因子StepTime=StepDis/Velocity !每个载荷步时间TotalTime=(Length+Radius*2)/Velocity !载荷持续时间(扫描一层) StepNum=(Length+Radius*2)/StepDis !载荷步数!===============================================================================================!建模BLOCK,0,Length,0,-0.0066,0,WidthClad,BLOCK,0,Length,0,-0.0066,WidthClad,0.0067BLOCK,0,Length,-0.0066,-HeightBase,0,WidthClad,K, ,0,0,WidthBase,K, ,Length,0,WidthBase,K, ,0,-HeightBase,WidthBase,K, ,Length,-HeightBase,WidthBase,V, 16, 13, 27, 25, 15, 14, 28, 26V, 24, 21, 27, 13, 23, 22, 28, 14BLOCK,0,Length,0,HeightDeposition,0,WidthClad,VGLUE,ALLNUMCMP,ALL!===============================================================================================!划分网格LSEL, S, LOC, Y, dt, HeightDeposition-dt, !高度方向LESIZE, ALL, , , Layer,LSEL,S,LOC,Y,-DT,-0.0066+DTLSEL,R,LOC,X,0LESIZE,ALL,,,4,2LSEL,S,LOC,Y,-DT,-0.0066+DTLSEL,R,LOC,X,LengthLESIZE,ALL,,,4,0.5LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DTLSEL,R, LOC, Z, 0, WidthBase-DT,LESIZE,ALL,,,4,2LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DTLSEL,R, LOC, Z, WidthBase,LESIZE,ALL,,,4,0.5LSEL, S, LOC, X, dt, Length-dt, !长度方向LESIZE, ALL, , , Length/StepDis,LSEL, S, LOC, Z, dt, WidthClad-dt, !宽度方向LESIZE, ALL, , , 1,LSEL,S,LOC,Z,WidthClad+DT,0.0067-DTLESIZE,ALL,,,4,LSEL,S,LOC,Z,0.0067+DT,WidthBase-DTLSEL,R,LOC,Y,-DT,-HeightBaseLESIZE,ALL,,,4,2LSEL,S,LOC,Z,0.0067+DT,WidthBase-DTLSEL,R,LOC,Y,0LESIZE,ALL,,,4,0.5VSEL,S,LOC,Y,0,HeightDeposition!网格划分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLVSEL,S,LOC,Y,-1,0 !网格划分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLALLSEL,ALL!============================================================================= ==================!基底边界条件、初始条件NSEL, S, LOC, Y, -HeightBase, 0 !基底初始温度IC,ALL,TEMP,InitTempNSEL, S, LOC, Y, -HeightBase, 0 !基底侧面,换热边界条件NSEL, R, LOC, Z, WidthBaseSF, ALL, CONV, CoffConv, TempALLSEL,ALLNSEL, S, LOC, Y, -HeightBase, 0 !基底左端面,换热边界条件NSEL, R, LOC, X, 0SF, ALL, CONV, CoffConv, TempALLSEL,ALLNSEL, S, LOC, Y, -HeightBase, 0 !基底右端面,换热边界条件NSEL, R, LOC, X, LengthSF, ALL, CONV, CoffConv, TempALLSEL,ALLNSEL, S, LOC, Y, 0 !基底上表面,换热边界条件NSEL, R, LOC, Z, WidthClad, WidthBaseSF, ALL, CONV, CoffConv, TempALLSEL,ALLFINISH/SOLU!===============================================================================================!瞬态分析参数设置ANTYPE,4 !分析类型:瞬态!*TRNOPT,FULL !求解方法:完全的N-R方法!*!DELTIM,0.01,0.001,0.05 !载荷子步(默认子步时间步长、最小、最大)——载荷步为0.333NSUBST, 4CNVTOL,HEAT, ,0.01,2,0.000001, !收敛准则:控制热流OUTRES,NSOL,LAST !结果输出:所有!===============================================================================================!杀死单元NSEL, S, LOC, Y, 0, HeightDeposition !杀死熔覆层单元ESLN, S, 1, ALLEKILL,ALLALLSEL,ALLESEL,S,LIVEEPLOTESEL,S,LIVE !激活单元的上表面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,0NSEL,R,LOC,Z,0,RadiusSF,ALL,CONV,CoffConv,TempALLSEL,ALL!===============================================================================================!预热*DO, i, 1, 2m=mod(i,2)*IF,m,EQ,1,THEN !如果为奇数层,向右扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i-1)+StepTime*k !载荷步结束时间KBC, 1 !载荷步内载荷随时间分布:常数LeftX=StepDis*(k-1)RightX=StepDis*kNSEL, S, LOC, Y, 0ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, RightX-2*Radius, RightXNSEL, R, LOC, Z, 0, Radius !激活单元的上表面,加热流密度ESLN, S, 1SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,ALLSEL, ALLSOLVESAVESFEDELE,ALL,4,HFLUX!删除热流密度载荷ALLSEL, ALLESEL,S,LIVE!激活单元的上表面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,0NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDisNSEL, R, LOC, Z, O, RadiusSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ELSE !如果为偶数层,向左扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i-1)+StepTime*k !载荷步结束时间KBC, 1 !载荷步内载荷随时间分布:常数LeftX=Length-StepDis*kRightX=Length-StepDis*(k-1)NSEL, S, LOC, Y, 0ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, RightX-2*Radius, RightXNSEL, R, LOC, Z, 0, Radius !激活单元的上表面,加热流密度ESLN, S, 1SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,ALLSEL, ALLSOLVESAVESFEDELE,ALL,4,HFLUX!激活单元的上表面,删除载荷ALLSEL, ALLESEL,S,LIVE!激活单元的上表面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,0NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDisNSEL, R, LOC, Z, O, RadiusSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ENDIF*ENDDOESEL,S,LIVEEPLOT!===============================================================================================!熔覆*DO, i, 1, Layer, 1m=mod(i,2)*IF,m,EQ,1,THEN !如果为奇数层,向右扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i+1)+StepTime*k !载荷步结束时间KBC, 1 !载荷步内载荷随时间分布:常数LeftX=StepDis*(k-1)RightX=StepDis*kNSEL, S, LOC, Y, HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALLALLSEL, ALLNSEL, S, LOC, Y, HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iNSEL, R, LOC, x, RightX-2*Radius, RightXNSEL, R, LOC, Z, 0, Radius !激活单元的上表面,加热流密度ESLN, S, 1SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,ALLSEL, ALLESEL,S,LIVE!激活单元的表面,如果包含左端面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,0NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活单元的表面,如果包含右端面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,LengthNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活单元的侧面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,Z,WidthCladNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLSOLVESAVESFEDELE,ALL,4,HFLUX!删除热流密度载荷ALLSEL, ALLESEL,S,LIVE!激活单元的上表面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDisSF,ALL,CONV,CoffConv,TempALLSEL,ALLNSEL,S,LOC,Y,HeightDeposition/Layer*(i-1) !激活单元的下表面,删除对流换热边条ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, LeftX-2*StepDis,RightX-2*StepDisNSEL, R, LOC, Z, 0, RadiusNSEL, U, LOC, Y, HeightDeposition/Layer*iESLN, S, 1NSEL, R, LOC, Y,HeightDeposition/Layer*(i-1)SFDELE, ALL, CONVALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ELSE !如果为偶数层,向左扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i+1)+StepTime*k !载荷步结束时间KBC, 1 !载荷步内载荷随时间分布:常数LeftX=Length-StepDis*kRightX=Length-StepDis*(k-1)NSEL, S, LOC, Y, HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALLALLSEL, ALLNSEL, S, LOC, Y, HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX, LeftX+2*RadiusNSEL, R, LOC, Z, 0, Radius !激活单元的上表面,加热流密度ESLN, S, 1SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,ALLSEL, ALLESEL,S,LIVE!激活单元的表面,如果包含左端面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,0NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活单元的表面,如果包含右端面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,LengthNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活单元的侧面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,Z,WidthCladNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLSOLVESAVESFEDELE,ALL,4,HFLUX!激活单元的上表面,删除载荷ALLSEL, ALLESEL,S,LIVE!激活单元的上表面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX+2*StepDis, RightX+2*StepDisSF,ALL,CONV,CoffConv,TempALLSEL,ALLNSEL,S,LOC,Y,HeightDeposition/Layer*(i-1) !激活单元的下表面,删除对流换热边条ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, LeftX+2*StepDis,RightX+2*StepDisNSEL, R, LOC, Z, 0, RadiusNSEL, U, LOC, Y, HeightDeposition/Layer*iESLN, S, 1NSEL, R, LOC, Y,HeightDeposition/Layer*(i-1)SFDELE, ALL, CONVALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ENDIF*ENDDOESEL,S,LIVEEPLOT!===============================================================================================!冷却!==============================================================================!~100s*DO, k, 1, 2, 1TIME,TotalTime*(Layer+2)+50*k!载荷步结束时间NSUBST, 5KBC, 1SOLVESAVE*ENDDO!==============================================================================!~1000s*DO, k, 1, 9, 1TIME,TotalTime*(Layer+2)+100+100*k!载荷步结束时间NSUBST, 5KBC, 1SOLVESAVE*ENDDO!==============================================================================!~3000s*DO, k, 1, 10, 1TIME,TotalTime*(Layer+2)+1000+200*k!载荷步结束时间NSUBST, 5KBC, 1SOLVESAVE*ENDDO!==============================================================================!~10000s*DO, k, 1, 14, 1TIME,TotalTime*(Layer+2)+3000+500*k!载荷步结束时间NSUBST, 5KBC, 1SOLVESAVE*ENDDO FINISH【我的硕士论文的一部分】求解应力场!修改速度、时间子步步长、载荷文件位置!如果修改基底的热膨胀系数,要修改宏文件!单位制:米、秒、摄氏度/CLEAR,START/FILNAME,stress,0/COM,ANSYS RELEASE 10.0 UP20050718 20:15:52 09/10/2007/CONFIG, NRES, 5000/PREP7/PAGE, 1000, , 1000,/VIEW,1,1,2,3/ANG,1/REP,FAST!*!===============================================================================================!指定单元ET,1,45!*!*!===============================================================================================!材料属性!=====================================================================!316LMPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,7850MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,ALPX,1,,1.75E-005MPTEMP,,,,,,,,MPTEMP,1,20MPTEMP,2,300MPTEMP,3,600MPTEMP,4,900MPTEMP,5,1300 MPDATA,EX,1,,2.0E+11MPDATA,EX,1,,1.7E+11 MPDATA,EX,1,,1.5E+11MPDATA,EX,1,,5.0E+10 MPDATA,EX,1,,0.4E+10MPDATA,PRXY,1,,0.25MPDATA,PRXY,1,,0.25MPDATA,PRXY,1,,0.25MPDATA,PRXY,1,,0.25MPDATA,PRXY,1,,0.25TB,KINH,1,5,4,0TBTEMP,20TBPT,,0,0TBPT,,7E-4,1.4E8TBPT,,0.0012,1.83E8TBPT,,0.1,2.16E9TBTEMP,300TBPT,,0,0TBPT,,5.5E-4,9.35E7TBPT,,0.0012,1.27E8TBPT,,0.1,1.84E9TBTEMP,600TBPT,,0,0TBPT,,3.2E-4,4.8E7TBPT,,0.0012,7.19E7TBPT,,0.1,1.54E9TBTEMP,900TBPT,,0,0TBPT,,2.5E-4,1.25E7TBPT,,0.0012,5.1E7TBPT,,0.1,5.45E8TBTEMP,1300TBPT,,0,0TBPT,,2.5E-4,1E6TBPT,,0.00375,1.13E7TBPT,,0.1,7.05E7!=====================================================================!A3ExpandCoeff=1.75E-005!structural->nonlinear->inelastic->rate independent->kinematic hardeningplasticity->mises plasticity->bilinear!===============================================================================================!定义常量WidthBase=0.025 !宽度HeightBase=0.02 !基底高度Length=0.09 !长度WidthClad=0.0015 !宽度HeightDeposition=0.00375 !覆层高度Layer=15 !层数HeightClad=HeightDeposition/layerdt=0.0001 !小量Temp=20 !环境温度InitTemp=300 !初始温度CoffConv=30 !对流换热系数!===============================================================================================!定义常量Velocity=0.003 !扫描速度StepDis=0.0015 !每个载荷步位移LaserPower=700 !激光功率Radius=0.0015 !激光光斑半径Area=3.14159265*(Radius**2) !激光光斑面积Factor=0.3 !吸收因子StepTime=StepDis/Velocity !每个载荷步时间TotalTime=(Length+Radius*2)/Velocity !载荷持续时间(扫描一层) StepNum=(Length+Radius*2)/StepDis !载荷步数!===============================================================================================!建模BLOCK,0,Length,0,-0.0066,0,WidthClad,BLOCK,0,Length,0,-0.0066,WidthClad,0.0067BLOCK,0,Length,-0.0066,-HeightBase,0,WidthClad,K, ,0,0,WidthBase,K, ,Length,0,WidthBase,K, ,0,-HeightBase,WidthBase,K, ,Length,-HeightBase,WidthBase,V, 16, 13, 27, 25, 15, 14, 28, 26V, 24, 21, 27, 13, 23, 22, 28, 14BLOCK,0,Length,0,HeightDeposition,0,WidthClad,VGLUE,ALLNUMCMP,ALL!===============================================================================================!划分网格LSEL, S, LOC, Y, dt, HeightDeposition-dt, !高度方向LESIZE, ALL, , , Layer,LSEL,S,LOC,Y,-DT,-0.0066+DTLSEL,R,LOC,X,0LESIZE,ALL,,,4,2LSEL,S,LOC,Y,-DT,-0.0066+DTLSEL,R,LOC,X,LengthLESIZE,ALL,,,4,0.5LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DTLSEL,R, LOC, Z, 0, WidthBase-DT,LESIZE,ALL,,,4,2LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DTLSEL,R, LOC, Z, WidthBase,LESIZE,ALL,,,4,0.5LSEL, S, LOC, X, dt, Length-dt, !长度方向LESIZE, ALL, , , Length/StepDis,LSEL, S, LOC, Z, dt, WidthClad-dt, !宽度方向LESIZE, ALL, , , 1,LSEL,S,LOC,Z,WidthClad+DT,0.0067-DTLESIZE,ALL,,,4,LSEL,S,LOC,Z,0.0067+DT,WidthBase-DTLSEL,R,LOC,Y,-DT,-HeightBaseLESIZE,ALL,,,4,2LSEL,S,LOC,Z,0.0067+DT,WidthBase-DTLSEL,R,LOC,Y,0LESIZE,ALL,,,4,0.5VSEL,S,LOC,Y,0,HeightDeposition!网格划分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLVSEL,S,LOC,Y,-1,0 !网格划分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLALLSEL,ALLVSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DT!删除热影响区外的单元VSEL,A,LOC,Z,0.0067+DT,WidthBase-DTVCLEAR,ALLVDELE,ALL, , ,1ALLSEL,ALLFINISH!===============================================================================================!瞬态分析参数设置/SOLANTYPE,4 !分析类型:瞬态TRNOPT,FULL !求解方法:对于材料非线性,这是唯一的方法NLGEOM,on !大变形分析LNSRCH, onNSUBST, 4NEQIT,30CNVTOL,U,,0.05,2,, !收敛准则:控制位移CNVTOL,F,,0.01,2, !收敛准则:控制力OUTRES,NSOL,LAST !结果输出:TREF, Temp!===============================================================================================!杀死单元NSEL, S, LOC, Y, 0, HeightDeposition !杀死熔覆层单元ESLN, S, 1EKILL,ALLALLSEL,ALLESEL,S,LIVEEPLOTNSEL,S,LOC,Z,0 !对称边条(相当于三个约束,UZ,ROTX,ROTY)D,ALL,UZ,0NSEL,S,LOC,Y,-0.0066 !固定中心点(增加两个约束,UX,UY)NSEL,R,LOC,Z,0NSEL,R,LOC,X,length/2D,ALL,ALL,0NSEL,S,LOC,Y,-0.0066,0 !固定中心线(增加一个约束,ROTZ)NSEL,R,LOC,Z,0NSEL,R,LOC,X,length/2D,ALL,UX,0ALLSEL,ALLSAVE!===============================================================================================!熔覆*DO,m,1,Layerk=mod(m,2)*IF,K,EQ,1,THEN*DO,n,1,StepNum,KBC,0TIME,TotalTime*(m-1)+StepTime*nLeftX=StepDis*(n-1)RightX=StepDis*nNSEL, S, LOC, Y, HeightDeposition/Layer*(m-1),HeightDeposition/Layer*m !熔覆层生长NSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALLALLSEL, ALLLDREAD,TEMP,StepNum*(m+1)+n,last,, ,'temp','rth','F:\temp\differentvelocity\3' !读取体载荷NSEL,S,BF,TEMP,1300,3000BF, ALL, TEMP, 1300ESEL,S,LIVE!显示生单元EPLOTALLSEL,ALLMyDBC!宏命令ALLSEL,ALLSOLVESAVE*ENDDO*ELSEIF,K,EQ,0,THEN*DO,n,1,StepNum,KBC,0TIME,TotalTime*(m-1)+StepTime*nLeftX=Length-StepDis*(n-1)RightX=Length-StepDis*nNSEL, S, LOC, Y, HeightDeposition/Layer*(m-1),HeightDeposition/Layer*mNSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALLALLSEL, ALLLDREAD,TEMP,StepNum*(m+1)+n,last,, ,'temp','rth','F:\temp\differentvelocity\3'NSEL,S,BF,TEMP,1300,3000BF, ALL, TEMP, 1300ESEL,S,LIVE!显示生单元EPLOTALLSEL,ALLMyDBC!宏命令ALLSEL,ALLSOLVESAVE*ENDDO*ENDIF*ENDDO!===============================================================================================!冷却*DO, k, 1, 35, 1TIME,TotalTime*Layer+0.5*k!载荷步结束时间NSUBST, 4KBC, 0 !载荷步内载荷随时间分布:常数LDREAD,TEMP,StepNum*(Layer+2)+K,last,, ,'temp','rth','F:\temp\different velocity\3'NSEL,S,BF,TEMP,1300,3000BF, ALL, TEMP, 1300ESEL,S,LIVE!显示生单元EPLOTALLSEL,ALLMyDBC!宏命令ALLSEL,ALLSOLVESAVE*ENDDOESEL,S,LIVEEPLOTFINISH【补充】其中的宏命令是加比较复杂的位移边条,可以去掉,然后加上你需要的边界条件即可;。
热处理过程流场-温度场-组织场-应力场耦合模拟研究

热处理过程流场-温度场-组织场-应力场耦合模拟研究摘要:本文研究了热处理过程中流场、温度场、组织场和应力场之间的耦合关系。
采用了ANSYS Fluent和ABAQUS有限元软件对工件进行了流场和热力学分析,并利用ABAQUS进行了热-组织-应力场耦合分析。
通过分析结果,得出了热处理工艺参数对工件性能的影响规律,为工件热处理过程中的优化设计提供了理论依据。
关键词:热处理;流场;温度场;组织场;应力场;耦合分析正文:1. 引言热处理是一种常用的工艺方法,它通过在一定的温度条件下改变材料的组织结构,从而使材料的力学性能得到提高。
在热处理过程中,流场、温度场、组织场和应力场之间存在着复杂的耦合关系,如何对这些场进行耦合分析,是优化热处理工艺设计的关键之一。
本文采用了ANSYS Fluent和ABAQUS有限元软件对工件进行了流场和热力学分析,并利用ABAQUS进行了热-组织-应力场耦合分析。
通过分析结果,得出了热处理工艺参数对工件性能的影响规律,为工件热处理过程中的优化设计提供了理论依据。
2. 流场分析热处理过程中,加热炉内气流的速度分布对工件表面的热传递有着重要的影响。
本文采用ANSYS Fluent软件对加热炉内气流进行了数值模拟,得到了炉内气流的速度场分布图。
图1为炉内气流的速度场分布图。
(插入图1)从图1可以看出,炉内气流的速度呈现出较大的不均匀性,气流速度较高的区域主要集中在加热炉内部的两侧,而中央区域的气流速度较低。
3. 温度场分析在热处理过程中,工件表面的温度分布对工件组织结构的形成以及机械性能的提高都有着重要的影响。
本文采用ANSYS Fluent软件对工件表面温度进行了数值模拟,得到了加热炉内的温度场分布图。
图2为炉内温度场分布图。
(插入图2)从图2可以看出,工件表面的温度分布呈现出明显的不均匀性。
整个工件表面温度的分布范围较大,在工件的上下部位温度较高,在中央区域温度较低。
4. 组织场分析在热处理过程中,材料的组织结构是影响材料力学性能的重要因素之一。
ABAQUS有限元软件温度应力模拟及其在机场运营中的应用前景

ABAQUS有限元软件温度应力模拟及其在机场运营中的应用前景发布时间:2021-06-29T11:00:26.487Z 来源:《基层建设》2021年第5期作者:王旭刘景文邢俊[导读] 摘要:温度应力的长期作用可造成水泥混凝土道面结构强度和刚度的下降,严重时会导致结构开裂,因此温度应力对水泥混凝土道面承载能力和使用寿命的影响不容忽视。
中铁北京工程局集团有限公司机场工程分公司北京 102300摘要:温度应力的长期作用可造成水泥混凝土道面结构强度和刚度的下降,严重时会导致结构开裂,因此温度应力对水泥混凝土道面承载能力和使用寿命的影响不容忽视。
目前,对混凝土道面温度应力的研究较多,理论计算已较为成熟和深入。
但由于计算基本假定与实际情况存在一定的偏差,因而开展现场道面板温度变形试验研究是极为必要的。
本人于2016年,在学校导师的支持下,选取天津滨海国际机场机坪区道面板为试验研究对象,应用光纤光栅传感器测量路面及路基内部的温度场与温度应变。
采用ABAQUS有限元软件对不同温度条件下道面结构温度变形响应开展仿真计算,并与实测结果进行对比,以获得有价值的分析结论。
研究结果表明:(1)试验实测的路面结构温度应变具有明显的滞回变化特征,单日最大变形水平达到120个微应变,滞回曲线规律与仿真结果较为一致;(2)对比板内相邻测点结果,距离板边接缝位置越远应变幅度越大,且在板跨中部达到最大值,两侧呈现近似对称分布特征;(3)采用二维平面分析模型可较好的模拟实际路面及路基内部的变温过程。
所得温度应力结果与理论计算有较好的一致性,可用于工程实际分析。
关键词:机场场道;ABAQUS有限元;温度应力;航空运输量管控引言水泥混凝土具有刚度大、抗压能力强;抗弯拉强度高、使用寿命长;耐久性、耐磨性较好,可不设路缘石;耐腐蚀性强,施工机械化程度较高等优点。
目前已被绝大多数高等级公路,机场所采用。
温度应力是影响混凝土道面使用寿命的重要参数之一,尤其是机场跑道板厚设计中并没有考虑到温度应力的影响,在我国北方地区,昼夜温差较大,温度应力很可能会产生较大影响。
利用ABAQUS分析混凝土构件的温度场

・
8 ・ 8
第3 3卷 第 3 1期 2 7年 11 0 0 月
山 西 建 筑
SHANXI ARCHI TECTURE
Vo . 3No. 1 13 3 NO . 2 o V o7
文 章 编 号 :096 2 (0 7 3 .0 80 10 .8 5 2 0 )10 8 —2
pa 、 c T I 十 T / \ T. /
的导热 系数 , /m・ ; 为火灾燃烧 时间 ,。 w ( K) t S
【 , l J
32个节点 , 图 1 5 见 。
其 中 , 为材料 的密度 ,g m3C为 比热 ,/ ・ ; 为材料 』 D k/ ; J ( K)
在火 灾情况下 , 件受 火面一般 同时存在 对流和辐射两种热 构 交换方 式 , : 即
( z+
…
十 十 一 ÷斗 H斗 H+n 十 ~ J } 卜
:j 十‘ : {H }菇: :㈠一# = { = f 簿 i H{ 量
噩潞 i 阵 ]W - j
’ _
一
1 计 算原理 与 建模
1 1 计算原 理 .
T=T +7 0 1 - ・53 ) 7 4 4 o 5 ( 一e 3 957 +1 0 1 7 7 /
.
() 3
其 中, T为炉 内气 体的温度 , ; o ℃ T 为试验开 始时炉 内温度 ; 对于梁柱构件来 说 , 般假 设沿构 件长 度方 向温度 相 同, 一 则 t 时间, 。 为 h 可简化 为二维温度场 。其热传导方程为 _ 3: 3 _ 由于对称性 , 选取 14建模 。共划 分为 3 5个 E-D / 1 L 4单元 , 2
∞2+0 e 2 60+0 0e 2 2 e十0 93 2 94 0 e+0 2 51 e+0 5 2 1 e+0 26 2 73 e+0 6 2
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
式。本文采用后者实现水化热的赋值。
3. ABAQUS 中徐变应力的计算
实验资料表明,当应力保持为常量时,混凝土的应变将随着时间而有所增加,这种现象 称为混凝土的徐变。在大体积混凝土的温度应力计算过程中,都要考虑徐变应力的计算。 ABAQUS 没有直接提供混凝土徐变应力计算这一功能,但是可以通过子程序 UMAT 来实现
ABAQUS 是功能强大的有限元软件,可以分析复杂的固体力学和结构力学系统,模拟 非常庞大复杂的模型,处理高度非线性的问题[1]。
为了了解 ABAQUS 在计算温度场和温度应力这一领域的精度,本文采用该有限元软件 对某水闸加固工程中的新浇筑底板进行了温度场与温度应力的计算,并且与工程现场的一些 实测数据进行了验证比较。
2. ABAQUS 对水泥水化热的模拟
水泥水化过程中放出的热量称为水泥的水化热。水泥的水化热是影响混凝土温度应力的
一个重要因素,水泥的水化热可以用以下三种表达式:
指数式:
Q(τ ) = Q0 (1 − e−mτ )
(1)
双曲线式:
Q(τ ) = Q0τ /(n + τ )
(2)
复合指数式:
Q(τ ) = 引言
在大体积混凝土结构中,温度变化不但可能引起裂缝,对结构的应力状态也具有重要影 响,有时温度应力在数值上可能超过其他的外荷载引起的应力。温度应力分析,温度控制和 防止裂缝的措施,是大体积混凝土结构设计与施工中十分重要的课题。对于实际工程中的温 度场和温度应力的求解一般采用有限元法来进行求解。
ABAQUS 在混凝土温度场和徐变温度应力场计算中 的应用
郑海远
河海大学工程力学系,南京(210098)
E-mail:zhenghaiyuan0412@
摘 要:本文采用了有限元软件 ABAQUS 对某一水闸加固工程中新浇筑底板的非稳定温度 场与徐变温度应力进行了计算,分析过程中采用了 ABAQUS 提供的用户子程序 DFLUX 和 UMAT。将计算结果与现场实测的温度场和温度应力的结果进行了比较后表明 ABAQUS 的 计算结果的精度较为理想,可以在混凝土温度场和温度应力等问题的计算中使用该软件。 关键词:温度场,徐变温度应力,二次开发,有限元
(3)
式中:Q(τ ) 为在龄期τ 时的累积水化热, kJ/kg;Q0 为τ → ∞ 时的最终水化热,kJ/kg; τ 为龄期,d;m、n、a、b等为常数[2]。
本文通过 ABAQUS 提供的 Body heat flux 对混凝土水化热进行赋值。Body heat flux 功
能中提供了两种赋值方式,一种为常量赋值,另外一种为根据户自定义的子程序 DFLUX 方
-1-
这一功能。
计算时刻 t 的徐变温度应力可按照下式计算:
n
∑ σ *(t) = �������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������������