SBFEM与非局部宏微观损伤模型相结合的准脆性材料开裂模拟

混凝土是一种被广泛应用于土木和水利工程中的准脆性材料, 在各种内外部因素的作用下, 开裂是混凝土结构最为普遍的破坏形式, 准确模拟结构的开裂过程, 对于结构的安全评估至关重要. 将比例边界有限元与非局部宏微观损伤模型相结合提出一种准脆性材料开裂模拟新方法. 以比例边界有限元子域的比例中心作为物质点, 通过两比例中心(物质点对)之间的物质键的正伸长率来定义微细观损伤, 将某点影响域内物质键的微细观损伤加权平均得到该点的宏观拓扑损伤. 再引入能量退化函数, 将宏观拓扑损伤嵌入到比例边界有限元的基本框架中. 充分利用比例边界有限元网格允许存在悬挂节点的优势, 采用四叉树网格离散技术进行快速、高质量的多级网格划分与过渡. 通过一个I型开裂与一个混合型开裂的两个典型算例, 验证了该方法可捕获结构裂纹扩展路径与荷载变形曲线. 与现有的方法相比, 本文的损伤模型可得到更准确的局部开裂损伤带, 结果更为合理, 且具有更高的计算精度和计算效率. 当损伤过程区网格尺寸小于影响域半径的1/5时, 计算结果不存在网格敏感性问题.      

一种模拟动边界绕流的锐利界面浸入边界法

When the structural wall moves over a fixed grid, the structure coverage will change, resulting in many dead and emerging elements. To avoid the influence of malformation and reconstruction of body-fitted grids on the calculation efficiency and accuracy of the fluid-structure interaction problems with coupled boundary movement on the fixed grid, an improved numerical method for describing the interaction between an immersed rigid body and fluid based on a sharp-interface is proposed. In this method, both the fluid and solid are regarded as pure fluid domains in the whole computational domain, and the solid boundary is divided into several Lagrangian grid points. The flow parameter or velocity is reconstructed by interpolation at the interface element, which is then directly used as the boundary condition of the flow field, thus reflecting the influence of the wall boundary conditions. The method constructs the calculation structure of “virtual point, force point and vertical foot point”, and the velocity of the virtual point is obtained by bilinear interpolation. Then, the velocity of the force point is calculated by forcing the solid boundary to meet the no-slip condition, and the equations of the coupling system based on the immersion boundary method are finally solved to realize the numerical simulation of the flow with a complex moving boundary. The numerical program for this immersed boundary method is established using C++, then the accuracy and reliability of the proposed method are validated by comparison with the literature and experimental results of the basic numerical example of flow around a cylinder. Furthermore, the effects of the structural shape and the angle of attack on the trailing vortex structure, the vortex shedding frequency, and the lift/ coefficient characteristics of the flow around the elliptical cylinder have been analyzed. The anti-symmetric S-type, “P+S” Ⅰ-type and “P+S” Ⅱ-type trailing vortex shedding modes, as well as the variation laws of the vortex structure size, vortex shedding frequency and lift-drag coefficients ratio with axis ratio and angle of attack, are captured. The critical angle of attack (25°) corresponding to the maximum lift-drag ratio is determined as 25°.     

集中载荷作用下不同边界条件的弯曲厚矩形板的受迫振动

运用边界积分法研究了四边简支、两对边固定另两对边简支、四边固定三种复杂边界条件下厚矩形板的受迫振动问题,求解过程清晰,从而给出了受迫振动控制方程和挠曲面方程。通过在Matlab平台上进行数值计算,得出了图表形式的计算结果,并与有限元模拟值进行对照。研究表明,边界积分法用于求解厚矩形板的受迫振动问题的准确性,本文推导的控制方程和挠曲面方程的正确性,进而对工程实际中的各种相关问题具有一定的现实意义,也为求解此类问题提供了一种新途径,可以直接运用到工程实际中。     

晶粒尺寸对Pr0.5(Ca0.47Ba0.03)MnO3中台阶效应的影响

