Based on the design of the multi-row sprocket with a new tooth profile,a cold semi-precision forging process for manufacturing 5052 aluminum alloy multi-row sprocket was presented.Through simulating the forging proces...Based on the design of the multi-row sprocket with a new tooth profile,a cold semi-precision forging process for manufacturing 5052 aluminum alloy multi-row sprocket was presented.Through simulating the forging process of 5052 aluminum alloy sprocket billet with 3D rigid-viscoplastic FEM,both the distributions of flow velocity field in axial(U_Z),radial(U_R) and circumferential(U_θ) directions and the curves of velocity component in different deformation regions were respectively obtained.By comparison and analysis of the velocity varying curves,the velocity component relation conditions for filling the die cavity were clarified.It shows that when the die cavity is almost fully filled,the circumferential velocity U_θ increases sharply,implying that U_θplays a key role in fully filling the die cavity.展开更多
In order to consider the influence of temperature and underground water movement, an elastoplastic model and a 2D FEM stress fields on the migration of radioactive nuclide with code for analysis of coupled thermo-hyd...In order to consider the influence of temperature and underground water movement, an elastoplastic model and a 2D FEM stress fields on the migration of radioactive nuclide with code for analysis of coupled thermo-hydro-mechanical (THM) processes in saturated and unsaturated porous media were extended and improved through introducing the percolation and migration equation, so that the code can be used for solving the temperature field, flow field, stress field and nuclide concentration field simultaneously. The states of temperatures, pore pressures and nuclide concentrations in the near field of a hypothetical nuclear waste repository were investigated. The influence of the half life of the radioactive nuclide on the temporal change of nuclide concentration was analyzed considering the thermo-hydro-mechanical-migratory coupling. The results show that, at the boundary of the vitrified waste, the concentration of radioactive nuclide with a half life of 10 a falls after a period of rising, with the maximum value of 0.182 mol/m3 and the minimum value of 0.181 mol/m^3 at the end of computation. For a half life of 1 000 a, the concentration of radioactive nuclide always increases with the increase of the time during the computation period; and the maximum value is 1.686 mol/m^3 at the end of the computation. Therefore, under the condition of THM coupling, the concentration of radioactive nuclide with a shorter half life will decrease more quickly with water flow; but for the radioactive nuclide with a longer half life, its concentration will keep at a higher level for a longer time in the migration process.展开更多
This paper deals with the mathematical modelling and 3D FEM study of the energy release rate(ERR)in the band crack’s front contained in the orthotropic thick rectangular plate which is stretched or compressed initial...This paper deals with the mathematical modelling and 3D FEM study of the energy release rate(ERR)in the band crack’s front contained in the orthotropic thick rectangular plate which is stretched or compressed initially before the loading of the crack's edge planes.The initial stretching or compressing of the plate causes uniformly distributed normal stress to appear acting in the direction which is parallel to the plane on which the band crack is located.After the appearance of the initial stress in the plate it is assumed that the crack's edge planes are loaded with additional uniformly distributed normal forces and the ERR caused with this additional loading is studied.The corresponding boundary value problem is formulated within the scope of the so-called 3D linearized theory of elasticity which allows the initial stress on the values of the ERR to be taken into consideration.Numerical results on the influence of the initial stress,anisotropy properties of the plate material,the crack’s length and its distance from the face planes of the plate on the values of the ERR,are presented and discussed.In particular,it is established that for the relatively greater length of the crack’s band,the initial stretching of the plate causes a decrease,but the initial compression causes an increase in the values of the ERR.展开更多
The pile-supported subgrade has been widely used in high-speed railway construction in China.To investigate the ground vibrations of such composite foundation subjected to moving loads induced by high-speed trains(HST...The pile-supported subgrade has been widely used in high-speed railway construction in China.To investigate the ground vibrations of such composite foundation subjected to moving loads induced by high-speed trains(HSTs),three-dimensional(3D)finite element method(FEM)models involving the pile,pile cap and cushion are established.Validation of the proposed model is conducted through comparison of model predictions with the field measurements.On this basis,ground vibrations generated by HSTs under different train speeds as well as the ground vibration attenuation with the distance away from the track centerline are investigated.In addition,the effects of piles and pile elastic modulus on ground vibrations are well studied.Results show that the pile-reinforcement of the subgrade could significantly contribute to the reduction of ground vibrations.In particular,the increase of elastic modulus of pile could lead to consistent reduction of ground vibrations.However,when the pile elastic modulus is beyond 10 GPa,this benefit of pile-reinforcement on vibration isolation can hardly be increased further.展开更多
The rapid development of high-speed railway networks requires advanced methods for analysing vibration and sound radiation characteristics of a fast rotating train wheel subject to a vertical harmonic wheel-rail force...The rapid development of high-speed railway networks requires advanced methods for analysing vibration and sound radiation characteristics of a fast rotating train wheel subject to a vertical harmonic wheel-rail force. In order to consider the rotation of the wheel and at the same time increase the computational efficiency, a procedure is adapted in this paper taking advantage of the axial symmetry of the wheel. In this procedure, a recently developed 2.5D finite element method, which can consider wheel rotation but only requires a 2D mesh over a cross section containing the wheel axis, is used to calculate the vibration response of the wheel. Then, the vibration response of the wheel is taken as acoustic boundary condition and the 2.5D acoustic boundary element method, which only requires a 1D mesh over the boundary of the above cross section, is utilised to calculate the sound radiation of the wheel. These 2.5D methods and relevant programs are validated by comparing results from this procedure with those from conventional 3D analyses using commercial software. The comparison also demonstrates that these 2.5D methods have a much higher computational efficiency. Using the 2.5D methods, we study the wheel rotation speed influences on the factors including the vertical receptance of the wheel at wheel-rail contact point, sound pressure level at a pre-defined standard measurement point, radiated sound power level, directivity of the radia- tion, and contribution of each part of the wheel. It can be concluded that the wheel rotation speed splits most peaks of the vertical receptance at the wheel-rail contact point, sound pressure levels at the field, and the sound power level of the wheel into two peaks. The directivity and power contribution of the wheel are also significantly changed by the wheel rotation speed. Therefore, the rotation of a train wheel should be taken into account when calculating its vibration and sound radiation.展开更多
The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress...The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress functions. The stress fields are initially assumed by means of the Lekhnitskii stress functions under the plane strain state. Applying the principle of complementary virtual work,the coupled ordinary differential equations are obtained in which the solutions can be obtained by solving a generalized eigenvalue problem. Then an iterative procedure is established to achieve convergent stress distributions. It should be noted that the stress function based extended Kantorovich method can satisfy both the traction-free and free edge stress boundary conditions during the iterative processes. The stress components near the free edges and in the interior regions are calculated and compared with those obtained results by finite element method(FEM). The convergent stresses have good agreements with those results obtained by three dimensional(3D) FEM. For generality, various layup configurations are considered for the numerical analysis. The results show that the proposed polynomial stress function based extended Kantorovich method is accurate and efficient in predicting the local stresses in composite laminates and computationally much more efficient than the 3D FEM.展开更多
Forced and natural vibrations of a rectangular pre-stressed orthotropic compositeplate containing two neighboring cylindrical cavities whose cross sections are rectangularwith rounded-off corners are investigated nume...Forced and natural vibrations of a rectangular pre-stressed orthotropic compositeplate containing two neighboring cylindrical cavities whose cross sections are rectangularwith rounded-off corners are investigated numerically. It is assumed that all the end surfacesof the rectangular pre-stressed composite plate are simply supported and subjected to auniformly distributed normal time-harmonic force on the upper face plane. The consideredproblem is formulated within the Three-Dimensional Linearized Theory of Elastic Waves inInitially Stressed Bodies (TDLTEWISB). The influence of mechanical and geometricalparameters as well as the initial stresses and the effect of cylindrical cavities on the dynamicalcharacteristics of the rectangular orthotropic composite plate are analyzed and discussed.展开更多
The blade precision forging process is a forming process with high temperature and large plastic deformation. Interaction of deformation and heat conduction leads to large uneven distribution of temperature. The uneve...The blade precision forging process is a forming process with high temperature and large plastic deformation. Interaction of deformation and heat conduction leads to large uneven distribution of temperature. The unevenness of temperature distribution has a great effect on mechanical properties and the microstracture of materials. So it is necessary to consider the influence of temperature on the precision forging process of blades. Taking a blade with a tenon into consideration, a 3D mechanical model in precision forging is built up. The distribution laws of temperature field and the influence of the temperature on the equivalem stress in the process are obtained by using 3-D coupled thermo-mechanical FEM code developed by the authors Theresuits obtained illustrate that the influence of the temperature field on the blade forging process is considerable. The achievements of predicting microstructure and mechanical properties for forged blades is significant.展开更多
In order to explore the mineral resources buried in sea mud,it is necessary to use seabed resistivity measuring equipment,which works closer to the sediments than ordinary ship-based geophysical measuring equipment. B...In order to explore the mineral resources buried in sea mud,it is necessary to use seabed resistivity measuring equipment,which works closer to the sediments than ordinary ship-based geophysical measuring equipment. Because of the harsh environment of seafloor,high pressure and highly conductive seawater,marine magnetotelluric method developed slowly. The sea floor environment is similar to the environment of logging, According to the design of dual lateral logging equipment,a new equipment for seafloor electrical resistivity measurement is designed. Four 3D FEM models that contain resistivity abnormal targets are built to test the ability of this equipment to locate different shape of shallow buried resistivity abnormal targets in sea mud. The authors propose the method to correct the response curve while the bottom surface of this equipment is suspended or not parallel to the seafloor. The resistivity of targets can be calculated accurately.展开更多
This study explored the seismic response of a house supporting base sides with a polymeric displacement control material and by magnetically levitating the foundation base. In this paper, we explore the possibility an...This study explored the seismic response of a house supporting base sides with a polymeric displacement control material and by magnetically levitating the foundation base. In this paper, we explore the possibility and efficacy of a seismic-isolated detached house as described above from both a shaking table experiment of model and three-dimensional finite element analysis. The seismic-isolated model showed stable response and its acceleration response was significantly reduced compared to the base-fixed model in the shaking table test. Three-dimensional finite element analysis was possible to simulate the experimental results. In the seismic response analysis of a full-scale detached house, the seismic-isolated model showed response reduction and its residual displacement was smaller than that of the sliding-base isolation model.展开更多
Springback is one of important factors influencing the forming quality of numerical control (NC) bending of thin-walled tube. In this paper, a numerical-analytic method for springback angle prediction of the process...Springback is one of important factors influencing the forming quality of numerical control (NC) bending of thin-walled tube. In this paper, a numerical-analytic method for springback angle prediction of the process was put forward. The method is based on springback angle model derived using analytic method and simulation results from three-dimensional (3D) rigid-plastic finite element method (FEM). The method is validated through comparison with experimental results. The features of the method are as follows: (1) The method is high in efficiency because it combines advantages of rigid-plastic FEM and analytic method. (2) The method is satisfactory in accuracy, since the field variables used in the model is resulting from 3D rigid-plastic FEM solution, and the effects both of axial force and strain neutral axis shift have been included. (3) Research on multi-factor effects can be carried out using the method due to its advantage inheriting from rigid-plastic FEM. The method described here is also of general significance to other bending processes.展开更多
The disk-pad brake used in automobile is divided i nt o two parts: the disk, geometrically axisymmetric, and the pad, of which the geo metry is three-dimensional. In the course of braking, all parameters of the pro ce...The disk-pad brake used in automobile is divided i nt o two parts: the disk, geometrically axisymmetric, and the pad, of which the geo metry is three-dimensional. In the course of braking, all parameters of the pro cesses (velocity, load, temperature, physicomechanical and tribological characte ristics of materials of the couple, and conditions of contacts) vary with the ti me. Considerable evidence has show that the contact temperature is an integral f actor reflecting the specific power friction influence at the combined effect of load, speed, friction coefficient, thermo physical and durability properties of materials of a frictional couple. Furthermore, the physic mechanical state of t he interface of the disk and pads is determined not only by the contact temperat ure but also by the nonstationary temperature field. Using the two-dimensional model for thermal analysis implies that the contact conditions and frictiona l heat flux transfer are independent of θ. This may lead to false thermal elast ic distortions and unrealistic contact conditions. An analytical model is presen ted in this paper for the determination of contact temperature distribution on t he working surface of a brake. To consider the effects of the moving heat source (the pad) with relative sliding speed variation, a transient finite element tec hnique is used to characterize the temperature fields of the solid rotor with ap propriate thermal boundary conditions. And the transient heat conduction problem can be solved as a nominal 3-D transient heat transfer problem with an immovab le heat source. Numerical results shows that the operating characteristics of th e brake exert an essentially influence on the surface temperature distribution a nd the maximal contact temperature. The temperature field presents a noaxisymmet ric characteristic (a function of θ) and proves to be strongly localized and po ssesses a sharp gradient in both axial and radial directions.展开更多
基金Projects(51175363,51274149)supported by the National Natural Science Foundation of China
文摘Based on the design of the multi-row sprocket with a new tooth profile,a cold semi-precision forging process for manufacturing 5052 aluminum alloy multi-row sprocket was presented.Through simulating the forging process of 5052 aluminum alloy sprocket billet with 3D rigid-viscoplastic FEM,both the distributions of flow velocity field in axial(U_Z),radial(U_R) and circumferential(U_θ) directions and the curves of velocity component in different deformation regions were respectively obtained.By comparison and analysis of the velocity varying curves,the velocity component relation conditions for filling the die cavity were clarified.It shows that when the die cavity is almost fully filled,the circumferential velocity U_θ increases sharply,implying that U_θplays a key role in fully filling the die cavity.
基金Project(2010CB732101) supported by the National Basic Research Program of China Project(SKLQ 008) supported by the Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering of China
文摘In order to consider the influence of temperature and underground water movement, an elastoplastic model and a 2D FEM stress fields on the migration of radioactive nuclide with code for analysis of coupled thermo-hydro-mechanical (THM) processes in saturated and unsaturated porous media were extended and improved through introducing the percolation and migration equation, so that the code can be used for solving the temperature field, flow field, stress field and nuclide concentration field simultaneously. The states of temperatures, pore pressures and nuclide concentrations in the near field of a hypothetical nuclear waste repository were investigated. The influence of the half life of the radioactive nuclide on the temporal change of nuclide concentration was analyzed considering the thermo-hydro-mechanical-migratory coupling. The results show that, at the boundary of the vitrified waste, the concentration of radioactive nuclide with a half life of 10 a falls after a period of rising, with the maximum value of 0.182 mol/m3 and the minimum value of 0.181 mol/m^3 at the end of computation. For a half life of 1 000 a, the concentration of radioactive nuclide always increases with the increase of the time during the computation period; and the maximum value is 1.686 mol/m^3 at the end of the computation. Therefore, under the condition of THM coupling, the concentration of radioactive nuclide with a shorter half life will decrease more quickly with water flow; but for the radioactive nuclide with a longer half life, its concentration will keep at a higher level for a longer time in the migration process.
文摘This paper deals with the mathematical modelling and 3D FEM study of the energy release rate(ERR)in the band crack’s front contained in the orthotropic thick rectangular plate which is stretched or compressed initially before the loading of the crack's edge planes.The initial stretching or compressing of the plate causes uniformly distributed normal stress to appear acting in the direction which is parallel to the plane on which the band crack is located.After the appearance of the initial stress in the plate it is assumed that the crack's edge planes are loaded with additional uniformly distributed normal forces and the ERR caused with this additional loading is studied.The corresponding boundary value problem is formulated within the scope of the so-called 3D linearized theory of elasticity which allows the initial stress on the values of the ERR to be taken into consideration.Numerical results on the influence of the initial stress,anisotropy properties of the plate material,the crack’s length and its distance from the face planes of the plate on the values of the ERR,are presented and discussed.In particular,it is established that for the relatively greater length of the crack’s band,the initial stretching of the plate causes a decrease,but the initial compression causes an increase in the values of the ERR.
基金Project(51978510)supported by the National Natural Science Foundation of China。
文摘The pile-supported subgrade has been widely used in high-speed railway construction in China.To investigate the ground vibrations of such composite foundation subjected to moving loads induced by high-speed trains(HSTs),three-dimensional(3D)finite element method(FEM)models involving the pile,pile cap and cushion are established.Validation of the proposed model is conducted through comparison of model predictions with the field measurements.On this basis,ground vibrations generated by HSTs under different train speeds as well as the ground vibration attenuation with the distance away from the track centerline are investigated.In addition,the effects of piles and pile elastic modulus on ground vibrations are well studied.Results show that the pile-reinforcement of the subgrade could significantly contribute to the reduction of ground vibrations.In particular,the increase of elastic modulus of pile could lead to consistent reduction of ground vibrations.However,when the pile elastic modulus is beyond 10 GPa,this benefit of pile-reinforcement on vibration isolation can hardly be increased further.
基金the support to this work from the Ministry of Science and Technology of China (2016YFE0205200)China Railway (2015Z003-B)Scientific Research Foundation of State Key Laboratory of Traction Power (2017TPL_T01)
文摘The rapid development of high-speed railway networks requires advanced methods for analysing vibration and sound radiation characteristics of a fast rotating train wheel subject to a vertical harmonic wheel-rail force. In order to consider the rotation of the wheel and at the same time increase the computational efficiency, a procedure is adapted in this paper taking advantage of the axial symmetry of the wheel. In this procedure, a recently developed 2.5D finite element method, which can consider wheel rotation but only requires a 2D mesh over a cross section containing the wheel axis, is used to calculate the vibration response of the wheel. Then, the vibration response of the wheel is taken as acoustic boundary condition and the 2.5D acoustic boundary element method, which only requires a 1D mesh over the boundary of the above cross section, is utilised to calculate the sound radiation of the wheel. These 2.5D methods and relevant programs are validated by comparing results from this procedure with those from conventional 3D analyses using commercial software. The comparison also demonstrates that these 2.5D methods have a much higher computational efficiency. Using the 2.5D methods, we study the wheel rotation speed influences on the factors including the vertical receptance of the wheel at wheel-rail contact point, sound pressure level at a pre-defined standard measurement point, radiated sound power level, directivity of the radia- tion, and contribution of each part of the wheel. It can be concluded that the wheel rotation speed splits most peaks of the vertical receptance at the wheel-rail contact point, sound pressure levels at the field, and the sound power level of the wheel into two peaks. The directivity and power contribution of the wheel are also significantly changed by the wheel rotation speed. Therefore, the rotation of a train wheel should be taken into account when calculating its vibration and sound radiation.
基金supported by the National Natural Science Foundation of China (Grants 11372145, 11372146, and 11272161)the State Key Laboratory of Mechanics and Control of Mechanical Structures (Nanjing University of Aeronautics and astronautics) (Grant MCMS-0516Y01)+1 种基金Zhejiang Provincial Top Key Discipline of Mechanics Open Foundation (Grant xklx1601)the K. C. Wong Magna Fund through Ningbo University
文摘The extended Kantorovich method is employed to study the local stress concentrations at the vicinity of free edges in symmetrically layered composite laminates subjected to uniaxial tensile load upon polynomial stress functions. The stress fields are initially assumed by means of the Lekhnitskii stress functions under the plane strain state. Applying the principle of complementary virtual work,the coupled ordinary differential equations are obtained in which the solutions can be obtained by solving a generalized eigenvalue problem. Then an iterative procedure is established to achieve convergent stress distributions. It should be noted that the stress function based extended Kantorovich method can satisfy both the traction-free and free edge stress boundary conditions during the iterative processes. The stress components near the free edges and in the interior regions are calculated and compared with those obtained results by finite element method(FEM). The convergent stresses have good agreements with those results obtained by three dimensional(3D) FEM. For generality, various layup configurations are considered for the numerical analysis. The results show that the proposed polynomial stress function based extended Kantorovich method is accurate and efficient in predicting the local stresses in composite laminates and computationally much more efficient than the 3D FEM.
文摘Forced and natural vibrations of a rectangular pre-stressed orthotropic compositeplate containing two neighboring cylindrical cavities whose cross sections are rectangularwith rounded-off corners are investigated numerically. It is assumed that all the end surfacesof the rectangular pre-stressed composite plate are simply supported and subjected to auniformly distributed normal time-harmonic force on the upper face plane. The consideredproblem is formulated within the Three-Dimensional Linearized Theory of Elastic Waves inInitially Stressed Bodies (TDLTEWISB). The influence of mechanical and geometricalparameters as well as the initial stresses and the effect of cylindrical cavities on the dynamicalcharacteristics of the rectangular orthotropic composite plate are analyzed and discussed.
基金supported by the Aeronautical Science Foundation of China(No.02H53061)the National Science Found of China for Distinguished Young Scholar(No.50225518)the Shaan'xi Provincial Natural Science Foundation of China(No.2001CS0401)
文摘The blade precision forging process is a forming process with high temperature and large plastic deformation. Interaction of deformation and heat conduction leads to large uneven distribution of temperature. The unevenness of temperature distribution has a great effect on mechanical properties and the microstracture of materials. So it is necessary to consider the influence of temperature on the precision forging process of blades. Taking a blade with a tenon into consideration, a 3D mechanical model in precision forging is built up. The distribution laws of temperature field and the influence of the temperature on the equivalem stress in the process are obtained by using 3-D coupled thermo-mechanical FEM code developed by the authors Theresuits obtained illustrate that the influence of the temperature field on the blade forging process is considerable. The achievements of predicting microstructure and mechanical properties for forged blades is significant.
文摘In order to explore the mineral resources buried in sea mud,it is necessary to use seabed resistivity measuring equipment,which works closer to the sediments than ordinary ship-based geophysical measuring equipment. Because of the harsh environment of seafloor,high pressure and highly conductive seawater,marine magnetotelluric method developed slowly. The sea floor environment is similar to the environment of logging, According to the design of dual lateral logging equipment,a new equipment for seafloor electrical resistivity measurement is designed. Four 3D FEM models that contain resistivity abnormal targets are built to test the ability of this equipment to locate different shape of shallow buried resistivity abnormal targets in sea mud. The authors propose the method to correct the response curve while the bottom surface of this equipment is suspended or not parallel to the seafloor. The resistivity of targets can be calculated accurately.
文摘This study explored the seismic response of a house supporting base sides with a polymeric displacement control material and by magnetically levitating the foundation base. In this paper, we explore the possibility and efficacy of a seismic-isolated detached house as described above from both a shaking table experiment of model and three-dimensional finite element analysis. The seismic-isolated model showed stable response and its acceleration response was significantly reduced compared to the base-fixed model in the shaking table test. Three-dimensional finite element analysis was possible to simulate the experimental results. In the seismic response analysis of a full-scale detached house, the seismic-isolated model showed response reduction and its residual displacement was smaller than that of the sliding-base isolation model.
基金This work was supported by the National Natural Science Foundation of China for Distinguished Young Scholars (Grant No. 50225518)the Teaching and Research Award Program for 0utstanding Young Teachers in Higher Education Institution of M0E, PRCthe Aeronautical Science Foundation of China (Grant No. 04H53057).
文摘Springback is one of important factors influencing the forming quality of numerical control (NC) bending of thin-walled tube. In this paper, a numerical-analytic method for springback angle prediction of the process was put forward. The method is based on springback angle model derived using analytic method and simulation results from three-dimensional (3D) rigid-plastic finite element method (FEM). The method is validated through comparison with experimental results. The features of the method are as follows: (1) The method is high in efficiency because it combines advantages of rigid-plastic FEM and analytic method. (2) The method is satisfactory in accuracy, since the field variables used in the model is resulting from 3D rigid-plastic FEM solution, and the effects both of axial force and strain neutral axis shift have been included. (3) Research on multi-factor effects can be carried out using the method due to its advantage inheriting from rigid-plastic FEM. The method described here is also of general significance to other bending processes.
文摘The disk-pad brake used in automobile is divided i nt o two parts: the disk, geometrically axisymmetric, and the pad, of which the geo metry is three-dimensional. In the course of braking, all parameters of the pro cesses (velocity, load, temperature, physicomechanical and tribological characte ristics of materials of the couple, and conditions of contacts) vary with the ti me. Considerable evidence has show that the contact temperature is an integral f actor reflecting the specific power friction influence at the combined effect of load, speed, friction coefficient, thermo physical and durability properties of materials of a frictional couple. Furthermore, the physic mechanical state of t he interface of the disk and pads is determined not only by the contact temperat ure but also by the nonstationary temperature field. Using the two-dimensional model for thermal analysis implies that the contact conditions and frictiona l heat flux transfer are independent of θ. This may lead to false thermal elast ic distortions and unrealistic contact conditions. An analytical model is presen ted in this paper for the determination of contact temperature distribution on t he working surface of a brake. To consider the effects of the moving heat source (the pad) with relative sliding speed variation, a transient finite element tec hnique is used to characterize the temperature fields of the solid rotor with ap propriate thermal boundary conditions. And the transient heat conduction problem can be solved as a nominal 3-D transient heat transfer problem with an immovab le heat source. Numerical results shows that the operating characteristics of th e brake exert an essentially influence on the surface temperature distribution a nd the maximal contact temperature. The temperature field presents a noaxisymmet ric characteristic (a function of θ) and proves to be strongly localized and po ssesses a sharp gradient in both axial and radial directions.