Tailoring particle shape for enhancing the homogeneity of powder mixtures:Experimental study and DEM modelling

The effect of particle shape modification on the segregation reduction of enzyme granules in laundry detergent powder mixtures was investigated,both experimentally and computationally using Deseret Element Method(DEM).The shape of modified enzyme particles was in such a way that the large and dense enzyme particles were layered by other fine particles in the detergent powder,by means of a process known in the literature as“seeded granulation”.It is found that the homogeneity of modified enzyme particles could be improved significantly comparing to the original spherical enzyme particles in powder mixtures.Overall,the results of this research demonstrated that the segregation-induced properties of the dense/spherical enzyme particles could be lowered by altering their shape,which could enable the enzyme particles to behave almost similar to other ingredients during the pile formation process.     

Improving the hydrodynamic performance of the SUBOFF bare hull model: a CFD approach

The main aims of this study are to investigate the hydrodynamic performance of an autonomous underwater vehicle(AUV),calculate its hydrodynamic coefficients,and consider the flow characteristics of underwater bodies.In addition,three important parts of the SUBOFF bare hull,namely the main body,nose,and tail,are modified and redesigned to improve its hydrodynamic performance.A three-dimensional(3D)simulation is carried out using the computational fluid dynamics(CFD)method.To simulate turbulence,the k-ωshear stress transport(SST)model is employed,due to its good prediction capability at reasonable computational cost.Considering the effects of the length-to-diameter ratio(LTDR)and the nose and tail shapes on the hydrodynamic coefficients,it is concluded that a hull shape with bullet nose and sharp tail with LTDR equal to 7.14 performs better than the SUBOFF model.The final proposed model shows lower drag by about 14.9%at u=1.5 m·s^-1.Moreover,it produces 8 times more lift than the SUBOFF model at u=6.1 m·s^-1.These effects are due to the attachment of the fluid flow at the tail area of the hull,which weakens the wake region.     

Analysis of a two-grid method for semiconductor device problem

The mathematical model of a semiconductor device is governed by a system of quasi-linear partial differential equations.The electric potential equation is approximated by a mixed finite element method,and the concentration equations are approximated by a standard Galerkin method.We estimate the error of the numerical solutions in the sense of the Lqnorm.To linearize the full discrete scheme of the problem,we present an efficient two-grid method based on the idea of Newton iteration.The main procedures are to solve the small scaled nonlinear equations on the coarse grid and then deal with the linear equations on the fine grid.Error estimation for the two-grid solutions is analyzed in detail.It is shown that this method still achieves asymptotically optimal approximations as long as a mesh size satisfies H=O(h^1/2).Numerical experiments are given to illustrate the efficiency of the two-grid method.     

Shape of free-fall arch in quasi-2D silo

In this study,the optical flow method is used to measure the velocity distribution of a granular flow in a rectangular quasi-two-dimensional silo.Using the velocity gradient,a free-fall arch(FFA)is obtained and its geometric characteristics are calculated.A parabola-shaped FFA structure is discovered above the orifice in the steady flow state.The shape of the FFA affects the flow rate through the orifice.Furthermore,as jamming begins to occur,the geometry of the FFA disappears gradually from both sides and then from the middle;finally,the FFA disappears completely in the state of jamming.As the boundary between finite-stress and stress-free regions,the FFA facilitates further studies regarding the discontinuity of the stress area above the orifice.     

A drag force formula for heterogeneous granular flow systems based on finite average statistical method

The existing drag models are mostly based on the assumption of homogenous fluidization.However,the use of a homogeneous drag model to predict a heterogeneous granular flow system will cause a deviation.In this study,we developed a drag force model based on the assumption of heterogeneous fluidization.To prevent weakening of the heterogeneous characteristics in the drag force formula,we propose a finite average statistical method to filter the information of the heterogeneous granular cluster.The filtered information was used to fit the modified drag formula,which can reflect the heterogeneity of the granular cluster considering different configurations.A comparison shows that the new proposed drag formula filtered by the finite average statistical method fits well with energy minimization multi-scale simulation results.     

A numerical investigation of CO2 dilution on the thermochemical characteristics of a swirl stabilized diffusion flame

