金属刚粘塑性变形的能量泛函与动力显式算法有限元分析

设计出刚粘塑性动态变形场的能量泛函，采用罚函数法和Ｌａｇｒａｎｇｉａｎ乘子法取消对运动容许速度场满足体积不可压缩的约束条件，运用虚功原理和广义变分原理，推导出刚粘塑性变形的动力分析显式算法有限元方程，以及速率形式的中心差分求解的时间积分显式，给出了两个金属变形的算例。     

功能梯度层与弹性层间的单退让接触问题

基于所有接触面间光滑的假设,研究了同时受压的功能梯度层与弹性层间的单退让平面接触问题.假设功能梯度层是各向同性的非均匀材料,其剪切模量按照指数函数形式变化.利用Fourier积分变换把问题转化为求解奇异积分方程.然后利用Gauss-Chebyshev求积公式和迭代法得到下层接触应力和退让接触半径的数值解.最后在数值算例中,分别讨论了两层间的厚度比值,功能梯度层的硬度参数,以及上层接触半径对退让接触半径与下层接触应力的影响.     

耦合GPU与PCG的EFG法并行计算及应用研究

针对迭代法求解无网格Galerkin法中线性方程组收敛速度慢的问题,提出了一种耦合GPU和预处理共轭梯度法的无网格Galerkin法并行算法,在对其总体刚度矩阵、总体惩罚刚度矩阵进行并行联合组装的同时即可得到对角预处理共轭矩阵,有效地节省了GPU的存储空间和计算时间;通过采用四面体积分背景网格,提高了所提算法对三维复杂几何形状问题的适应性。通过2个三维算例验证了所提算法的可行性,且预处理共轭梯度法与共轭梯度法相比,其迭代次数最大可减少1686倍,最大的迭代时间可节省1003倍;同时探讨了加速比与线程数和节点个数之间的关系,当线程数为64时其加速比可达到最大,且预处理共轭梯度法的加速比与共轭梯度法相比可增大4.5倍,预处理共轭梯度法的加速比最大达到了88.5倍。     

接触-碰撞算法研究进展

接触-碰撞广泛存在于实际工程问题中,是影响数值计算效率与计算精度的重要因素。本文针对变形体间接触-碰撞问题的显式有限元计算,介绍接触-碰撞算法近30年来取得的主要进展。首先,简要介绍接触-碰撞问题的界面离散模型;然后,从全局接触搜索、局部搜索、接触约束施加以及接触计算的并行化四方面详细阐述目前主要算法的基本思想与特点,并分析其优势与不足;最后,对接触-碰撞算法相关研究方向给出建议。     

????????????????????????????

合隐式和显式时间积分技术,对结构非线性动力反应分析提出一种并行混合时间积分算 法. 该算法采用区域分解技术. 将并发性引入到算法中,即利用显式时间积分技术进行界面 节点积分而利用隐式算法求解局部子区域. 为实现并行混合时间积分算法,设计了灵活的并 行数据信息流. 编写了该算法的程序,在工作站机群实现了数值算例,验证了算法的精度和 性能. 计算结果表明该算法具有良好的并行性能,优于隐式算法.     

广义Hamilton系统的保结构算法

基于广义Hamilton系统微分方程解析解理论。给出了构造保持系统真解典则性的高阶显式积分格式的方法,并说明其可推广到广义Hamilton控制系统。该方法保持了原系统的几何定性特征,因而是稳定的。数值例子说明了算法的有效性。     

岩质圆形隧洞围岩应力场弹塑性新解

针对动态接触问题的有限元并行计算，提出了一种新的接触算法. 新算法引入局部拉氏 乘子技术来计算接触力. 由于同时考虑了无穿透的接触约束条件和相邻接触对的相互影响， 较之广泛使用的罚参数法，新算法使接触约束条件和系统平衡方程得到更充分的满足. 虽然 为提高接触计算精度而在局部采用了迭代技术，但算法仍然具有较高的效率，且与显式时间 积分方案完全相容. 此外，通过构造专门的区域分解方案，实现了将现有为串行程序开发的 搜索算法平滑移植到并行环境的目标. 数值算例表明，所提出的接触算法具有很好的并行性， 在保证了接触问题并行计算精度的同时，取得了满意的并行效率.     

三维板料成形过程的显式有限元分析

本文采用基于随动坐标系的假定应变域壳单元及显式有限元格式求解三维板料成形问题。板材料服从Ｈｉｌ各向异性弹塑性准则，板料与模具之间的接触界面由主仆面接触搜寻法处理，接触力由罚参数法计算。文中给出了几个三维成形过程的计算实例。数值算例表明，本文方法具有较高的计算精度和计算效率，可在微机上分析中等复杂程度的成形过程     

