Reliability of elasto-plastic structure using finite element method

A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases. The project supported by the Research Grant Council of Hong Kong (HKUST 722196E, 6039197E), the National Natural Science Foundation of China(59809003) and the Foundation of University Key Teacher by the Chinese Ministry of Education     

直杆碰撞刚性壁弹塑性动力后屈曲有限元分析

利用显式动力学有限元方法对直杆弹塑性动力后屈曲进行了分析,模拟了直杆轴向碰撞动力屈曲的变形及发展过程。分析中在直杆碰撞端局部临界屈曲长度范围内引入半正弦波形式的初始缺陷,计算结果与文献中的实验数据获得了很好的一致。分析结果表明,随着碰撞过程中所产生的应力波逐渐向前传播,后屈曲变形过程中所呈现的多个半波形式的高阶屈曲模态由初始具有单个半波形式的简单屈曲模态迅速演变而成。分析结果同时也揭示了直杆动力屈曲变形发展的机理,以及轴向应力波和屈曲变形的相互作用规律。     

Pulse shaping techniques for testing elastic-plastic materials with a split Hopkinson pressure bar

We present pulse shaping techniques to obtain compressive stress-strain data for elastic-plastic materials with a split Hopkinson pressure bar. The conventional split Hopkinson pressure bar apparatus is modified by placing a combination of copper and steel pulse shapers on the impact surface of the incident bar. After impact by the striker bar, the copper-steel pulse shaper deforms plastically and spreads the pulse in the incident bar so that the sample is nearly in dynamic stress equilibrium and has a nearly constant strain rate in the plastic response region. We present analytical models and data that show a broad range of incident strain pulses can be obtained by varying the pulse shaper geometry and striking velocity. For an application, we present compressive stress-strain data for 4340 R_c 43 steel.     

Exact integration of constitutive equations of kinematic hardening materials and its extended applications

In this paper, an exact formula for the integration of the constitutive equations of kinematic hardening material is presented. Its algorithms are simple and clear. For isotropic hardening or mixed hardening material, the formula is still an exact solution for the case of radial loading, and it is an approximate solution with reasonable accuracy for the case of non-radial loading. The computation results show that the procedure proposed in this paper improves both accuracy and efficiency of numerical integration schemes adopted widely in elastic-plastic finite element analysis.     

岩石弹塑性破裂过程的数值模拟研究

分析了对于非均质、非连续且材料参数分布具有分散性和随机性的岩体介质，采用传统有限元法解决其工程问题的局限性．阐述了用岩石力学性质的分布特性来描述岩石材料不均匀性，并用蒙特卡洛模拟方法，将符合Weibull分布规律的材料参数赋予岩体结构中的各个单元体，进而用计算机进行数值模拟的研究方法．应用了在Windows95环境下利用FORTRAN Power Station4．0平台开发的具有可视性的“岩石弹塑性破裂过程分析”软件——REPFPA软件，对典型的岩石试样和巷道模型的弹塑性破裂过程进行数值模拟研究，得到了有价值的实验结果．证明把材料参数随机赋值方法引入到岩土类材料的弹塑性破裂过程分析中去，并用计算机进行数值模拟是研究岩石弹塑性破裂过程的一种有效方法。     

内爆加载金属界面不稳定性的数值分析

采用自行研制的多介质弹塑性流体力学欧拉程序,对柱形内爆加载金属界面不稳定性进行了数值研究,数值模拟结果与文献实验数据吻合较好。数值结果表明:材料强度对界面不稳定性发展有不可忽略的抑制作用;材料屈服强度对较高模数不稳定性增长的抑制较强,而剪切模量对不稳定性发展的影响相似但敏感性相对较弱;金属界面不稳定性增长存在最不稳定模数,最不稳定模数随屈服强度增加而减小,并近似与屈服强度的对数呈线性关系;随着壳的厚度减小,扰动增长加快。     

基于新弹塑性模型全程模拟金属疲劳失效

