Computation of elastic properties of 3D digital cores from the Longmaxi shale

The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Longmaxi black shale samples using X-ray CT, we built detailed 3D digital images of cores with porosity properties and mineral contents. Next, we used finite-element （FE） methods to derive the elastic properties of the samples. The FE method can accurately model the shale mineralogy. Particular attention is paid to the derived elastic properties and their dependence on porosity and kerogen. The elastic moduli generally decrease with increasing porosity and kerogen, and there is a critical porosity （0.75） and kerogen content （ca. ≤3%） over which the elastic moduli decrease rapidly and slowly, respectively. The derived elastic moduli of gas- and oil-saturated digital cores differ little probably because of the low porosity （4.5%） of the Longmaxi black shale. Clearly, the numerical experiments demonstrated the feasibility of combining microstructure images of shale samples with elastic moduli calculations to predict shale properties.

Exceptional Performance of 3D Additive Manufactured NiFe Phosphite Oxyhydroxide Hollow Tubular Lattice Plastic Electrode for Large-Current-Density Water Oxidization

In this article,we report a 3D NiFe phosphite oxyhydroxide plastic electrode using high-resolution digital light processing(DLP)3D-printing technology via induced chemical deposition method.The as-prepared 3D plastic electrode exhibits no template requirement,freedom design,low-cost,robust,anticorrosion,lightweight,and micro-nano porous characteristics.It can be drawn to the conclusion that highly oriented open-porous 3D geometry structure will be beneficial for improving surface catalytic active area,wetting performance,and reaction–diffusion dynamics of plastic electrodes for oxygen evolution reaction(OER)catalysis process.Density functional theory(DFT)calculation interprets the origin of high activity of NiFe(PO_(3))O(OH)and demonstrates that the implantation of the–PO_(3)can effectively bind the 3d orbital of Ni in NiFe(PO_(3))O(OH),lead to the weak adsorption of intermediate,make electron more active to improve the conductivity,thereby lowing the transform free energy of*O to*OOH.The water oxidization performance of as-prepared 3D NiFe(PO_(3))O(OH)hollow tubular(HT)lattice plastic electrode has almost reached the state-of-the-art level compared with the as-reported large-current-density catalysts or 3D additive manufactured plastic/metal-based electrodes,especially for high current OER electrodes.This work breaks through the bottleneck that plagues the performance improvement of low-cost high-current electrodes.

ELASTIC BEHAVIOR ANALYSIS OF 3D ANGLE-INTERLOCK WOVEN CERAMIC COMPOSITES

A micromechanical model for elastic behavior analysis of angle-interlock woven ceramic composites is proposed in this paper. This model takes into account the actual fabric structure by considering the fiber undulation and continuity in space, the cavities between adjacent yarns and the actual cross-section geometry of the yarn. Based on the laminate theory, the elastic properties of 3D angle-interlock woven ceramic composites are predicted. Different numbers of interlaced wefts have almost the same elastic moduli. The thickness of ceramic matrix has little effect on elastic moduli. When the undulation ratio increases longitudinal modulus decreases and the other Young＇s moduli increase. Good agreement between theoretical predictions and experimental results demonstrates the feasibility of the proposed model in analyzing the elastic properties of 3D angle-interlock woven ceramic composites. The results of this paper verify the fact that the method of analyzing polyester matrix composites is suitable for woven ceramic composites.

3D elastic waveform modeling with an optimized equivalent staggered-grid finite-difference method

Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.

Stability of High-Order Staggered-Grid Schemes for 3D Elastic Wave Equation in Heterogeneous Media

In this paper, we firstly derive the stability conditions of high-order staggered-grid schemes for the three-dimensional (3D) elastic wave equation in heterogeneous media based on the energy method. Moreover, the plane wave analysis yields a sufficient and necessary stability condition by the von Neumann criterion in homogeneous case. Numerical computations for 3D wave simulation with point source excitation are given.

3D bioprinting of complex biological structures with tunable elastic modulus and porosity using freeform reversible embedding of suspended hydrogels

Three-dimensional(3D)bioprinting has been used widely for the construction of hard tissues such as bone and cartilage.However,constructing soft tissues with complex structures remains a challenge.In this study,complex structures characterized by both tunable elastic modulus and porosity were printed using freeform reversible embedding of suspended hydrogels(FRESHs)printing methods.A mixture of alginate and gelatin was used as the main functional component of the bioink.Rheological analysis showed that this bioink possesses shear thinning and shear recovery properties,supporting both cryogenic and FRESH printing methods.Potential printing capabilities and limitations of cryogenic and FRESH printing were then analyzed by printability tests.A series of complex structures were printed by FRESH printing methods which could not be realized using conventional approaches.Mechanical tests and scanning electron microscopy analysis showed that the printed structure is of excellent flexibility and could be applied in various conditions by adjusting its mechanical modulus and porosity.L929 fibroblast cells maintained cell viability in cell-laden-printed structures,and the addition of collagen further improved the hydrogels’biocompatibility.Overall,all results provided useful insight into the building of human soft tissue organ blocks.

