基于改进位移模式的一维C1有限元超收敛算法

提出了基于改进位移模式的一维C1有限元超收敛算法。利用单元内部需满足平衡方程的条件,推导了超收敛计算解析公式的显式,即将高阶有限元解的位移模式用常规有限元解的位移模式表示。用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式。采用积分形式推导了单元刚度矩阵。该算法在前处理阶段使用了超收敛计算公式,在常规试函数的基础上,增加了高阶试函数,使得单元内平衡方程的残差减少,从而达到提高精度的目标。对于Hermite单元,本文的结点和单元的位移、导数都达到了h4阶的超收敛精度。     

基于整体局部1,3位移模式的层合板层间应力研究

为提高层合板层间应力计算的准确性,对Reddy型高阶剪切理论的基本位移模式进行改进,提出整体-局部1,3高阶位移模式．在满足层间位移连续,层间剪切应力连续,以及上下表面自由的条件下,与前人提出的整体-局部1,2-3位移模式相比,层合板板结构每个节点的独立变量由13缩减到11,并且不随层数的增加而变化．将整体-局部1,3高阶位移模式位移和应力的数值解与解析解进行对比,验证了整体-局部1,3位移模式的准确性,可应用于复合材料层合板的位移和应力分析．     

复合材料面层夹层板中转动一致有效理论

采用分层位移模式，引入夹层板层间和表面剪应力协调条件，导出位移场的修正形式。根据能量误差一致原则，应用最小位能原理建立了复合材料层夹层板的基本方程和边界条件。     

单位载荷法中虚位移限制条件辨析

单位载荷法是一种求解位移的通用方法,在材料力学能量方法中占据重要地位.使用虚位移模式表述单位载荷方法时,满足位移边界条件是虚位移的主要限制条件.在静不定问题、特别是带有自由边界条件单位载荷系统构造时,对于该限制条件的理解决定了使用单位载荷方法的灵活性.本文使用几个经典例题剖析该限制条件的作用,为教学活动中解释相关问题提...     

精化直接刚度法的不协调位移模式与收敛性分析

基于精化直接刚度法的不协调位移模式及有限元收敛理论中的ＩＰＴ条件，建立了Ｃ￣０类可分离位移型精化元的收敛性充分条件——强ＩＰＴ条件。进一步推出一族平面及空间Ｃ￣０类精化元并验证了其收敛性。     

新型整体-局部高阶剪切变形理论及层合板精化三角形单元

基于整体-局部位移方法,建立了一种高阶剪切变形理论。整体位移部分采用的是Reddy理论的位移模式（1984）,局部位移为LIXY等（1997）建立的1,2-3理论的局部函数。这一理论使满足自由表面条件的Red@理论进一步满足层间位移、应力连续,同时有效减少了1,2—3理论的未知数个数。基于此理论深入开展了有限元法研究,建立了满足C^1连续条件的精化三节点三角形单元（每个节点参数为9个）。计算结果表明：建立的精化单元能准确计算整体位移和层间应力。     

基于改进位移模式的一维有限元超收敛算法

将多尺度方法的思想与超收敛计算的解析公式结合起来,提出了改进有限元位移模式的算法。利用超收敛计算的解析公式,将高阶有限元解的位移模式用常规有限元解的位移模式表示。用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式,采用积分形式推导了单元刚度矩阵。该算法在前处理和后处理两个阶段都使用超收敛计算公式,在常规试函数的基础上,增加了高阶试函数,使得单元内平衡方程的残差减少,从而达到提高精度的目标。对于线性单元,本文结点和单元的位移、导数都达到了h4阶的超收敛精度。     

Reissner型平板弯曲断裂问题分析

本文在Reissner型平板裂纹尖端位移场展开式基础上,采用高阶奇异元计算中厚板弯曲应力强度因子。本文在基本公式中考虑剪切变形影响,计算分析了有限尺寸板在不同厚度、不同宽度以及不同支承条件下应力强度因子及其变化,并对奇异元位移模式项数的选择、奇异元最佳尺寸的选取等问题进行了分析讨论。     

橡胶圆柱体的平面外剪切变形分析

利用Gao Y C给出的一类不可压缩的应变能函数,分析了不均匀橡胶圆柱体受平面外剪切变形问题.通过平面外剪切变形问题的变形模式,并结合应力和位移的边界条件以及界面连续条件得到了位移场和应力场在瞬时构型中的解析解,并讨论了材料的不均匀性对位移场及应力场分布的影响.结果表明:不均匀模型可以方便地退化成均匀模型;橡胶圆柱体的不均匀性对位移场和应力场的分布有着重要影响.     

位移模式下PMI泡沫夹层结构疲劳损伤试验研究

