Multi-scale Modeling and Finite Element Analyses of Thermal Conductivity of 3D C/SiC Composites Fabricating by Flexible-Oriented Woven Process

Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.

Contribution to the Full 3D Finite Element Modelling of a Hybrid Stepping Motor with and without Current in the Coils

The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.

3-D finite element modeling for evolution of stress field and interaction among strong earthquakes in Sichuan-Yunnan region

Based on the latest achievement about activities of geological structure, a 3-D finite-element model containing four layers of upper crust, lower crust (two layers) and upper mantle is established in the paper. By repeated tests and revisions, the boundary conditions of the model are determined. And then the background stress field, the stress field caused by fault creep and the stress field triggered by strong earthquake in Sichuan-Yunnan region, as well as their dynamic variations are calculated. The results indicate that the latter earthquake often occurs in the area with positive Coulomb rupture stress change associated with the former one, the former earthquake has a triggering effect on the latter one to a certain extent, and strong earthquake often occur in groups under the background of high stress, which is of great significance for distinguishing seismic anomalies, as well as for improving the level of earthquake prediction.

Grain refinement of Mg-Al alloys by optimization of process parameters based on three-dimensional finite element modeling of roll casting

To study the influence of roll casting process parameters on temperature and thermal-stress fields for the AZ31 magnesium alloy sheets,three-dimensional geometric and 3D finite element models for roll casting were established based on the symmetry of roll casting by ANSYS software.Meshing method and smart-sizing algorithm were used to divide finite element mesh in ANSYS software.A series of researches on the temperature and stress distributions during solidification process with different process parameters were done by 3D finite element method.The temperatures of both the liquid-solid two-phase zone and liquid phase zone were elevated with increasing pouring temperature.With the heat transfer coefficient increasing,the two-phase region for liquid-solid becomes smaller.With the pouring temperature increasing and the increase of casting speed,the length of two-phase zone rises.The optimized of process parameters（casting speed 2 m/min,pouring temperature 640 ℃ and heat transfer coefficient 15 kW/（m2·℃） with the water pouring at roller exit was used to produce magnesium alloy AZ31 sheet,and equiaxed grains with the average grain size of 50 μm were achieved after roll casting.The simulation results give better understanding of the temperature variation in phase transformation zone and the formation mechanism of hot cracks in plates during roll casting and help to design the optimized process parameters of roll casting for Mg alloy.

3D finite element analysis on pile-soil interaction of passive pile group

The interaction between pile and soft soil of the passive pile group subjected to soil movement was analyzed with three-dimensional finite element model by using ANSYS software. The soil was assumed to be elastic-plastic complying with the Drucker-Prager yield criterion in the analysis. The large displacement of soil was considered and contact elements were used to evaluate the interaction between pile and soil. The influences of soil depth of layer and number of piles on the lateral pressure of the pile were investigated, and the lateral pressure distributions on the (2×1) pile group and on the (2×2) pile group were compared. The results show that the adjacent surcharge may result in significant lateral movement of the soft soil and considerable pressure on the pile. The pressure acting on the row near the surcharge is higher than that on the other row, due to the "barrier" and arching effects in pile groups. The passive load and its distribution should be taken into account in the design of the passive piles.

A Numerical Algorithm Based on Quadratic Finite Element for Two-Dimensional Nonlinear Time Fractional Thermal Diffusion Model

In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+
t2)is obtained,which is verified by the numerical results.

Three-dimensional Modeling of Ore-forming Elements and Mineralization Prognosis for the Yechangping Mo Deposit,Henan Province,China

