Fourth-order phase-field modeling for brittle fracture in piezoelectric materials

Failure analyses of piezoelectric structures and devices are of engineering and scientific significance.In this paper,a fourth-order phase-field fracture model for piezoelectric solids is developed based on the Hamilton principle.Three typical electric boundary conditions are involved in the present model to characterize the fracture behaviors in various physical situations.A staggered algorithm is used to simulate the crack propagation.The polynomial splines over hierarchical T-meshes(PHT-splines)are adopted as the basis function,which owns the C1continuity.Systematic numerical simulations are performed to study the influence of the electric boundary conditions and the applied electric field on the fracture behaviors of piezoelectric materials.The electric boundary conditions may influence crack paths and fracture loads significantly.The present research may be helpful for the reliability evaluation of the piezoelectric structure in the future applications.

Research on Three-Dimensional Simulation of the Internal Arc Gear Skiving

Aiming at the problems that the simulation accuracy which is reduced due to the simplification of the model,a three-dimensional simulation method based on solid modeling is being proposed.By analyzing the motion relationship and positional relationship between the caries knife and the workpiece,the coordinate system of the caries machining was established.With the MATLAB software,the cutting edge model and the blade sweeping surface model of the boring cutter are sequentially established.Boolean operation is performed on the blade swept surface formed by the tooth cutter teeth with time t and the workpiece tooth geometry as well as the undeformed three-dimensional chip geometry model and the instantaneous cogging geometry model are obtained at different times.Through the compare between gear end face simulation tooth profile and the theoretical inner arc tooth profile,we verified the accuracy and rationality of the proposed method.

Model Prediction and Optimal Control of Gas Oxygen Content for A Municipal Solid Waste Incineration Process

In the municipal solid waste incineration process,it is difficult to effectively control the gas oxygen content by setting the air flow according to artificial experience.To address this problem,this paper proposes an optimization control method of gas oxygen content based on model predictive control.First,a stochastic configuration network is utilized to establish a prediction model of gas oxygen content.Second,an improved differential evolution algorithm that is based on parameter adaptive and t-distribution strategy is employed to address the set value of air flow.Finally,model predictive control is combined with the event triggering strategy to reduce the amount of computation and the controller's frequent actions.The experimental results show that the optimization control method proposed in this paper obtains a smaller degree of fluctuation in the air flow set value,which can ensure the tracking control performance of the gas oxygen content while reducing the amount of calculation.

Modeling the dynamic process of tsunami earthquake by liquid-solid coupling model

Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock absorber.The acceleration at the sea shore is about 1/10 as large as at the epicenter.The maximum vertical velocity at the epicenter is 1.4 times that on the fault.The maximum vertical displacement at the fault is less than that at the epicenter.The difference between them is the amplitude of the tsunami at the epicenter.The time of the maximum displacement to occur on the fault is not at the beginning of the fault slipping but retards 23 s.

Application of 3-D Geoscience Modeling Technology for the Estimation of Solid Mineral Reserves

Applying new approaches, methods, and technologies for the estimation of reserves can effectively improve the efficiency and accuracy of assessments of solid mineral resources. After analyzing the development of 3-D geoscience modeling technology （3-D GMT）, this paper discusses the application of 3-D GMT for the estimation of solid mineral reserves, emphatically introducing its workflow and two key technologies, 3-D orebody surface modeling, and property modeling. Moreover, the paper analyzes the limitations of traditional methods, such as the section method and geological block method, and points out the advantages of 3-D GMT： building more accurate 3-D orebody models, expressing the internal inhomogeneous attributes of an orebody, reducing the potential for errors in the estimation of reserves, and implementing dynamic estimations of reserves.

Three-Dimensional Tidal Model and Its Application to Numerical Simulation of Water Quality in Coastal Waters

