Fatigue analysis of closed-cell aluminium foam using different material models

The fatigue analyses of AlSi7 closed-cell aluminium foam were performed using a real porous model and three different homogenised material models:crushable foam model,isotropic hardening model and kinematic hardening model.The numerical analysis using all three homogenised material models is based on the available experimental results previously determined from fatigue tests under oscillating tensile loading with the stress ratio R=0.1.The obtained computational results have shown that both isotropic and kinematic hardening models are suitable to analyse the fatigue behaviour of closed-cell aluminium foam.Besides,the kinematic hardening material model has demonstrated significantly shorter simulation time if compared to the isotropic hardening material model.On the other hand,the crushable foam model is recognized as an inappropriate approach for the fatigue analyses under tension loading conditions.

A powder-scale multiphysics framework for powder bed fusion of fiber-reinforced polymer composites

Additive manufacturing of fiber-reinforced polymer composites has garnered great interest due to its potential in fabricating functional products with lightweight characteristics and unique material properties.However,the major concern in polymer composites remains the presence of pore defects,as a thorough understanding of pore formation is insufficient.In this study,a powder-scale multiphysics framework has been developed to simulate the printing process of fiber-reinforced polymer composites in powder bed fusion additive manufacturing.This numerical framework involves various multiphysics phenomena such as particle flow dynamics of fiber-reinforced polymer composite powder,infrared laser–particle interaction,heat transfer,and multiphase fluid flow dynamics.The melt depths of one-layer glass fiber–reinforced polyamide 12 composite parts fabricated by selective laser sintering are measured to validate modelling predictions.The numerical framework is employed to conduct an in-depth investigation of pore formation mechanisms within printed composites.Our simulation results suggest that an increasing fiber weight fraction would lead to a lower densification rate,larger porosity,and lower pore sphericity in the composites.

Abrasive Waterjet Machining Simulation by Coupling Smoothed Particle Hydrodynamics /Finite Element Method

In dealing with abrasive waterjet machining（AWJM） simulation,most literatures apply finite element method（FEM） to build pure waterjet models or single abrasive particle erosion models.To overcome the mesh distortion caused by large deformation using FEM and to consider the effects of both water and abrasive,the smoothed particle hydrodynamics（SPH） coupled FEM modeling for AWJM simulation is presented,in which the abrasive waterjet is modeled by SPH particles and the target material is modeled by FEM.The two parts interact through contact algorithm.Utilizing this model,abrasive waterjet with high velocity penetrating the target materials is simulated and the mechanism of erosion is depicted.The relationships between the depth of penetration and jet parameters,including water pressure and traverse speed,etc,are analyzed based on the simulation.The simulation results agree well with the existed experimental data.The mixing multi-materials SPH particles,which contain abrasive and water,are adopted by means of the randomized algorithm and material model for the abrasive is presented.The study will not only provide a new powerful tool for the simulation of abrasive waterjet machining,but also be beneficial to understand its cutting mechanism and optimize the operating parameters.

A structural topological optimization method for multi-displacement constraints and any initial topology configuration

In density-based topological design, one expects that the final result consists of elements either black （solid material） or white （void）, without any grey areas. Moreover, one also expects that the optimal topology can be obtained by starting from any initial topology configuration. An improved structural topological optimization method for multidisplacement constraints is proposed in this paper. In the proposed method, the whole optimization process is divided into two optimization adjustment phases and a phase transferring step. Firstly, an optimization model is built to deal with the varied displacement limits, design space adjustments, and reasonable relations between the element stiffness matrix and mass and its element topology variable. Secondly, a procedure is proposed to solve the optimization problem formulated in the first optimization adjustment phase, by starting with a small design space and advancing to a larger deign space. The design space adjustments are automatic when the design domain needs expansions, in which the convergence of the proposed method will not be affected. The final topology obtained by the proposed procedure in the first optimization phase, can approach to the vicinity of the optimum topology. Then, a heuristic algorithm is given to improve the efficiency and make the designed structural topology black/white in both the phase transferring step and the second optimization adjustment phase. And the optimum topology can finally be obtained by the second phase optimization adjustments. Two examples are presented to show that the topologies obtained by the proposed method are of very good 0/1 design distribution property, and the computational efficiency is enhanced by reducing the element number of the design structural finite model during two optimization adjustment phases. And the examples also show that this method is robust and practicable.

