PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论,通过解析分析,导出了一种分析有限压电板裂纹问题的解析数值方法. 首先,计算了含中心裂纹有限板的断裂参数,与Woo和Wang的解析数值法(Int J Fract, 1993, 62: 203$\sim$218)相比较,表明该方法具有很高的精度和很好的计算效率. 随后,采用该方法和有限元法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数,发现各断裂参数的临界值分散性很大,不能作为压电材料的单参数断裂准则. 进而,针对试样真实的裂隙形状,采用有限元法计算了裂隙尖端的应力、电位移场,比较了裂隙内介质的介电性能对裂隙尖端场的影响,计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.      

高边坡结构可靠度的二次处理有限元分析

本文依据强度折减理论,利用MIDAS/GTS有限元软件,分析计算了高边坡结构的安全系数K,找到边坡滑裂带的位置。在此基础上,对有限元输入数据和输出结果进行二次处理,建立基本随机变量c,f与滑裂带中单元的最大(和最小)主应力σ1(和σ3)的拟合关系f1(和f3),将其代入高边坡结构的功能函数Z中,使Z由隐式形式变为显式。基于该显式表示的Z,利用Monte Carlo法计算滑裂带中所有失效单元的可靠指标β1,并将其单元面积A1作为权重系数,经过加权平均得到边坡结构的整体可靠指标β。上述方法使得结合有限元软件计算边坡结构的整体可靠度得以简化。经实例分析可知,本文提出的方法是合理可行的,可使边坡结构整体可靠性分析得以简化,也可为高边坡结构整体可靠性分析提供理论参考。     

PARALLEL FINITE ELEMENT ANALYSIS OF HIGH FREQUENCY VIBRATIONS OF QUARTZ CRYSTAL RESONATORS ON LINUX CLUSTER

Quartz crystal resonators are typical piezoelectric acoustic wave devices for frequency control applications with mechanical vibration frequency at the radio-frequency （RF） range. Precise analyses of the vibration and deformation are generally required in the resonator design and improvement process. The considerations include the presence of electrodes, mountings, bias fields such as temperature, initial stresses, and acceleration. Naturally, the finite element method is the only effective tool for such a coupled problem with multi-physics nature. The main challenge is the extremely large size of resulted linear equations. For this reason, we have been employing the Mindlin plate equations to reduce the computational difficulty. In addition, we have to utilize the parallel computing techniques on Linux clusters, which are widely available for academic and industrial applications nowadays, to improve the computing efficiency. The general principle of our research is to use open source software components and public domain technology to reduce cost for developers and users on a Linux cluster. We start with a mesh generator specifically for quartz crystal resonators of rectangular and circular types, and the Mindlin plate equations are implemented for the finite element analysis. Computing techniques like parallel processing, sparse matrix handling, and the latest eigenvalue extraction package are integrated into the program. It is clear from our computation that the combination of these algorithms and methods on a cluster can meet the memory requirement and reduce computing time significantly.     

????Timoshenko????????????

将理性有限元法引入到Timoshenko梁问题中,提出了一种理性Timoshenko梁单元,克服了 剪切锁死现象. 在推导控制方程时,与传统有限元方法采用Lagrange插值不同, 理性有限元法用Timoshenko梁弯曲问题的基本解逼近单元内部场. 运用该梁单元分析 Timoshenko梁时,无需缩减积分,就能避免剪切锁死,并且极大地提高了计算精度,说明 理性有限元法具有广泛的应用前景.     

基于动态本构模型的整流罩铰链系统自冲击分析

整流罩铰链系统结构分析是大型整流罩分离分析的重要环节,也是国内新一代运载火箭研制过程中的关键技术之一。基于30CrMnSiA和ZTC4-1的力学性能试验研究,得到材料的动态本构模型,并建立整流罩及动态本构模型下的铰链系统有限元模型。通过非线性显式动力学分析方法,对铰链系统自冲击过程进行仿真,得到冲击作用下铰链系统的应力和应变情况,并通过内力分析得到自冲击过程中铰链系统的主要内力形式。最后,通过比较铰链系统动态本构与准静态本构下的自冲击响应,验证并说明了动态本构模型对铰链系统自冲击仿真的影响。计算结果表明,自冲击作用会使铰链系统发生塑性变形,且动态本构模型有利于反映真实结果,为新一代运载火箭的研制提供一定的参考价值。     

Post-divergence and post-flutter analysis of sub-tangentially loaded and damped Beck column using a nonlinear FE procedure

Both post-divergence and post-flutter behaviors of damped Beck columns subjected to a sub-tangentially follower force are rigorously explored using an exact co-rotational frame element. First a linear stability theory of the damped Beck column is summarized using a stability map. A geometrically nonlinear frame element based on the co-rotational formulation is then formulated including mass matrix, Rayleigh damping matrix, and load-correction stiffness matrix due to circulatory forces. The dynamic FE analysis is performed using Newmark integration method. Finally a Beck column model is parametrically analyzed in order to investigate non-linear stability characteristics of the internally and externally damped non-conservative system. In particular, interesting nonlinear behaviors of Beck column quite different from those predicted by the linear theory are reported through static and dynamic nonlinear analysis with variation of sub-tangentiality of the follower force.     

Finite element analysis and experiment on high pressure apparatus with split cylinder

Ultra-high pressure belt-type die was designed with a large sample volume prism cavity and a split cylinder which was divided into eight segments to eliminate circumferential stress. The cylinder of this type die has no cambered surface on inner wall, and the inner hole is a hexagonal prism-type cavity. The divided bodies squeeze with each other, providing the massive support and lateral support effect of the cylinder. Simulation results indicate that the split cylinder with the prism cavity possesses much smaller stress and more uniform stress distribution. The split cylinder with the prism cavity has been shown to bear larger compressive stresses in radial, circumferential and axial directions due to its structure, and tungsten carbide is most effective in pure compression so this type cylinder could bear higher pressure. Experimental results prove that the high pressure apparatus with a prism-type cavity could bear higher pressure. The apparatus with a prism cavity could bear 52.2% more pressure than the belt-type die.     

“吕”字形巷道布置的卸压机理分析

利用有限元法分析了综放采场瓦斯巷和回风巷或运输巷采用“吕”字形布置后的卸压机理，得到了上部瓦斯巷宽度对下部回风巷或运输巷卸压效果的定量关系，为改革综放采场巷道布置并提高瓦斯巷抽放效果提供了新的设计思路和可靠依据。     

基于多孔介质理论的混凝土徐变力学行为的有限元分析

水分是影响混凝土徐变的最重要因素之一,但通常为了计算的简便往往忽略了水分的影响,造成计算结果有很大的误差.文章根据混合物理论结合罚参数法、Galerkin加权残值法以及徐变应力分析的隐式解法,推导得到了混凝土徐变两相多孔介质模型.建立了有限元计算模型,给出了求解算法并编制了有限元程序.计算和试验比较分析表明,该理论和计算方法得到的结果遵循徐变相关性质,说明了该理论和计算模型的合理性.最后,指出了模型中存在的不足以及进一步的研究方向.     

An FEM approximation for a fourth-order variational inequality of second kind

A fourth-order variational inequality of the second kind arising in a plate frictional bending problem is considered. By using regularization method, the original problem can be formulated as a differentiable variational equation, and the corresponding discrete FEM variational equation is presented afterwards. Abstract error estimates and error estimates of the approximation are derived in terms of energy norm and L^2-norm.