The turbulence mechanism plays an important part in the mixing process and momentum transfer of turbulence. A three-dimensional Prandtl mixing length tidal model has been developed to simulate tidal flows and water quality. The eddy viscosities and diffusivities are computed from the Prandtl mixing length model. In order to model the water quality of an estuary or coastal area many interdependent processes need to be simulated. These may be conveniently separated into three main groups: transport and mixing processes, biochemical interaction of water quality variables and the utilization and re-cycling of nutrients by living matter. The model simulates full oxygen and nutrient balance, primary productivity and the transport, reaction mechanism and fate of pollutants over tidal time-scales. The model is applied to numerical simulation of tidal flows and water quality in Dalian Bay. The model has been calibrated against a limited data set of historical water quality observations and in general demonstrates excellent agreement with all available data.

The three-dimension model for the rock-breaking mechanism of disc cutter andanalysis of rock-breaking forces

To study the rock deformation with three- dimensional model under rolling forces of disc cutter, by car- rying out the circular-grooving test with disc cutter rolling around on the rock, the rock mechanical behavior under rolling disc cutter is studied, the mechanical model of disc cutter rolling around the groove is established, and the the- ory of single-point and double-angle variables is proposed. Based on this theory, the physics equations and geometric equations of rock mechanical behavior under disc cutters of tunnel boring machine （TBM） are studied, and then the bal- ance equations of interactive forces between disc cutter and rock are established. Accordingly, formulas about normal force, rolling force and side force of a disc cutter are de- rived, and their validity is studied by tests. Therefore, a new method and theory is proposed to study rock- breaking mech- anism of disc cutters.

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.

Modeling for Collapsing Cavitation Bubble near Rough Solid Wall by Mulit-Relaxation-Time Pseudopotential Lattice Boltzmann Model

Cavitation bubble collapse near rough solid wall is modeled by the multi-relaxation-time (MRT) pseudopotential lattice Boltzmann (LB) model. The modified forcing scheme, which can achieve LB model’s thermodynamic consistency by tuning a parameter related with the particle interaction range, is adopted to achieve desired stability and density ratio. The bubble collapse near rough solid wall was simulated by the improved MRT pseudopotential LB model. The mechanism of bubble collapse is studied by investigating the bubble profiles, pressure field and velocity field evolution. The eroding effects of collapsing bubble are analyzed in details. It is found that the process and the effect of the interaction between bubble collapse and rough solid wall are affected seriously by the geometry of solid boundary. At the same time, it demonstrates that the MRT pseudopotential LB model is a potential tool for the investigation of the interaction mechanism between the collapsing bubble and complex geometry boundary.

Modeling and Control of the Modular Multilevel Converter (MMC) based Solid State Transformer (SST) with Magnetic Integration

Solid state transformer(SST)can provide more advanced functionalities compared with conventional transformer,and has great potential in smart grid application.Recently,the SST with medium frequency(MF)isolation link and magnetic integration feature has been proposed,which can reduce the system volume and thus increase the power density.However,the magnetic integration also introduces strong coupling between the line frequency(LF)and MF variables,which poses a great challenge on modeling and control issues.This paper proposes a modeling and control method for an SST with magnetic integration and mixed-frequency modulation.A mathematical model based on dual d-q references is deduced,and then a cascaded control system is designed according to the model.Parameters of the controller for the variables at one frequency are properly designed to avoid disturbance from the variables at the other frequency.The simulation and experimental results show good decoupling effect and satisfactory dynamics performance of the proposed control system.

A new numerical approach of coupled modeling for solid deformation and gas leak flow in multi-coal-seams

From the viewpoint of interaction mechanics for solid and gas, a coupled mathematical model was presented for solid coal/rock deformation and gas leak flow in parallel deformable coal seams. Numerical solutions using the SIP (Strong Implicit Proce- dure) method to the coupled mathematical model for double parallel coal seams were also developed in detail. Numerical simulations for the prediction of the safety range using protection layer mining were performed with experimental data from a mine with potential danger of coal/gas outbursts. Analyses show that the numerical simulation results are consistent with the measured data in situ.

Modeling of gas-solid flow in a CFB riser based on computational particle fluid dynamics

