PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

Stress Distributions Created by Short and Regular Implants Placed in the Anterior Maxilla at Different Angles: A Finite Element Analysis

Although short implants are seen as alternative treatments that require additional surgical techniques in posterior region, they can be applied to anterior maxilla and various studies are required on this subject. The purpose of this study was to examine and compare the peak von Mises stress distributions in the crown, implant and abutment by using finite element analysis (FEA). Besides, a comparison of the implant-abutment connection types in the short implant with the FEA method was established. A short implant (4 × 5 mm) with a taper-lock connection and a regular implant (4 × 9 mm) with a screw connection were used in maxillary central incisor tooth area. Three different titanium abutments with 0?, 15? and 25? angles were used for abutments. In addition, in order to determine whether the stress change in short implants is due to the length of the implant-abutment connection, a screw was designed for a short implant and it was also evaluated in the same three angles. A total of three groups and nine models were generated. 114.6N load was applied to the cingulum area of the crown at an angle of 135? to the long axis of the crowns. A torque load of 25 Ncm was applied to the regular and short implant screw. Von Mises stress distributions of implants, abutments and crowns were evaluated by using FEA. Increased angle in implants increased von Mises stress values of implant, abutment and crown. Screw connection was found higher at all angles in short implants. Close values were found at different angles in taper-lock short implant crowns. The length and the angle in the bone of implant with the type of implant-abutment connection results in the accumulated stress values. Clinical Implications Taper implant-abutment connection system was found to be more promising in terms of stress accumulation in crowns. Although the amount of stress on the abutment increased due to the length of the implant in short implants, taper implant-abutment connection system slightly reduced related to this increase.

The finite element analysis of cemented long and short stem prosthetic replacement in aged patients with comminuted intertrochanteric fracture

Objective To investigate the stress distribution of the femur after cemented prosthetic replacement in aged patients with comminuted intertrochanteric fracture and to analyze the difference of stress distribution between cemented long

Damage prediction for magnesium matrix composites formed by liquid-solid extrusion process based on finite element simulation

A damage prediction method based on FE simulation was proposed to predict the occurrence of hot shortness crocks and surface cracks in liquid-solid extrusion process. This method integrated the critical temperature criterion and Cockcroft ＆ Latham ductile damage model, which were used to predict the initiation of hot shortness cracks and surface cracks of products, respectively. A coupling simulation of deformation with heat transfer as well as ductile damage was carried out to investigate the effect of extrusion temperature and extrusion speed on the damage behavior of Csf/AZ91D composites. It is concluded that the semisolid zone moves gradually toward deformation zone with the punch descending. The amplitude of the temperature rise at the exit of die from the initial billet temperature increases with the increase of extrusion speed during steady-state extrusion at a given punch displacement. In order to prevent the surface temperature of products beyond the incipient melting temperature of composites, the critical extrusion speed is decreased with the increase of extrusion temperature, otherwise the hot shortness cracks will occur. The maximum damage values increase with increasing extrusion speed or extrusion temperature. Theoretical results obtained by the Deform^TM-2D simulation agree well with the experiments.

Parametric investigation of railway fastenings into the formation and mitigation of short pitch corrugation

Short pitch corrugation has been a problem for railways worldwide over one century.In this paper,a parametric investigation of fastenings is conducted to understand the corrugation formation mechanism and gain insights into corrugation mitigation.A three-dimensional finite element vehicle-track dynamic interaction model is employed,which considers the coupling between the structural dynamics and the contact mechanics,while the damage mechanism is assumed to be differential wear.Various fastening models with different configurations,boundary conditions,and parameters of stiffness and damping are built up and analysed.These models may represent different service stages of fastenings in the field.Besides,the effect of train speeds on corrugation features is studied.The results indicate:(1)Fastening parameters and modelling play an important role in corrugation formation.(2)The fastening longitudinal constraint to the rail is the major factor that determines the corrugation formation.The fastening vertical and lateral constraints influence corrugation features in terms of spatial distribution and wavelength components.(3)The strengthening of fastening constraints in the longitudinal dimension helps to mitigate corrugation.Meanwhile,the inner fastening constraint in the lateral direction is necessary for corrugation alleviation.(4)The increase in fastening longitudinal stiffness and damping can reduce the vibration amplitudes of longitudinal compression modes and thus reduce the track corrugation propensity.The simulation in this work can well explain the field corrugation in terms of the occurrence possibility and major wavelength components.It can also explain the field data with respect to the small variation between the corrugation wavelength and train speed,which is caused by frequency selection and jump between rail longitudinal compression modes.

