TWO-SCALE FEM FOR ELLIPTIC MIXED BOUNDARY VALUE PROBLEMS WITH SMALL PERIODIC COEFFICIENTS

Presents a study which proposed a dual approximate expression of the exact solution for mixed boundary value problems of second order elliptic partial differential equation with small periodic coefficients. Dual approximate error estimates; Main results.

Study on vibration reduction of two-scale system coupled with dynamic vibration absorber

The dynamic vibration absorber with inerter and grounded stiffness(IGDVA)is used to control a two-scale system subject to a weak periodic perturbation.The vibration suppression effect is remarkable.The amplitude of the main system coupled with absorber is significantly reduced,and the high frequency vibration completely disappears.First,through the slow-fast analysis and stability theory,it is found that the stability of the autonomous system exerts a notable regulating effect on the vibration response of the non-autonomous system.After adding the dynamic vibrator absorber,the center in the autonomous system changes to an asymptotically stable focus,consequently suppressing the vibration in the non-autonomous system.Further research reveals that the parameters of the absorber affect the real parts of the eigenvalues of the autonomous system,thereby regulating the stability of the system.Transitioning from a qualitative standpoint to a quantitative approach,a comparison of the solutions before and after the introduction of the dynamic absorber reveals that,when the grounded stiffness ratio and the mass ratio of the dynamic absorber are not equal,the high-frequency part in the analytical solution disappears.As a result,this leads to a reduction in the amplitude of the trajectory,achieving a vibration reduction effect.

Acid-rock reaction kinetics in a two-scale model based on reaction order correction

The reaction order plays a crucial role in evaluating the response rate of acid-rock.However,the conventional two-scale model typically assumes that the reaction order is constant as one,which can lead to significant deviations from reality.To address this issue,this study proposes a novel multi-order dynamic model for acid-rock reaction by combining rotating disk experimental data with theoretical derivation.Through numerical simulations,this model allows for the investigation of the impact of acidification conditions on different orders of reaction,thereby providing valuable insights for on-site construction.The analysis reveals that higher response orders require higher optimal acid liquid flow rates,and lower optimal H+diffusion coefficients,and demonstrate no significant correlation with acid concentration.Consequently,it is recommended to increase the displacement and use high-viscosity acid for reservoirs with high calcite content,while reducing the displacement and using low-viscosity acid for reservoirs with high dolomite content.

Concurrent Two-Scale Topology Optimization of Thermoelastic Structures Using a M-VCUT Level Set Based Model of Microstructures

By analyzing the results of compliance minimization of thermoelastic structures,we observed that microstructures play an important role in this optimization problem.Then,we propose to use a multiple variable cutting(M-VCUT)level set-based model of microstructures to solve the concurrent two-scale topology optimization of thermoelastic structures.A microstructure is obtained by combining multiple virtual microstructures that are derived respectively from multiple microstructure prototypes,thus giving more diversity of microstructure and more flexibility in design optimization.The effective mechanical properties of microstructures are computed in an off-line phase by using the homogenization method,and then a mapping relationship between the design variables and the effective properties is established,which gives a data-driven model of microstructure.In the online phase,the data-driven model is used in the finite element analysis to improve the computational efficiency.The compliance minimization problem is considered,and the results of numerical examples prove that the proposed method is effective.

基于FEM-SPG耦合的盾构滚刀破岩仿真研究

隧道施工过程中岩石在盾构滚刀破岩力的作用下会发生破碎,传统的有限元方法(FEM)无法探究这一现象。为解决使用FEM面临网格畸变及无法同时分析岩石破碎的问题,引入光滑粒子伽辽金(SPG)分析方法,通过FEM和SPG耦合的方式建立了滚刀回转切削岩石模型。结果表明:FEM-SPG方法能够有效地展示岩石的密实核破碎过程及岩渣飞溅的现象。通过与CSM模型进行对比验证,发现仿真得到的滚刀垂直力与CSM模型相比误差为1.8%,验证了模型的可靠性。该方法与模型能够为滚刀精准破岩研究及分析岩石破碎机理提供新思路。

基于声学FEM-BEM的户内变电站噪声场计算