采用溶胶-凝胶法制备了不同颗粒尺寸的Pr0.5(Ca0.47Ba0.03)MnO3多晶样品.对样品进行电阻测量发现,随颗粒边界的减小,其电阻率随磁场的变化即ρ(H)曲线上台阶数目和位置发生变化;实验同时发现,热循环对样品低温下的基态也有很大的影响.实验结果表明,颗粒尺寸大小,可以有效地调节反铁磁相(AFM)和铁磁相(FM)之间的相互作用,从而对台阶效应产生影响.进一步证明了台阶效应与相应力有关.     

氢键在分子自组装中的花样

介绍氢键的基本花样类型。由氢键组装的复杂结构大都是由4种基本类型构成。对醇类、羧酸、酰胺、α-氨基酸分别研究,阐述了不同剩余集团对氢键花样的影响。简单介绍了氢键作用在超分子化学中的作用。     

单个球形生物质颗粒热解过程的数值模拟

本文对单个球形生物质颗粒的热解过程进行数值模拟．利用双倒易边界元法和四阶龙格-库塔方法分别对非线性的导热方程和化学反应动力学方程组进行求解．讨论了生物质颗粒的大小与环境温度对热解时间和热解产物中相关组分质量分数的影响．     

DSMC方法下给定热流边界条件的实施

本文详细讨论了DSMC方法中流体温度、能量及边界热流的统计方法,发展了一种从边界热流求得与壁面碰撞分子反射速度的方法。该方法被称为逆温度抽样算法(ITS,Inverse Temperature Sampling)方法。在此基础上,本文发展了DSMC方法中壁面处给定热流边界条件的实施方法。计算结果表明: ITS方法能准确抽样反射分子的特征温度,进而求得分子反射速度。基于该方法的给定热流边界条件可以准确求得壁面处温度分布,以及流场内的压力、速度。     

一类新的求解双边障碍问题的数值方法

针对双边障碍问题的离散互补形式,提出了一类新的格式将其等价转化为方程组的形式,并采用牛顿迭代法进行求解.实验结果显示所提算法能快速,有效地计算出数值解和接触集.     

孔洞裂纹问题的一种数值方法

研究孔洞与裂纹的相互作用问题,通过把适于单一裂纹的Bueckner原理扩充到含有多孔洞多裂纹的一般体系,将原问题分解为承受远处载荷不含裂纹不含孔洞的均匀问题,和在远处不承受载荷但在裂纹面上和孔洞表面上承受面力的多孔洞多裂纹问题.于是,以应力强度因子作为参量的问题可以通过考虑后者来解决,而利用笔者提出的杂交位移不连续法,这种多孔洞多裂纹问题是容易数值求解的.算例说明该数值方法对分析平面弹性介质中孔洞与裂纹的相互作用既简单又有效.     

平面波输入下海水-海床-结构动力相互作用分析

海洋工程结构的地震反应分析是保证海洋工程结构地震安全的重要环节.由于其所处的复杂环境, 该问题涉及到流固耦合和土-结相互作用.本文基于海水、饱和海床、基岩流固耦合统一计算框架,采用Davidenkov模型和修正的Masing准则考虑饱和海床的非线性,在脉冲SV波垂直入射下, 进行了海域场地和海洋工程结构的动力响应分析. 首先,对比分析了线性自由场和非线性自由场输入情形的海域场地非线性反应,结果表明线性自由场输入时反应不合理,自由场分析和场地分析应该采用相一致的本构模型. 然后,对比分析了海床分别为线性和非线性情形时,海域场地以及海水-海床-结构体系的反应特征. 与线性海床情形相比,非线性对海床反应的影响主要由如下两方面因素控制: 一方面,非线性导致饱和海床模量减小, 饱和海床与基岩间的波阻抗比减小,由基岩到饱和海床间的反射系数和透射系数增加, 导致反应增大; 另一方面,非线性导致阻尼加大, 使海床反应减小. 对于本文算例而言,阻尼对非线性海床结果的影响占主导作用.