双侧下颌支矢状劈开截骨术后患者可施加合适咬合力的数值研究

双侧下颌支矢状劈开截骨术在临床手术中广泛用于治疗缩颌、凸颌和颌面不对称患者的颌面畸形症状。本文用数值方法研究双侧下颌支矢状劈开截骨术并前移下颌骨后患者可施加的合适咬合力。首先,用缩颌患者的CT断层扫面数据建立下颌骨、关节盘、关节窝和相关软组织的高度几何相似性的几何模型;然后,用有限元模型模拟双侧下颌支矢状劈开截骨术后患者的咬合过程。通过分析发现,术后关节盘上的位移处于正常的范围,而应力却远远超过正常范围,因此,将应力因素作为确定术后患者可施加的合适咬合力的标准。当施加正常咬合力的25%时,关节盘上的应力分布和术前没有显著性差异,同时下颌骨上的应力水平也处于正常范围。     

Crowe IV型发育性髋关节发育不良在全髋关节置换术中横行截骨最佳位置的有限元分析

Crowe IV型髋关节发育不良（developmental dysplasia of the hip,DDH）的全髋关节置换术（total hip arthroplasty, THA）中，多数患者需要截骨保证复位，目前横行截骨是运用最广泛的术式，但截骨的最佳高度仍存在争议。本文通过建立Crowe IV型DDH患者的有限元模型，分析了在不同载荷下，不同截骨高度模型的应力分布与位移情况。结果表明：在多项指标下，小转子下0 cm的截骨高度均表现出了较优性能。在不影响假体安放和断端对合情况允许的情况下，推荐在手术中使用0 cm的截骨高度。     

骨的力–电效应研究方法与进展

骨组织在受到应力作用(正常的生理活动)变形后在骨内产生电位的现象称为骨的力–电效应,它主要包括压电效应和动电效应.研究骨在动态过程中产生的电位幅值和分布特点,不仅是了解电刺激骨生长机理的必要步骤,也是实现骨治疗和重建的生理基础.它一方面是用数学方法来描述外力作用下其电位大小与应力、应变、应变率及加载速率的关系,另一方面是考察生理环境(pH值、离子浓度、温度、湿度等)对电位的影响.首先对力–电理论进行了简单的介绍,重点总结了其研究方法,包括理论模型和分离式霍普金生杆冲击、弯曲变形及缓冲液中的动态测试等实验方法.此外,对骨替代材料和牙本质领域的力–电效应研究也进行了一定的综述.     

腰椎椎间植骨融合有限元分析

为了评估植骨融合术对腰椎运动和应力分布的影响, 利用影像诊断技术(CT扫描) 和CAD三维重建技术获得人体腰椎三维模型, 椎间盘纤维环采用各向异 性纤维加强超弹性本构模型, 建立了包括椎体、椎间盘和韧带的正常有限元模型以 及L3-L4融合的腰椎L2-L4节段有限元模型, 并通过对比模拟和实验结果验证其有效 性. 利用商业有限元软件ABAQUS/Standard 进行前屈、后伸和轴向旋转载荷情况下 的模拟分析, 对比两种模型在不同状态下的运动范围及应力分布情况. 模拟结果显 示: 在相同的载荷情况下, 融合模型的运动范围和正常模型相比明显偏小; 相邻椎 体的应力分布与正常模型明显不同, 但植骨融合对相邻椎间盘的影响相对较小. 此 外, 小关节对维持脊椎正常生理功能起着重要作用, 小关节功能丧失会使相应节 段的椎间盘髓核压力增大.     

ACL单束/双束重建对多屈曲角位姿的胫股关节力学特性影响

进行膝关节的轴移仿真, 在证实了单例个体膝关节的重建模型能反映不同重建术后膝关节力学特性有效性的基础上, 应用SONATA MAESTRO 1.5T扫描得到的屈曲角位姿为0°/25°/60°/80° 的膝关节MRI图像数据, 建立了对应屈曲角的正常/ACL(anterior cruciate ligament)单束/双束重建胫股关节3D模型, 通过对该系列关节模型施加轴向压力与压扭载荷, 来分析多屈曲角度位姿的ACL单束/双束重建法对胫股关节软骨、半月板和韧带应力分布影响及韧带张力特性影响. 结果表明: (1) ACL单束/双束重建关节的软骨与半月板上的应力分布改变跟屈曲角度相关, 某些角度位姿下软骨应力分布改变显著; (2) 单束重建关节的软骨和半月板上的最大等效应力有较明显的增大, 最大增幅达40%左右; 双束重建关节在各屈曲角位姿下软骨和半月板上的最大等效应力值更接近于正常关节; (3) ACL单束重建虽降低了PCL(posterior cruciate ligament)的最大等效应力值, 却使高屈曲角位姿的内侧/外侧副韧带等效应力明显增大; 双束重建后MCL(medial collateral ligament)上的最大等效应力随屈曲角的变化明显, 但LCL(lateral collateral ligament)和PCL的最大等效应力随屈曲角度的变化趋势与正常关节相一致性; (4) 单束/双束重建后MCL上平均张力要高于另外两条韧带, 而双束重建的韧带张力特性比单束重建更接近于正常关节. 总之, 综合ACL重建后软骨、半月板和3条韧带随屈曲角度的等效应力分布、张力变化等多种特性表明: ACL双束重建的胫股关节力学特性比单束重建更接近于正常关节, 且无论ACL单束还是双束重建, 都引起术后关节软骨与韧带上应力分布变化与等效应力峰值增大, 这将是诱发术后关节慢性退变与膝关节并发症的根源.     

