New Sixth-Order Compact Schemes for Poisson/Helmholtz Equations

Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two- and three-dimensions are developed and analyzed. Different from a few sixth-order compact finite difference schemes in the literature, the finite difference and weight coefficients of the new methods have analytic simple expressions. One of the new ideas is to use a weighted combination of the source term at staggered grid points which is important for grid points near the boundary and avoids partial derivatives of the source term. Furthermore, the new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6. The coefficient matrices of the new schemes are $M$-matrices for Helmholtz equations with wave number $K≤0,$ which guarantee the discrete maximum principle and lead to the convergence of the new sixth-order compact schemes. Numerical examples in both 2D and 3D are presented to verify the effectiveness of the proposed schemes.     

八节点等参元滑动接触算法的改进

分析认为采用八节点等参元计算滑动接触问题时,往往存在较大的接触应力计算误差的主要原因在于这种单元具有两种不同的节点。为此,本文提出了双模拟接触边界法,即分别采用两种不同的节点,形成两个模拟接触边界,并限制“从接触边界”上的节点只能和“主接触边界”上其同类节点形成的模拟接触边界形成接触关系,从而有效的减小了采用八节点等参元计算滑动接触问题时的误差。     

材料的动态损伤和失效

本文对动载荷下材料内部损伤演化过程的研究进行了较为全面、系统的介绍和总结,并对现有的各种动态损伤模型(包括笔者所提出的模型)进行了必要的比较、讨论和评述。     

弹塑性岩土介质中渗透固结问题的参数二次规划法

岩土类介质在承受荷载时,不仅产生弹塑性变形,还伴随着渗流和固结,是一个关于时间的动态过程。本文建立起处理该问题的参变量变分原理以及相应的有限元方法,这样将原问题化为求解带约束条件(本构状态方程)的二次规划问题.文中讨论了单元的选取形式及具体的实施过程,还给出了一个实例.     

影响无网格方法求解精度的因素分析

基于移动最小二乘法的无网格方法的计算精度除受到节点的分布密度和基底函数的阶次影响外，还受到其它因素的影响，其中权函数的选取、权函数影响域的大小及位移边界条件的引入对计算精度影响较大。本文分析了几种常用权函数在数值计算时的特点，包括计算精度、收敛情况、计算效率等，同时分析了影响域大小及边界条件的引入对计算精度的影响。通过分析给出了确定权函数及其影响域大小的方法。当受约束的自由度较多时，通过配点法引入位移边界条件会引起计算结果的振荡，通过施加稳定项可以消除振荡现象，通过对带孔方板的受力分析证明了其可行性。应用以上结论对J23—10曲柄压力机机身进行了受力分析，应力集中部位的计算结果得到了较高的精度。     

SOIL AGGREGATE AND ITS RESPONSE TO LAND MANAGEMENT PRACTICES

This paper provides a broad review of the existing study on soil aggregate and its responses to land management practices. Soil aggregate is used for structural unit, which is a group of primary soil particles that cohere to each other more strongly than other surrounding particles. The mechanism of soil particle aggregation may be expressed by a hierarchical model, which is based upon the hypothesis that macroaggregates （〉250μm） are collections of smaller microaggregates （〈250μm） held together with organic binding agents. Primary particles form microaggregates and then macroaggregates. Carbon （C）-rich young plant residues form and stabilize macroaggregates, whereas old organic C is occluded in the microaggregates. The interaction of aggregate dynamics with soil organic carbon （SOC） is complex and embraces a range of spatial and temporal processes within macroaggregates and microaggregates. The nature and properties of aggregates are determined by the quantity and quality of coarse residues and humic compounds and by the degree of their interaction with soil particles. The mechanisms resulting in the binding of primary soil particles into stable aggregates vary with soil parent material, climate, vegetation, and land management practices. Land management practices, including tillage methods, residue management, amendments, and soil fertility management, enhance soil aggregation. However, there is still much uncertainty in the dynamics of organic matter in macroaggregation and microaggregation, and research is still needed to understand further the mechanisms of aggregate formation and its responses to human activities.     

参数拟二次规划法及其在复合材料失效过程分析中的应用

本文介绍关于不连续线弹性静力学问题的多变量变分原理与相应的参数拟二次规划算法在复合材料层压板强度分析及纤维增强脆性基体材料单元细观失效过程分析中的应用。     

ADAPTIVE MESHLESS METHOD BASED ON LOCAL FIT TECHNOLOGY

An h-adaptive meshless method is proposed in this paper. The error estimation is based on local fit technology, usually confined to Voronoi Cells. The error is achieved by comparison of the computational results with smoothed ones, which are projected with Taylor series. Voronoi Cells are introduced not only for integration of potential energy but also for guidance of refinement.New nodes are placed within those cells with high estimated error. At the end of the paper, two numerical examples with severe stress gradient are analyzed. Through adaptive analysis accurate results are obtained at critical subdomains, which validates the efficiency of the method.     

考虑晶界效应的多晶体有限变形分析

将晶界及其影响区综合考虑，建立了考虑晶界效应的力学模型，结合晶体塑性理论，利用有限变形有限元对多晶体进行数值模拟，数值结果显示了细观层次下晶粒变形场的特点，理论计算同实验定性一致。     

A screw dislocation near an elliptic hole in non-local elasticity

Summary The solutions are given for the stress field and elastic energy of a screw dislocation near an elliptic hole in non-local elasticity by a conformal mapping technique. The numerical results are discussed for the stress along thex-axis. Allowing one axis of the ellipse to dwindle to zero, the solution of the finite length crack is obtained.

---

Schraubenversetzung neben einem elliptischen Loch bei nichtlokaler Elastizität

Übersicht Durch eine konforme Abbildungstechnik werden die Lösungen für das Spannungsfeld und die elastische Energie der Schraubenversetzung neben einem elliptischen Loch in nichtlokaler Elastizität gegeben. Die numerischen Ergebnisse werden für die Spannung entlang derx-Achse diskutiert. Läßt man eine Ellipsenachse gegen Null gehen, so gewinnt man die Lösung für den Riß endlicher Länge.