Dynamic ModeⅡfracture behavior of rocks under hydrostatic pressure using the short core in compression(SCC)method

The shear failure of rocks under both a static triaxial stress and a dynamic disturbance is common in deep underground engineering and it is therefore essential for the design of underground engineering to quantitively estimate the dynamic ModeⅡfracture toughness KⅡCof rocks under a triaxial stress state.However,the method for determining the dynamic KⅡCof rocks under a triaxial stress has not been developed yet.With an optimal sample preparation,the short core in compression(SCC)method was designed and verified in this study to measure the dynamic KⅡCof Fangshan marble(FM)subjected to different hydrostatic pressures through a triaxial dynamic testing system.The formula for calculating the dynamic KⅡCof the rock SCC specimen under hydrostatic pressures was obtained by using the finite element method in combination with secondary cracks.The experimental results indicate that the failure mode of the rock SCC specimen under a hydrostatic pressure is the shear fracture and the KⅡCof FM increases as the loading rate.In addition,at a given loading rate the dynamic rock KⅡCis barely affected by hydrostatic pressures.Another important observation is that the dynamic fracture energy of FM enhances with loading rates and hydrostatic pressures.

Experimental Study on the Axial Compression Behavior of Short Columns of Steel-Fiber-Reinforced Recycled Aggregate Concrete

In order to study the axial compression performances of short columns made of recycled aggregate concrete,four samples were designed with different recycled aggregate replacement rates and carbon fibre reinforced plastics(CFRP)sheets.Then,monotonic loading was implemented to assess the variation trends of their axial compression properties.The ABAQUS finite element software was also used to determined the compression performances.Good agreement between experimental and numerical results has been found for the different parameters being considered.As shown by the results,recycled coarse aggregates result in improved ductility and better deformation performance of the specimens.The failure of specimens caused by pre damage can be compensated by using CFRP sheets,by which both the resistance to deformation and the axial carrying capacity of columns can be increased.

The Young's moduli prediction of random distributed short-fiber-reinforced polypropylene foams using finite element method

The elastic moduli of short-fiber-reinforced foams depend critically on the fiber content and fiber length, as well as on the fiber orientation distribution. Based on periodic tetrakaidecahedrons, the finite element models with short-fiber reinforcement were proposed in this paper to examine the effects of the fiber content and fiber length on Young's modulus. The fiber length distribution and fiber orientation distribution were also considered. The proposed models featured in a three-dimensional diorama with random short-fiber distribution within or on the surfaces of the walls and edges of the closed-cells of polypropylene (PP) foams. The fiber length/orientation distributions were modeled by Gaussian prob-ability density functions. Different fiber volume fractions, different lengths, and different distributions were investigated. The predicted Young's moduli of the PP foams with short-glass-fiber or short-carbon-fiber reinforcement were compared with other theoretic and experimental results, and the agreement was found to be satisfactory. The proposed finite element models were proved to be acceptable to predict the Young's moduli of the grafted closed-cell PP foams with short-fiber reinforcement.

COMPUTER SIMULATION OF CREEP DAMAGE AT CRACK TIP IN SHORT FIBRE COMPOSITES

Creep damage at crack tip in short fibre composites has been sim- ulated by using the finite element method(FEM).The well-known Schapery non- linear viscoelastic constitutive relationship was used to characterize time-dependent behaviour of the material.A modified recurrence equation was adopted to accelerate the iteration.Kachanov-Rabotnov's damage evolution law was employed.The growth of the damage zone with time around the crack tip was calculated and the results were shown with the so-called 'digit photo',which was produced by the printer.