Material Removal Model Considering Influence of Curvature Radius in Bonnet Polishing Convex Surface

The bonnet tool polishing is a novel, advanced and ultra-precise polishing process, by which the freeform surface can be polished. However, during the past few years, not only the key technology of calculating the dwell time and controlling the surface form in the bonnet polishing has been little reported so far, but also little attention has been paid to research the material removal function of the convex surface based on the geometry model considering the influence of the curvature radius. Firstly in this paper, for realizing the control of the freeform surface automatically by the bonnet polishing, on the basis of the simplified geometric model of convex surface, the calculation expression of the polishing contact spot on the convex surface considering the influence of the curvature radius is deduced, and the calculation model of the pressure distribution considering the influence of the curvature radius on the convex surface is derived by the coordinate transformation. Then the velocity distribution model is built in the bonnet polishing the convex surface. On the basis of the above research and the semi-experimental modified Preston equation obtained from the combination method of experimental and theoretical derivation, the material removal model of the convex surface considering the influence of the curvature radius in the bonnet polishing is established. Finally, the validity of the model through the simulation method has been validated. This research presents an effective prediction model and the calculation method of material removal for convex surface in bonnet polishing and prepares for the bonnet polishing the free surface numerically and automatically.

Numerical Simulation of Warm Forming Behavior of High Strength Aluminum Alloy 7075

Numerical analysis is critically important to understanding the complex deformation mechanics that occur during sheet forming processes.It has been widely used in simulation of sheet metal forming processes at room temperature in the automotive industry.However,material at elevated temperature behaves more differently than at room temperature and specific material parameters and models need to be developed for the simulation of warm forming.Based on the experimental investigation of material behavior of high strength aluminum alloy 7075(AA7075),constitutive equations with strain rate sensitivity at 140,180 and 220 ℃ are developed.Anisotropic yield criterion Barlat 89 is used in the simulation.Warm forming of limit dome height tests and limit drawing ratio tests of AA7075 at 140,180 and 220℃are performed.Forming limit diagrams developed from experiment at several elevated temperatures in the previous study are used to predict the failure in the simulation results.Punch force and displacement predicted from simulation are compared with the experimental data.Simulation results agree with experimental results,so the developed material model can be used to accurately predict material behavior during isothermal warm forming of the AA7075-T6 alloy.

Investigation of flue gas water-alternating gas (flue gas–WAG) injection for enhanced oil recovery and multicomponent flue gas storage in the post-waterflooding reservoir

Flue gas fooding is one of the important technologies to improve oil recovery and achieve greenhouse gas storage.In order to study multicomponent fue gas storage capacity and enhanced oil recovery(EOR)performance of fue gas water-alternating gas(fue gas-WAG)injection after continuous waterfooding in an oil reservoir,a long core fooding system was built.The experimental results showed that the oil recovery factor of fue gas-WAG fooding was increased by 21.25%after continuous waterfooding and fue gas-WAG fooding could further enhance oil recovery and reduce water cut signifcantly.A novel material balance model based on storage mechanism was developed to estimate the multicomponent fue gas storage capacity and storage capacity of each component of fue gas in reservoir oil,water and as free gas in the post-waterfooding reservoir.The ultimate storage ratio of fue gas is 16%in the fue gas-WAG fooding process.The calculation results of fue gas storage capacity showed that the injection gas storage capacity mainly consists of N_(2) and CO_(2),only N_(2) exists as free gas phase in cores,and other components of injection gas are dissolved in oil and water.Finally,injection strategies from three perspectives for fue gas storage,EOR,and combination of fue gas storage and EOR were proposed,respectively.

Mesoscopic Modeling Approach and Application for Steel Fiber Reinforced Concrete under Dynamic Loading:A Review

