Characterization for elastic constants of fused deposition modelling-fabricated materials based on the virtual fields method and digital image correlation

Fused deposition modelling(FDM), a widely used rapid prototyping process, is a promising technique in manufacturing engineering. In this work, a method for characterizing elastic constants of FDM-fabricated materials is proposed. First of all, according to the manufacturing process of FDM, orthotropic constitutive model is used to describe the mechanical behavior. Then the virtual fields method(VFM) is applied to characterize all the mechanical parameters(Q, Q, Q, Q) using the full-field strain,which is measured by digital image correlation(DIC). Since the principal axis of the FDM-fabricated structure is sometimes unknown due to the complexity of the manufacturing process, a disk in diametrical compression is used as the load configuration so that the loading angle can be changed conveniently. To verify the feasibility of the proposed method, finite element method(FEM) simulation is conducted to obtain the strain field of the disk. The simulation results show that higher accuracy can be achieved when the loading angle is close to 30?. Finally, a disk fabricated by FDM was used for the experiment. By rotating the disk, several tests with different loading angles were conducted. To determine the position of the principal axis in each test, two groups of parameters(Q, Q, Q, Q) are calculated by two different groups of virtual fields. Then the corresponding loading angle can be determined by minimizing the deviation between two groups of the parameters. After that, the four constants(Q, Q, Q, Q) were determined from the test with an angle of 27?.

Volumetric lattice Boltzmann method for pore-scale mass diffusionadvection process in geopolymer porous structures

Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.

Characterization of material mechanical properties using strain correlation method combined with virtual fields method

The isotropic and anisotropic behaviors are considered as the important formats of the constitutive behaviors,and can also be called the global properties.To improve the identification ability of virtual fields method(VFM)when the global properties are unknown,this paper proposes the strain correlation method(SCM)to determine the global properties before the parameter identification using the VFM.Firstly,the basic principle of SCM is described in detail.Then,the feasibility and accuracy of SCM are verified through the numerical experiments based on the three-point bending configuration and the real experiment of polymethyl methacrylate(PMMA).The influence of the additive Gaussian white noise,local errors in the strain fields,and missing data at the specimen edges on the characterization results are evaluated.The results show that the SCM has good noise immunity and lower accuracy requirements for the strain fields.As an application,the mechanical properties of Ti-6A1-4V alloys fabricated by selective laser melting(SLM)are characterized by the SCM.The results show that the alloys are isotropic,and the isotropic VFM is utilized to determine the mechanical parameters.By using the SCM,the accuracy of identification results can be improved for the isotropic or bidirectional reinforced orthotropic materials when using VFM.

Design and analysis of an advanced thermal management system for the solar close observations and proximity experiments spacecraft

In this paper,the mission and the thermal environment of the Solar Close Observations and Proximity Experiments(SCOPE)spacecraft are analyzed,and an advanced thermal management system(ATMS)is designed for it.The relationship and functions of the integrated database,the intelligent thermal control system and the efficient liquid cooling system in the ATMS are elaborated upon.For the complex thermal field regulation system and extreme space thermal environment,a modular simulation and thermal field planning method are proposed,and the feasibility of the planning algorithm is verified by numerical simulation.A solar array liquid cooling system is developed,and the system simulation results indicate that the temperatures of the solar arrays meet the requirements as the spacecraft flies by perihelion and aphelion.The advanced thermal management study supports the development of the SCOPE program and provides a reference for the thermal management in other deep-space exploration programs.

Rolling Force and Rolling Moment in Spline Cold Rolling Using Slip-line Field Method

Rolling force and rolling moment are prime process parameter of external spline cold rolling. However, the precise theoretical formulae of rolling force and rolling moment are still very fewer, and the determination of them depends on experience. In the present study, the mathematical models of rolling force and rolling moment are established based on stress field theory of slip-line. And the isotropic hardening is used to improve the yield criterion. Based on MATLAB program language environment, calculation program is developed according to mathematical models established. The rolling force and rolling moment could be predicted quickly via the calculation program, and then the reliability of the models is validated by FEM. Within the range of module of spline m=0.5-1.5 mm, pressure angle of reference circle α=30.0°-45.0°, and number of spline teeth Z=19-54, the rolling force and rolling moment in rolling process （finishing rolling is excluded） are researched by means of virtualizing orthogonal experiment design. The results of the present study indicate that： the influences of module and number of spline teeth on the maximum rolling force and rolling moment in the process are remarkable; in the case of pressure angle of reference circle is little, module of spline is great, and number of spline teeth is little, the peak value of rolling force in rolling process may appear in the midst of the process; the peak value of rolling moment in rolling process appears in the midst of the process, and then oscillator weaken to a stable value. The results of the present study may provide guidelines for the determination of power of the motor and the design of hydraulic system of special machine, and provide basis for the farther researches on the precise forming process of external spline cold rolling.

