A High Order Bound Preserving Finite Difference Linear Scheme for Incompressible Flows

We propose a high order finite difference linear scheme combined with a high order bound preserving maximum-principle-preserving (MPP) flux limiter to solve the incompressible flow system. For such problem with highly oscillatory structure but not strong shocks, our approach seems to be less dissipative and much less costly than a WENO type scheme, and has high resolution due to a Hermite reconstruction. Spurious numerical oscillations can be controlled by the weak MPP flux limiter. Numerical tests are performed for the Vlasov-Poisson system, the 2D guiding-center model and the incompressible Euler system. The comparison between the linear and WENO type schemes, with and without the MPP flux limiter, will demonstrate the good performance of our proposed approach.     

Simulation of transcranial magnetic stimulation in head model with morphologically-realistic cortical neurons

BackgroundTranscranial magnetic stimulation (TMS) enables non-invasive modulation of brain activity with both clinical and research applications, but fundamental questions remain about the neural types and elements TMS activates and how stimulation parameters affect the neural response.ObjectiveTo develop a multi-scale computational model to quantify the effect of TMS parameters on the direct response of individual neurons.MethodsWe integrated morphologically-realistic neuronal models with TMS-induced electric fields computed in a finite element model of a human head to quantify the cortical response to TMS with several combinations of pulse waveforms and current directions.ResultsTMS activated with lowest intensity intracortical axonal terminations in the superficial gyral crown and lip regions. Layer 5 pyramidal cells had the lowest thresholds, but layer 2/3 pyramidal cells and inhibitory basket cells were also activated at most intensities. Direct activation of layers 1 and 6 was unlikely. Neural activation was largely driven by the field magnitude, rather than the field component normal to the cortical surface. Varying the induced current direction caused a waveform-dependent shift in the activation site and provided a potential mechanism for experimentally observed differences in thresholds and latencies of muscle responses.ConclusionsThis biophysically-based simulation provides a novel method to elucidate mechanisms and inform parameter selection of TMS and other cortical stimulation modalities. It also serves as a foundation for more detailed network models of the response to TMS, which may include endogenous activity, synaptic connectivity, inputs from intrinsic and extrinsic axonal projections, and corticofugal axons in white matter.     

Three-Dimensional Soft Tissue Prediction Using Finite Elements

Abstract Background and Aim: The prediction of soft tissue esthetics is important for achieving an optimal esthetic outcome in orthodontic treatment planning. Applicable procedures have so far been restricted to two-dimensional profile predictions that have not proven to be very reliable. The goal of this investigation was therefore to develop a novel finite element-based procedure that allows a three-dimensional, easily visualized, quantitative analysis and prediction of soft tissue behavior for the clinician. The procedure to be developed should be easy to handle and not entail any additional radiation exposure for the patient. Material and Methods: Using a three-dimensional scanner, the facial surfaces of 20 probands were digitalized and individual FEM models were generated. Results: After reduction of data redundancy via several conversion steps, a patient-specific simulation model was prepared consisting of 20,000 to 40,000 individual elements to which specific physical properties could be assigned. The average time required for generating a virtual model was 50 minutes. Problems occurring during model generation were rare (mainly shadowing phenomena and movement artifacts). Conclusion: The procedure outlined herein makes the reliable generation of patient-specific simulation models possible for facial soft tissue prediction in orthodontics.     

经Lister结节穿针多针阶梯形排列髓内弹性固定的3D有限元分析