Three-dimensional geochemical modeling of ore-forming elements is crucial for predicting deep mineralization.This approach provides key information for the quantitative prediction of deep mineral localization,three-dimensional fine interpolation,analysis of spatial distribution patterns,and extraction of quantitative mineral-seeking markers.The Yechangping molybdenum(Mo)deposit is a significant and extensive porphyry-skarn deposit in the East Qinling-Dabie Mo polymetallic metallogenic belt at the southern margin of the North China Block.Abundant borehole data on oreforming elements underpin deep geochemical predictions.The methodology includes the following steps:(1)Threedimensional geological modeling of the deposit was established.(2)Correlation,cluster,and factor analyses post delineation of mineralization stages and determination of mineral generation sequence to identify(Cu,Pb,Zn,Ag)and(Mo,W,mfe)assemblages.(3)A three-dimensional geochemical block model was constructed for Mo,W,mfe,Cu,Zn,Pb,and Ag using the ordinary kriging method,and the variational function was developed.(4)Spatial distribution and enrichment characteristics analysis of ore-forming elements are performed to extract geological information,employing the variogram and w(Cu+Pb+Zn+Ag)/w(Mo+W)as predictive indicators.(5)Identifying the western,northwestern,and southwestern areas of the mine with limited mineralization potential,contrasted by the northeastern and southeastern areas favorable for mineral exploration.

3-D Finite-Element Numerical Simulation of Centrifugal Induction Electrosalg Casting Solidification Process

A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series of corresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of temperature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R-1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESC casting and structure or casting is fine and compact.

3D Finite Element Simulation of Tunnel Boring Machine Construction Processes in Deep Water Conveyance Tunnel

Applying stiffness migration method,a 3D finite element mechanical model is established to simulate the excavation and advance processes.By using 3D nonlinear finite element method,the tunnel boring machine(TBM) excavation process is dynamically simulated to analyze the stress and strain field status of surrounding rock and segment.The maximum tensile stress of segment ring caused by tunnel construction mainly lies in arch bottom and presents zonal distribution.The stress increases slightly and limitedly in the course of excavation.The maximum and minimum displacements of segment,manifesting as zonal distribution,distribute in arch bottom and vault respectively.The displacements slightly increase with the advance of TBM and gradually tend to stability.

Propagations of Rayleigh and Love waves in ZnO films/glass substrates analyzed by three-dimensional finite element method

Propagation characteristics of surface acoustic waves（SAWs） in ZnO films/glass substrates are theoretically investigated by the three-dimensional（3D） finite element method. At first, for（11ˉ20） ZnO films/glass substrates, the simulation results confirm that the Rayleigh waves along the [0001] direction and Love waves along the [1ˉ100] direction are successfully excited in the multilayered structures. Next, the crystal orientations of the ZnO films are rotated, and the influences of ZnO films with different crystal orientations on SAW characterizations, including the phase velocity, electromechanical coupling coefficient, and temperature coefficient of frequency, are investigated. The results show that at appropriate h/λ, Rayleigh wave has a maximum k^2 of 2.4% in（90°, 56.5°, 0°） ZnO film/glass substrate structure; Love wave has a maximum k^2 of 3.81% in（56°, 90°, 0°） ZnO film/glass substrate structure. Meantime, for Rayleigh wave and Love wave devices, zero temperature coefficient of frequency（TCF） can be achieved at appropriate ratio of film thickness to SAW wavelength. These results show that SAW devices with higher k^2 or lower TCF can be fabricated by flexibly selecting the crystal orientations of ZnO films on glass substrates.

A multiscale 3D finite element analysis of fluid/solute transport in mechanically loaded bone

The transport of fluid, nutrients, and signaling molecules in the bone lacunar-canalicular system （LCS） is critical for osteocyte survival and function. We have applied the fluorescence recovery after photobleaching （FRAP） approach to quantify load-induced fluid and solute transport in the LCS in situ, but the measurements were limited to cortical regions 30-50 μm underneath the periosteum due to the constrains of laser penetration. With this work, we aimed to expand our understanding of load-induced fluid and solute transport in both trabecular and cortical bone using a multiscaled image-based finite element analysis （FEA） approach. An intact murine tibia was first re-constructed from microCT images into a three-dimensional （3D） linear elastic FEA model, and the matrix deformations at various locations were calculated under axial loading. A segment of the above 3D model was then imported to the biphasic poroelasticity analysis platform （FEBio） to predict load-induced fluid pressure fields, and interstitial solute/fluid flows through LCS in both cortical and trabecular regions. Further, secondary flow effects such as the shear stress and/or drag force acting on osteocytes, the presumed mechano-sensors in bone, were derived using the previously developed ultrastructural model of Brinkman flow in the canaliculi. The material properties assumed in the FEA models were validated against previously obtained strain and FRAP transport data measured on the cortical cortex. Our results demonstrated the feasibility of this computational approach in estimating the fluid flux in the LCS and the cellular stimulation forces （shear and drag forces） for osteocytes in any cortical and trabecular bone locations, allowing further studies of how the activation of osteocytes correlates with in vivo functional bone formation. The study provides a promising platform to reveal potential cellular mechanisms underlying the anabolic power of exercises and physical activities in treating patients with skeletal deficiencies.