On the field method in non-holonomic mechanics

This paper deals with the generalization of the fieldmethod to non-holonomic systems whose motion is subject toeither non-linear constraints or those of a higher order,whiletheir motion is modeled by the generalized Lagrange equa-tions of the second kind.Two examples are given to illustratethe theory.

Anisotropic growth of multigrain in equiaxial solidification simulated with the phase field method

The phase field method has been mainly used to simulate the growth of a single crystal in the past. But polycrystalline materials predominate in engineering. In this work, a phase field model for multigrain solidification is developed, which takes into account the random crystallographic orientations of crystallites and preserves the rotational invariance of the free energy. The morphological evolution of equiaxial multigrain solidification is predicted and the effect of composition on transformation kinetics is studied. The numerical results indicate that due to the soft impingement of grains the Avrami exponent varies with the initial melt composition and the solidification fraction.

Average vector field methods for the coupled Schrdinger KdV equations

The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.

Analysis on high-temperature oxidation and growth stress of iron-based alloy using phase field method

High-temperature oxidation is an important property to evaluate thermal protection materials. However, since oxidation is a complex process involving microstructure evolution, its quantitative analysis has always been a challenge. In this work, a phase field method （PFM） based on the thermodynamics theory is developed to simulate the oxidation behavior and oxidation induced growth stress. It involves microstructure evolution and solves the problem of quantitatively computational analysis for the oxidation behavior and growth stress. Employing this method, the diffusion process, oxidation performance, and stress evolution axe predicted for Fe-Cr-A1-Y alloys. The numerical results agree well with the experimental data. The linear relationship between the maximum growth stress and the environment oxygen concentration is found. PFM provides a powerful tool to investigate high-temperature oxidation in complex environments.

Interval finite difference method for steady-state temperature field prediction with interval parameters

A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.

A FIELD METHOD FOR INTEGRATING THE EQUATIONS OF MOTION OF NONHOLONOMIC CONTROLLABLE SYSTEMS

This paper presents a field method for integrating the equations of motion of nonholonomic controllable systems. An example is given to illustrate the application of the method.

Effects of Method for Returning Straw to Field on Soil Properties,Straw Decomposition and Nutrient Release

[Objectives]To alleviate the influence of meteorological conditions on soil environment(temperature and water content)and maintain high and stable grain yield.[Methods]Taking Sunzhen Experimental Station of Weinan Academy of Agricultural Sciences as the experimental base,the effects of returning double-crop wheat and corn straw to field(Twm),returning single-crop corn straw to field(Tm),returning single-crop wheat straw to field(Tw)on soil temperature,water content,straw decomposition rate and nutrient release,soil organic matter and bulk density were studied systematically.[Results]Twm treatment could effectively alleviate the effects of meteorological conditions on soil temperature and water content.The decomposition rate of straw treated with Twm was 4.7%higher than that of Tm treatment,3.8%higher than that of Tw treatment,10.5%higher than that of Tm treatment,and the decomposition rate of straw showed a trend of"first fast,then slow and then fast".The release of nitrogen from straw was basically similar to that of straw decay,and the release of potassium and phosphorus increased at first and then remained basically unchanged.The release rate of potassium was the highest,followed by phosphorus and nitrogen.The content of soil organic matter in Twm treatment increased by 11.67%annually,an annual average of 0.998 g/kg.The soil bulk density of Twm treatment decreased by 0.058 g/cm^(3) annually,an annual average of 4.29%.The fundamental reason is that Twm treatment provides conditions(temperature,water content,nutrition)for microbial growth,reproduction,enzyme production and biochemical reaction,and increases the exchange capacity of soil and external water,heat,gas and fertilizer.[Conclusions]It is expected is to help people change their understanding of returning straw to field from"quick harvest"to"fertilizer transformation".

Novel energy dissipative method on the adaptive spatial discretization for the Allen–Cahn equation