Steel fiber reinforced concrete(SFRC)has drawn extensive attention in recent years for its superior mechanical response to dynamic and impact loadings.Based on the existing test results,the highstrength steel fibers embedded in a concrete matrix usually play a strong bridging effect to enhance the bonding force between fiber and the matrix,and directly contribute to the improvement of the post-cracking behavior and residual strength of SFRC.To gain a better understanding of the action behavior of steel fibers in matrix and further capture the failure mechanism of SFRC under dynamic loads,the mesoscopic modeling approach that assumes SFRC to be composed of different mesoscale phases(i.e.,steel fibers,coarse aggregates,mortar matrix,and interfacial transition zone(ITZ))has been widely employed to simulate the dynamic responses of SFRC material and structural members.This paper presents a comprehensive review of the state-of-the-art mesoscopic models and simulations for SFRC under dynamic loading.Generation approaches for the SFRC mesoscale model in the simulation works,including steel fiber,coarse aggregate,and the ITZ between them,are reviewed and compared systematically.The material models for different phases and the interaction relationship between fiber and concrete matrix are summarized comprehensively.Additionally,some example applications for SFRC under dynamic loads(i.e.,compression,tension,and contact blast)simulated using the general mesoscale models are given.Finally,some critical analysis on the current shortcomings of the mesoscale modeling of SFRC is highlighted,which is of great significance for the future investigation and development of SFRC.

Deformation Stability of GH4033 Superalloy in the Hot Continuous Rolling Process Based on Dynamic Material Model and Finite Element Model

The flow stress behavior of GH4033 superalloy was determined by the hot compression tests at the temperatures of 1223-1473 K and the total strains of 0.6 with the strain rates of 0.001-30.0 s^(-1) by using cylindrical samples.The processing maps based on the dynamic material model(DMM)combined with the corresponding microstructure observations indicate the reasonable processing domain locating at the strain rates of 0.1-1.0 s^(-1) and the deformation temperature of 1273-1423 K.Meanwhile,the numerical simulation based on finite element model(FEM)described the variation of the effective strain,effective strain rate and the temperature for the core node,and unveiled the influence of the hot rolling parameters considering the initial temperature(T_(0))range of 1223-1473 K and the first-stand biting velocity(v_(0))range of 0.15-0.35 m·s^(-1).Furthermore,the deformation stability of GH4033 superalloy in the round rod hot continuous rolling(HCR)process is described and analyzed by coupling the three-dimensional(3-D)processing map,and the spatial trajectory lines were determined by the numerically simulated temperatures,the strains and the strain rates.Finally,the results show that the hot deformation stability of GH4033 can be achieved by the rolling process parameters located at T_(0)=1423 K and v_(0)=0.25 m·s^(-1).Additionally,the practical HCR processes as T_(0)=1423 K and v_(0)=0.15,0.25,0.35 m·s^(-1) were operated to verify the influence of the hot rolling parameters on the hot deformation stability by the microstructure observation of the final products.

Correlation between Single Factor Material Constitutive Model Parameters and Temperature for Ti6Al4V Alloy

Accurate material constitutive model is considered highly necessary to perform finite element simulation and analysis.However,it is difficult to establish the material constitutive model because of uncertainty of mathematical relationship and constraint of existing experimental condition.At present,there exists considerable gap between finite element simulation result and actual cutting process.Particular emphases were put on investigating the correlation between ＂single factor＂ material constitutive model parameters and temperature for Ti6Al4V alloy,and also establishment of material constitutive model for this kind of material.Theoretical analyses based on dislocation theory and material functional relations showed that material model was deeply affected by variation temperature.By the least squares best fit to the available quasi-static and high-speed impact compression experiment data,material parameters at various temperatures were found.Experimental curves analyses and material parameters comparison showed that the ＂single factor＂ material constitutive model parameters were temperature dependent.Using the mathematical mapping between material parameters and temperature,＂single factor＂ material constitutive model of Ti6Al4V alloy was established,which was proven to be right by comparing with experimental measurements.This work makes clear that the ＂single factor＂ material constitutive model parameters of Ti6Al4V alloy are temperature dependent.At the same time,an accurate material constitutive model is established,which helps to optimize cutting process and control machining distortion for Ti6Al4V alloy aerospace parts.

Elasto-plastic analysis of masonry with anisotropic plastic material model

This paper establishes an anisotropic plastic material model to analyze the elasto-plastic behavior of masonry in plane stress state.Being an anisotropic material,masonry has different constitutive relation and fracture energies along each orthotropic axes.Considering the unique material properties of masonry,a new yield criterion for masonry is proposed combining the Hill's yield criterion and the Rankine's yield criterion.The new yield criterion not only introduces compression friction coefficient of shear but also considers yield functions for independent stress state along two material axes of tension.To solve the involved nonlinear equations in numerical analysis,several nonlinear methods are implemented,including Newton-Raphson method for nonlinear equations and Implicit Euler backward mapping algorithm to update stresses.To verify the proposed material model of masonry,a series of tests are operated.The simulation results show that the new developed material model implements successfully.Compared with isotropic material model,the proposed model performs better in elasto-plastic analysis of masonry in plane stress state.The proposed anisotropic model is capable of simulating elasto-plastic behavior of masonry and can be used in related applications.