Tunable three-dimensional nonreciprocal transmission in a layered nonlinear elastic wave metamaterial by initial stresses

In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and second harmonic with the quasilongitudinal(qP)and quasi-shear(qS_(1) and qS_(2))modes are derived.Based on the transfer and stiffness matrices,band gaps with initial stresses are obtained by the Bloch theorem.The transmission coefficients are calculated to support the band gap property,and the tunability of the nonreciprocal transmission by the initial stress is discussed.This work is expected to provide a way to tune the nonreciprocal transmission with vector characteristics.

Hybrid absorbing boundary condition for threedimensional elastic wave modeling

Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.

3D Printing of Polylactic Acid Bioplastic–Carbon Fibres and Twisted Kevlar Composites Through Coextrusion Using Fused Deposition Modeling

Polylactic acid(PLA)bioplastic is a common material used in Fused Deposition Modeling(FDM)3D printing.It is biodegradable and environmentally friendly biopolymer which made out of corn.However,it exhibits weak mechanical properties which reduced its usability as a functional prototype in a real-world application.In the present study,two PLA composites are created through coextruded with 3K carbon fibres and twisted Kevlar string(as core fibre)to form a fibre reinforced parts(FRP).The mechanical strength of printed parts was examined using ASTM D638 standard with a strain rate of 1 mm/min.It has been demonstrated that the FRPs coextruded with 3K carbon fibres had achieved significant improvement in Young’s modulus(+180.6%,9.205 GPa),ultimate tensile strength(+175.3%,103 MPa)and maximum tensile strain(+21.6%,1.833%).Although the Young’s modulus of Kevlar FRP was found to be similar to as compared to unreinforced PLA(~3.29 GPa),it has gained significant increment in terms of maximum tensile strain(+179.7%,104.64 MPa),and maximum tensile strain(+257%,5.384%).Thus,this study revealed two unique composite materials,in which the 3K carbon FRP can offer stiff and high strength structure while Kevlar FRP offers similar strength but at a higher elasticity.

ANALYSIS OF SLAB EDGING BY A 3-D RIGID VISCO-PLASTIC FINITE ELEMENT METHOD

3-D rigid visco-plastic finite element method (FEM) is used in the analysisof metal forming processes, including strip and plate rolling, shape rolling, slab edging, specialstrip rolling. The shifted incomplete Cholesky decomposition of the stiffness matrix with thesolution of the equations for velocity increment by the conjugate gradient method is combined. Thistechnique, termed the shifted ICCG method, is then employed to solve the slab edging problem. Theperformance of this algorithm in terms of the number of iterations, friction variation, shiftedparameter psi and the results of simulation for processing parameters are analysed. Numerical testsand application of this technique verify the efficiency and stability of the shifted ICCG method inthe analysis of slab edging.

Elastic Modulus Prediction of Three-dimension-4 Directional Braided C_f/SiC Composite Based on Double-scale Model

Double-scale model for three-dimension-4 directional（3D-4d） braided C/SiC composites has been proposed to investigate its elastic properties. The double-scale model involves micro-scale that takes fiber/ matrix/porosity in fibers tows into consideration with unit cell which considers the 3D-4d braiding structure. Micro-optical photographs of composites have been taken to study the braided structure. Then a parameterized finite element model that reflects the structure of 3D-4d braided composites is proposed. Double-scale elastic modulus prediction model is developed to predict the elastic properties of 3D-4d braided C/SiC composites. Stiffness and eompliance-averaging method and energy method are adopted to predict the elastic properties of composites. Static-tension experiments have been conducted to investigate the elastic modulus of 3D-4d braided C/SiC composites. Finally, the effect of micro-porosity in fibers tows on the elastic modulus of 3D-4d braided C/SiC composites has been studied. According to the conclusion of this thesis, elastic modulus predicted by energy method and stiffness-averaging method both find good agreement with the experimental values, when taking the micro-porosity in fibers tows into consideration. Differences between the theoretical and experimental values become smaller.

前陆盆地挠曲沉降和沉积过程3-D模型研究

前陆盆地是在造山带负荷作用下岩石圈发生挠曲沉降而形成的,并且被主要从造山带搬运的沉积物所充填。为了更好地理解和认识前陆盆地的形成演化机制,特别是受控于周缘多个造山带活动所形成的前陆盆地的演化机制,本文通过建立前陆盆地挠曲沉降与沉积过程的3-D模型,模拟展示了造山带逆冲推覆作用、岩石圈挠曲沉降响应及在山盆体系中由于动力地形变化而导致的河流体系的发育变化及其产生的剥蚀和沉积过程。模型的建立和实验完整体现了逆冲推覆、弹性挠曲沉降和沉积物搬运这三者之间的耦合机制,为全面深入研究前陆盆地动力学提供了理论依据和方法。