A three-dimensional model for gas-solid flow in a circulating fluidized bed（CFB） riser was developed based on computational particle fluid dynamics（CPFD）.The model was used to simulate the gas-solid flow behavior inside a circulating fluidized bed riser operating at various superficial gas velocities and solids mass fluxes in two fluidization regimes,a dilute phase transport（DPT） regime and a fast fluidization（FF） regime.The simulation results were evaluated based on comparison with experimental data of solids velocity and holdup,obtained from non-invasive automated radioactive particle tracking and gamma-ray tomography techniques,respectively.The agreement of the predicted solids velocity and holdup with experimental data validated the CPFD model for the CFB riser.The model predicted the main features of the gas-solid flows in the two regimes;the uniform dilute phase in the DPT regime,and the coexistence of the dilute phase in the upper region and the dense phase in the lower region in the FF regime.The clustering and solids back mixing in the FF regime were stronger than those in the DPT regime.

Finding Intersection Points of Curves and Surfaces for Solid Modeling in Curved-Surface Environment

Intersecting is an important factor which influences the effociency androbustness of Boolean algorithms in solid modeling based on surved-surfaces,andintersecting algorithms are closely related to geometric representations of curved-surfaces.Although surfaces can be commonly represented with NURBS,unnecessary complexitiesare caused in the intersecting of quadric surfaces.Quadrics are frequently used to des-cribe geometric features of shafts,holes and grooves etc.in mechanical part designing,therefore;their intersection algorithms are required to have higher accuracy,higher efficiency and higher robustness.For this reason,a practical representation ofquadric surfaces is studied in detail,and on the basis of that,algorithms of intersectingpoints are developed between quadric suraces and their boundaies,i.e.,conics,quarticnonplanar space curves.

Modeling transportation of suspended solids in Zhujiang River estuary, South China

A three-dimensional transportation model for suspended solids (SS) in Zhujiang (Pearl) River estuary, South China, was developed by coupling with a three-dimensional hydrodynamic model. The model was validated using hourly measured data of sediment contents during 25–26, July 1999. The results showed that modeled contents matched well with measured ones and that the modeled top layer distribution agreed with the remotely sensed image of suspended solids in summer. The modeled results showed clearly the layers of sus- pended solids in depth, with larger sediment contents in lower layers though in the interface between salt water and freshwater the lowest contents appeared in middle layer. In overall, the suspended solids inflow from 8 rivers, transport southwestward, and carried by strong coastal flow in Zhujiang River estuary. Contours of sediment contents in the estuary spread further to the open sea during ebb tide rather than flood tide which reflects that the suspended solids in the estuary are land sourced.

Three-dimensional Simulation of Gas/Solid Flow in Spout-fluid Beds with Kinetic Theory of Granular Flow

A three-dimensional Eulerian multiphase model, with closure law according to the kinetic theory of granular flow, was used to study the gas/solid flow behaviors in spout-fluid beds. The influences of the coefficient of restitution due to non-ideal particle collisions on the simulated results were tested. It is demonstrated that the simulated result is strongly affected by the coefficient of restitution. Comparison of simulations with experiments in a small spout-fluid bed showed that an appropriate coefficient of restitution of 0.93 was necessary to simulate the flow characteristics in an underdesigned large size of spout-fluid bed coal gasifier with diameter of 1m and height of 6m. The internal jet and gas/solid flow patterns at different operating conditions were obtained. The simulations show that an optimal gas/solid flow pattern for coal gasification is found when the spouting gas flow rate is equal to the fluidizing gas flow rate and the total of them is two and a half times the minimum fluidizing gas flow rate. Be-sides, the radial distributions of particle velocity and gas velocity show similar tendencies; the radial distributions of particle phase pressure due to particle collisions and the particle pseudo-temperature corresponding to the macro-scopic kinetic energy of the random particle motion also show similar tendencies. These indicate that both gas drag force and particle collisions dominate the movement of particles.

Trimming self-intersections in swept volume solid modeling