The Influence of Thrust Fault Structure on Cross-fault Short-leveling Survey

Aiming at different variation patterns of cross-fault short-leveling before earthquakes,the paper establishes the 2-D finite-element models with different crustal stratification and fault occurrence perpendicular to Longmenshan fault zone. By contact analysis and viscoelastic finite element method,the influence of fault structure on cross-fault short-leveling is obtained under the same constraint conditions,the results show that:with the increase of the horizontal projection distance of fault,the cumulative displacements of surface increase gradually in the models with fixed dip angles of the fault plane(model 1). However,when the horizontal projection distance exceeds 20 km,the influence of fault?s dip angle on the cumulative displacements of surface short-leveling will decrease significantly,and the cumulative displacements are maintained at about 1. 5 m. However,in the listric fault models(model 2),when the projection distance is less than 20 km,the listric fault structure impedes the sliding of fault. The short-leveling variation rate is only half of model 1;as a result,the ability to reflect the regional stress enhancement by cross-fault short-leveling is further weakened. But when the horizontal projection distance exceeds25 km,the cumulative displacements significantly increase,with the maximum displacement reaching 1. 75 m. The results of equivalent stress show that the listric fault structure causes a sudden increasement in equivalent stress when the horizontal projection distance is 10 km,higher equivalent stress values are accumulated between the projection distance of 5-20 km,and then high-low stress difference zones are formed at the bottom of the fault plane and near the transition zone of low-high dip angle.

Calculation of Eddy Current Loss and Short-circuit Force in SSZ11-50000/110 Power Transformer

The ?method is used in this paper to calculate the leakage magnetic field of SSZ11-50000/110 Power transformer, and by which the structures’ influences to the main leakage flux are analyzed. Through the combination of the product and TEAM Problem 21B, the surface impedance method shows its great advantage in the calculation of eddy current loss.

Short-Crested Waves Interaction with A Concentric Porous Cylinder System with Partially Porous Outer Cylinder

In this paper, based on the linear wave theory, the interaction of short-crested waves with a concentric dual cylindrical system with a partially porous outer cylinder is studied by using the scaled boundary finite element method （SBFEM）, which is a novel semi-analytical method with the advantages of combining the finite element method （FEM） with the boundary element method （BEM）. The whole solution domain is divided into one unbounded sub-domain and one bounded sub-domain by the exterior cylinder. By weakening the governing differential equation in the circumferential direction, the SBFEM equations for both domains can be solved analytically in the radial direction. Only the boundary on the circumference of the exterior porous cylinder is discretized with curved surface finite elements. Meanwhile, by introducing a variable porous-effect parameter G, non-homogeneous materials caused by the complex configuration of the exterior cylinder are modeled without additional efforts. Comparisons clearly demonstrate the excellent accuracy and computational efficiency associated with the present SBFEM. The effects of the wide range wave parameters and the structure configuration are examined. This parametric study will help determine the various hydrodynamic effects of the concentric porous cylindrical structure.

新型同步调相机匝间短路故障机理分析

随着新型同步调相机大规模投入运行,调相机的安全问题逐渐凸显。绕组匝间短路故障是调相机1种典型故障,对调相机长期运行和电网稳定都有严重影响。该文从两种故障发生后对气隙磁场磁动势的影响出发,建立转子总电流、定子绕组相电流、支路电流及电磁转矩的统一数学表征,从幅值、相位、谐波含量等多方面对两种故障进行深入的对比分析,揭示新型同步调相机绕组匝间短路故障机理。该文基于已投入使用的TTS-300-2型新型同步调相机,开发有限元软件的仿真程序,搭建调相机模拟机的实验平台。理论分析结果分别得到了仿真与实验结果的验证。

更高速下轮轨瞬态“滚-滑-跳”接触行为及中/短波不平顺临界限值研究