进化算法与确定性算法在优化控制问题中的收敛性对比

对比了进化算法(基因算法)与确定性算法(共轭梯度法)在优化控制问题中的优化效率．两种方法都与分散武优化策略-Nash对策进行了结合，并成功地应用于优化控制问题。计算模型采用绕NACA0012翼型的位流流场．区域分裂技术的引用使得全局流场被分裂为多个带有重叠区的子流场，使用4种不同的方法进行当地流场解的耦合，这些算法可以通过当地的流场解求得全局流场解。数值计算结果的对比表明．进化算法可以得到与共轭梯度法相同的计算结果．并且进化算法的不依赖梯度信息的特性使其在复杂问题及非线性问题中具有广泛的应用前景。     

黏性阻尼系统时域响应灵敏度及其一致性研究

时域响应灵敏度分析是时域梯度优化算法的基础. 灵敏度分析通常只涉及对设计变量的微分运算, 但时域响应灵敏度问题还涉及时间域的离散化. 因此, 微分和离散的先后顺序可能对时域响应灵敏度结果产生影响. 针对黏性阻尼系统时域响应灵敏度求解问题, 基于改进精细积分方法, 分别推导了先微分后离散和先离散后微分两种伴随变量方法. 其中, 先微分后离散法首先对由伴随变量构造的增广函数微分, 再利用改进精细积分方法在各离散时间点求解时域响应灵敏度; 而先离散后微分方法则首先在各离散时间点引入残值方程构造增广函数, 再对各增广函数进行微分以求解时域响应灵敏度. 通过数值算例验证了所提出方法的有效性和准确性, 并与传统基于Newmark的方法进行比较. 结果表明, 积分方案、数值离散误差以及离散和微分的先后顺序共同影响灵敏度的一致性误差. 综合考虑精度、效率和一致性问题, 基于改进精细积分的先微分后离散伴随变量法表现更优, 最适合应用于黏性阻尼系统时域梯度优化算法.      

粗糙表面接触研究进展

工程中的任何接触表面都不会是绝对光滑的,对于某些特殊接触问题必须考虑表面形貌的接触效应.粗糙表面接触研究是一门年轻的交叉学科.它主要研究粗糙表面接触力学,接触导电和导热,接触密封,接触振动和噪声,接触刚度,骨关节接触等问题,涉及固体力学,流体力学,随机过程及概率统计,计量学,电学和传热学,振动理论以及生物学等学科知识.本文综述了粗糙表面接触的研究进展.      

The wetting problem of fluids on solid surfaces. Part 2: the contact angle hysteresis

In part 1 (Gouin, [13]), we proposed a model of dynamics of wetting for slow movements near a contact line formed at the interface of two immiscible fluids and a solid when viscous dissipation remains bounded. The contact line is not a material line and a Young-Dupré equation for the apparent dynamic contact angle taking into account the line celerity was proposed. In this paper we consider a form of the interfacial energy of a solid surface in which many small oscillations are superposed on a slowly varying function. For a capillary tube, a scaling analysis of the microscopic law associated with the Young-Dupré dynamic equation yields a macroscopic equation for the motion of the contact line. The value of the deduced apparent dynamic contact angle yields for the average response of the line motion a phenomenon akin to the stick-slip motion of the contact line on the solid wall. The contact angle hysteresis phenomenon and the modelling of experimentally well-known results expressing the dependence of the apparent dynamic contact angle on the celerity of the line are obtained. Furthermore, a qualitative explanation of the maximum speed of wetting (and dewetting) can be given.Received: 5 June 2001, Accepted: 24 May 2003, Published online: 29 July 2003PACS: 02.90, 47.50, 66.20, 68.03, 68.08     

表面织构对界面接触参数影响的实验研究

机械接触表面的接触刚度和接触阻尼是直接影响界面振动水平的两个重要动力学参数。本文介绍了一种测试机械连接界面接触刚度和接触阻尼的实验装置和方法,并针对用两种不同加工方法(激光加工及线切割)加工出的具有不同尺寸织构的试样与光滑平面相接触的界面进行了测量,初步探索了平面接触刚度和接触阻尼的表面织构效应,并得出结论：织构及润滑对接触刚度会产生一定影响,而实验中均未观察到这两者对接触阻尼有明显影响。     

Contact radius of stamps in reversible adhesion