Swept volume solid modeling has been applied to many areas such as NC machining simulation and verification, robot workspace analysis, collision detection, and CAD. But self-intersections continue to be a challenging problem in the boundary representation of swept volume solids. A novel algorithm is presented in this paper to trim self-intersection regions in swept volume solids modeling. This trimming algorithm consists of two major steps: (1) roughly detecting self-intersection regions by checking intersections or overlapping of the envelop profiles; (2) splitting the whole envelop surfaces of the swept volume solid into separate non-self-intersecting patches to trim global self-intersections, and to trim local self-intersections, dividing local self-intersecting regions into patches and replacing self-intersecting patches with non-self-intersecting ones. Examples show that our algorithm is efficient and robust.

One-dimensional Dynamic Modeling and Simulation of a Planar Direct Internal Reforming Solid Oxide Fuel Cell

This article aims to investigate the transient behavior of a planar direct internal reforming solid oxide fuel cell(DIR-SOFC) comprehensively.A one-dimensional dynamic model of a planar DIR-SOFC is first developed based on mass and energy balances,and electrochemical principles.Further,a solution strategy is presented to solve the model,and the International Energy Agency(IEA) benchmark test is used to validate the model.Then,through model-based simulations,the steady-state performance of a co-flow planar DIR-SOFC under specified initial operating conditions and its dynamic response to introduced operating parameter disturbances are studied.The dynamic responses of important SOFC variables,such as cell temperature,current density,and cell voltage are all investigated when the SOFC is subjected to the step-changes in various operating parameters including both the load current and the inlet fuel and air flow rates.The results indicate that the rapid dynamics of the current density and the cell voltage are mainly influenced by the gas composition,particularly the H2 molar fraction in anode gas channels,while their slow dynamics are both dominated by the SOLID(including the PEN and interconnects) temperature.As the load current increases,the SOLID temperature and the maximum SOLID temperature gradient both increase,and thereby,the cell breakdown is apt to occur because of excessive thermal stresses.Changing the inlet fuel flow rate might lead to the change in the anode gas composition and the consequent change in the current density distribution and cell voltage.The inlet air flow rate has a great impact on the cell temperature distribution along the cell,and thus,is a suitable manipulated variable to control the cell temperature.

Flower solid modeling based on sketches

In this paper we propose a method to model flowers of solid shape. Based on (Ijiri et al., 2005)’s method, we separate individual flower modeling and inflorescence modeling procedures into structure and geometry modeling. We incorporate interactive editing gestures to allow the user to edit structure parameters freely onto structure diagram. Furthermore, we use free-hand sketching techniques to allow users to create and edit 3D geometrical elements freely and easily. The final step is to automatically merge all independent 3D geometrical elements into a single waterproof mesh. Our experiments show that this solid modeling approach is promising. Using our approach, novice users can create vivid flower models easily and freely. The generated flower model is waterproof. It can have applications in visualization, animation, gaming, and toys and decorations if printed out on 3D rapid prototyping devices.

Application of Solid Thinking Evolve in Product Modeling Design

In view of application characteristics of Solid Thinking Evolve in product modeling design, and the practical design of cosmetics packaging, this paper introduced the unique Construction TreeTM, fast, convenient and easy features of Solid Thinking Evolve, and its application prospects in product development.

Experimental measurement and thermodynamic modeling of binary and ternary solid–liquid phase equilibrium for the systems formed by L-arabinose,D-xylose and water

Solid–liquid phase equilibrium data for binary(L-arabinose–water) and(D-xylose–water) systems at temperatures from(269.85–298.05) K and ternary(L-arabinose–D-xylose–water) system at temperatures of 273.85 K,278.85 K and 284.45 K were measured at atmospheric pressure.The ternary phase diagrams of the systems were constructed on the base of the measured solubility.Two pure solid phases were formed at given temperatures,including pure L-arabinose and pure D-xylose,which were con firmed and determined by the method of Schreinemakers' wet residue.At the same temperature,the crystallization region of L-arabinose was larger than D-xylose's.The acquired solubility data were then correlated using the NRTL model,Wilson model and Xu model.The calculated solubility with the three models agreed well with the experimental values.