The turbulent combustion flow modeling is performed to study the effects of CO_2 addition to the fuel and oxidizer streams on the thermochemical characteristics of a swirl stabilized diffusion flame. A flamelet approach along with three well-known turbulence models is utilized to model the turbulent combustion flow field. The k-ω shear stress transport(SST) model shows the best agreement with the experimental measurements compared with other models. Therefore, the k-ω SST model is used to study the effects of CO_2 dilution on the flame structure and strength, temperature distribution, and CO concentration. To determine the chemical effects of CO_2 dilution, a fictitious species is replaced with the regular CO_2 in both the fuel stream and the oxidizer stream. The results indicate that the flame temperature decreases when CO_2 is added to either the fuel or the oxidizer stream. The flame length reduction is observed at all levels of CO_2 dilution. The H radical concentration indicating the flame strength decreases, following by the thermochemical effects of CO_2 dilution processes. In comparison with the fictitious species dilution, the chemical effects of CO_2 addition enhance the CO mass fraction. The numerical simulations show that when the dilution level is higher, the rate of the flame length reduction is more significant at low swirl numbers.     

Comparative study of two fluid model and dense discrete phase model for simulations of gas-solid hydrodynamics in circulating fluidized beds

Computational fluid dynamics(CFD)has become a valuable tool to study the complex gas-solid hydrodynamics in the circulating fluidized bed(CFB).Based on the two fluid model(TFM)under the Eulerian-Eulerian framework and the dense discrete phase model(DDPM)under the Eulerian-Lagrangian framework,this work conducts the comparative study of the gas-solid hydrodynamics in a CFB riser by these two different models.Results show that DDPM could be used to predict gas-solid hydrodynamics in the circulating fluidized bed,and there are differences between TFM and DDPM,especially in the radial distribution profiles of solid phase.Sensitivity analysis results show that the gas-solid drag model exhibits significant effects on the results for both the two models.The specularity coefficient and the restitution coefficient in the TFM,as well as the reflection coefficient and the parcel number in the DDPM,exhibit less impact on the simulated results.     

Static output feedback stabilization for second-order singular systems using model reduction methods

In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.     

Geometric and material nonlinearities of sandwich beams under static loads

A new theory developed from extended high-order sandwich panel theory(EHSAPT)is set up to assess the static response of sandwich panels by considering the geometrical and material nonlinearities simultaneously.The geometrical nonlinearity is considered by adopting the Green-Lagrange-type strain for the face sheets and core.The material nonlinearity is included as a piecewise function matched to the experimental stress-strain curve using a polynomial fitting technique.A Ritz technique is applied to solve the governing equations.The results show that the stress stiffening feature is well captured in the geometric nonlinear analysis.The effect of the geometric nonlinearity in the face sheets on the displacement response is more significant when the stiffness ratio of the face sheets to the core is large.The geometric nonlinearity decreases the shear stress and increases the normal stress in the sandwich core.By comparison with open literature and finite element simulations,the present nonlinear EHSAPT is shown to be sufficiently precise for estimating the nonlinear static response of sandwich beams by considering the geometric and material nonlinearities simultaneously.     

第七届全国周培源大学生力学竞赛试题

一、小球在高脚玻璃杯中的运动(20分) 一半球形高脚玻璃杯,半径r=5cm,其质量m1=0.3kg,杯底座半径R=5cm,厚度不计,杯脚高度h=10cm.如果有一个质量m2=0.1kg的光滑小球自杯子的边缘由静止释放后沿杯的内侧滑下,小球的半径忽略不计.     

平面Cosserat模型有限元分析的4和8节点单元与分片试验研究

经典连续体理论不包括物质内部尺度，当考虑应变软化问题时，有限元结果对网格具有很强的依赖性。与经典连续介质力学理论不同，Cosserat连续体模型在传统平动自由度的基础上添加了一独立的旋转自由度，在本构模型中引入了内尺度参数。本文研究了基于Cosserat理论的平面4和8节点等参元以及8（4）节点线、角位移混合插值等参单元，给出Cosserat单元分片试验的实施过程。最后将单元运用到小孔应力集中问题的分析当中，通过计算结果与理论解的比较，表明了4和8节点以及8（4）节点等参元的适用性，为问题的非线性分析打下基础。     

基于第二类四边形面积坐标的弹性力学纯弯解及MacNeal局限定理的破解