A New Finite Element Model with Manufactured Error for Additive Manufacturing

Additive manufacturing(AM),adding materials layer by layer,can be used to produce objects of almost any shape or geometry.However,AM techniques cannot accurately build parts with large overhangs,especially for the large features close to horizontal,hanging over the void.The overhangs will make the manufactured model deviate from the design model,which will result in the performance of the manufactured model that cannot satisfy the design requirements.In this paper,we will propose a new finite element(FE)analysis model that includes the manufacturing errors by mimicking the AM layer by layer construction process.In such FE model,an overhang coefficient is introduced to each FE,which is defined by the support elements in the lower layer.By mimicking the AM process from the bottom layer to the top layer,all the FE properties are updated based on their overhang coefficients,which makes the computational model be able to predict the manufactured model with manufacturing errors.The proposed model can be used to predict the performance of the AM objects in the design stage,which will help the designers to improve their design by the simulation results.

STUDY ON THREE-DIMENSIONAL FINITE BODIES CONTAINING CRACKS USING THE FINITE ELEMENT METHOD OF LINES

The three-dimensional finite element method of lines is presented, and the basic processing description of 3D FEMOL in cracking questions is given in detail. Applications to 3D bodies with cracks indicate that good accuracy can be obtained with relatively coarse girds. In particular, application to the tension specimen shows very good agreement with the evaluation of stress intensity factors, which is better than the results of other methods. This implies a considerable potential for using this method in the 3D analysis of finite geometry solids and suggests a possible extension of this technique to nonlinear material behavior.

3D Finite Elements Technique for Collapse Causes of the Pylons in Egyptian Temples: A Study of the Great Pylon of Ramesseum Temple, Luxor, Egypt

This research presents damage causes of the pylons in the ancient Egyptian temples based on 3D finite elements analysis. The main purpose of the research determines the failure causes of the first pylon of the Ramessium temple, which is situated in Upper Egypt, at Luxor “Thebes” on the west bank of the Nile River. The first pylon of Ramessium temple subjected to seismic activity effects on long term, combined with several structural damage factors such as the defects resulting from the construction technique, where the builder used the poor quality of stones in foundations of the pylon, the building materials residue was used as filler for the core of the pylon walls, and it lacked vertical joints between the courses. In addition to it founded on alluvial soil that is vulnerable to contaminated water, it is still suffering damage factors and urban trespasses at the moment. All of the former factors helped the pylon to be affected by the earthquakes loads that occurred on it. The structural behavior of the pylon under self-weight and earthquakes loads were carried out by Numerical analysis to find out the loads and stresses which caused collapsing of the pylon. Results of the study indicated that the pylon subjected to a horizontal displacement due to old earthquakes force, led to collapse of the pylon. Finally, the study represents use of modern technique to study the structural behavior of the most important architectural units in ancient Egyptian temples to identify the causes of its collapse.

A study on Parallel Computation Based on Finite Element Forward Modeling of 2D Magnetotelluric

Magnetotelluric sounding method based on the difference of the rock’s resistivity is an exploration method about doing research in earth’s resistivity and phase using the native electromagnetic field. The paper adopts 2D finite element method as the magnetotelluric forward method and calculates the total field by primary field (also named background field) plus secondary field. We can?get more accurate forward result through the finite element method and we can get the result effected by the dense degree of grid slightly by the total field. But the method is not effective?enough when the model is divided into relative big grid. When the frequency changes, program solves relevant equation separately. According to the feature of the algorithm, we apply MPI parallel method in the algorithm. Every process solves relevant equation. The account of frequency?that a process needs to solve in parallel computation is less than the account that the process?needs to solve in serial algorithm. We can see that the forward result is the same with the serial algorithm and proves the correctness of algorithm. We do statistics about the efficiency of the parallel algorithm. When the account of processes is from 2 to 8, the speedup is from 1.63 to 2.64. It proves the effectiveness of the parallel algorithm.