金属材料疲劳失效问题是工程中一类主要问题,选择合理的材料本构关系对于准确确定疲劳失效强度至关重要．本文提出一个简单形式的自由光滑新弹塑性模型,该模型在不涉及通常的屈服条件以及加卸载条件的意义上完全自由,避免了模量间断问题,且比经典模型更简单、更符合实际．结果表明,该模型可直接模拟金属材料在循环加载下直至疲劳破坏的全过程,不涉及任何损伤变量以及人为假设的失效判据,特别地,可直接模拟金属材料的高周、低周疲劳直至最终失效破坏行为．应用P92钢的相关实验数据,给出了循环加载下的数值模拟结果,结果表明,在各种加载-卸载循环情形下,随着应力幅值的增大,疲劳失效循环数减少,这与实际材料行为相一致．     

基于数值子结构方法的结构弹塑性分析

结合震害调研及数值分析可知,结构最终失效可能仅由部分关键构件破坏引起,大部分构件仍处于弹性或小变形状态。因此为提高计算效率,在结构全过程分析中一致采用非线性单元建模并非必要,同时为准确考虑关键构件的非线性响应,本文提出一种新的数值子结构建模策略。进入弹塑性状态后,针对一般钢构件或钢筋混凝土构件采用动态替换子结构方法在单元或截面层次将其替换成非线性单元或非线性截面,并基于OpenSees平台开发了两类新单元予以实现;针对可能发生严重损伤的关键构件,采用隔离子结构方法将其隔离并建立精细化分析模型,考虑主、子结构间不同尺度边界耦合,并推导了切线刚度的传递关系,采用Client/Server技术在OpenSees平台开发了一类新的接口单元予以实现主、子结构之间的信息传递。为验证新开发单元的合理性,分别以钢及钢筋混凝土平面框架结构为例,采用纤维单元、动态替换子结构方法以及隔离子结构方法建模进行静、动力分析。计算结果表明,采用本文提出的动态替换子结构方法与常规建模方法的计算结果完全吻合并且可大幅缩短计算耗时,随着荷载水平的增大,结构中受到动态替换的构件比例急剧增大,计算效率提高程度略有降低,但仍远高于常规模型;采用本文提出的接口单元可准确传递主、子结构间的界面信息,为隔离数值子结构方法在结构弹塑性分析中的应用提供了基础。     

Elastic-plastic deformation of sandwich rod on elastic basis

Sandwich composite material possesses advantages of both light weight and high strength. Although the mechanical behaviors of sandwich composite material with the influence of single external environment have been intensively studied, little work has been done in the study of mechanical property, in view of the nonlinear behavior of sandwich composites in the complicated external environments. In this paper, the problem about the bending of the three-layer elastic-plastic rod located on the elastic base, with a compressibly physical nonlinear core, has been studied. The mechanical response of the designed three-layer elements consisting of two bearing layers and a core has been examined. The complicated problem about curving of the three-layer rod located on the elastic base has been solved. The convergence of the proposed method of elastic solutions is examined to convince that the solution is acceptable. The calculated results indicate that the plasticity and physical nonlinearity of materials have a great influence on the deformation of the sandwich rod on the elastic basis.     

Improved nonlinear fluid model in rotating flow

The pseudoplastic circular Couette flow (CCF) in annuli is investigated. The viscosity is dependent on the shear rate that directly affects the conservation equations solved by the spectral method in the present study. The pseudoplastic model adopted here is shown to be the suitable representative of nonlinear fluids. Unlike the previous studies, where only the square of the shear rate term in the viscosity expression is considered to ease the numerical manipulations, in the present study, the term containing the quadratic power is also taken into account. The curved streamlines of the CCF can cause the centrifugal instability leading to toroidal vortices, known as the Taylor vortices. It is further found that the critical Taylor number becomes lower as the pseudoplastic effect increases. The comparison with the existing measurements on the pseudoplastic CCF results in good agreement.