收藏本站
《山东大学》 2005年
收藏 | 手机打开
二维码
手机客户端打开本文

淬火过程有限元模拟关键技术及工艺参数优化的研究

李辉平  
【摘要】:针对淬火过程温度和组织变化的特点,提出了处理相变潜热的计算方法。使用集中热容矩阵、网格细化和动态时间步长的方法,编写了淬火过程有限元模拟程序,并提出了一种新的计算时间步长的方法—基于最大和最小温差的自适应步长,该方法是用前一步模拟得到的温度场与当前步模拟所得的温度场的差值控制模拟时间步长。选择有精确解或解析解的算例,用其结果与本文中有限元模拟程序求得的结果进行比较,从模拟结果的对比中可以看出,本文提出的方法可以有效地避免数据振荡,而且模拟精度较高。 自行设计了新的端部淬火实验工装和基于ISA总线的数据高速采集系统,对P20试样进行端部淬火实验,并对试样相应部分进行定量金相分析及硬度检验,得到了准确的实验数值。建立了P20钢端部淬火的有限元模型,用本文编写的程序对P20钢端部淬火实验进行有限元模拟,得到相应实验位置的各新生成相的含量及相应位置的硬度计算值,将这些数据与实验数据进行比较,检验编写的有限元模拟程序进行检验,验证了本文的有限元程序是准确性和可靠性。 提出了一种新的求解随温度变化的换热系数方法,该方法把有限元方法引入到反向热传导问题中,结合使用最优化方法中的进退法和试探法。在计算过程中,利用有限元法可以方便地计算出各个单元在整个过程的相变情况,并能得到各单元在相应时间段所产生相变潜热,从而在计算随温度变化的换热系数时,将各单元的相变潜热与单元温度场进行耦合计算。为了使进退法适用于该类反传热问题,对其算法进行了改进,并用其确定优化的搜索区间,然后用试探法(黄金分割法)在搜索区间内找到换热系数的最佳值。建立了淬火过程求解换热系数的有限元模型,利用有限元法编写了淬火过程换热系数的求解程序。利用文献中的换热系数数据和热物性参数,用大型有限元分析软件MARC及二次开发的计算相变量的用户子程序,计算得到淬火零件相应位置的温度变化曲线,然后根据温度变化曲线,用本文开发的基于有限元和最优化方法的反传热分析程序求解随温度变化的换热系数,将求得的换热系数与文献中的换热系数比较,结果表明:用本文中的方法及有限元求解程序求解的换热系数精度较好,而且计算的迭代收敛速度较快。利用自行设计基于ISA总线的数据高速采集系统采集了P20钢在20℃和60℃水
【学位授予单位】:山东大学
【学位级别】:博士
【学位授予年份】:2005
【分类号】:TG156.3