Multi-Material Topology Optimization of Structures Using an Ordered Ersatz Material Model

This paper proposes a new element-based multi-material topology optimization algorithm using a single variable for minimizing compliance subject to a mass constraint.A single variable based on the normalized elemental density is used to overcome the occurrence of meaningless design variables and save computational cost.Different from the traditional material penalization scheme,the algorithm is established on the ordered ersatz material model,which linearly interpolates Young’s modulus for relaxed design variables.To achieve a multi-material design,the multiple floating projection constraints are adopted to gradually push elemental design variables to multiple discrete values.For the convergent element-based solution,the multiple level-set functions are constructed to tentatively extract the smooth interface between two adjacent materials.Some 2D and 3D numerical examples are presented to demonstrate the effectiveness of the proposed algorithm and the possible advantage of the multi-material designs over the traditional solid-void designs.

Finite Element Simulation of Flexible Roll Forming with Supplemented Material Data and the Experimental Verification

Flexible roll forming is a promising manufacturing method for the production of variable cross section products. Considering the large plastic strain in this forming process which is much larger than that of uniform deformation phase of uniaxial tensile test, the widely adopted method of simulating the forming processes with non-supplemented material data from uniaxial tensile test will certainly lead to large error. To reduce this error, the material data is supplemented based on three constitutive models. Then a finite element model of a six passes flexible roll forming process is established based on the supplemented material data and the original material data from the uniaxial tensile test. The flexible roll forming experiment of a B pillar reinforcing plate is carried out to verify the proposed method. Final cross section shapes of the experimental and the simulated results are compared. It is shown that the simulation calculated with supplemented material data based on Swift model agrees well with the experimental results, while the simulation based on original material data could not predict the actual deformation accurately. The results indicate that this material supplement method is reliable and indispensible, and the simulation model can well reflect the real metal forming process. Detailed analysis of the distribution and history of plastic strain at different positions are performed. A new material data supplement method is proposed to tackle the problem which is ignored in other roll forming simulations, and thus the forming process simulation accuracy can be greatly improved.

Development of Mathematical Model on Preparation of Functionally Graded Material by Co-sedimentation

From the process of sedimentation the mathematical relationships among deposition Volume and powder properties as well as sedimentation parameters were deduced. Based on the formula a mathematical model was set up and simulated through the computer. At last the validity of mathematical model was supported by the representative experiment on Ti-Mo system FGM prepared by co-sedimentation.

Experimental study on the movement law of overlying rock non-pillar coal overhead mining

Pillarless coal mining technology is a new practical technology.Based on the compensating mechanical behavior of the Negative Poisson’s Ratio(NPR)anchor cable on the roof,the roadway was successfully retained by the top cutting and pressure relief technology.This study utilizes the Digital Speckle Monitoring(DIC monitoring),stress-strain monitoring,and infrared thermal imaging systems to conduct physical model experiment of similar materials from the displacement,stress-strain,and temperature fields to investigate in depth the fracture change law of the overlying rock.In addition,it uses FLAC3D numerical simulation to invert the surface displacement settlement.The results show that the non-pillar overhead mining under the 110 mining method has little influence on the rock crack in the middle of the coal seam,and the crack development area is mainly concentrated in the overlying rock mass of the upward coal seam.The compensatory mechanical behavior of NPR anchor cable and the dilatation characteristics of rock mass have a good effect of retaining roadway along goaf,and can also reduce surface settlement.The 110 mining method provides a scientific basis for ecological environment protection and the development of other kilometer deep soft rock high ground stress underground projects.

Hybrid graded element model for transient heat conduction in functionally graded materials

This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials （FGMs）. First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest＇s algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.

MODELING OF LASER MACHINING ON POLYMETHYL METHACRYLATE TO FABRICATE MICROFLUIDIC CHIP