THE EFFECT OF THE ACL SINGLE-BUNDLE/DOUBLE-BUNDLE RECONSTRUCTION ON THE BIOMECHANICS OF TIBIO-FEMORAL JOINT AT DIFFERENT FLEXION ANGLES 1）

进行膝关节的轴移仿真,在证实了单例个体膝关节的重建模型能反映不同重建术后膝关节力学特性有效性的基础上,应用SONATA MAESTRO 1.5T扫描得到的屈曲角位姿为0°/25°/60°/80°的膝关节MRI图像数据,建立了对应屈曲角的正常/ACL（anterior cruciate ligament）单束/双束重建胫股关节3D模型,通过对该系列关节模型施加轴向压力与压扭载荷,来分析多屈曲角度位姿的ACL单束/双束重建法对胫股关节软骨、半月板和韧带应力分布影响及韧带张力特性影响.结果表明：（1）ACL单束/双束重建关节的软骨与半月板上的应力分布改变跟屈曲角度相关,某些角度位姿下软骨应力分布改变显著;（2）单束重建关节的软骨和半月板上的最大等效应力有较明显的增大,最大增幅达40%左右;双束重建关节在各屈曲角位姿下软骨和半月板上的最大等效应力值更接近于正常关节;（3）ACL单束重建虽降低了PCL（posterior cruciate ligament）的最大等效应力值,却使高屈曲角位姿的内侧/外侧副韧带等效应力明显增大;双束重建后MCL（medial collateralligament）上的最大等效应力随屈曲角的变化明显,但LCL（lateral collateral ligament）和PCL的最大等效应力随屈曲角度的变化趋势与正常关节相一致性;（4）单束/双束重建后MCL上平均张力要高于另外两条韧带,而双束重建的韧带张力特性比单束重建更接近于正常关节.总之,综合ACL重建后软骨、半月板和3条韧带随屈曲角度的等效应力分布、张力变化等多种特性表明：ACL双束重建的胫股关节力学特性比单束重建更接近于正常关节,且无论ACL单束还是双束重建,都引起术后关节软骨与韧带上应力分布变化与等效应力峰值增大,这将是诱发术后关节慢性退变与膝关节并发症的根源.     

浅谈动脉流体动力学分析

就我们所熟知,绝大部分正常动脉流,其血液的流动特性是属于层流范围,但随着弯曲和分 支部分会产生血液流之二次回流区,进而形成所谓近似非稳态流及紊流. 因此动脉流体的特 性会随动脉外形及条件的改变而改变. 在某些情形下,异常动脉的血液动力特性会造成动脉 的病变. 因此,近年来动脉血液流体的特性的研究,常着重于异常动脉的血液动力特性所形 成剪应力和病变部位动脉粥状硬化关系的探讨. 动脉血液流动经常包含分离流或二次回流运动,而这是流体力学的分析或数值模拟最困难的 部分. 有关分离流或二次回流的研究包括正常血管流和异常血管流,藉由二次回流的模拟与 测量可以观察血管病变的形成与演变,其中最受注目探讨题目是窄化血管如粥状斑块相关的 血液流动分析. 将回顾二维和三维、稳态、非稳态之动脉血流与窄化血管相关的几何外形作模拟研究和 实验. 并提供对血液动力学的研究方向,以作为未来医疗诊断与发展相关器材之参考.     

用于截瘫FES行走的三维上肢力动态测量系统