准确计算户内变电站大型、复杂的噪声场分布,进而评价可采用降噪措施的减噪效果,是解决户内变电站噪声污染的关键问题。为此,综合声学有限元法(finite element method,FEM)求解复杂声场收敛性好及精度高的优点,及声学边界元法(boundary element method,BEM)降维求解大型声场的优势,提出了一种基于声学FEM-BEM的户内变电站噪声场求解算法。首先,建立变电站内部声源声固耦合模型,采用声学FEM求解混响噪声作用下的声固耦合响应;然后,基于声学FEM-BEM耦合理论,求解内、外耦合边界处结构单元受声固耦合激励产生的位移及应力载荷;最后,根据声压及应力载荷激发的外场声波扩散模型,基于常规Gauss数值积分法,建立外部空间声域2维BEM声学积分方程,求解外部声场。该算法在湖南某110 kV户内变电站噪声场的求解分析中得到了成功应用,与实测值的相对误差为3.61%~4.87%。

基于DEM-FEM耦合方法的碎屑流对桥墩最大冲击力的影响因素研究

碎屑流是我国山区最危险的地质灾害之一,山区桥墩常受到碎屑流冲击而开裂、倾斜甚至倒塌,给山区桥梁建设、运营带来严重的安全隐患。采用离散元方法(discrete element method,DEM)和有限元方法(finite element method,FEM)耦合的三维数值模拟方法模拟了碎屑流对双柱式桥墩的冲击效应,并结合斜槽试验,验证了耦合方法的准确性,进一步分析了碎屑流冲击坡度、距离和体积密度对桥墩冲击力的影响规律。结果表明,最大冲击力与碎屑流冲击坡度、距离和体积密度分别呈幂函数(指数大于1)、幂函数(指数小于1)和线性正相关。冲击坡度、距离和体积密度对最大冲击力的敏感度值分别为3.012、0.202、0.804,在桥梁碎屑流灾害防治时需重视冲击坡度和体积密度的影响。将冲击力的数值模拟值与流体动力学模型预测值对比分析表明,流体动力学模型理论公式能较好地预测桥墩所受的最大冲击力,最大预测误差低于23.6%。相关研究结果可为山区桥梁碎屑流灾害防治与设计提供一定的参考依据。

A two-scale method for identifying mechanical parameters of composite materials with periodic configuration

In this paper, a two-scale method （TSM） is presented for identifying the mechanics parameters such as stiffness and strength of composite materials with small periodic configuration. Firstly, a formulation is briefly given for two-scale analysis （TSA） of the composite materials. And then a two-scale computation formulation of strains and stresses is developed by displacement solution with orthotropic material coefficients for three kinds of such composites structures, i.e., the tension column with a square cross section, the bending cantilever with a rectangular cross section and the torsion column with a circle cross section. The strength formulas for the three kinds of structures are derived and the TSM procedure is discussed. Finally the numerical results of stiffness and strength are presented and compared with experimental data. It shows that the TSM method in this paper is feasible and valid for predicting both the stiffness and the strength of the composite materials with periodic configuration.

考虑强度速率衰减效应的地震滑坡SPH-FEM模拟

基于SPH和FEM耦合的数值计算方法,引入滑动面摩擦强度的速率衰减模型,提出了一种能够模拟地震诱发滑坡破坏过程的数值模拟方法。基于所提数值方法模拟了唐家山地震滑坡,模拟结果与现场勘查结果及室内试验现象较为一致。基于模拟的滑动面上摩擦系数的演化过程,将唐家山滑坡的发生分为4个阶段:启动阶段、摩擦衰减阶段、低摩擦滑移阶段和逐步稳定阶段。模拟结果表明速度增加和摩擦强度衰减的相互促进,是触发滑体的高速运动的根本原因。提出采用滑体上作用的动摩擦力fd和动下滑力Td的比值R作为判别指标用于判断大型滑坡的启动,当首次出现R小于1时认为滑动面发生整体贯通并出现失稳启动。基于滑动面不同位置摩擦系数的演化,揭示了滑坡启动中滑动面摩擦强度衰减和滑动面的渐进贯通过程,解释了地震作用与滑动面摩擦参数速率衰减效应共同作用触发大型滑坡发生破坏的内在机理。

A nonlinear k_p-ε_p particle two-scale turbulence model and its application