简述轮轨滚动接触行为预测方法现状的基础上,选择最适于400 km/h及以上更高速和中/短波不平顺激励下轮轨“滚-滑-跳”接触行为分析的瞬态滚动接触显式有限元建模方法,建立更高速轮轨三维瞬态滚动接触时域分析模型,数值模拟了轮对在典型曲线轨道上的瞬态曲线通过行为。根据现营高速轮轨的中/短波不平顺调研,以波长30~210 mm钢轨波磨为典型不平顺,考虑波磨激励下左右两侧轮轨动力作用相互影响,分析波磨几何等因素对速度500 km/h及以下轮轨瞬态“滚-滑-跳”接触行为的影响。考虑低、中频动力减载等因素的前提下,仅从轮轨接触脱离角度,提出轮轨中/短波不平顺的临界管理限值建议,讨论其重要性、合理性和应用局限,并与《高速铁路钢轨打磨管理办法》规定的波磨整治限度进行对比。

基于端部漏磁场分析的集成式电机驱动系统在线匝间短路故障诊断

基于宽禁带器件的变频器具有高功率密度和耐高温工作的优势,通过小型化和轻薄化设计可以将变频器置于电机端盖与端绕组之间,构成结构紧凑的集成式电机驱动系统。该系统发生匝间短路故障时,短路绕组电流增大不仅会使电机内部温度升高导致电机绝缘性能降低,还将引起磁场变化影响变频器的正常工作。以一台1.5 kW的交流感应电机为例,搭建基于宽禁带器件的集成式电机驱动系统三维物理模型,利用数学解析法分析并计算匝间短路故障前后的端部漏磁场,并提取电机内端部印制电路板(printed circuit board,PCB)表面与轴心距离相等、机械角度相差45°的两点处磁感应强度在直角坐标系下椭圆曲线的离心率反映端部漏磁场的对称性,以此为特征指标判断电机是否出现匝间短路故障。仿真与实验结果验证了所提方法的正确性和有效性。

高强型钢超高性能混凝土短柱轴压性能试验及其有限元分析

为研究高强型钢超高性能混凝土短柱的轴压性能,对试件进行了轴心受压试验,设计参数包括配钢率、箍筋间距和箍筋形式。观察了试件的轴压破坏过程,获得了破坏形态及轴向荷载-位移曲线,分析了承载能力、变形能力、刚度等轴压性能指标,以及超高性能混凝土、纵筋、箍筋和型钢的应变发展规律。在试验研究基础上,采用ABAQUS建立了高强型钢超高性能混凝土短柱轴压性能有限元分析模型,计算结果与试验结果吻合较好,进而进行了参数分析。结果表明:试件发生的是典型轴压破坏,中部表面出现“锯齿形”裂缝;轴向荷载-位移曲线会出现2次荷载峰值;箍筋间距减小时,试件的承载能力和变形能力均提高,但刚度退化速率没有显著变化;配钢率增大时,试件的承载能力提高,变形能力先增大后减小,刚度退化变缓;型钢强度增大时,试件的承载能力和变形能力均提高,提高速度先快后慢;超高性能混凝土抗压强度增大时,试件的承载能力提高,变形能力降低;超高性能混凝土受拉能够达到峰值应变;纵筋在极限点之前发生受压屈服;箍筋在极限点时应变不足屈服应变的1/3,但在破坏点前达到受拉屈服;型钢翼缘在极限点前后发生屈服,腹板屈服晚于翼缘。

圆钢管超高性能混凝土短柱偏压力学性能研究

为研究钢管超高性能混凝土短柱的偏压性能,以荷载偏心率和钢管径厚比为变化参数,设计了12个圆钢管UHPC短柱试件并对其进行偏心受压加载试验,分析了该类构件的破坏模式、荷载-挠度曲线、钢管应变和变形系数等,考察了主要因素对短柱偏压性能的影响规律。结果表明:试件的破坏特征为钢管屈服和核心混凝土压坏;荷载-挠度曲线有较明显的峰值点,偏心率越大和径厚比越小,曲线的下降段越平缓;达到60%峰值荷载时,钢管开始产生明显的约束作用,达到90%峰值荷载时,截面变形不再符合平截面假定;偏心率增大使试件的承载力和刚度下降,而径厚比减小可降低这种不利影响。在试验基础上,结合数值模拟对短柱的N-M曲线进行分析,建立了临界偏心率和套箍系数的关系表达式,并基于此提出短柱偏压承载力实用计算方法,理论与试验结果吻合较好。

