PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

Conversion between solid and beam element solutions of finite element method based on meta-modeling theory:development and application to a ramp tunnel structure

In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.

Treatment of discontinuous interface in liquid-solid forming with extended finite element method

Extended finite element method(XFEM) is proposed to simulate the discontinuous interface in the liquid-solid forming process.The discontinuous interface is an important phenomenon happening in the liquid-solid forming processes and it is difficult to be simulated accurately with conventional finite element method(CFEM) because it involves solid phase and liquid phase simultaneously.XFEM is becoming more and more popular with the need of solving the discontinuous problem happening in engineering field.The implementation method of XFEM is proposed on Abaqus code by using UEL(user element) with the flowchart.The key is to modify the element stiffness in the proposed method by using UEL on the platform of Abaqus code.In contrast to XFEM used in the simulation of solidification,the geometrical and physical properties of elements were modified at the same time in our method that is beneficial to getting smooth interface transition and precise analysis results.The analysis is simplified significantly with XFEM.

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

Coupling of discrete-element method and smoothed particle hydrodynamics for liquid-solid flows

Particle based methods can be used for both the simulations of solid and fluid phases in multiphase medium, such as the discrete-element method for solid phase and the smoothed particle hydrodynamics for fluid phase. This paper presents a computational method combining these two methods for solid-liquid medium. The two phases are coupled by using an improved model from a reported Lagrangian-Eulerian method. The technique is verified by simulating liquid-solid flows in a two-dimensional lid-driven cavity.

A boundary element method for the simulation of non-spherical bubbles and their interactions near a free surface

The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix are eliminated using coordinate transformation and so-called 4π rule. The solid angle for the open surface is treated in direct method based on its definition. Several kinds of configurations for the bubbles and free surface have been investigated. The pressure contours during the evolution of bubbles are obtained in our model and can better illuminate the mechanism underlying the motions of bubbles and free surface. The bubble dynamics and their interactions have close relation with the standoff distances, buoyancy parameters and initial sizes of bubbles. Completely different bubble shapes, free surface motions, jetting patterns and pressure distributions under different parameters can be observed in our model, as demon- strated in our calculation results.

Approach of technical decision-making by element flow analysis and Monte-Carlo simulation of municipal solid waste stream

This paper deals with the procedure and methodology which can be used to select the optimal treatment and disposal technology of municipal solid waste （MSW）, and to provide practical and effective technical support to policy-making, on the basis of study on solid waste management status and development trend in China and abroad. Focusing on various treatment and disposal technologies and processes of MSW, this study established a Monte-Carlo mathematical model of cost minimization for MSW handling subjected to environmental constraints. A new method of element stream （such as C, H, O, N, S） analysis in combination with economic stream analysis of MSW was developed. By following the streams of different treatment processes consisting of various techniques from generation, separation, transfer, transport, treatment, recycling and disposal of the wastes, the element constitution as well as its economic distribution in terms of possibility functions was identified. Every technique step was evaluated economically. The Mont-Carlo method was then conducted for model calibration. Sensitivity analysis was also carried out to identify the most sensitive factors. Model calibration indicated that landfill with power generation of landfill gas was economically the optimal technology at the present stage under the condition of more than 58% of C, H, O, N, S going to landfill. Whether or not to generate electricity was the most sensitive factor. If landfilling cost increases, MSW separation treatment was recommended by screening first followed with incinerating partially and composting partially with residue landfilling. The possibility of incineration model selection as the optimal technology was affected by the city scale. For big cities and metropolitans with large MSW generation, possibility for constructing large-scale incineration facilities increases, whereas, for middle and small cities, the effectiveness of incinerating waste decreases.

Cup-drawing formation of steel sheet with nickel coating by finite element method

The finite element method(FEM) simulation of deep-drawing of steel sheet with nickel coating based on the solid element and dynamic explicit method was reported. Penalty function method was used to treat the contact algorithm. The friction between the punch and coating sheet was based on a Coulomb formulation. The combination of coating and substrate was defined as tied with failure contact. The results of the simulation illustrate that the steel sheet and the nickel coating do not delaminate at the interface. The stress field of the nickel coating is more complicated than that of the steel substrate. Furthermore,it is found that the punch-nose radius is the most unsubstantial part for the intensity of the entire deep-drawing part and the thinnest part,it is a dangerous zone for the break. At this zone,the thickness thinning of the steel sheet and the nickel coating are up to 4.8% and 6.7%,respectively. Meanwhile,it is found that the curve of the variable blankholder force(VBHF) designed can improve the formability of sheet.

A method for 3D simulation of internal gas effects on thermal-mechanical behaviors in nuclear fuel elements

A new method for three-dimensional simulation of the interaction between the gas and the solid around is developed.The effects of the gas on the thermal-mechanical behaviors within the surrounded solid are performed by replacing the internal gas with an equivalent solid in the modeling,which can make it convenient to simulate the thermal-mechanical coupling effects in the solid research objects with gases in them.The applied thermal expansion coefficient,Young's modulus and Poisson's ratio of the equivalent solid material are derived.A series of tests have been conducted;and the proposed equivalent solid method to simulate the gas effects is validated.