A particle nonlinear two-scale kp-εp turbulence model is proposed for simulating the anisotropic turbulent two-phase flow. The particle kinetic energy equation for two-scale fluctuation, particle energy transfer rate equation for large-scale fluctuation, and particle turbulent kinetic energy dissipation rate equation for small-scale fluctuation are derived and closed. This model is used to simulate gas-particle flows in a sudden-expansion chamber. The simulation is com- pared with the experiment and with those obtained by using another two kinds of tow-phase turbulence model, such as the single-scale k-ε two-phase turbulence model and the particle two-scale second-order moment （USM） two-phase turbulence model. It is shown that the present model gives simulation in much better agreement with the experiment than the single-scale k-ε two-phase turbulence model does and is almost as good as the particle two-scale USM turbu-lence model.

Improved Microwave Ocean Emissivity and Reflectivity Models Derived from Two-Scale Roughness Theory

The Geometrical Optics(GO)approach and the FAST Emissivity Model(FASTEM)are widely used to estimate the surface radiative components in atmospheric radiative transfer simulations,but their applications are limited in specific conditions.In this study,a two-scale reflectivity model(TSRM)and a two-scale emissivity model(TSEM)are developed from the two-scale roughness theory.Unlike GO which only computes six non-zero elements in the reflectivity matrix,The TSRM includes 16 elements of Stokes reflectivity matrix which are important for improving radiative transfer simulation accuracy in a scattering atmosphere.It covers the frequency range from L-to W-bands.The dependences of all TSRM elements on zenith angle,wind speed,and frequency are derived and analyzed in details.For a set of downwelling radiances in microwave frequencies,the reflected upwelling brightness temperature(BTs)are calculated from both TSRM and GO and compared for analyzing their discrepancies.The TSRM not only includes the effects of GO but also accounts for the small-scale Bragg scattering effect in an order of several degrees in Kelvins in brightness temperature.Also,the third and fourth components of the Stokes vector can only be produced from the TSRM.For the emitted radiation,BT differences in vertical polarization between a TSEM and FASTEM are generally less than 5 K when the satellite zenith angle is less than 40°,whereas those for the horizontal component can be quite significant,greater than 20 K.

Two-scale finite element method for piezoelectric problem in periodic structure

The prediction of the mechanical and electric properties of piezoelectric fibre composites has become an active research area in recent years. By means of introducing a boundary layer problem, some new kinds of two-scale finite element methods for solutions to the electric potential and the displacement for composite material in periodic struc- ture under the coupled piezoelectricity are derived. The coupled two-scale relation of the electric potential and the displacement is set up, and some finite element approximate estimates and numerical examples which show the effectiveness of the method are presented.

饱和地基水泥土复合桩近场主动隔振的BEM-FEM耦合分析

人工振动污染是当前城市发展过程中面临的一个难题,单排或多排桩是人工振动污染防治的一种重要方式。传统基桩施工容易造成泥浆污染、振动等各种负面环境问题,而水泥土复合桩是一种低振动、低噪声、无挤土、少排泥浆的新型基桩,特别适合在城区构建排桩屏障。针对饱和地基中单排水泥土复合桩的近场主动隔振问题,建立了动力机器基础环境振动影响的半解析BEM方程;在此基础上,分别采用半解析BEM对饱和地基、FEM对水泥土复合桩进行建模,根据饱和地基-水泥土桩接触面的平衡和相容条件,建立了水泥土复合桩近场主动隔振的半解析BEM-FEM耦合方程,给出了耦合方程稀疏矩阵的存储策略,并对单排水泥土复合桩的近场主动隔振效果进行了计算分析。研究结果表明:单排水泥土复合桩能够有效地对动力机器基础引起环境振动进行隔振;等长芯桩复合桩的隔振效果要优于短芯桩。单排水泥土复合桩的隔振效果随着桩数的增加而逐渐提高;较小桩间距并不一定取得更好的隔振效果,建议相邻桩桩间净距取2.0~2.5λ_(R)(Rayleigh波波长)。随着单排水泥土复合桩距振源距离的增加,屏障隔振效果逐渐降低,但降低幅度相对较小。此外,对水泥土复合桩而言,当内插预制桩的外轮廓尺寸相同时,预制芯桩型式对水泥土复合桩隔振效果影响相对较小,可根据实际工程条件选择合适的芯桩类型。

A two-scale second-order moment two-phase turbulence model for simulating dense gas-particle flows

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision and through a unified treatment of these two kinds of fluctuations. The proposed model is used to simulate gas-particle flows in a channel and in a downer. Simulation results are in agreement with the experimental results reported in references and are near the results obtained using the sin- gle-scale second-order moment two-phase turbulence model superposed with a particle collision model （USM-θ model） in most regions.