手机知网App
【引证文献】
中国期刊全文数据库 前2条
1 李辉平;赵国群;栾贻国;牛山廷;;高压气体淬火工艺有限元模拟及其工艺参数评估[J];金属热处理;2006年S1期
2 李辉平;赵国群;黄传真;牛山廷;;气体淬火工艺参数多目标优化方法[J];机械工程学报;2007年01期
中国博士学位论文全文数据库 前1条
1 贺连芳;硼钢B1500HS的热冲压关键参数测试及其淬火性能研究[D];山东大学;2012年
中国硕士学位论文全文数据库 前8条
1 韩涛;高铬锻钢支承辊淬冷过程有限元分析[D];燕山大学;2011年
2 金珍;等离子弧表面淬火过程的温度场数值模拟[D];中国科学技术大学;2011年
3 商松;Cr5锻钢支承辊喷雾淬火过程的数值模拟与工艺优化[D];燕山大学;2012年
4 巩再艳;轴类高频淬火机床设计及其关键技术研究[D];辽宁工业大学;2012年
5 张茜;磨削淬硬层残余应力的有限元分析[D];山东大学;2008年
6 唐焱;数控火焰淬火机控制系统的设计与应用研究[D];太原理工大学;2012年
7 谷磊;履带板在连续热处理生产线中的变形控制研究[D];机械科学研究总院;2012年
8 张时锋;喷雾淬火设备冷却特性的研究[D];昆明理工大学;2012年
【参考文献】
中国期刊全文数据库 前10条
1 吴景之;;大锻件加热[J];大型铸锻件;1992年03期
2 吴景之;;温度场计算机模拟在国外大锻件生产中的应用[J];大型铸锻件;1993年02期
3 叶健松,李勇军,潘健生,胡明娟;大型支承辊热处理过程的数值模拟[J];机械工程材料;2002年06期
4 程眉,许小忠,程军,原少勤,刘国柱;淬火工件/介质边界换热条件研究[J];华北工学院学报;1997年02期
5 张丽萍;均匀设计和最优化方法在热处理中的应用[J];航天工艺;1994年04期
6 管一弘,陈铁力,陈君若,王洪纲;激光淬火温度场及材料性能的数值模拟[J];中国激光;1999年03期
7 邹壮辉;高守义;张立文;林立;;淬冷热应力场研究和应用进展[J];金属热处理学报;1993年03期
8 李勇军,潘健生,胡明娟,田东,张伟民;淬冷过程非线性有限元分析的单元选择[J];金属热处理学报;1999年04期
9 姚善长;T.Ericsson;;淬火过程的计算机模拟[J];金属热处理;1987年08期
10 李勇军,潘健生,顾剑锋,胡明娟,张幸,于文平;70Cr3Mo钢大型支承辊淬火加热计算机模拟[J];金属热处理;2000年09期
【共引文献】
中国期刊全文数据库 前10条
1 崔晓龙,万妮丽;大型锻件热处理过程的数值模拟研究[J];鞍钢技术;2004年06期
2 孙慧元;丁宁;;有限元模拟技术在工字钢控冷过程中的应用[J];鞍钢技术;2010年01期
3 宋立国,胡林,邢桂菊,张玉华;圆钢控轧时温度场分析[J];鞍山钢铁学院学报;2002年02期
4 艾家和,董建新,胡尧和,谢锡善;4Cr13线材轧制后冷却速率对表面宏观应力和微观组织影响的研究[J];兵器材料科学与工程;2001年03期
5 胡锐,朱冠勇,毕晓勤,李金山;液相浸渗用多孔预制体制备的研究现状[J];兵器材料科学与工程;2004年05期
6 雷声;赵雅丽;桂中祥;陈从升;吴跃波;;轴承滚道激光淬火的温度场数值模拟[J];兵器材料科学与工程;2010年01期
7 钤立妙,李强;柴油机曲轴淬火温度场模拟计算及淬硬层深度预测[J];北方交通大学学报;2002年04期
8 付宇明;郑丽娟;任兴利;;含半埋藏共线双裂纹的Cr12合金钢试件电磁热止裂分析[J];兵工学报;2008年08期
9 汪苏,陈佳,李建明;复杂陶瓷件烧结过程温度场数值模拟及仿真[J];北京航空航天大学学报;2004年08期
10 王军,贾成厂,郭志猛,殷声;逐步熔融凝固启动阶段电磁场和温度场的数值计算[J];北京科技大学学报;2000年04期
中国重要会议论文全文数据库 前10条
1 徐小青;韩庆;周石光;;新型热轧带钢层流冷却工艺控制系统[A];全国冶金自动化信息网2012年年会论文集[C];2012年
2 胡平;马宁;郭威;申国哲;武文华;;超高强度汽车结构件热冲压技术研究进展[A];2010’力学与工程应用学术研讨会论文集[C];2010年
3 宇慧平;隋允康;张轩;龙连春;;响应面方法的改进[A];北京力学会第13届学术年会论文集[C];2007年
4 陈明;龙连春;;激光半径对金属材料热场及应力场的影响[A];北京力学会第15届学术年会论文摘要集[C];2009年
5 ;Development and Application of Computer Simulation in Heat Treatment Process[A];2005年全国计算材料、模拟与图像分析学术会议论文集[C];2005年
6 胡跃新;;某型活塞式航空发动机气缸筒加工工艺研究[A];2010年“航空航天先进制造技术”学术交流论文集[C];2010年
7 何静波;黄因慧;田宗军;刘志东;王玮琳;刘朝晖;;激光重熔过程流场及温度场的数值模拟研究现状和发展[A];2005年中国机械工程学会年会论文集第11届全国特种加工学术会议专辑[C];2005年
8 陈炜;王冬颖;郝建国;;基于Ansys Workbench离心压缩机合金结构钢法兰淬火过程温度场的数值模拟[A];第八届沈阳科学学术年会论文集[C];2011年
9 李晓涛;陈国需;郭小川;;界面摩擦渗透的数学模型[A];中国石化协会、润滑脂专业委员会一届二次全体成员大会暨全国第十五届润滑脂技术交流会论文集[C];2009年
10 孙会魁;饶良斌;张金柱;;数值模拟在金属材料中的应用[A];2008年全国冶金物理化学学术会议专辑(上册)[C];2008年
中国博士学位论文全文数据库 前10条
1 杨秀芝;低合金高强钢双丝焊热过程数值模拟及接头组织转变行为研究[D];华中科技大学;2010年
2 班丽丽;中碳Si-Mn系高强度TRIP钢高周疲劳破坏行为研究[D];昆明理工大学;2008年
3 张洪兵;结合电弧喷涂与快速原形技术快速制造钢基模具的研究[D];沈阳工业大学;2010年
4 马宁;高强度钢板热成形技术若干研究[D];大连理工大学;2011年
5 张恒良;动力机械应用功能梯度材料寿命评估与优化研究[D];华中科技大学;2011年
6 廖凯;铝合金厚板淬火—预拉伸内应力形成机理及其测试方法研究[D];中南大学;2010年
7 李立英;ASTM4130钢焊接热影响区组织与性能研究[D];中国石油大学;2011年
8 丁敏;10Ni5CrMoV钢旋转电弧NG-GMAW稳定性与物理冶金研究[D];上海交通大学;2011年
9 安俊超;多因素下乙烯裂解炉管焊接接头寿命预测研究[D];天津大学;2011年
10 张建军;07MnNiCrMoVDR钢焊接性的物理和数值模拟[D];天津大学;2011年
中国硕士学位论文全文数据库 前10条
1 张玉莲;带参非线优化问题的上下界估计[D];郑州大学;2010年
2 李勇;热处理对模具钢组织性能的影响[D];郑州大学;2010年
3 冯媛媛;热处理零件组织性能分布的计算机预测及其结构设计[D];山东农业大学;2010年
4 纪文文;超重力燃烧合成YAG陶瓷热力学计算及过程仿真[D];山东农业大学;2010年
5 何烜坤;奥氏体化温度对Fe-0.88C-1.45Si-0.95Cr-0.37Mn钢贝氏体相变的影响[D];天津理工大学;2010年
6 武淑珍;水溶液在热钢材表面运动状态的分子动力学研究[D];昆明理工大学;2009年
7 孙建阁;PQ600低合金高强钢焊接HAZ组织与性能研究[D];昆明理工大学;2008年
8 王吉满;低碳贝氏体钢组织及强韧性机制研究[D];昆明理工大学;2009年
9 张文东;H型钢热轧温度场及轧制力模拟[D];昆明理工大学;2010年
10 田艳青;中密度纤维板热压控制参数优化系统的研究[D];东北林业大学;2011年
【同被引文献】
中国期刊全文数据库 前10条
1 邹有武;;大型锻件的热处理[J];鞍钢技术;1993年10期
2 曹忠孝;李文斌;庞兆夫;李桂艳;李阳;;微量硼对低碳贝氏体钢过冷奥氏体转变的影响[J];鞍钢技术;2009年01期
3 田欣利,徐燕申,彭泽民,林彬;结构陶瓷磨削表面残余应力的产生机理[J];兵工学报;1998年04期
4 秦子然;刘国权;宋月鹏;尹江辉;;45~#钢零件感应热处理过程温度场的有限元模拟及组织研究[J];北京科技大学学报;2006年07期
5 张少军;边智慧;杨春彦;单艳华;;喷雾冷却中射流出口高度和高温钢板表面温度对换热系数的影响[J];北京科技大学学报;2007年01期
6 刘文今 ,钟敏霖;钢和铸铁激光淬硬的工艺原理探讨[J];中国表面工程;2001年01期
7 黄晨;5CrMnMo钢冲剪机上模的等离子弧表面硬化工艺研究[J];表面技术;1999年06期
8 武淑珍;陈敬超;冯晶;杜晔平;洪振军;;钢铁淬火过程模拟的研究进展[J];材料导报;2008年06期
9 邢忠文;包军;杨玉英;钟朝廷;;可淬火硼钢板热冲压成形实验研究[J];材料科学与工艺;2008年02期
10 黄鹏;魏兴钊;;用三维有限元求解金属淬火过程的换热系数[J];材料科学与工艺;2008年02期
中国重要会议论文全文数据库 前1条
1 陈乃录;张伟民;;数字化淬火冷却控制技术的应用[A];第九次全国热处理大会论文集(一)[C];2007年
中国博士学位论文全文数据库 前9条
1 李强;淬火过程的计算机模拟与试验研究[D];燕山大学;2003年
2 吕成;数值模拟技术在燃气轮机零部件锻造及热处理过程中的应用[D];大连理工大学;2007年
3 胡庆贤;穿孔等离子弧焊接温度场的有限元分析[D];山东大学;2007年
4 牛山廷;淬火冷却过程三维有限元模拟及工艺参数优化的研究[D];山东大学;2007年
5 彭勇刚;模糊控制工程应用若干问题研究[D];浙江大学;2008年
6 王婧;流场计算机模拟在热处理炉设计中的应用[D];上海交通大学;2008年
7 郭永献;喷雾液膜流动理论及电子器件喷雾冷却实验研究[D];西安电子科技大学;2009年
8 左训伟;水-空交替控时淬火冷却设备的研究与应用[D];上海交通大学;2010年
9 王亚青;喷雾冷却无沸腾区换热特性研究[D];中国科学技术大学;2010年
中国硕士学位论文全文数据库 前10条
1 郑国昌;23MnB履带板热处理质量控制的研究[D];江西理工大学;2010年
2 尚建军;变速箱齿轮热处理的变形与控制研究[D];吉林大学;2011年
3 王劲;高强精轧螺纹钢筋连续感应热处理及控制变形研究[D];机械科学研究总院;2011年
4 曾志;高效的喷雾冷却技术的实验与数值研究[D];哈尔滨工业大学;2010年
5 王晓燕;低压转子淬火工艺过程数值模拟的研究[D];上海交通大学;2011年
6 李贤君;Φ1700×7000mm井式炉的温度场模拟研究和优化设计[D];机械科学研究院;2002年
7 许静;感应加热电源数字控制技术的研究[D];浙江大学;2002年
8 钱滨;45钢热处理过程组织转变量预测的试验研究[D];清华大学;2002年
9 万爱霞;高线穿水冷却过程温度场的数值模拟[D];燕山大学;2004年
10 王同珍;淬火过程的计算机模拟[D];燕山大学;2004年
【二级引证文献】
中国期刊全文数据库 前5条
1 刘明贺;韦建很;张修铭;修世超;;小切深条件下磨削加工表面的温度场及表面完整性[J];东北大学学报(自然科学版);2012年10期
2 范俊锴;杜凤山;刘才;黄华贵;;大锻件热处理残余应力及其对微孔隙周围氢聚集的影响[J];燕山大学学报;2013年04期
3 李辉平;贺连芳;赵国群;;硼钢B1500HS界面传热系数与压力关系的研究[J];机械工程学报;2013年16期
4 夏永辉;陈岩;顾敏;苏兴武;宋志强;;计算流体力学及其在热处理中的应用研究[J];热处理;2010年01期
5 张雪姣;侯云成;张作贵;殷福星;;支承辊热处理残余应力研究现状[J];一重技术;2013年02期
中国博士学位论文全文数据库 前3条
1 张莹;单程平面磨削淬硬工件变形及淬硬层深度均匀性研究[D];山东大学;2013年
2 张磊;高强度钢板热冲压过程的模具温度控制与数值模拟技术研究[D];山东大学;2013年
3 盈亮;高强度钢热冲压关键工艺试验研究与应用[D];大连理工大学;2013年
中国硕士学位论文全文数据库 前8条
1 庄杰真;磨削淬硬层及其均匀性研究[D];集美大学;2011年
2 张汝波;齿条成形磨削淬硬耦合应力场研究[D];山东大学;2011年
3 龚玉霞;曲轴电磁感应淬火残余应力模拟及其对圆角应力影响的研究[D];中北大学;2012年
4 张振果;磨削淬硬温度场分析及其工艺参数确定方法研究[D];山东大学;2010年
5 王伟;ADI热处理工艺参数对其力学性能影响机理的研究[D];苏州大学;2010年
6 孟召焕;高强度热冲压钢板强韧性工艺优化及试验分析[D];大连理工大学;2012年
7 刘祥;大型半开放式铜模压铸的温度参数优化[D];昆明理工大学;2012年
8 刘冰琪;高强钢板热成形局部硬化的数值模拟与实验研究[D];吉林大学;2014年
【二级参考文献】
中国期刊全文数据库 前8条
1 郭可讱,金俊泽,高钦,黄琪瑞,徐安顺,王卫东;大型铸件凝固进程的数字模拟[J];大连工学院学报;1980年02期
2 雷晓燕,杜庆华;三维非线性有限元与弹性边界元耦合数值方法[J];固体力学学报;1991年01期
3 陈红桔,刘清友;模具钢P20的动态断裂韧性[J];钢铁研究学报;1998年02期
4 胡明娟,潘健生,李兵,田东;界面条件剧变的淬火过程三维温度场的计算机模拟[J];金属热处理学报;1996年S1期
5 李廷威;;渗氮变形规律板块模型的研究[J];金属热处理;1989年11期
6 潘健生,胡明娟,张伟民,阮冬,田东;我国热处理计算机模拟研究与应用的进展[J];热处理;1999年01期
7 田东,胡明娟,潘健生;T8 钢淬火过程三维温度场计算及实验[J];上海交通大学学报;1998年02期
8 顾剑锋,潘健生,胡明娟;淬火介质换热系数的计算机测算[J];热加工工艺;1998年05期
【相似文献】
中国期刊全文数据库 前10条
1 李鹏志;齐红梅;;45号钢轴类零件淬火裂纹形成原因及预防措施[J];河北工业科技;2008年03期
2 孙朝阳;方刚;雷丽萍;曾攀;;温度和相变影响的弹塑性本构关系及应用[J];材料热处理学报;2008年01期
3 蔡长景;;防止床身导轨淬火时的拉弧现象[J];金属热处理;1988年02期
4 安敬竹;;在淬火过程中回火对高速钢性能的影响[J];热处理技术与装备;1990年04期
5 朱守星;钢板淬火过程中表面传热系数的确定(一)[J];国外金属热处理;2001年01期
6 原思宇,张立文,张国梁,江国栋,裴继斌;大型锻件淬冷过程数值模拟与实验验证[J];大连理工大学学报;2005年04期
7 陈洪;刘勇;杨贤镛;;森吉米尔中间辊淬火温度场的有限元模拟[J];材料热处理学报;2006年04期
8 魏启武;宛农;周悠;;45Cr4NiMoV支承辊淬冷温度场仿真与淬硬深度预测[J];轻工机械;2007年03期
9 韩斌慧;阎献国;;基于ANSYS的丝锥淬火冷却过程计算机模拟[J];制造技术与机床;2008年06期
10 李海升;杨川;高国庆;郑周;;45Mn2螺栓非正常拉伸断裂原因分析[J];热加工工艺;2010年04期
中国重要会议论文全文数据库 前10条
1 陈炜;王冬颖;郝建国;;基于Ansys Workbench离心压缩机合金结构钢法兰淬火过程温度场的数值模拟[A];第八届沈阳科学学术年会论文集[C];2011年
2 乔馨;于峰;刘源;;薄规格钢板淬火过程中的板形控制[A];2011年全国中厚板生产技术交流会论文集[C];2011年
3 石伟;武世勇;张向琨;刘庄;;淬火过程的三维数值模拟研究与软件开发[A];首届中国热处理活动周论文集[C];2002年
4 孔繁宇;;轧机类重载圆柱齿轮(轴)渗碳淬火过程中变形问题的探讨[A];第十五届华东地区热处理年会暨华东地区热处理年会三十周年纪念活动论文摘要集[C];2006年
5 王同珍;蔡大勇;李强;廖波;肖立春;;淬火过程组织转变孕育期的模拟计算[A];第三届中国热处理活动周暨第六次全国热处理生产技术改造会议论文专辑[C];2005年
6 周宏军;鞠颂东;贾国平;王晓慧;王磊;;圆柱工件淬火冷却过程温度场的数学模拟[A];数学·物理·力学·高新技术研究进展——1998(7)卷——中国数学力学物理学高新技术交叉研究会第7届学术研讨会论文集[C];1998年
7 陈洪;刘勇;邓雪莲;;基于PLC的森吉米尔中间辊淬火控制系统的设计与实现[A];第十一届全国自动化应用技术学术交流会论文集[C];2006年
8 周谦永;张伟民;陈乃录;王鸣华;宋冬利;;特种端淬实验方法的研究[A];第三届中国热处理活动周暨第六次全国热处理生产技术改造会议论文专辑[C];2005年
9 谢琰军;聂宝华;刘培英;;7050铝合金淬火析出模型[A];2011中国材料研讨会论文摘要集[C];2011年
10 吕成;张立文;牟正君;邰清安;郑渠英;杨钢;;重型燃机压气机盘淬火过程的数值模拟与工艺优化[A];第三届中国热处理活动周暨第六次全国热处理生产技术改造会议论文专辑[C];2005年
中国重要报纸全文数据库 前9条
1 记者 林庆雄;一条手链“淬”出遂昌旅游“黄金成色”[N];丽水日报;2009年
2 先越蓉;非晶态硅钢的组织变化[N];世界金属导报;2002年
3 王祝堂;大挤压铝材:中国走进世界前列(二)[N];中国有色金属报;2009年
4 耿晓稚 ;淬火炼精兵[N];人民公安报·消防周刊;2005年
5 袁国;新钢高强碳素中厚板辊式淬火机研制[N];世界金属导报;2011年
6 谢倩 袁国;新余钢铁公司辊式淬火机实现薄规格板材淬火稳定生产[N];世界金属导报;2009年
7 ;从民间幽默到文化品牌[N];山西日报;2003年
8 胡向泽 朱慧;万荣笑话:群众自己奶大的民间幽默[N];山西日报;2004年
9 王超;华菱涟钢实现国内最薄规格辊式淬火板连续稳定生产[N];世界金属导报;2011年
中国博士学位论文全文数据库 前10条
1 李强;淬火过程的计算机模拟与试验研究[D];燕山大学;2003年
2 李辉平;淬火过程有限元模拟关键技术及工艺参数优化的研究[D];山东大学;2005年
3 牛山廷;淬火冷却过程三维有限元模拟及工艺参数优化的研究[D];山东大学;2007年
4 王秋成;航空铝合金残余应力消除及评估技术研究[D];浙江大学;2003年
5 吕成;数值模拟技术在燃气轮机零部件锻造及热处理过程中的应用[D];大连理工大学;2007年
6 林高用;高性能7×75系铝合金厚板加工技术相关基础研究[D];中南大学;2006年
7 徐雪霞;模拟方法在贝氏体钢研究中的应用[D];清华大学;2008年
8 杨夏炜;铝合金大型复杂构件热处理过程的多场耦合模型与变形预报[D];哈尔滨工业大学;2013年
9 刘文军;Al-Zn-Mg-Cu铝合金淬火析出行为及淬火敏感性研究[D];中南大学;2011年
10 明平剑;基于非结构化网格气液两相流数值方法及并行计算研究与软件开发[D];哈尔滨工程大学;2008年
中国硕士学位论文全文数据库 前10条
1 华望;中厚板辊式淬火过程数值模拟研究[D];东北大学;2011年
2 丁冬芳;9SiCr常压气雾淬火过程的数值模拟与实验研究[D];昆明理工大学;2011年
3 刘海江;30CrMnMo无缝钢管淬火过程温度场、应力场数值模拟[D];内蒙古科技大学;2010年
4 徐斯强;重轨淬火过程的有限元分析[D];武汉科技大学;2013年
5 郭明清;中厚板辊式淬火过程温度场及应力场数值模拟[D];内蒙古科技大学;2012年
6 韩毅;中厚板辊式淬火机淬火过程的温度场及板形控制研究[D];东北大学;2008年
7 丁立勇;强约束条件下伞齿轮淬火工艺仿真及其变形控制[D];哈尔滨工业大学;2012年
8 冯媛媛;热处理零件组织性能分布的计算机预测及其结构设计[D];山东农业大学;2010年
9 郑泽花;钢淬火换热系数计算及温度与微观组织模拟[D];大连理工大学;2010年
10 王同珍;淬火过程的计算机模拟[D];燕山大学;2004年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026