20米自動(dòng)伸縮門設(shè)計(jì)[含CAD高清圖紙和文檔資料]

【溫馨提示】 dwg后綴的文件為CAD圖，可編輯，無水印，高清圖，壓縮包內(nèi)文檔可直接點(diǎn)開預(yù)覽，需要原稿請自助充值下載，請見壓縮包內(nèi)的文件及預(yù)覽，所見才能所得，請細(xì)心查看有疑問可以咨詢QQ：414951605或1304139763 外文翻譯 Artificial hip joint學(xué)生姓名 班 級 學(xué)院名稱 專業(yè)名稱 指導(dǎo)教師 1. IntroductionIt has been recognised by a good number of researchers that the computation of the pressure distribution and contact area of artificial hip joints during daily activities can play a key role in predicting prosthetic implant wear 1, 2, 3 and 4. The Hertzian contact theory has been considered to evaluate the contact parameters, namely the maximum contact pressure and contact area by using the finite element method 1 and 2. Mak and his co-workers 1 studied the contact mechanics in ceramic-on-ceramic (CoC) hip implants subjected to micro-separation and it was shown that contact stress increased due to edge loading and it was mainly dependent on the magnitude of cup-liner separation, the radial clearance and the cup inclination angle 3 and 4. In fact, Hertzian contact theory can captured slope and curvature trends associated with contact patch geometry subjected to the applied load to predict the contact dimensions accurately in edge-loaded ceramic-on-ceramic hips 5. Although the finite element analysis is a popular approach for investigating contact mechanics, discrete element technique has also been employed to predict contact pressure in hip joints 6. As computational instability can occur when the contact nodes move near the edges of the contact elements, a contact smoothing approach by applying Gregory patches was suggested 7. Moreover, the contributions of individual muscles and the effect of different gait patterns on hip contact forces are of interest, which can be determined by using optimisation techniques and inverse dynamic analyses 8 and 9. In addition, contact stress and local temperature at the contact region of dry-sliding couples during wear tests of CoC femoral heads can experimentally be assessed by applying fluorescence microprobe spectroscopy 10. The contact pressure distribution on the joint bearing surfaces can be used to determine the heat generated by friction and the volumetric wear of artificial hip joints 11 and 12. Artificial hip joint moment due to friction and the kinetics of hip implant components may cause prosthetic implant components to loosen, which is one of the main causes of failure of hip replacements. Knee and hip joints moment values during stair up and sit-to-stand motions can be evaluated computationally 13. The effect of both body-weight-support level and walking speed was investigated on mean peak internal joint moments at ankle, knee and hip 14. However, in-vivo study of the friction moments acting on the hip demands more research in order to assess whether those findings could be generalised was carried out 15.The hypothesis of the present study is that friction-induced vibration and stick/slip friction could affect maximum contact pressure and moment of artificial hip joints. This desideratum is achieved by developing a multibody dynamic model that is able to cope with the usual difficulties of available models due to the presence of muscles, tendons and ligaments, proposing a simple dynamic body diagram of hip implant. For this purpose, a cross section through the interface of ball, stem and lateral soft and stiff tissues is considered to provide the free body diagram of the hip joint. In this approach, the ball is moving, while the cup is considered to be stationary. Furthermore, the multibody dynamic motion of the ball is formulated, taking the friction-induced vibration and the contact forces developed during the interaction with cup surface. In this study, the model utilises available information of forces acting at the ball centre, as well as angular rotation of the ball as functions of time during a normal walking cycle. Since the rotation angle of the femoral head and their first and second derivatives are known, the equation of angular momentum could be solved to compute external joint moment acting at the ball centre. The nonlinear governing equations of motion are solved by employing the adaptive RungeKuttaFehlberg method, which allows for the discretisation of the time interval of interest. The influence of initial position of ball with respect to cup centre on both maximum contact pressure and the corresponding ball trajectory of hip implants during a normal walking cycle are investigated. Moreover, the effects of clearance size, initial conditions and friction on the system dynamic response are analysed and discussed throughout this work.2. Multibody dynamic model of the artificial hip jointThe multibody dynamic model originaly proposed by Askari et al. 16 has been considered here to address the problem of evaluating the contact pressure and moment of hip implants. A cross section A-A of a generic configuration of a hip joint is depicted in the diagram of in Fig. 1, which represents a total hip replacement. Fig. 1 also shows the head and cup placed inside of the pelvis and separated from stem and neck. The forces developed along the interface of the ball and stem are considered to act in such a way that leads to a reaction moment, M. This moment can be determined by satisfying the angular motion of the ball centre during a walking cycle. The available data reported by Bergmann et al. 17 is used to define the forces that act at the ball centre. This data was experimentally obtained by employing a force transducer located inside the hip neck of a live patient. The information provided deals with the angular rotation and forces developed at the hip joint. Thus, the necessary angular velocities and accelerations can be obtained by time differentiating the angular rotation. Besides the 3D nature of the global motion of the hip joint, in the present work a simple 2D approach is presented, which takes into account the most significant hip action, i.e. the flexion-extension motion. With regard to Fig. 2 the translational and rotational equation of motion of the head, for both free flight mode and contact mode, can be written by employing the NewtonEulers equations 18 and 19, yieldingequation(1)MOk=Ik,MO=Mk(Rj)nFPjt0Mk0equation(2)FX=mx,FX=fx+(FPjt+FPjn)i0fx0equation(3)FY=my,FY=fy+(FPjn+FPjt)jmg0fy0where FPjn and FPjt denote the normal and tangential contact forces developed during the contact between the ball and cup, as it is represented in the diagram of Fig. 3. In Eqs. (1), (2) and (3), x, y and are the generalised coordinates used to define the systems configuration. In turn, variable m and I are the mass and moment of inertia of ball, respectively. The external generalised forces are denoted by fx, fy and M and they act at the centre of the ball as it is shown in Fig. 3. The gravitational acceleration is represented by parameter g, Rj is the ball radius and represents relative penetration depth between the ball and cup surfaces.Fig. 1. A schematic of the artificial hip implant with the cross section A-A (Left figure), and the head and cup separated from the neck and stem through the cross section A-A (Right figure).Figure options Download full-size image Download high-quality image (142 K) Download as PowerPoint slideFig. 2. A schematic of the head and cup interaction observed in the Sagittal plane.Figure options Download full-size image Download high-quality image (87 K) Download as PowerPoint slideFig. 3. Free body diagram of ball and corresponding external, internal and body forces and moment.Figure options Download full-size image Download high-quality image (87 K) Download as PowerPoint slideTheDownload as PowerPoint slideThe penetration depth can be expressed as 20equation(4)=r(RbRj)=r(RbRj)in which Rb denotes the cup radius and (RbRj) represents the joint radial clearance, which is a parameter specified by user.In the present study, the cup is assumed to be stationary, while the head describes the global motion. With regard to Fig. 2, it can be observed that Oj and Ob denote the head and cup centres, respectively. While Pj and Pb represent the contact points on the head and cup, respectively. The magnitude and orientation of the clearance vector are denoted by r and , respectively. In general, r and can be expressed as functions of the generalised coordinates used to describe the configuration of multibody mechanical system. The normal and tangential unit vectors at the contact point can be written asequation(5)n=cosi+sinjequation(6)t=sini+cosjIn order to compute the normal contact and tangential forces, it is first necessary to evaluate the relative tangential and normal velocities at the contact points, which can be obtained as follows: 17equation(7)vpj/pb=rn+(r+Rjj)t=vnn+vttwhere vnvn and vtvt are module of the normal and tangential velocities, respectively. Thus, Eqs. (2) and (3) can be re-written as follows:equation(8)m00mxy=FXFYUsing now the concept of the state space representation, the second order equations of motion (8), can be expressed as a first order equation set asequation(9)z=H(z)where z=z1z2z3z4=xyxy and H(z)H(z) is expressed as follows:equation(10)z=z1z2z3z4=z3z4FX(z)FY(z)It must be mentioned that the r, and their time derivatives can be obtained with respect to state space parameters as follows:equation(11)=atan(z2z1)equation(12)r=z12+z22equation(13)=z2z3+z1z4z12+z22z12equation(14)r=z1z3+z2z4z12+z22It is known that the evaluation of the contact forces developed during an impact event plays a crucial role in the dynamic analysis of mechanical systems 21, 22 and 23. The contact forces must be computed by using a suitable constitutive law that takes into account material properties of the contacting bodies, the geometric characteristics of impacting surfaces and impact velocity. Additionally, the numerical approach for the calculation of the contact forces should be stable in order to allow for the integration of the mechanical systems equations of motion 24. Different constitutive laws are suggested in the literature, being one of the more prominent proposed by Hertz 25. However, this law is purely elastic in nature and cannot explain the energy loss during the impact process. Thus, Lankarani and Nikravesh 26 overcame this difficulty by separating the contact force into elastic and dissipative components asequation(15)Fpjn=(Kn+D)nRegarding Lankarani and Nikravesh model, normal contact force on the head is expressed asequation(16)Fpjn=K3/2(1+3(1ce2)4()nwhere and () are the relative penetration velocity and the initial impact velocity, respectively, and ce is the coefficient of restitution. The generalised stiffness parameter K depends on the geometry and physical properties of the contacting surfaces, which for two internal spherical contacting bodies with radii Ri and Rj can be expressed as 25equation(17)K=43(i+j)(RiRjRiRj)2in which the material properties ii and jj are given byequation(18)z=1z2EzAt this stage, it must be said the use of Eq. (15) is limited by Loves criterion, that is, it is only valid for impact velocities lower than the propagation velocity of elastic waves across the solids 27.It is known that the way in which the friction phenomena are modelled, plays a key role in the systems behaviour 28. In the present study, the tangential friction force are evaluated by using a modified Coulomb friction law, which can be expressed as 29 and 30equation(19)Fpjt=(vt)|Fpjn|vt|vt|tThe friction force is described in the sense of Coulombs approach, and is proportional to the magnitude of the normal force developed at the contact points, where the ratio is the coefficient of friction, , which is dependent on the relative tangential velocity. The model considered in reference 17 is employed here for the purpose of evaluating the coefficient of friction, which can be written asequation(20)(vt)=(cfv02(|vt|v0)2+cf)sgn(vt),|vt|v0(cd+(cfcd)exp(|vt|v0)sgn(vt),|vt|v0where is a slop parameter. The first part of friction coefficient function is exhibits a continuous behaviour when the function is close to zero, in order to avoid the numerical instabilities associated with null tangential velocity. Fig. 4 shows the plot of Eq. (20), which represents the Stribeck effect. Thus, it can be stated that with this approach, the stick/slip effect can be taken into account. In fact, the model can represent the dry friction behaviour and address stick/slip phenomenon in relative low velocity case accurately. Moreover, this modified Coulombs friction model can avoid computational instability associated with the change of velocity direction.Fig. 4. Stribeck characteristic for dry friction.Figure options Download full-size image Download high-quality image (102 K) Download as PowerPoint slideNormal and tangential forces described above are present if the system isin contact situation, which means detecting impact or contact is one important step. Moreover, to compute the contact force, the initial impact velocity has to be calculated as an initial condition for following regimes, which could be in contact or in free flight, the following condition should be checked during the solution process by progressing time. Therefore, a contact event is detected when the following condition is verifiedequation(21)(ti)0Indeed, the precise instant i