基于Solidworks的激光去溢料机固定架有限元分析

通过Solidworks软件建立固定架重件的三维模型,依据物料抓取对机械手的运动状态,对其进行静态有限元分析,发现零件结构设计存在的缺陷并捉出有效的改进措施。通过模拟机械手动作过程中的振动情况,找出固定架零件在设备运行中可能发生破坏的频率。实践表明,采用有限元分析的方法进行零件结构强度和振动分析,有效地提高了机械设计产品的准确性和稳定性。

Numerical Analysis of Hydrodynamic Pressure Induced by Fluid-Solid Impact

As a further development of the authors' work (Huang and Qian, 1993), in this paper a new numerical method based on the time domain boundary element technique is proposed for solving fluid-solid coupling problems, in which a rigid body impacts normally on the calm surface of a half-space fluid. A fundamental solution to the half-space potential flow problem is first derived with the method of images. Then, an equivalent boundary integral equation in the Laplace transform domain is established by means of Green's second identity. Through the inverse Laplace transform and discretization in both time and boundary of the fluid region, the numerical calculation for the problem under consideration has been carried out. Several examples demonstrate that the present method is more efficient than existing ones, from which it is also seen that the shape of the impacting body has a considerable effect on the total impact force.

The effect of agglomerate on micro-structural evolution in solid-state sintering

Discrete element method （DEM） is used in the present paper to simulate the microstructural evolution of a planar layer of copper particles during sintering. Formation of agglomerates and the effect of their rearrangement on densification are mainly focused on. Comparing to the existing experimental observations, we find that agglomerate can form spontaneously in sintering and its rearrangement could accelerate the densification of compacts. Snapshots of numerical simulations agree qualitatively well with experimental observations. The method could be readily extended to investigate the effect of agglomerate on sintering in a three- dimensional model, which should be very useful for understanding the evolution of microstructure of sintering systems.

Fluid-solid Interaction Model for Hydraulic Reciprocating O-ring Seals

Elastohydrodynamic lubrication characteristics of hydraulic reciprocating seals have significant effects on sealing and tribology performances of hydraulic actuators, especially in high parameter hydraulic systems. Only elastic deformations of hydraulic reciprocating seals were discussed, and hydrodynamic effects were neglected in many studies. The physical process of the fluid-solid interaction effect did not be clearly presented in the existing fluid-solid interaction models for hydraulic reciprocating O-ring seals, and few of these models had been simultaneously validated through experiments. By exploring the physical process of the fluid-solid interaction effect of the hydraulic reciprocating O-ring seal, a numerical fluid-solid interaction model consisting of fluid lubrication, contact mechanics, asperity contact and elastic deformation analyses is constructed with an iterative procedure. With the SRV friction and wear tester, the experiments are performed to investigate the elastohydrodynamic lubrication characteristics of the O-ring seal. The regularity of the friction coefficient varying with the speed of reciprocating motion is obtained in the mixed lubrication condition. The experimental result is used to validate the fluid-solid interaction model. Based on the model, The elastohydrodynamic lubrication characteristics of the hydraulic reciprocating O-ring seal are presented respectively in the dry friction, mixed lubrication and full film lubrication conditions, including of the contact pressure, film thickness, friction coefficient, liquid film pressure and viscous shear stress in the sealing zone. The proposed numerical fluid-solid interaction model can be effectively used to analyze the operation characteristics of the hydraulic reciprocating O-ring seal, and can also be widely used to study other hydraulic reciprocating seals.

RESEARCH ON SOLID-LIQUID COUPLING DYNAMICSOF PIPE CONVEYING FLUID

On the basis of Hamilton principle. the equation of sonlid-liquid coupling vibration of pipe conveying fluid is deduced. An asymmetrical sonlid-liquid coupling damp matrix and a symmetrical solid-liquid coupling Stiffness matrix are obtained. Using QR method , pipe’s nature frequencies are calculated. The curves of the first four orders of natural frequency-flow velocity of pipe waw given .The influence of flowing velocity ,pressure, solid-liquid coupling damp and solid-liquid coupling stiffness on natural frequency are discussed respectively.The dynamic respondence of the pipes for stepload with different flow velocity are calculated by Newmark method .It is found that,with the flow velocity increased, the nature frequency of the pipes reduced, increased,reduced again and so on.

LIQUID-SOLID COUPLED SYSTEM OF MICROPUMP

This paper employs the integral-averaged method of thickness to approximate the periodical flows in a piezoelectric micropump, with a shallow water equation including nonlinearity and viscous damp presented to characterize the flows in micropump. The finite element method is used to obtain a matrix equation of fluid pressure. The fluid pressure equation is combined with the vibration equation of a silicon diaphragm to construct a liquid-solid coupled equation for reflecting the interaction between solid diaphragm and fluid motion in a micropump. Numerical results of a mode analysis of the coupled system indicate that the natural frequencies of the coupled system are much lower than those of the non-coupled system. The influence of additional mass and viscous damp of fluid on the natural frequencies of the coupled system is more significant as the pump thickness is small. It is found that the vibration shape functions of silicon diaphragm of the coupled system are almost the same as those of the non-coupled system. This paper also gives the first-order amplitude-frequency relationship of the silicon diaphragm, which is necessary for the flow-rate-frequency analysis of a micropump.