功能性电刺激(FES)技术已被广泛用于截瘫病人行走能力的恢复,但作为评估截瘫FES行走效率的重要指标,三维上肢力状况的研究至今仍十分有限,亟待深入开展。本文研制出了一套可用于截瘫FES行走的三维上肢力动态测量系统,该系统在一个标准步行器框架上装备有 12个应力应变片电桥,通过专用的隔离放大及A/D转换装置与计算机相连,经过标定后,截瘫FES行走时的上肢支撑力数据就可以被系统实时获得并显示。测力系统在应变片位置的确定以及“冗余-优化”方法的采用等多个方面保证了测量的可靠性和准确性,后续的校准实验显示系统的非线性、交互干扰和力学精度误差分别小于 2. 901%, 3. 188%和 1. 01%。研究中,一例截瘫患者作为志愿者按照正常的FES行走训练程序进行了相关测试,以验证新系统在行走能力评估方面的临床前景,结果表明其可用于: 1)截瘫FES行走效率的评估; 2)病人自我训练的控制指导; 3)作为反馈信号选择一种有效的FES模式和阵列。     

Photoelastic Analysis of the Effects of Tooth Position on Apical Stress

The anatomical form of the teeth and the density of the base bone differ greatly at different tooth positions. Thus, it is necessary to investigate the effect of tooth position variation on the apical occlusal stress because it would be helpful to understand the mechanism of occusal force transduction along the mandible. Epoxy resin models of seven different occlusal types set up with identical plastic teeth were employed to determine the stress directions and magnitudes in the apical area under 2.0 Kg of vertical load. Both the values and directions of the main stress were measured. Our results indicated that variation of the tooth position had no effect on the directions of principle stress around the apex area in the symmetrical occlusal type. For the asymmetrical occlusal type, the mesio-distal and the lingual-buccal directions of anterior teeth were found to be affected by variations in tooth position. The magnitude of the apical main stress was also apparently altered with changes of tooth position, but the differences were insignificant with occlusal types. Therefore, we conclude that variation in the tooth position could affect the direction of the main apical stress when the symmetry of the occlusion type was also involved. In contrast, the apical stress magnitude was much more related to tooth position than to occlusal symmetry.     

Quasi-Linear Viscoelastic Model of the Articular Disc of the Temporomandibular Joint

A precise characterization of the articular disc of the temporomandibular joint (TMJ) is essential to study the masticatory biomechanics. The disc is responsible for the load distribution over the articular surface and for absorbing impacts during mastication. The main objective of this work is to characterize the mechanical behaviour of the articular disc under compression, the usual stress state during mastication. A quasi-linear viscoelastic (QLV) model, with a hyperelastic response for the elastic function, is proposed to describe the mechanical behaviour of the articular disc. The validity of that simplified model relies on the independence of their constants with the strain level and strain rate. The independence of the strain level was proved in a previous work. In this paper, different loading rates were tested to fully confirm the validity of the model in the physiological range of loads. Moreover, the strong non-linearity of the stress-strain relation made the exponential strain energy function the most suitable of the different models tried to represent the elastic response of the QLV model.     

Interface Forces in Laterally Impacted Steel Tubes

In engineering mechanisms, solid bodies frequently come to a contact with a variety of geometric and kinematic situations. There has been a trend to express the interaction of the contacted bodies in the form of equivalent interface forces. The interface force represents the level of load transferred from one body to another. In a static or quasi-static contact problem, the interface forces could be sensibly evaluated by integrating the normal and the shear stresses over the common interfaces. In a dynamic contact, the interface stresses and subsequently the interface forces happen to be more complicated. They are, furthermore, affected by other parameters such as inertia forces, stress waves propagation, the material strain rate dependency and damping. This paper reports interface forces recorded in a series of experimental impact tests on axially pre-compressed steel tubes along with those from the numerical simulations of the tests. To monitor the so called impact loads, two small steel rings as load cells have been built in the striker. Based on the experimental and numerical results, the concept of equivalent interface forces in the impacted tubes has been verified. It has been highlighted that the interface forces (or the impact loads) may vary on the striker or the specimen themselves, depending on the measuring locations. Effects of axial compressions on the interface forces picked up by the striker load cells, impact loads imparted to the specimen supports and on the impact duration have been discussed. It has been reported that the initial compression applied to the tubes does not remain constant during and after the impact event. The amount of variations also depends on the initial level of the tube compression.     

髓核固定/滑动型人工颈椎间盘应力特征对比及关节面摩擦学优化设计

对比了髓核固定型和滑动型人工颈椎间盘应力和变形特征,根据应力和λ值(最小膜厚与等效粗糙度之比),采用正交试验法进行关节面摩擦学优化设计. 结果表明:两款球窝假体应力集中均发生在相对运动关节面和组件连接处,且髓核上表面应力最大值呈环形分布特征. 滑动型相对固定型,优势是运动模式更丰富、应力传递更均匀,劣势是复杂运动会加大髓核应力(从1.36 MPa到2.43 MPa). 固定型的髓核变形量仅是滑动型的22.5%,但总变形均集中在前后两端. 关节面参数对应力影响权重是球缺高度>球头半径>球窝间隙,且高度和半径越大、间隙越小越好;对润滑状态影响权重是球窝间隙>球头半径,且间隙越小、半径越大,润滑状态越佳. 半径、间隙和高度的最佳组合为16.00 mm、10.0 μm和1.2 mm,在此组合下滑动型假体最大应力数值最小(2.79 MPa),且润滑状态最佳(λ值0.573).      