3D numerical simulation in acoustic-elastic coupled media with staggered-grid finite-difference method

Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.

3D Finite Element Analysis of a Man Hip Joint Femur under Impact Loads

Objective： The biomechanical characters of the bone fracture of the man femoral hip joint under impact loads are explored. Methods ：A biosystem model of the man femoral hip joint by using the GE （ General Electric） lightspeed multi-lay spiral CT is conducted. A 3D finite element model is established by employing the finite element software ANSYS. The FE analysis mainly concentrates on the effects of the impact directions arising from intense movements and the parenchyma on the femoral hip joint on the stress distributions of the proximal femur. Results：The parenchyma on the hip joint has relatively large relaxation effect on the impact loads. Conclusion：Effects of the angle δ of the impact load to the anterior direction and the angle γ of the impact load to the femur shaft on the bone fracture are given;δ has larger effect on the stress and strain distributions than the angle γ,which mainly represents the fracture of the upper femur including the femoral neck fracture when the posterolateral femur is impacted, consistent with the clinical resuits.

Spatially-graded 3D-printed viscoelastic truss metamaterials for impact trajectory control and energy absorption

This study demonstrates that two-and three-dimensional spatially graded,truss-based polymeric-material metamaterials can be designed for beneficial impact mitigation and energy absorption capabilities.Through a combination of numerical and experimental techniques,we highlight the broad property space of periodic viscoelastic trusses,realized using 3D printing via selective laser sintering.Extending beyond periodic designs,we investigate the impact response of spatially variant viscoelastic lattices in both two and three dimensions.Our result reveal that introducing spatial variations in lattice topology allows for redirecting of the impact trajectory,opening new opportunities for engineering and tailoring lightweight materials with target impact functionality.This is achieved through the combined selection of base material and metamaterial design.

Numerical analysis of the failure process of soil-rock mixtures through computed tomography and PFC3D models

Soil-rock mixture （SRM） is a unique type of geomaterial characterized by a heterogeneous composition and a complicated structure. It is intractable for the continuum-based soil and rock mechanics theories to accurately characterize and predict the SRM＇s mechanical properties. This study reports a novel numerical method incorporating microfocus computed tomography and PFC3D codes to probe the deformation and failure processes of SRM. The three-dimensional （3D） PFC models that represent the SRM＇s complex structures were built. By simulating the entire failure process in PFC3D, the SRM＇s strength, elastic modulus and crack growth were obtained. The influence of rock ratios on the SRM＇s strength, deformation and failure processes, as well as its internal mesoscale mechanism, were analyzed. By comparing simulation results with experimental data, it was verified that the 3D PFC models were in good agreement with SRM＇s real structure and the SRM＇s compression process, deformation and failure patterns; its intrinsic mesomechanism can be effectively analyzed based on such 3D PFC models.

3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing

A grain-based distinct element model featuring three-dimensional （3D） Voronoi tessellations （randompoly-crystals） is proposed for simulation of crack damage development in brittle rocks. The grainboundaries in poly-crystal structure produced by Voronoi tessellations can represent flaws in intact rockand allow for numerical replication of crack damage progression through initiation and propagation ofmicro-fractures along grain boundaries. The Voronoi modelling scheme has been used widely in the pastfor brittle fracture simulation of rock materials. However the difficulty of generating 3D Voronoi modelshas limited its application to two-dimensional （2D） codes. The proposed approach is implemented inNeper, an open-source engine for generation of 3D Voronoi grains, to generate block geometry files thatcan be read directly into 3DEC. A series of Unconfined Compressive Strength （UCS） tests are simulated in3DEC to verify the proposed methodology for 3D simulation of brittle fractures and to investigate therelationship between each micro-parameter and the model＇s macro-response. The possibility of numericalreplication of the classical U-shape strength curve for anisotropic rocks is also investigated innumerical UCS tests by using complex-shaped （elongated） grains that are cemented to one another alongtheir adjoining sides. A micro-parameter calibration procedure is established for 3D Voronoi models foraccurate replication of the mechanical behaviour of isotropic and anisotropic （containing a fabric） rocks. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.