Safety analysis on truing process of solid propellant

In order to investigate safety character of solid propellant, the truing process was analyzed by a thermalmechanical coupled finite element method. Based on the viscoelastic integral constitutive relation, the temperature response in truing process was calculated. The result showed that the maximum temperature in solid propellant truing process could attain up to 65 ℃ when the truing velocity was 10 mm/s and the truing thickness was 1 mm, and there was no remarkable temperature change under the truing thickness of 1.0 mm, 1.5 mm and 2.5 mm. The hazard in truing process of solid propellant was mainly close to the edge of cutting area.

CONSERVATION LAWS IN FINITE MICROCRACKING BRITTLE SOLIDS

This paper addresses the conservation laws in finite brittle solids with microcracks. The discussion is limited to the 2-D cases. First, after considering the combination of the Pseudo-Traction Method and the indirect Boundary Element Method, a versatile method for solving multi-crack interacting problems in finite plane solids is proposed, by which the fracture parameters （SIF and path-independent integrals） can be calculated with a desirable accuracy. Second, with the aid of the method proposed, the roles the conservation laws play in the fracture analysis for finite microcracking solids are studied. It is concluded that the conservation laws do play important roles in not only the fracture analysis but also the analysis of damage and stability for the finite microcracking system. Finally, the physical interpretation of the M-integral is discussed further. An explicit relation between the M-integral and the crack face area, i.e., M = GS, has been discovered using the analytical method, which can shed some light on the Damage Mechanics issues from a different perspective.

基于微观流固耦合的超深层致密砂岩气藏应力敏感性分析

超深层致密砂岩气藏具有强应力敏感性,目前常用的研究方法包括压力脉冲实验法和实时在线CT扫描法两种,压力脉冲实验法不能揭示产生应力敏感性的微观机理,而实时在线CT扫描法也无法模拟深部地层高压、高应力的条件。为解决应力敏感性实验研究的不足,基于离散单元法与管道网络模型建立了微观流固耦合算法,编制了模拟器,并对模拟器力学计算和流固耦合模块的正确性进行了验证,分析了应力大小、加载方向和孔隙压力对岩心渗透率的影响,最后从微观上揭示了超深层致密砂岩气藏的应力敏感性机理。研究结果表明:①应力通过增加与之垂直方向上喉道两侧的法向压力,减小喉道的水力半径,进而降低储层的渗透率;②较高的孔隙压力能够阻碍岩石颗粒在应力作用下的移动,从而减缓了孔隙和喉道的变形,使模型保持较高的渗透率;③致密砂岩气藏的渗透率受到应力和地层压力的共同控制,并且具有各向异性,在垂直于最小主应力方向上形成渗透率较大的优势通道;④异常高压阻碍了地应力的压实作用,有利于保护储层孔隙,使地层有较好的储集性能和较高的渗透率。结论认为,根据离散元法结合孔隙网络模型建立的流固耦合方法可为理解超深层致密砂岩应力敏感性提供理论参考,并为超深层致密砂岩气藏的科学高效开发提供指导。

Multi-scale Simulation on Bonding Mechanism of Solid-Liquid Cast-Rolling of Cu/Al Cladding Strip based on FEM and MD

To explore the complex thermal-mechanical-chemical behavior in the solid-liquid cast-roll bonding(SLCRB) of Cu/Al cladding strip, numerical simulations were conducted from both macro and micro scales. In macro-scale, with birth and death element method, a thermo-mechanical coupled finite element model(FEM) was set up to explore the temperature and contact pressure distribution at the Cu/Al bonding interface in the SLCRB process. Taking these macro-scale simulation results as boundary conditions, we simulated the atom diffusion law of the bonding interface by molecular dynamics(MD) in micro-scale. The results indicate that the temperature in Cu/Al bonding interface deceases from 700 to 320 ℃ from the entrance to the exit of caster, and the peak of contact pressure reaches up to 140 MPa. The interfacial diffusion thickness depends on temperature and rolling reduction, higher temperature results in larger thickness, and the rolling reduction below kiss point leads to significant elongation deformation of cladding strip which yields more newborn interface with fresh metal and make the diffusion layer thinner. The surface roughness of Cu strip was found to be benefit to atoms diffusion in the Cu/Al bonding interface. Meanwhile, combined with the SEM-EDS observation on the microstructure and composition in the bonding interface of the experimental samples acquired from the castrolling bite, it is revealed that the rolling reduction and severe elongation deformation in the solid-solid contact zone below kiss point guarantee the satisfactory metallurgical bonding with thin and smooth diffusion layer. The bonding mechanisms of reactive diffusion, mechanical interlocking and crack bonding are proved to coexist in the SLCRB process.