作用在岩体裂隙网络中的渗透力分析

从岩体裂隙网络渗流的特点出发, 以单裂隙渗透力分析为基础, 分析了二维及三维情况下岩体裂隙网络渗流对岩体裂隙壁施加的两种作用力: 垂直于裂隙壁使裂隙产生扩容的法向渗透静水压力以及平行于裂隙壁和裂隙水流方向一致的切向动水压力, 推导出二维及三维情况下裂隙单元因这两种作用力而产生的等效结点力, 并应用算例定量分析了岩体裂隙网络渗透静水压力和动水压力共同作用对岩体应力的影响, 结果显示: (a)渗透力作用下裂隙上部岩体压应力减小, 而裂隙下部岩体压应力增大, 最大压应力增大 10 .5 3%; (b)渗透力作用下裂隙岩体拉应力增大, 最大拉应力增大 9.0 9%; (c)裂隙渗透力使岩体剪应力增大, 最大值达 2 3.75 %。     

Probability density functions of contact forces for cohesionless frictional granular materials

A theory is developed for the probability density functions of contact forces for cohesionless, frictional granular materials in quasi-static equilibrium. This theory is based on a maximum information entropy principle, with an expression for information entropy that is appropriate for granular materials. Entropy is maximized under the constraints of a prescribed stress and that the normal component of the contact force is compressive and that the tangential component of the contact force is limited by Coulomb friction. The theory results in a dependence of the probability density function for the tangential contact forces on the friction coefficient. The theoretical predictions are compared with results from discrete element simulations on isotropic, two-dimensional assemblies under hydrostatic stress. Good qualitative agreement is found for means and standard deviations of contact forces and the shape of the probability density functions, while the quantitative agreement is fairly good. Discrepancies between theory and simulations, such as the difference in shape of the probability density function for the normal force and the observed dependence on elastic properties of the exponential decay rate of tangential forces, are attributed to the fact that the method does not take into account any kinematics, which are essential in relation to elastic effects.     

Contact forces in anisotropic frictional granular materials

The probability density functions (PDFs) of contact forces in anisotropic, cohesionless and frictional granular materials are studied numerically and theoretically. Using discrete element simulations of biaxial deformation of a large two-dimensional assembly consisting of 200,000 disks, it is observed that the PDFs for the normal and tangential components of the contact forces depend significantly on contact orientation. The PDFs exhibit exponential decay and the PDF for the tangential component of the contact forces is not always symmetrical with respect to zero tangential force. The shape of the PDF for the normal component of the contact forces changes with shear strain. A qualitative explanation for this change is given that is related to the biaxial deformation mechanism in which the disrupted contacts are predominantly oriented in the direction of the minor principle stress.A maximum entropy method is employed to study these PDFs theoretically, using a prescribed stress tensor as constraint. It is found that the theoretical results correspond qualitatively to many of the results obtained from the discrete element simulations. Discrepancies between theory and simulations are attributed to the fact that the kinematics have not been taken into account in the theory.     

A General Formulation of Stress Distribution in Cylinders Subjected to Non-Uniform External Pressure

This work presents a two-dimensional stress analysis for elastic solid cylinders subjected to combined loading. The loading is generally formed with a number of concentrated and partially distributed forces all applied radially on the outer surface. The distributed forces cause pressures with non-uniform intensity along the circumferential direction. The cylinder is assumed to be long so that a state of plane-strain is valid. To obtain the stress distribution for the problem of partially distributed forces a new approach is followed first introduced in this paper. It is based on the expressions formed after using the theory of simple radial stress distribution when point-forces are applied on the cylinder and leads to the solution after direct integration. The total stresses due to both concentrated and distributed forces are obtained using the method of superposition. Apart from its simplified formulation, this general solution is always preferable since it proved to have a great advantage. As a result of not containing Fourier series, it eliminates some problems of convergence of the series at the boundaries that appear due to the Gibbs phenomena when the boundary conditions are a discontinuous function. Numerical results are presented for some interesting cases of loading conditions. This revised version was published online in August 2006 with corrections to the Cover Date.     

集中载荷下含偏心裂纹椭圆盘的应力强度因子计算

本文采用边界配置方法计算了受集中载荷作用时椭圆盘中偏心裂纹两端的应力强度因子,其中包括椭圆两半轴不同比值。不同裂纹长度和不同偏心程度的情况,在其特例椭圆盘中心裂纹情形,本文结果与Isida一致;在圆盘偏心裂纹情形,本文给出比Rooke等人更好的结果。