论文选取磁悬浮列车车身所用的由铝面板与聚甲基苯丙酰亚胺(PMI)聚合物泡沫芯层所组成的轻质夹层复合材料为研究对象,对夹层结构在室温下进行静态强度测试和进行以位移为控制变量的疲劳损伤演化试验,探讨分析了PMI泡沫夹层结构在交变位移控制下的疲劳性能和破坏行为.给出了PMI泡沫夹层结构在静态载荷作用下的力学性能参数.参考静态试验结果加载适当的位移载荷进行疲劳试验,发现在位移控制模式下,夹层结构的疲劳损伤过程和破坏模式明显区别于载荷控制模式.当最大控制位移较小时破坏形式为面板与泡沫层脱离;位移较大时为面板断裂和泡沫芯层塌陷.通过引入载荷的变化作为损伤参量建立了位移控制模式下损伤演变公式,并对两种模式的破坏行为进行了较好的预测.     

A NEW HYBRID QUADRILATERAL FINITE ELEMENT FOR MINDLIN PLATE

ＡＮＥＷＨＹＢＲＩＤＱＵＡＤＲＩＬＡＴＥＲＡＬＦＩＮＩＴＥＥＬＥＭＥＮＴＦＯＲＭＩＮＤＬＩＮＰＬＡＴＥＣｈｉｎＹｉ（秦奕）（ＴｉａｎｊｉｎＡｒｃｈｉｔｅｃｔｕｒａｌＤｅｓｉｇｎＩｎｓｔｉｔｕｔｅ，Ｔｉａｎｊｉｎ）ＺｈａｎｇＪｉｎｇ－ｙｕ（张敬宇）（Ｉ...     

Alternative state space formulations for magnetoelectric thermoelasticity with transverse isotropy and the application to bending analysis of nonhomogeneous plates

By introducing two displacement functions and two stress functions, the governing equations of the linear theory of magneto-electro-thermo-elasticity with transverse isotropy are simplified. On selecting certain physical quantities as the basic unknowns, two new state equations are established. Each of them is order reduced when compared with the one reported recently in literature, leading to a higher numerical efficiency. The material inhomogeneity along the axis of symmetry (z-direction) can be taken into account and an approximate laminate model is employed to facilitate deriving analytical solutions. The validity of new formulations is examined by considering a laminated magneto-electro-elastic rectangular plate and good agreement is obtained with existent results. A plate with a functionally graded property is then analyzed. The effect of magnetoelectric coupling in a BaTiO₃–CoFe₂O₄ composite predicted from the micromechanics simulation is studied quantitatively.     

An experimental and numerical investigation of air side heat transfer and flow characteristics on finned plate configuration

In this paper, a new type of finned plate heat exchanger (FPHE) is presented to recover the waste heat from exhaust flue gases. A finned plate configuration causes low pressure drop and it is especially appropriate for heat transfer at the flue gas side. Meanwhile, this paper presents a detailed experimental and numerical study of convection heat transfer and pressure drop of the new structure. Three-dimensional numerical simulation results using the CFD code FLUENT6.3 were compared with experimental data to select the best model. The heat transfer and pressure drop with different geometry pattern was then studied numerically using the selected model. And the velocity field and temperature distribution of air flow in the finned plate channel are presented with different geometry patterns. These results provide insight into improved designs of FPHEs.     

Fundamental-solution-based finite element model for plane orthotropic elastic bodies

A new hybrid finite element formulation is presented for solving two-dimensional orthotropic elasticity problems. A linear combination of fundamental solutions is used to approximate the intra-element displacement fields and conventional shape functions are employed to construct elementary boundary fields, which are independent of the intra-element fields. To establish a linkage between the two independent fields and produce the final displacement-force equations, a hybrid variational functional containing integrals along the elemental boundary only is developed. Results are presented for four numerical examples including a cantilever plate, a square plate under uniform tension, a plate with a circular hole, and a plate with a central crack, respectively, and are assessed by comparing them with solutions from ABAQUS and other available results.     

A multi-entropy-level lattice Boltzmann model for the one-dimensional compressible Euler equations

In this paper, we propose a new lattice Boltzmann model for the one-dimensional compressible Euler equations. The new model is based on a three-entropy-level and three-speed lattice Boltzmann equation by using a method of higher-order moments of the equilibrium distribution functions. In order to obtain the second-order accuracy model, we employ the ghost field distribution functions to remove the non-physical dissipation terms in the Euler equations. We also use the conditions of the higher-order moments of the ghost field equilibrium distribution functions to obtain the equilibrium distribution functions. The numerical examples show that the numerical results can be compared with those classical methods.     

弹性薄板弯曲问题的边界轮廓法

