NUMERICAL SIMULATIONS OF WHOLE FIELD BEHAVIOR FOR ROUND BUOYANT JET

ＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮＳＯＦＷＨＯＬＥＦＩＥＬＤＢＥＨＡＶＩＯＲＦＯＲＲＯＵＮＤＢＵＯＹＡＮＴＪＥＴＬｉＷｅｉ，ＨｕａｉＷｅｎ－ｘｉｎ（ＷｕｈａｎＵｎｉｖｅｒｓｉｔｙｏｆＨｙｄｒａｕｌｉｃａｎｄＥｌｅｃｔｒｉｃＥｎｇｉｎｅｅｒｉｎｇ，...     

An eight‐node hybrid‐stress solid‐shell element for geometric non‐linear analysis of elastic shells

This paper presents eight‐node solid‐shell elements for geometric non‐linear analysis of elastic shells. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc modified generalized laminate stiffness matrix are employed. A selectively reduced integrated element is formulated with its membrane and bending shear strain components taken to be constant and equal to the ones evaluated at the element centroid. With the generalized stresses arising from the modified generalized laminate stiffness matrix assumed to be independent from the ones obtained from the displacement, an extended Hellinger–Reissner functional can be derived. By choosing the assumed generalized stresses similar to the assumed stresses of a previous solid element, a hybrid‐stress solid‐shell element is formulated. Commonly employed geometric non‐linear homogeneous and laminated shell problems are attempted and our results are close to those of other state‐of‐the‐art elements. Moreover, the hybrid‐stress element converges more readily than the selectively reduced integrated element in all benchmark problems. Copyright © 2002 John Wiley & Sons, Ltd.     

A 3D contact smoothing method using Gregory patches

In this work, a method is developed for smoothing three‐dimensional contact surfaces. The method can be applied to both regular and irregular meshes. The algorithm employs Gregory patches to interpolate finite element nodes and provide tangent plane continuity between adjacent patches. The resulting surface interpolation is used to calculate gaps and contact forces, in a variationally consistent way, such that contact forces due to normal and frictional contact vary smoothly as slave nodes transition from one patch to the next. This eliminates the ‘chatter’ which typically occurs in a standard contact algorithm when a slave node is situated near a master facet edge. The elimination of this chatter provides a significant improvement in convergence behaviour, which is illustrated by a number of numerical examples. Furthermore, smoothed surfaces also provide a more accurate representation of the actual surface, such that resulting stresses and forces can be more accurately computed with coarse meshes in many problems. This fact is also demonstrated by the numerical examples. Published in 2002 by John Wiley & Sons, Ltd.     

多约束的桥梁结构拓扑优化

基于ESO (Evolutionary Structural Optimization)的拓扑优化方法,本文提出了一种适应于桥梁结构的拓扑优化方法。引进了两种性能指标公式来确定最佳拓扑设计。另外,为了更有效地尝试改进最终设计的细节,而又不进行更精细有限元网格的完整分析,这里提出了一种精细网络设计方案。再者,结合细啃技术,完成了考虑应力、位移、频率约束的斜拉桥优化设计。通过几种桥梁的设计优化,表明该方法的可应用性、简洁性和有效性。     

随机徐变应力影响下重力坝时变可靠度初探

由于重力坝的随机温度场以及随机徐变应力场均为非平稳随机过程，由此所涉及的结构可靠度是非平稳向量随机过程非线性组合的时变可靠度．本文计入重力坝徐变应力场和混凝土强度随机性的影响，将时变可靠度的数值解法与随机有限元法相结合，模拟施工过程，对一典型重力坝的时变可靠度进行了计算，得出了一些有益的结论．     

双材料挤压成形的数值分析

基于非线性有限元理论和方法，并引入摩擦单元技术，对双材料的挤压成形进行了数值分析，揭示了挤压过程中双材料的流动规律。     

各向加权异性有限元法在不可压粘性流动计算中的应用

本文利用各向异性张力有限单元，考虑到其独特的性质，提出了一种各向加权异性的算法，即在不同方向采用不同连续性的插值函数，然后再在原连续性较弱的方向以连续性较强的插值函数计算，最后把两次计算的结果加权处理。这种方法可以适应两个方向连续性要求不同或相同的各种情况，即克服了一般各向异性张力有限单元的局限性，又大大减少了存贮量和计算量。     

The effect of particle distribution on damage formation in particulate reinforced metal matrix composites deformed in compression

Image analysis results are reported on the generation of damage in particulate reinforced metal matrix composites during compressive deformation. The technique allows the automated collection of data on the incidence of particle fracture and void formation in the matrix as a function of important microstructural parameters such as local particle volume fraction and particle size. There is a strong relationship between damage and the local volume fraction of the reinforcement proving that damage formation is accentuated in regions of particle clustering. With the SiC reinforced materials examined, there was observed to be a change in dominance of damage mechanism from particle fracture at low local volume fractions to void formation in the matrix within strongly clustered regions. The results are compared with finite element (FE) modelling of the compressive deformation of clustered particles using a simple cluster of equi-spaced particles. The FE results suggest that plastic flow is generally inhibited in clustered regions. In certain highly clustered configurations shielding is such that flow does not occur in the heart of the cluster even at high levels of average plastic strain. The modelling suggests that the change in dominance of damage mechanism is related to the dramatic increase in tensile hydrostatic stresses in the matrix with higher levels of particle clustering.     

Pi-coefficient analysis of update algorithms for adaptive systems

This paper finds the appropriate pi-coefficients for a parameter estimation adaptive system and uses them to analyze the stability of two estimation algorithms. The estimation error dynamics of the system are modeled by a linear time-invariant subsystem and a nonlinear time-varying update law in a feedback loop. Then the so-called max-p problems are formulated and solved to obtain the pi-coefficients for the linear subsystem and nonlinear update low. For the investigated system, the quantitative results show that the least-squares update algorithm has larger stability range than that of the gradient algorithm, and the σ-modification scheme gives larger stability ranges for both algorithms.