预应力薄壁圆钢管约束混凝土短柱轴压性能初探

为提高薄壁钢管混凝土柱轴心受压性能,提出了预应力薄壁圆钢管约束混凝土柱,即采用对拉螺栓对薄壁圆钢管施加环向预应力,形成对内部混凝土初始约束围压,以提高内部混凝土的轴心抗压强度及承载力性能。为探究其轴向受压性能提升效果,进行了预应力薄壁圆钢管约束混凝土短柱轴心受压试验;并采用ABAQUS有限元计算软件建立了相应的非线性有限元数值模型,基于试验所获得荷载-位移数据,验证了模型的可靠性。在合理可靠试验模型的基础上,进行了不同预应力度(k=0、0.1、0.2、0.3、0.4和0.5)和不同钢管厚度(t=1.5mm、2.0mm、2.5mm和3.0mm)情况下的轴心受压性能影响分析。结果表明:采用对拉螺栓对薄壁圆钢管混凝土柱施加环向预应力,可有效提高其轴心抗压强度、刚度和延性;在合适预应力范围内,预应力度(k=0~0.4)越大,预应力薄壁钢约束混凝土短柱轴心受压性能提升越明显;但较大的预应力度(k=0.5)会让构件在达到峰值荷载后,其承载力降低较为明显,从而使构件延性降低;当预应力度k不变,构件的轴心受压性能随钢管厚度t的增大而增大。

两种方案设计下支架材料对无牙上颌种植固定修复影响的有限元分析

背景:在上颌无牙颌种植体支持式固定修复治疗中,上部支架结构材料的选择及不同远端种植体植入方式的设计密切影响着全口种植修复体的长期稳定性。目的:采用三维有限元方法,综合研究3种材料的上部支架及2种固定种植设计方案对上颌种植固定修复生物力学的影响。方法:根据1名健康成人正常颌的锥形束CT资料,采用Mimics软件分离上颌骨及上颌牙列三维实体模型,利用Geomagic Studio软件结合具体的模型参数构建无牙上颌全牙弓种植固定义齿的三维有限元模型。依据上颌后牙区远中种植体植入方式设计的不同,建立两种方案模型:方案1,按照临床常规应用的“All-on-4”设计,即2个种植体垂直植入上颌骨的双侧侧切牙区,另2个种植体则倾斜30°植入到双侧第二前磨牙区;方案2,即2个种植体垂直植入上颌骨的侧切牙区,后牙区采用2个短植体垂直植入双侧第二前磨牙区;在两种方案设计中使用3种材料(钛、氧化锆及聚醚醚酮)分别对上部支架结构进行赋值,得到6个模型,随后对比分析在斜向加载力方向下种植体、周围骨组织及上部支架结构的生物力学效果。结果与结论:①无论采用何种材料的上部支架,在两种无牙颌种植固定修复方案设计下,每个模型的最大应力峰值均分布于后侧种植体周围的颈部区域和皮质骨处;②与方案2中采用短种植体垂直植入的替代设计相比,方案1在上颌骨上显示出更为理想的应力分布;③由聚醚醚酮材料构建的支架模型,向上颌骨接近载荷区的种植体和周围骨组织传递了更高的应力,其次是钛、氧化锆;对于支架本身,聚醚醚酮上部支架内部应力峰值均明显小于氧化锆和钛支架;在3种上部支架材料中,选用氧化锆支架更有利于分散种植体及骨组织的应力分布;④结果提示,两种无牙颌种植固定修复设计方案都可以应用于临床实践,但从生物力学角度上考虑,方案1设计中的种植体、周围骨组织及上部支架的应力分散较为理性,更有利于无牙颌患者固定种植修复的远期预后;上部支架材料对种植体-骨界面的应力分布有着一定的影响。