导出了弹性薄板弯曲问题边界积分方程的另一种形式,基于这种方程,提出了平板弯曲问题的边界轮廓法,讨论了三次边界单元边界轮廓法的计算列式,并给出了计算内力的边界轮廓法方程。该法无需进行数值积分计算,完全避免了角点问题和奇异积分计算。给出的算例,与解析解相比较,证实该方法的有效性。     

新修正偶应力层合板理论中的几个基本假设对自由振动计算的影响

基于新修正偶应力理论,在对微细观尺度的复合材料层合梁/板进行力学响应计算时,往往采用一系列假设来简化模型。现有文献都全部或部分应用了这些假设,但对这些假设是否会对计算结果造成影响尚未进行充分讨论分析。本文建立了未经简化的新修正偶应力Reddy层合板模型,并对其自由振动进行了分析。通过数值算例的对比,讨论了常用的几个简化假设对微细观复合材料四边简支方板自振频率的影响以及适用范围。算例结果表明,常用的几个简化假设对于微尺度层合薄板自由振动的影响很小,对于厚板的低阶频率影响也很小,但对厚板的高阶频率影响显著。     

A non-classical Kirchhoff plate model incorporating microstructure,surface energy and foundation effects

A new non-classical Kirchhoff plate model is developed using a modified couple stress theory, a surface elasticity theory and a two-parameter elastic foundation model. A variational formulation based on Hamilton’s principle is employed, which leads to the simultaneous determination of the equations of motion and the complete boundary conditions and provides a unified treatment of the microstructure, surface energy and foundation effects. The new plate model contains a material length scale parameter to account for the microstructure effect, three surface elastic constants to describe the surface energy effect, and two foundation moduli to represent the foundation effect. The current non-classical plate model reduces to its classical elasticity-based counterpart when the microstructure, surface energy and foundation effects are all suppressed. In addition, the newly developed plate model includes the models considering the microstructure dependence or the surface energy effect or the foundation influence alone as special cases and recovers the Bernoulli–Euler beam model incorporating the microstructure, surface energy and foundation effects. To illustrate the new model, the static bending and free vibration problems of a simply supported rectangular plate are analytically solved by directly applying the general formulas derived. For the static bending problem, the numerical results reveal that the deflection of the simply supported plate with or without the elastic foundation predicted by the current model is smaller than that predicted by the classical model. Also, it is observed that the difference in the deflection predicted by the new and classical plate models is very large when the plate thickness is sufficiently small, but it is diminishing with the increase of the plate thickness. For the free vibration problem, it is found that the natural frequency predicted by the new plate model with or without the elastic foundation is higher than that predicted by the classical plate model, and the difference is significant for very thin plates. These predicted trends of the size effect at the micron scale agree with those observed experimentally. In addition, it is shown both analytically and numerically that the presence of the elastic foundation reduces the plate deflection and increases the plate natural frequency, as expected.     

Bending analysis of magnetoelectroelastic nanoplates resting on Pasternak elastic foundation based on nonlocal theory

Based on the nonlocal theory and Mindlin plate theory, the governing equations (i.e., a system of partial differential equations (PDEs) for bending problem) of magnetoelectroelastic (MEE) nanoplates resting on the Pasternak elastic foundation are first derived by the variational principle. The polynomial particular solutions corresponding to the established model are then obtained and further employed as basis functions with the method of particular solutions (MPS) to solve the governing equations numerically. It is confirmed that for the present bending model, the new solution strategy possesses more general applicability and superior flexibility in the selection of collocation points. The effects of different boundary conditions, applied loads, and geometrical shapes on the bending properties of MEE nanoplates are evaluated by using the developed method. Some important conclusions are drawn, which should be helpful for the design and applications of electromagnetic nanoplate structures.     

Computations of modes I and II stress intensity factors of sharp notched plates under in-plane shear and bending loading by the fractal-like finite element method

The fractal-like finite element method (FFEM) is used to compute the stress intensity factors (SIFs) for different configurations of cracked/notched plates subject to in-plane shear and bending loading conditions. In the FFEM, the large number of unknown variables in the singular region around a notch tip is reduced to a small set of generalised co-ordinates by performing a fractal transformation using global interpolation functions. The use of exact analytical solutions of the displacement field around a notch tip as the global interpolation functions reduces the computational cost significantly and neither post-processing technique to extract SIFs nor special singular elements to model the singular region are required. The results of numerical examples of various configurations of cracked/notched plates are presented and validated via published data. Also, new results for cracked/notched plate problems are presented. These results demonstrate the accuracy and efficiency of the FFEM to compute the SIFs for notch problems under in-plane shear and bending loading conditions.