We propose a novel energy dissipative method for the Allen–Cahn equation on nonuniform grids.For spatial discretization,the classical central difference method is utilized,while the average vector field method is applied for time discretization.Compared with the average vector field method on the uniform mesh,the proposed method can involve fewer grid points and achieve better numerical performance over long time simulation.This is due to the moving mesh method,which can concentrate the grid points more densely where the solution changes drastically.Numerical experiments are provided to illustrate the advantages of the proposed concrete adaptive energy dissipative scheme under large time and space steps over a long time.

A field method for integrating the equations of motion of mechanico-electrical coupling dynamical systems

A field method for integrating the equations of motion for mechanico-electrical coupling dynamical systems is studied. Two examples in mechanico-electrical engineering are given to illustrate this method.

A differential wide field electromagnetic method and its application in alkaline-surfactant-polymer(ASP) flooding monitoring

To improve effectiveness of ASP flooding, it is necessary to establish a reliable parameter design and tracking adjustment method to monitor the process of oil displacement. A differential wide field electromagnetic method was proposed and applied to the ASP displacement monitoring test in a block of the Daqing Oilfield. In the process of ASP flooding, the electromagnetic field was measured many times. The data acquired before the ASP flooding were set as the background field, and the resistivity model was obtained by inversion. Then, the resistivity data were calibrated by logging data and the resistivity model was established. Finally, the range and front of ASP flooding were deduced with the residual gradient from the spatial domain first-order difference of the resistivity model. Production data of well groups in this block have proved that this method can work out the range and front of ASP flooding accurately, providing support for optimization of ASP flooding parameters.

New Method for 3D Transient Eddy Current Field Calculation and Its Application in Magneto-Acoustic Tomography

A new method of 3D transient eddy current field calculation is proposed. The Maxwell equations with time component elimination （METCE） are derived under the assumption of magnetic quasi static approximation, especially for the sample of low conductivity. Based on METCE, we deduce a more efficient reconstruction algorithm of a 3D transient eddy current field. The computational burden is greatly reduced through the new algorithm, and the computational efficiency is improved. This new algorithm decompounds the space-time variables into two individual variables. The idea is to solve the spatial vector component firstly, and then multiply it by the corresponded time component. The iterative methods based on METCE are introduced to recover the distribution of conductivity in magneto-acoustic tomography. The reconstructed images of conductivity are consistent with the original distribution, which validate the new method.

Simulation of Temperature and Flow Field by Three Dimension Finite ElementMethod for Castex Process of AS Wire

Castex of AS wire is a new technology of near net shape. To study the variation of temperature and velocity of liquid (or semisolid) aluminum during dynamic solidification the numerical simulation was carried out with the theory of heat-transfer and hydrodynamics by means of 3-dimensional finite element method. From simulation results, it is found that the variation of temperature and velocityis mainly influenced by the casting temperature of aluminum, rotating speed of Castex wheel and flow of cooling water. Among theseinfluencing factors, the casting temperature distributes most to the length of liquid phase metal. Moreover, the faster the metal solidifies,the higher the metal there moves with the overall trend of descending from the bottom of the wheel to the shoe wall as well as from sidewalls to the center of wheel groove. In comparison with the practical value, the simulation is reliable.

Thermo-hydro-mechanical-air coupling finite element method and its application to multi-phase problems

In this paper, a finite element method （FEM）-based multi-phase problem based on a newly proposed thermal elastoplastic constitutive model for saturated/unsaturated geomaterial is discussed. A program of FEM named as SOFT, adopting unified field equations for thermo-hydro-mechanical-air （THMA） behavior of geomaterial and using finite element-finite difference （FE-FD） scheme for so/l-water-air three-phase coupling problem, is used in the numerical simulation. As an application of the newly proposed numerical method, two engineering problems, one for slope failure in unsaturated model ground and another for in situ heating test related to deep geological repository of high-level radioactive waste （HLRW）, are simulated. The model tests on slope failure in unsaturated Shirasu ground, carried out by Kitamura et al. （2007）, is simulated in the framework of soil-water-air three-phase coupling under the condition of constant temperature. While the in situ heating test reported by Munoz （2006） is simulated in the same framework under the conditions of variable temperature hut constant air pressure.

Finite Element Method Based on Equivalent Magnetic Energy Method for Computation of 2D Nonlinear Eddy Current Field

In this paper, the finite element method using vector potential in applications to 2D nonlinear eddy current field is discussed. The authors use the equivalent magnetic energy method to deal with magnetization curve of ferromagnetic material,and present the formulation of 2D nonlinear eddy current field.With this method the authors analyze the eddy current field in an induction ladle furnace and the force distribution in the charge (molten metal),and plot the corresponding curves.