A mixed differential quadrature and finite element free vibration and buckling analysis of thick beams on two-parameter elastic foundations

As a first endeavor, a mixed differential quadrature (DQ) and finite element (FE) method for boundary value structural problems in the context of free vibration and buckling analysis of thick beams supported on two-parameter elastic foundations is presented. The formulations are based on the two-dimensional theory of elasticity. The problem domain along axial direction is discretized using finite elements. The resulting system of equations and the related boundary conditions are discretized in the thickness direction and in strong-form using DQM. The method benefits from low computational efforts of the DQ in conjunction with the effectiveness of the FE method in general geometry and systematic boundary treatment resulting in highly accurate and fast convergence behavior solution. The boundary conditions at the top and bottom surface of the beams are implemented accurately. The presented formulations provide an effective analysis tool for beams free of shear locking. Comparisons are made with results from elasticity solutions as well as higher-order beam theory.     

Topology optimization of microfluidic mixers

This paper demonstrates the application of the topology optimization method as a general and systematic approach for microfluidic mixer design. The mixing process is modeled as convection dominated transport in low Reynolds number incompressible flow. The mixer performance is maximized by altering the layout of flow/non‐flow regions subject to a constraint on the pressure drop between inlet and outlet. For a square cross‐sectioned pipe the mixing is increased by 70% compared with a straight pipe at the cost of a 2.5 fold increase in pressure drop. Another example where only the bottom profile of the channel is a design domain results in intricate herring bone patterns that confirm findings from the literature. Copyright © 2008 John Wiley & Sons, Ltd.     

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

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

某边坡变形破坏机理的数值分析

某花岗岩矿岩土体主要有填土、粉质粘土、淤泥及淤泥质土、含泥中粗砂、粘土、残积砂质粘土,其下为花岗岩基岩。边坡最大高度75m,其中残积土层以上的土质边坡高度约50m。本文运用有限元强度折减法对该边坡的稳定性进行了研究,得出了潜在的破坏模式及各个滑动部分的安全系数。分析结果表明,本边坡存在3个潜在滑动面,其中两个为浅层滑动面,另一个为深层滑动面。两个浅层滑动面的稳定系数在1.13～1.17,处于临界状态。深层滑动面稳定性系数在1.29,目前尚处在相对稳定状态。根据分析结果,提出了如下加固措施:（1）对滑坡体进行削坡处理;（2）在其滑动面下方设置抗滑台阶;（3）采用预应力锚索框架梁结合肋柱进行加固处理。     

非线性动力学常微分方程组高精度数值积分方法

建立了一种求解非线性动力学常微分方程组初值问题的新方法．若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解．引入一般自治或非自治非线性系统的首次近似Jacobi矩阵，不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel卷积分解析表达式，在时间步长内进行数值积分选代求解，在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解．积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法．数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率．算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题．基于以上理论和数值方法，计算了线性非线性算例并进行了分析，验证了本方法简捷而有广泛的适应性，可以有足够的精确性．     

压电传感器探测层合板分层的数值分析

用修正的准三维有限元法分析了受单轴拉伸表面覆盖有压电薄膜的复合材料层合板。文中首先讨论了压电薄膜层对层合板层间应力的影响，然后在存在分层裂纹的情况下，计算了压电层的电压分布。结果表明将压电薄膜作为传感器来进行分层探测是可行的。文中最后考虑了多种因素对传感灵敏度的影响。     

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.     

三维非线性有限元与弹性边界元耦合数值方法

本文系统地讨论了以下三个问题:(1) 有限元与边界元耦合中的几个数值问题,其中包括:边界积分方程的凝聚、等效刚度矩阵的对称化及面力不连续的处理;(2) 弹塑性有限元与弹性边界元的耦合;(3) 弹粘塑性有限元与弹性边界元的耦合及数值计算稳定性条件。     

微重力环境下Marangoni对流的有限元数值模拟

本文以悬浮区为背景研究液桥中气液交界面上由表面张力所驱动的流体对流。我们采用有限元方法对轴对称、定常运动方程,能量方程,自由面切向、法向应力平衡条件迭代求解,首次考虑了边界形状的影响,得到了自洽的流区位形和流场、温度场、表面压强分布。一般说来,流区自由面呈弯月形。结果表明,只要表面张力数 C_a<1,静态平衡的流区位形就是动态情况的良好近似。本文还分析了C_a数及G_r 数对流区位形的影响,得到了不同 M_a 数及散热条件下的温度场和流函数分布。