A mechanics model is developed for the contact radius of stamps with pyramid tips in transfer printing. This is important to the realization of reversible control of adhesion, which has many important applications, such as climbing robots, medical tapes, and transfer printing of electronics. The contact radius is shown to scale linearly with the work of adhesion between the stamp and the contacting surface, and inversely with the plane-strain modulus of the stamp. It also depends on the cone angle and tip radius of the stamp, but is essentially independent of details of the tip geometry.     

Separated steady state solutions for two thermoelastic half-planes in contact with out-of-plane sliding

When two materials are placed in contact along an interface, thermoelastic effects can separate the surfaces and create “hot spots” when there is sufficient frictional heating fVp generated at the interface, even if the two surfaces are nominally flat. Additionally, heat can flow because the bodies are generally at different temperatures, and this is an independent cause of separation, generally when heat flows into the less distortive material. These two effects have been considered separately, and here we consider the case with interaction of the two effects, showing possible non-existence, multiplicity and instability of solutions. Approximate Hertzian solutions for the separated contact regime are very limited, particularly for the frictional heating case. Hence, a new efficient full numerical solution is developed, and compared with direct FEM results, the latter permitting also the assessment of stability in the transient regime. Connection to previous results for simple rod models is made. The case of heat flow into the more distortive material is discussed.     

Contact angle calculations from the contact/maximum diameter of sessile drops

This paper presents two algorithms for computing the contact angle of sessile liquid drops given data about the drops. The first yields the contact angle given the volume, surface tension and maximum diameter (or contact diameter) of a single drop. This algorithm is an extension of existing algorithms based on knowledge of the maximum diameter or of the contact diameter of a drop. A sensitivity analysis is included for this algorithm, allowing estimates to be made of the error in computed contact angle caused by errors in the measurement of the volume and/or diameter. The second algorithm requires only the volume and maximum or contact diameter of two different drops as input, and it produces both the contact angle and surface tension as output. Both algorithms are based on Newton's method applied to a function whose value is computed by solving a system of ordinary differential equations obtained from the Laplace equation of capillarity. The techniques are applicable to both hydrophobic and hydrophilic surfaces. Copyright © 2000 John Wiley & Sons, Ltd.     

A model for computational investigation of elasto-plastic normal and tangential contact considering adhesion effect

With the rapid development of Micro-Electro-Mechanical System (MEMS), we enter a field in which the surface effects have dominated many of the micro-scale phenomena, and the adhesive contact is one of the focuses. In this paper, a feasible model for finite element computation is presented via a macroscopic and microscopic combination approach, in which the adhesive forces are simulated by some non-linear spring elements considering the softening stage. Two basic problems concerning the adhesion effect were considered; through specific quantitative analysis, the results show a consistency with the current elastic continuum theories of adhesion and a brief investigation into the effects of adhesion on plastic deformation and tangential contact will be carried out as well. The project supported by the National Natural Science Foundation of China (10172050, 90205022) and Key Grant Project of Chinese MoE (0306)     

线接触弹性接触变形的解析算法

以一般光滑性体接触理论为基础,结合有限长弹性体接触的特点,求出线接触弹性接蟹变形的解析公式,并发现其解析解与数值解具有很好的一致性,所得公式可以对赫兹线接触理论加以补充,与经验公式相比,它能够确切反映材料、载荷以及曲率半径等对接触变形的影响,为工程中的精确计算提供了方便。     

A contact searching algorithm for contact-impact problems

A new contact searching algorithm for contact-impact systems is proposed in this paper. In terms of the cell structure and the linked-list, this algorithm solves the problem of sorting and searching contacts in three dimensions by transforming it to a retrieving process from two one-dimensional arrays, and binary searching is no longer required. Using this algorithm, the cost of contact searching is reduced to the order ofO(N) instead ofO(Nlog₂ N) for traditional ones, whereN is the node number in the system. Moreover, this algorithm can handle contact systems with arbitrary mesh layouts. Due to the simplicity of this algorithm it can be easily implemented in a dynamic explicit finite element program. Our numerical experimental result shows that this algorithm is reliable and efficient for contact searching of three dimensional systems. The project supported by the National Natural Science Foundation of China (59875045), and the State Key Laboratory of Automobile Safety and Energy Saving (K9705)     

Indentation of a Punch with a Fine-Grained Base into an Elastic Foundation

The linear contact problem for a system of small punches located periodically on a part of the boundary of an elastic foundation is studied. An averaged contact problem is derived using the Marchenko–Khruslov averaging theory. An asymptotic formula is obtained for the translational capacity of a smooth punch with a fine-grained flat base.