基于DEM-FEM耦合的环锤式破碎机反击板可靠性分析

由于环锤式破碎机破碎物料的过程是一个极其复杂多变的过程,且与反击板所受的碰撞力是不连续、无规则的,因此,无法使用公式推导或使用单一软件对反击板受力进行计算,为了能够更好地确定破碎机破碎过程中反击板的载荷大小、分布和寿命,文中提出了一种基于离散元(DEM)和有限元(FEM)耦合计算的环锤式破碎机反击板可靠性分析方法,通过EDEM建立环锤式破碎机破碎物料仿真过程,并提取出破碎过程中反击板的载荷分布,与ANSYS Workbench建立耦合计算通道,计算得到在正常工作参数下反击板的受力大小和位置,在此基础上采用n Code Design Life疲劳分析软件对反击板进行寿命分析,得到反击板的受力云图和损伤云图,研究转子转速、反击板倾角、转子与筛板距离等不同因素下对反击板可靠性的影响,结果表明:转子转速为1000 r∕min时反击板的最大受力值为247.9 MPa、反击板角度为60°时应力最大值为206.48 MPa、转子与筛板距离为530 mm时受力最大值为268.37 MPa。

EFFECT OF EMPIRICAL COEFFICIENTS ON SIMULATION IN TWO-SCALE SECOND-ORDER MOMENT PARTICLE-PHASE TURBULENCE MODEL

A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concept of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision. The proposed model is used to simulate gas-particle downer reactor flows. The computational results of both particle volume fraction and mean velocity are in agreement with the experimental results. After analyzing effects of empirical coefficient on prediction results, we can come to a conclusion that, inside the limit range of empirical coefficient, the predictions do not reveal a large sensitivity to the empirical coefficient in the downer reactor, but a relatively great change of the constants has important effect on the prediction.

Derivation of a second-order model for Reynolds stress using renormalization group analysis and the two-scale expansion technique

With the two-scale expansion technique proposed by Yoshizawa,the turbulent fluctuating field is expanded around the isotropic field.At a low-order two-scale expansion,applying the mode coupling approximation in the Yakhot-Orszag renormalization group method to analyze the fluctuating field,the Reynolds-average terms in the Reynolds stress transport equation,such as the convective term,the pressure-gradient-velocity correlation term and the dissipation term,are modeled.Two numerical examples：turbulent flow past a backward-facing step and the fully developed flow in a rotating channel,are presented for testing the efficiency of the proposed second-order model.For these two numerical examples,the proposed model performs as well as the Gibson-Launder （GL） model,giving better prediction than the standard k-ε model,especially in the abilities to calculate the secondary flow in the backward-facing step flow and to capture the asymmetric turbulent structure caused by frame rotation.

Second-order two-scale computational method for ageing linear viscoelastic problem in composite materials with periodic structure

The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspon- dence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale （SOTS） computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calcu- lating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement ap- proximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.

Second-order two-scale computations for conductive radiative heat transfer problem in periodic porous materials

In this paper, a kind of second-order two-scale （SOTS） computation is developed for conductive-radiative heat trans- fer problem in periodic porous materials. First of all, by the asymptotic expansion of the temperature field, the cell problem, homogenization problem, and second-order correctors are obtained successively. Then, the corresponding finite element al- gorithms are proposed. Finally, some numerical results are presented and compared with theoretical results. The numerical results of the proposed algorithm conform with those of the FE algorithm well, demonstrating the accuracy of the present method and its potential applications in thermal engineering of porous materials.

Polarimetric entropy of the ocean surface with a two-scale scattering model

The relationships among an ocean wave spectrum,a fully polarimetric coherence matrix,and radar parameters are deduced with an electromagnetic wave theory.Furthermore,the relationship between the polarimetric entropy and ocean wave spectrum is established based on the definition of entropy and a twoscale scattering model of the ocean surface.It is the first time that the polarimetric entropy of the ocean surface is presented in theory.Meanwhile,the relationships among the fully polarimetric entropy and the parameters related to radar and ocean are discussed.The study is the basis of further monitoring targets on the ocean surface and deriving oceanic information with the entropy from the ocean surface.The contrast enhancement between human-made targets and the ocean surface with the entropy is presented with quad-pol airborne synthetic aperture radar（AIRSAR） data.