Uniaxial experimental study of the acoustic emission and deformation behavior of composite rock based on 3D digital image correlation(DIC)

In this paper, uniaxial compression tests were carried out on a series of composite rock specimens with different dip angles, which were made from two types of rock-like material with different strength. The acoustic emission technique was used to monitor the acoustic signal characteristics of composite rock specimens during the entire loading process. At the same time, an optical non-contact 3 D digital image correlation technique was used to study the evolution of axial strain field and the maximal strain field before and after the peak strength at different stress levels during the loading process. The effect of bedding plane inclination on the deformation and strength during uniaxial loading was analyzed. The methods of solving the elastic constants of hard and weak rock were described. The damage evolution process, deformation and failure mechanism, and failure mode during uniaxial loading were fully determined. The experimental results show that the θ = 0?–45?specimens had obvious plastic deformation during loading, and the brittleness of the θ = 60?–90?specimens gradually increased during the loading process. When the anisotropic angle θincreased from 0?to 90?, the peak strength, peak strain,and apparent elastic modulus all decreased initially and then increased. The failure mode of the composite rock specimen during uniaxial loading can be divided into three categories:tensile fracture across the discontinuities(θ = 0?–30?), slid-ing failure along the discontinuities(θ = 45?–75?), and tensile-split along the discontinuities(θ = 90?). The axial strain of the weak and hard rock layers in the composite rock specimen during the loading process was significantly different from that of the θ = 0?–45?specimens and was almost the same as that of the θ = 60?–90?specimens. As for the strain localization highlighted in the maximum principal strain field, the θ = 0?–30?specimens appeared in the rock matrix approximately parallel to the loading direction,while in the θ = 45?–90?specimens it appeared at the hard and weak rock layer interface.

A Method of Obtaining Deformation Information by Tracing Deformation from Sections for 3D FE Simulation

With the development of computer technology and finite element method, the priority research area of plasticforming has focused on 3D FE simulation of forming processes for components with complicated geometrical shape.These processes have complex deforming mechanism, and different sections have different deforming characteristics.Therefore, for making a simple, convenient, and practical analysis of its deforming law, how to obtain deformationinformation of key sections from the results of 3D FE simulation has become one of problems urgently to be solved.So, a method of obtaining deformation information by tracing deformation from sections for 3D FE simulation hasbeen proposed. From the deformation information got by this method, the deformation law of key locations and thewhole deforming body can be obtained. This method can also help to compare the result from FE simulation withthat from physical modeling. Key procedures of this method have been presented in detail, and it has been testedby applying to 3D FE simulation of precision forging of the blade with a damper platform. The result shows thatthe method is practicable and reliable, and it can also be applied to 3D FE simulation of plastic forming processes ofother components.

Development and Properties of Glass Fiber Reinforced Plastics Geogrid

Glass fiber reinforced plastics geogrid has a wide application in the field of soil reinforcement because of its high strength, good toughness, and resistance to environmental stress, creep resistance and strong stability. In order to get high-powered glass fiber reinforced plastics geogrid and its mechanical characteristics, the properties and physical mechanical index of geogrid have been got through the study of its raw material, production process and important quality index. The analysis and study have been made to the geogrid＇s mechanical properties with loading speed, three-axial compression, temperature tensile test and FLAC3D numerical simulation, thus obtain the mechanical parameters of its displacement time curve, breaking strength and elongation at break. Some conclusions can be drawn as follows： （a） Using glass fiber materials, knurling and coated projection process, the f^acture strength and corrosion resistance of geogrid are greatly improved and the interlocking bite capability of soil is enhanced. （b） The fracture strength of geogrid is related to temperature and loading rate. When the surrounding rock pressure is fixed, the strength and anti-deformation ability of reinforced soil are significantly enhanced with increasing reinforced layers. （c） The pullout test shows the positive correlation between geogrid displacement and action time. （d） As a new reinforced material, the glass fiber reinforced plastics geogrid is not mature enough in theoretical research and practical experience, so it has become an urgent problem both in theoretical study and practical innovation.

Fiber Line Optimization in Single Ply for 3D Printed Composites