为了提高有限元的性能,弹性力学的解析解(齐次方程的通解)常常可用作有限元的试探函数。然而单元自由度数与完备的直角坐标解析解个数并不匹配,不完备的试函数会导致单元有方向依赖性。利用新型局部自然坐标——第二类四边形面积坐标QACM-II(S,T),给出了平面问题对应任意方向纯弯曲状态的应力函数解析解,即S3和T3的线性组合,并推导出了这两组应力函数对应的应力、应变和位移解析解。之后,利用QACM-II表示的解析解构造了非对称的平面4节点8自由度单元USQ4,该单元可以同时通过常应力/应变分片检验和纯弯测试,从而破解了MacNeal局限定理对平面低阶单元的限制。     

非协调四结点平面等参位移元新列式方法

本文导出一种等参协调元位移函数的新的表示方法,在此基础上建立起了构造等参非协调元的新方法。作为实例,构造出两个可以给出单刚显式的四结点平面非协调新单元。     

A penalty-hybrid finite element analysis of stokes flow

A type of penalty-hybrid variational principle is suggested for the analysis of Stokesian flow. On such a basis, a finite element model is formulated featuring, among others, a priori satisfaction of the deviatoric stress and hydrostatic pressure on linear momentum balance equations. Also in the present scheme the hydrostatic pressure is successfully eliminated at the element level, leaving only nodal velocities as solution unknowns. A series of 4-node and 8-node quadrilateral elements are derived and examined. Numerical examples demonstrating their characteristic behaviors are also included.The project Supported by National Natural Science Foundation of China     

Introducing new cracked finite elements and a method for SIF calculation of cracks

Two different types of 8-node cracked quadrilateral finite element are presented for fracture applications. The first element contains a central crack and the other one includes an edge crack. The introduced elements are applicable in 2D problems. The crack is not physically modeled within the element, but instead, its effects on the stiffness matrix are taken into account by utilizing linear fracture mechanics laws. Furthermore, a simple and practical procedure is proposed for calculation of stress intensity factor (SIF) by employing proposed cracked elements. Several numerical examples are presented to evaluate the capabilities of the proposed elements and procedure.     

基于高斯过程回归的有限元应力解的改善研究

本文应用高斯过程回归方法对有限元应力解进行了改善研究．考题是一简化为平面应力问题的各向同性且受均布载荷的等截面悬臂深梁，应力考察量取Mises 应力，高斯积分点为样本点，单元角结点为改善点．4结点单元有限元模型和8 结点单元有限元模型的计算结果表明：(1)改善点的总体误差比样本点的总体误差都小，且4 结点明显、8 结点不明显；(2)边界结点的改善效果均较传统整体应力修匀的效果显著；(3)改善点应力具有置信区间；(4)较传统分片应力修匀方法，高斯过程回归方法可将所选取区域内的所有角结点的应力同时给予改善，且边界角结点改善效果好．     

平面八节点四边形理性元

推导平面八节点四边形理性有限元列式，采用具有直到四次多项式的平面问题的微分方程的解作为插值函数，算例结果表明平面八节点四边形理性有限元的有效性。     

一类新的具有无外力圆柱表面的杂交应力元

根据Hellinger-Reissner原理,建立了进一步改进的具有一个无外力圆柱表面三维杂交应力元.元内假定应力场满足以柱坐标表示的平衡方程,及圆柱面上的无外力边界条件.当退化为二维时,也满足协调方程.数值算例表明,当分析带圆弧的槽孔板、块时,在稀疏的有限元网格下,这类单元即可提供较以前各类特殊元、一般假定位移元及一般假定应力元远为准确的三维及二维应力分布.     

Finite elements based on a first-order shear deformation moderate rotation shell theory with applications to the analysis of composite structures

The first-order shear deformation moderate rotation shell theory of Schmidt and Reddy [R. Schmidt and J. N. Reddy, J. Appl. Mech. 55, 611–617 (1988)] is used as a basis for the development of finite element models for the analysis of the static, geometrically non-linear response of anisotropic and laminated structures. The incremental, total Lagrangian formulation of the theory is developed, and numerical solutions are obtained by using the isoparametric Lagrangian 9-node and Serendipity 8-node shell finite elements. Various integration schemes (full, selective reduced, and uniformly reduced integration) are applied in order to detect and to overcome the effects of shear and membrane locking on the predicted structural response. A number of sample problems of isotropic, orthotropic, and multi-layered structures are presented to show the accuracy of the present theory. The von Kármán-type first-order shear deformation shell theory and continuum 2D theory are used for comparative analyses.