目的：考察经Lister结节穿针多针阶梯形排列髓内弹性固定的力学效应。方法：采用CT扫描与计算机图象处理相结合的建模方法，进行尺桡骨受力的3D有限元分析，对比无克氏针固定、多针固定、传统针固定3种状态下，尺桡骨干在压缩和弯曲两种外载荷条件下的受力状态。结果：受压缩载荷作用时，桡骨在产生压缩变形的同时，横向（X轴方向）也产生一定的位移，且随着压缩载荷的增大，横向位移（绝对值）增加；尺骨压缩时横向位移很小。未植入克氏针与植入不同数量克氏针的尺桡骨压缩时，横向位移差异明显。传统针固定和多针固定时，受压缩载荷的尺骨受力较均匀，且等效应力低；而桡骨的等效应力分布不均匀，在桡曲的中间段内外侧应力明显集中；植入克氏针后骨干受力有所变化，尤其是尺骨的等效应力值明显下降。在受到横向的弯曲载荷时，尺桡骨同时产生横向弯曲，其中间部分的最大桡曲位移量随着弯曲载荷的增加而增大。在相同外载荷作用下，无克氏针固定和植入克氏针固定的尺桡骨的最大横向位移有一定差别，并且弯曲载荷越大，无克氏针固定、多针固定、传统针固定的差别越明显。克氏针固定后，其最大等效应力下降明显，尤以多针固定下降最明显。随着横向弯曲载荷的增加，尺桡骨横向弯曲变形越来越严重，外侧表现出纵向拉应力。植入克氏针尤其是多针后，桡骨内侧最大的压应力（绝对值）和外侧的最大拉应力值均明显降低，桡骨承载能力明显提高。结论：经Lister结节穿针多针阶梯形排列髓内弹性固定能有效分担外来载荷对尺桡骨的力学影响，提高固定的稳定性，维持桡骨的生理曲度，从而提高疗效。     

铰链式冠外附着体的结构强度分析

目的探讨铰链式冠外附着体设计的合理性.方法采用有限元应力分析的实验方法,分析铰链式冠外附着体阳性结构L形球体的应力分布情况,对其结构强度进行测试.结果在附着体阳性结构L形球体的转角处应力最大,应力值为321.847MPa.结论该铰链式附着体的结构强度符合临床设计要求.     

Elemental composition and water content of myelinated axons and glial cells in rat central nervous system

The distribution of elements (e.g. Na, Cl, K) and water in CNS cells is unknown. Therefore, electron probe X-ray microanalysis (EPMA) was used to measure water content and concentrations (mmol/kg dry or wet weight) of Na, Mg, P, S, Cl, K and Ca in morphological compartments of myelinated axons and glial cells from rat optic nerve and cervical spinal cord white matter. Axons in both CNS regions exhibited similar water content ( 90%), and relatively high concentrations (wet and dry weight) of K with low Na and Ca levels. The K content of axons was related to diameter, i.e. small axons in spinal cord and optic nerve had significantly less (25–50%) K than larger diameter axons from the same CNS region. The elemental composition of spinal cord mitochondria was similar to corresponding axoplasm, whereas the water content (75%) of these organelles was substantially lower than that of axoplasm. In glial cell cytoplasm of both CNS areas, P and K (wet and dry weight) were the most abundant elements and water content was approximately 75%. CNS myelin had predominantly high P levels and the lowest water content (33–55%) of any compartment measured. The results of this study demonstrate that each morphological compartment of CNS axons and glia exhibits a characteristic elemental composition and water content which might be related to the structure and function of that neuronal region.     

Experimentelle und klinische Studie zur Stabilit?t einer Osteosyntheseplatte für den atrophen Unterkiefer