In conventional manufacturing processes of composites, Carbon Fibre Reinforced Plastic (CFRP) laminates have been made by stacking unidirectional or woven prepreg sheets. Recently, as a manufacturing process of CFRP, 3D printing of CFRP composites has been developed. The 3D printing process of CFRP composites enables us to fabricate CFRP laminates with arbitrary curvilinear fibre plies. This indicates that the optimization of the in-plane curved carbon fibre placement in a planar ply is strongly required to realize superior 3D printed composites. In the present paper, in-plane curved carbon fibre alignment of a ply with an open hole is optimized in terms of maximization of the fracture strength. For the optimization process, a genetic algorithm is adopted. To describe curved carbon fibre alignments in a planar ply, stream lines of perfect flow is employed. By using the stream lines of the perfect flow, number of optimization parameters is significantly reduced. After the optimization, the fracture strength of CFRP laminate is compared with the results of unidirectional CFRP ply. The curved fibre placement in a planar ply shows superior fracture improvement.

Scaling laws governing the elastic properties of 3D graphenes

In this study,we comprehensively investigated the scaling law for elastic properties of three-dimensional honeycomb-like graphenes(3D graphenes)using hybrid neural network potential-based molecular dynamics simulations and theoretical analyses.The elastic constants were obtained as functions of honeycomb hole size,denoted by the graphene wall length L.All five independent elastic constants in the large-L limit are proportional to L^(-1).The associated coefficients are combinations of elastic constants of two-dimensional graphene.High-order terms including L^(-2)and L^(-3)emerge for finite L values.They have three origins,the distorted areas close to the joint lines of 3D graphenes,the variation in solid angles between graphene plates,and the bending distortion of graphene plates.Significantly,the chirality becomes essential with decreasing L because the joint line structures are different between the armchair and zigzag-type 3D graphenes.Our findings provide insights into the elastic properties of graphene-based superstructures and can be used for further studies on graphene-based materials.

Recipe Development and Mechanical Characterization of Carbon Fibre Reinforced Recycled Polypropylene 3D Printing Filament

Recycled polypropylene filaments for fused filament fabrication were investigated with and without 14 wt% short fibre carbon reinforcements. The microstructure and mechanical properties of the filaments and 3D printed specimens were characterized using scanning electron microscopy and standard tensile testing. It was observed that recycled polypropylene filaments with 14 wt% short carbon fibre reinforcement contained pores that were dispersed throughout the microstructure of the filament. A two-stage filament extrusion process was observed to improve the spatial distribution of carbon fibre reinforcement but did not reduce the pores. Recycled polypropylene filaments without reinforcement extruded at high screw speeds above 20 rpm contained a centreline cavity but no spatially distributed pores. However, this cavity is eliminated when extrusion is carried out at screw speeds below 20 rpm. For 3D printed specimens, interlayer cavities were observed larger for specimens printed from 14 wt% carbon fibre reinforced recycled polypropylene than those printed from unreinforced filaments. The values of tensile strength for the filaments were 21.82 MPa and 24.22 MPa, which reduced to 19.72 MPa and 22.70 MPa, respectively, for 3D printed samples using the filaments. Likewise, the young’s modulus of the filaments was 1208.6 MPa and 1412.7 MPa, which reduced to 961.5 MPa and 1352.3 MPa, respectively, for the 3D printed samples. The percentage elongation at failure for the recycled polypropylene filament was 9.83% but reduced to 3.84% for the samples printed with 14 wt% carbon fiber reinforced polypropylene filaments whose elongation to failure was 6.58%. The SEM observations on the fractured tensile test samples showed interlayer gaps between the printed and the adjacent raster layers. These gaps accounted for the reduction in the mechanical properties of the printed parts.

Investigation of semi-solid aluminum alloy filling-plastic flowing in thixoforging

The semi-solid filling-plastic flowing integrated forging process of semi-solid 6061 Al alloy was simulated by commercial finite element software DEFORM-3D.Temperature,fluid and stress-strain fields were considered in numerical simulation.The simulation results show that the plastic deformation of billet of the ends is higher than that of billet in the straight cylinder.The value of plastic deformation varies with loading mode and plastic deformation fields at the stage of increasing pressure to constant value.When the thixoforging experiments were performed at 590 ℃,15 mm/s of punch velocity and 46 MPa of pressure side urn,it gets the filling wholly and dense internal organization of semi-solid thixoforging parts is gotten.Finite element analysis results are compatible with experimental ones.

Preliminary Report on the Theory of Elastic Circular Plate with no Kirchhoff-Love Assumptions

In this paper, the theory of elastic circular plate with no classical Kirchhoff-Love assumptions is established on the basis of a previous paper. In this theory, no classical Kirchhoff-Love assumptions are pre-assumed and the axial symmetrical analytic solution of fixed circular plate under the action of uniform pressure is obtained. Comparison of this solution and the known classical solution shows that this new solution agrees better than classical solution with the experiment measurement.This gives also the quantitative effect of the thickness on the deflection of circular plate with moderate thickness.