点石文库
全部分类
  • 学术论文 >
    学术论文
    机械毕业设计 电气毕业论文 土木工程论文 医学论文 法学论文 管理论文 人力资源论文 计算机论文 软件工程 英语学论文 教育学论文 开题报告 冶金学论文 任务书 通信工程论文 生物学论文 毕业设计 毕业论文 文献综述 外文翻译 答辩PPT 参考文献 课程设计 期刊论文 减速器设计 其他
  • 机械毕业设计精选 >
    机械毕业设计精选
    含3D结构设计 结构设计/毕业设计 注塑模具设计(含模 注塑模具设计(含数 注塑模具设计 夹具设计类 夹具设计(机床类课 夹具设计(含三维模 数控工艺类 冲压模具(单工序) 冲压模具(复合模) 冲压模具(级进模) 压铸及其他模具类 matlab仿真类 ansys仿真类 plc控制类 单片机控制类 其他控制类 减速器课程类
  • 图纸模型 >
    图纸模型
    机械设备 零部件模型 交通运输 电子产品 生活用品 电子电工 五金工具 军工模型 建筑模型 模具图纸 钣金图纸 设计方案 机器人模型 航空航天 海洋船舶 艺术品/工艺品 CAD建筑图纸 文体用品 科幻模型 其他模型
  • 土木建筑 >
    土木建筑
    工程造价 室内装修 建筑图纸 规划方案 市政工程 园林工程 结构设计 环保行业 建筑设计 水电图 建筑标准 安全施工 建筑材料 技术标书 其他 施工组织
  • 行业资料 >
    行业资料
    机械类 仿真类 国家标准 企业标准 机械标准 金融经济 医学类 交通电力 全国省级标准 环保消防 机械行业标准 铁路行业标准 化工行业标准 建筑行业标准 城建市政标准 其他类
  • 办公文档 >
    办公文档
    PPT模板 工作汇报/总结 演讲稿/致辞 工作计划 活动策划 读后感/观后感 调研文书 招标投标 产品使用说明 会议纪要 商业策划 合同协议 商务礼仪 财务报表 广告营销 通知/申请 制度体系 个人简历 其他
  • 认证考试 >
    认证考试
    财会类 学历类 公务员/事业单位类 教师资格考试 技工职业考试 司法考试 网络工程师考试 质量工程师 成考/自考/函授 建筑类 外语类 资格类 外贸类 医药类 计算机类 其他
  • 教育辅导 >
    教育辅导
    幼儿教育 小学资料 中学资料 高中资料 成人高考 大学教育 研究生教育 自考/成人/函授 考试试卷 高中政治 高中生物 高中地理 高中数学 高中语文 高中外语 高中历史 其他
  • 生活休闲 >
    生活休闲
    运动健康 养生知识 服装配饰 科普知识 时政新闻 游戏攻略 旅游攻略 两性情感 美食烹饪 摄影摄像 其他 党团政务
  • 首页 点石文库 > 资源分类 > PDF文档下载
     

    (有限元分析与设计)磨削轴承钢52100引起相变和残余应力的数值模拟.pdf

    • 资源ID:3311       资源大小:1.13MB        全文页数:11页
    • 资源格式: PDF        下载权限:游客/注册会员/VIP会员    下载费用:1
    换一换
    游客快捷下载 游客一键下载
    会员登录下载
    下载资源需要1

    邮箱/手机:
    温馨提示:
    支付成功后,系统会根据您填写的邮箱或者手机号作为您下次登录的用户名和密码(如填写的是手机,那登陆用户名和密码就是手机号),方便下次登录下载和查询订单;
    特别说明:
    请自助下载,系统不会自动发送文件的哦;
    支付方式: 微信支付    支付宝   
    验证码:   换一换

          加入VIP,下载共享资源
     
    友情提示
    2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,既可以正常下载了。
    3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
    4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰   

    (有限元分析与设计)磨削轴承钢52100引起相变和残余应力的数值模拟.pdf

    dKeywordsResidual stressesGrindingstressandmaterial.The thermo-metallurgical and mechanical analysis has been performed using the commercial finiteelement software Abaquss/Standard with various user subroutines developed to model the thermal,processes,unitthe energytovarying internal stress/strain field. The evolution of temperaturesmaindis-usedtemperature. The FE model developed in this work takes intoContents lists available at SciVerse ScienceDirectFinite Elements in AnalysFinite Elements in Analysis and Design 61 2012 1–11Peclet Pe number, non-dimensional heat transfer coefficient HE-mail address michel.coretinsa-lyon.fr M. Coret.during grinding can be well predicted using models presented by account the transformation strains associated with martensitictransformation along with the temperature dependent materialproperties. The variations of the residual stresses and strains atintegration points have been considered, and the effects of the0168-874X/-see front matter fax 33 4 72438913.occurs when the effective stress exceeds the yield strength. Allthese factors interact with each other and eventually lead to asteels in engineering components such as rolling bearings. More-over, its properties are easy to acquire over a wide range of[1–4], and therefore some variations of the material propertieslike hardness, etc.. These phenomena may in fine play a role onthe surface integrity of the component. Phase transformations inmost steels introduce volumetric changes, transformation plasti-city and changes in mechanical properties. Local plastic flowpart of these phenomena. However, the influence of metallurgicaltransformations remains a problem still widely open. Theobjective of this work is to investigate the internal stresstribution and their evolution that occurs during grinding processof AISI-52100 steel which is one of the most commonlyuniform high temperatures, which may then result in solid-statephase transformations. The consequences of this rapid heatingand cooling include the formation of a heat affected zone HAZ,the generation of residual stresses, possible shrinkage or crackingof the material, often chemical modifications of the materialstresses in the near-surface layer.Numerical simulation of such a problem requires modeling ofthree different types of phenomena thermal, metallurgical, andmechanical, which are, mostly, fully coupled. Many researchers[10–14] have proposed various models which account for all or1. IntroductionCompared with other machiningan extremely high energy input persurface layer of the material. Most ofheat which is concentrated in the grindinginteracts with the workpiece. This leadsmetallurgical and mechanical behavior of the material. The heat generated during grinding process wasassumed as a moving heat flux with elliptical distribution. The effects of the Peclet number and heattransfer coefficient on the phase transformations and residual stresses have been analyzed. It was foundthat an optimal combination of grinding conditions could produce the desired magnitude ofcompressive residual stresses at the surface of the machined workpiece. It is also shown that omittingphase transformations could lead to a strong difference in the prediction of residual stresses.8Finally the parameters kiT and niT are obtained by solvingEq. 8kiTC01tEniTln1C0099zeq9niT1ln tStEC0C1lnln1C0001ln1C0099zeqC20C2110The martensitic transformation is a displacive transformationonly controlled by the temperature. The kinetics is given by theKoistinen–Marburger equation [25]zMzg1C0exp C0bMSC0TC2C3C8C911where, zMand zgare the martensitic and austenitic phaseproportions, respectively, b is a material dependent coefficient;MSis the martensitic transformation start temperature, and T thetemperature.2.4. Mechanical modelAfter the thermo-metallurgical computations, temperaturesand volume phase fraction become inputs for the mechanicalsimulation. The temperature variations and the phase transfor-mations involve dilatational strains into the solid. The majorremaining difficulty is to obtain the mechanical behavior of themixture of phases. As a first assumption, the macroscopic beha-vior is supposed to follow an isotropic hardening with the vonMises criterion where the yield stress is obtained by a mixture ofthe yield stress of each phase. The flow stresses which take intoaccount hardening, thermal and viscous effects follow a Johnson–Cook model as initially proposed by Umbrello et al. [26] andadapted to multiphase materials.In the following, the total strain rate tensor is divided into arecoverable elastic part _eeijand an irrecoverable plastic one _epij_eij _eeij _epij122.4.1. Elastic strainIn the case of coupling of mechanical field with temperatureand phase change, the relation for the elastic strain can beexpressed aseeij1vEsijC0vEskkdijEthmdij13EthmT,za1C0zaethgTzaethaT 14ethaTaaTC0TrefC2C315ethgTagTC0TrefC2C3C0De25 1Cag16where E is Young’s modulus and v Poisson’s ratio.2.4.2. Yield function and plastic strain rateThe yield function criterion is expressed in stress space asF sij,ep,Hm,T,ziC16C170 17ijwith sayand sbythe yield stresses of the ferrite and austenitephases, respectively.The equivalent plastic strain is given byeplZt0_epldt 21with_eplffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi23_eplij_eplijr222.4.3. Flow stressesFor an accurate simulation of grinding, the mechanical beha-vior must take into account the hardening, thermal, and viscouseffects. A usual way to do that is the use of the Johnson–Cookrelation, well adapted for severe loadings. In [26], the authorspropose a slight modification of Johnson–Cook relation fornumerical simulation of hard machining of AISI 52100 bearingsteel. The stresses are a product of three terms, representing thehardening curve at room temperature, the influence of tempera-ture and the last one the visco-plastic part Eq. 23. In [26], thetemperature influence function zfactis interpolated as a 5th orderpolynomial function whose coefficients are given in Table 1. Thestrain hardening multiplier needs two other constants m, A, alsogiven in Table 1.sgeq CenC2C3zfactTC2C31ln_eplmC0Ahi23withzfactexp aT5bT4cT3dT2eTfC16C1724The previous single phase model must be adapted for multiphasematerials. For sake of simplicity, the following two assumptions aremade1 the flow curve for each phase follows an exponential lawsziCziepC0C1nzi;Table 1Parameters used in Eqs. 23 and 24 for AISI52100 steel.Parameters Valuea 3.81C210–15b C04.29C210–12c C06.91C210–9d 5.50C210–6e C01.60C210–3f 2.44C210C02m 0.1259A 0.0567sgyand calculated by1C0zasgyepl,T,ziC16C17zasayepl,T,ziC16C17202 the mixture is calculated as a simple linear rule sPi 1,nzziszihiwhich is consistent with a Voigt hypothesison the strains. Table 2The stress–strain curves at room temperatures for each phasesare presented in Appendix A.Finally, the relation used issgeqe, _e,T,ziXi 1,nzziszi“zfactTC2C31ln_eplmC0Ahi25with sziCzienzi .3. Finite element simulationthe AISI 52100 steel is listed in Table 2 and the key physical andmechanical properties of the material are given in Appendix A3.3. Initial and boundary conditionsThe initial temperature considered for the workpiece is theroom temperature i.e. Tt020 1C.On lateral faces, heat flux is imposed as linear convectivetransfer lawqconvhconvTC0T026where T and T0are the temperature of the semi-infinite solid andthe ambient temperature, respectively, and hconvW/K m2 is theconvective heat transfer coefficient of the cooling media. Heat lossfrom the bottom surface was assumed to be zero i.e. q0. Thethermal boundary conditions are schematically shown in Fig. 6.lemS.M. Shah et al. / Finite Elements in Analysis and Design 61 2012 1–11 50.1m 200 eThe principle of numerical simulation of grinding entails theknow-how of a comprehensive database with reference to geo-metry, thermo-mechanical properties, initial conditions, bound-ary and loading conditions. A brief description of the mainfeatures of the FE model is given below.3.1. Finite element meshThe mesh density is generally defined by the applied loadingand/or boundary conditions. Since grinding processes involvehigh temperature gradient in and near the grinding zone, a veryfine mesh is required to capture the temperature distribution inthe contact area. As the temperature gradient becomes low faraway from the grinding zone, a relatively coarser mesh is theresufficient for the analysis.In this study the workpiece is considered as a 2D semi-infiniteplate of 0.1 m length and 0.03 m width. The finite element FEmesh Fig. 5 consists of CPE4T 4-node plane strain thermallycoupled quadrilateral bilinear displacement and temperaturetype elements totaling over 3216 nodes and 3000 elements withthe smallest element in the mesh as [5C210C048C210C05] m2.3.2. Material bearing steel AISI 52100AISI 52100 also known as 100Cr6 in Europe is a high carbon–chrome–manganese alloy steel which finds its applications inseveral rotating parts like anti-friction bearings, cams, crank shaft,etc. for its good resistance to corrosion and fatigue [27]. Com-pared to low-carbon steels, high-carbon steels can carry highercontact stresses, such as those encountered in point contactloading in rolling bearings [28]. The chemical composition ofFig. 5. Finite element3.4. Imposed heat sourcesThe thermal loading consists in applying a surface heat fluxthrough a moving heat source. Jaeger [29] and Carslaw and Jaeger[5] have presented solutions for uniform moving rectangular heatsources and a uniform stationary heat source using the heatsource method. The temperature distribution in a sliding contactwas then estimated by several authors based on Jaeger’s theory[30–34]. There are differing views among researchers on whichdistribution of heat flux is best to use for grinding. Some [7,35,36]have used a rectangular uniform distribution, so as to simplifysubsequent calculations. However, due to the localized ‘‘spike’’temperatures during a very short time, others [4,37,38] haveargued that the assumption of a uniform heat flux field may notlead to accurate predictions. Keeping in view the contact origin ofthe heat source, theoretically the pressure and correspondingheat flux distribution – if one assumes a uniform friction coeffi-cient in the contact area – should be modeled according to asliding/rolling contact approach. Also by recalling that the Hertzcontact pressure distribution between a cylinder and a plane iselliptical in shape, it seems reasonable to assume an ellipticaldistribution of the heat flux. Fig. 6 shows the presence of aschematic heat source moving with velocity Vwon the top surfaceof the FE model. Here, the length of the heat source is equal to thecontact length 2aLc between the grinding wheel and theworkpiece. The heat flux distribution entering the work piece istherefore given byqx2Qpaffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi1C0x2a2C18C19s27where Q is the total heat per unit length in W/m.ents mesh.Fig. 7. Temperature distribution at the surface of the workpiece for different heatsource shapes.0200400600800100012000Distance 2x/LcTemperature [C]2y/Lc02y/Lc0.1582y/Lc0.3422y/Lc0.5562y/Lc1.092y/Lc1.42Q 400 W/mLc 0.003mH 0.01Pe 1.0t΄33.61xyNodal TemperatureAc1Ac35 10152025Fig. 8. Temperature profile along the x-axis at various depths.S.M. Shah et al. / Finite Elements in Analysis and Design 61 2012 1–116In Abaquss/Standard, the moving heat source is integrated inthe finite element model through a FORTRAN subroutine, calledDFLUX. The heat source is moving along the horizontal x- axis.Before carrying out a complete thermo-mechanical analysis, afew heat-transfer simulations were run to perform a sensitivityanalysis of various types of heat source distributions. An elliptical,a triangular, and a uniform heat source were used. The resultingtemperature distributions can be compared in Fig. 7. It was foundthat the peak surface temperatures in all three cases are veryclose, however, the distribution of temperatures over the surfacevaries to some extent. An interesting observation was that theflux distribution from elliptical source lies almost midwaybetween the triangular and the uniform heat sources.3.5. Implementation of the mechanical behavior of a multiphasematerialAs shown in Fig. 3, the mechanical behavior depends on thephase proportion but this coupling is not directly available withInsulation q0WORKPIECE Fig. 6. Thermal loading and boundary conditions.Table 2Chemical composition of AISI 52100 100Cr6 bearing steel.Element C Si MnMass 0.95–1.10 0.15–0.35 0.20–0.40Convective heat transfer q Heat source VAbaquss/Standard. Therefore a UMAT subroutine has been devel-oped to go beyond this difficulty. The UMAT subroutine calls threeother subroutines PHASE, PROP and UEXPAN. For a given tem-perature field, PHASE compute the austenitic phase proportionformed during heating, and the martensitic phase proportionformed during cooling see Section 2.3. Knowing the phaseproportion, UEXPAN subroutine gives the expansion coefficientfor the mixed material at a given temperature see Eqs. 15 and16. Finally, PROP subroutine compute the material properties ata given temperature based on the fraction of phases. A linearmixture rule was used for the identification of the multiphasematerial properties.4. Results and discussion4.1. Temperature distributionThe temperature variation in time or profile along thehorizontal axis calculated at a given time instance t033.61and at various dimensionless depth y/a2y/Lc is shown in Fig. 8for a specific set of parameters Q, Lc, H and Pe. From the values ofpeak temperatures, it is found that up to a certain dimensionlessdepth here 2y/Lc0.30 the temperature goes beyond Ac1andS Cr Mo Pr0.025 1.35–1.60 r0.10 r0.03003006009001200-20Distance related to start of Flux [mm]Max. Temperature [C] Elliptical Heat distributionTriangualr Heat DistributionUniform Heat DistributionQ400W/m Lc0.003m Pe1.0-15 -10 -5 0 5 10Ac3750 1C and 815 1C [39] for AISI 52100, respectively. Duringcooling the temperature at these points will quickly drops belowMs250 1C [39]. It means that at high cooling the transformationof austenite to martensite will occur at the top most surface.In Fig. 9 the maximum surface temperature as a function of thePeclet number is plotted for a given set of grinding parameters Qand Lcas specified. It is shown that the peak temperaturedecreases when increasing the Peclet number. A comparison withthe analytical solution of Blok [34] is also provided and a goodagreement is found, which validates the numerical model. Theeffect of the variation of the dimensionless heat transfer coeffi-cient, H, or the contact length, Lc, is illustrated in Fig. 10. It can beobserved that an increase of the heat transfer coefficientdecreases the maximum temperature as an increase in thecontact ar

    注意事项

    本文((有限元分析与设计)磨削轴承钢52100引起相变和残余应力的数值模拟.pdf)为本站会员(王牌秘书)主动上传,点石文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知点石文库(发送邮件至3339525602@qq.com或直接QQ联系客服),我们立即给予删除!

    温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




    关于点石文库 - 投资与合作 - 会员权益 - 联系我们 - 声明 - 人才招聘

    本站资源为会员上传分享,如有侵犯您的版权,请联系我们删除

    网站客服QQ:3339525602  网上上传投稿QQ群862612017

      copyright@ 2016-2020  dswenku.com 网站版权所有   

    经营许可证编号:湘ICP备18013834 


    收起
    展开