ZusammenfassungHintergrund Der Anteil an älteren Patienten steigt ständig und damit auch die Zahl traumatischer altersbedingter Verletzungen wie Frakturen des Unterkiefers. Die Frakturversorgung bei älteren Menschen stellt spezielle Anforderungen. Durch das Design einer neu entwickelten Osteosyntheseplatte sollte versucht werden, diese speziellen Gesichtspunkte zu erfüllen.Material und Methoden Im Gegensatz zu den 2.0-Miniplatten (Medartis AG, Basel) besitzt die Pencilbone-2.0-Platte, die aus diesen Miniplatten entwickelt wurde, einen oval geformten Mittelsteg zur Stabilisierung des frakturnahen Knochens und zwei sphärische Gleitlöcher jenseits der Fraktur. An den verstärkten Teil der Osteosyntheseplatte schließen sich jeweils 2 bzw. 3 normale unverstärkte Löcher an, die sich sehr leicht an den frakturfernen Knochen adaptieren lassen. Dies ermöglicht dem Operateur den intraoralen Zugangsweg und eine Handhabung, welche er von der Miniplattenosteosynthese gewohnt ist, bietet aber gleichzeitig eine höhere Stabilität im Vergleich mit den Standard-2.0.-Miniplatten.Ergebnisse und Diskussion Nach experimentellen und biomechanischen Untersuchungen, die alle positive Ergebnisse zeigten, wurde die neue Platte zwischen Oktober 2000 und November 2001 in zwei Kliniken an 16 Frakturen des atrophischen Unterkiefers bei 14 Patienten erfolgreich angewendet. 15 Frakturen heilten primär, lediglich bei einer Fraktur wurde eine neue operative Versorgung wegen Knochendislokation nach erneutem Sturz notwendig.     

额、枕部直接冲击导致中脑损伤的力学机制研究

目的分析额、枕部直接冲击造成中脑应力响应,探讨颅脑直接冲击导致脑干损伤的力学机制。方法构建头颅三维有限元模型,模拟计算额、枕部直接冲击损伤,分析中脑上、下交界面节点应力分布,总结中脑应力响应。结果中脑-间脑界面在冲击初始阶段,2ms即出现应力增高,中脑-桥脑界面应力升高开始于冲击后期(10ms),两交界面节点应力在冲击时程中均持续升高,额部冲击比枕部冲击在中-间脑界面造成的应力峰值高。结论中脑在颅脑直接冲击伤的致伤机制中占有突出地位,中脑的应力改变可解释脑震荡的意识障碍及部分原发脑干伤的临床表现。     

冲击载荷作用下骨盆三维有限元分析及其生物力学意义

目的:探讨骨盆受到冲击载荷作用的力学行为特征,为临床分析及判断骨盆骨折类型、力学分布、冲击载荷影响提供力学基础。方法:采用计算机仿真模拟方法,将所构建骨盆三维实体模型导入三维有限元分析软件AN-SYS7.0,计算单侧髂前上棘和单侧髂骨正后方部位冲击载荷作用下骨盆的力学行为表现,冲击载荷为0～8000N,峰值8000N,时间自0～40ms,分析主应力值、应力分布情况以及主应力方向上骨盆单元的位移。结果:冲击载荷作用于单侧髂前上棘时,0～10ms主要的应力沿髂后上棘到髂前上棘分布,在冲击后期10ms以后,骶髂关节、髋臼和耻骨支都会产生较大的应力分布;作用于单侧髂骨正后方部位时,0～20ms主要的应力沿髂骨纵行传导,应力分布并没有在冲击峰值10ms时达到最大,而是从20ms之后,应力开始向骶髂关节、双侧耻骨支、坐骨支以及髋臼等部位传导并可见到明显的应力分布。结论:分析冲击载荷作用下骨盆各部位应力分布以及骨盆各个单元在应力作用下的位移变化,有助于临床上进行骨盆损伤内固定力点的选择以及进一步明确骨盆内在应力值分布。     

足骨三维有限元模型对足跖骨缺损重建的指导意义

目的用足骨有限元模型模拟、分析多跖骨缺损对足弓的影响,同时,对临床采用带血管自体骨重建后足的功能进行评价。方法在足骼骨三维模型上模拟出胫侧组、腓侧组跖骨正常和缺损不同情况,以有限元法计算足弓应力和位移的改变;并结合临床评价采用自体骨重建后足位移及应力恢复数据和结果。结果多跖骨缺损后,其有限元模型显示:足弓应力和位移发生明显改变,对足功能影响较大,髂骨重建后,足弓及足支撑点恢复,有限元模拟结果及临床随访效果满意。结论通过足骨有限元模型分析,我们认为:前足损伤致足多跖骨缺损临床上应予重建修复,恢复足的三点支撑和足弓。采用带血管游离自体骨加皮瓣组合移植是较为理想的手术方式。