The use of a CO2 laser system for fabrication of microfluidic chip on polymethyl methacrylate （PMMA） is presented to reduce fabrication cost and time of chip. The grooving process of the laser system and a model for the depth of microchannels are investigated. The relations between the depth of laser-cut channels and the laser beam power, velocity or the number of passes of the beam along the same channel are evaluated. In the experiments, the laser beam power varies from 0 to 50 W, the laser beam scanning velocity varies from 0 to 1 000 mm/s and the passes vary in the range of 1 to 10 times. Based on the principle of conservation of energy, the influence of the laser beam velocity, the laser power and the number of groove passes are examine. Considering the grooving interval energy loss, a modified mathematical model has been obtained and experimental data show good agreement with the theoretical model. This approach provides a simple way of predicting groove depths. The system provides a cost alternative of the other methods and it is especially useful on research work of rnicrofluidic prototyping due to the short cycle time of production.

Selective ensemble modeling based on nonlinear frequency spectral feature extraction for predicting load parameter in ball mills

Strong mechanical vibration and acoustical signals of grinding process contain useful information related to load parameters in ball mills. It is a challenge to extract latent features and construct soft sensor model with high dimensional frequency spectra of these signals. This paper aims to develop a selective ensemble modeling approach based on nonlinear latent frequency spectral feature extraction for accurate measurement of material to ball volume ratio. Latent features are first extracted from different vibrations and acoustic spectral segments by kernel partial least squares. Algorithms of bootstrap and least squares support vector machines are employed to produce candidate sub-models using these latent features as inputs. Ensemble sub-models are selected based on genetic algorithm optimization toolbox. Partial least squares regression is used to combine these sub-models to eliminate collinearity among their prediction outputs. Results indicate that the proposed modeling approach has better prediction performance than previous ones.

Thermal-stress induced phenomena in two-component material:part I

The paper deals with analytical fracture mechanics to consider elastic thermal stresses acting in an isotropic multi-particle-matrix system. The multi-particle-matrix system consists of periodically distributed spherical particles in an infinite matrix. The thermal stresses originate during a cooling process as a consequence of the difference αm - αp in thermal expansion coefficients between the matrix and the particle, αm and αp, respectively. The multi-particle-matrix system thus represents a model system applicable to a real two-component material of a precipitation-matrix type. The infinite matrix is imaginarily divided into identical cubic cells. Each of the cubic cells with the dimension d contains a central spherical particle with the radius R, where d thus corresponds to inter-particle distance. The parameters R, d along with the particle volume fraction v = v（R, d） as a function of R, d represent microstructural characteristics of a twocomponent material. The thermal stresses are investigated within the cubic cell, and accordingly are functions of the microstructural characteristics. The analytical fracture mechanics includes an analytical analysis of the crack initiation and consequently the crack propagation both considered for the spherical particle （q = p） and the cell matrix （q = m）. The analytical analysis is based on the determination of the curve integral Wcq of the thermal-stress induced elastic energy density Wq. The crack initiation is represented by the determination of the critical particle radius Rqc = Rqc（V）. Formulae for Rqc are valid for any two-component mate- rial of a precipitate-matrix type. The crack propagation for R 〉 Rqc is represented by the determination of the function fq describing a shape of the crack in a plane perpendicular

An Elastic-Plastic Iceberg Material Model Considering Temperature Gradient Effects and its Application to Numerical Study

To simulate the FPSO-iceberg collision process more accurately, an elastic-plastic iceberg material model considering temperature gradient effects is proposed and applied. The model behaves linearly elastic until it reaches the ‘Tsai-Wu’ yield surfaces, which are a series of concentric elliptical curves of different sizes. Decreasing temperature results in a large yield surface. Failure criteria, based on the influence of accumulated plastic strain and hydrostatic pressure, are built into the model. Based on published experimental data on the relationship between depth and temperature in icebergs, three typical iceberg temperature profiles are proposed. According to these, ice elements located at different depths have different temperatures. The model is incorporated into LS-DYNA using a user-defined subroutine and applied to a simulation of FPSO collisions with different types of icebergs. Simulated area-pressure curves are compared with design codes to validate the iceberg model. The influence of iceberg shape and temperature on the collision process is analyzed. It is indicated that FPSO structural damage not only depends on the relative strength between the iceberg and the structure, but also depends on the local shape of the iceberg.