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.

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 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.

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.

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.

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.

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 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.

Numerical Study of the Biomechanical Behavior of a 3D Printed Polymer Esophageal Stent in the Esophagus by BP Neural Network Algorithm

Esophageal disease is a common disorder of the digestive system that can severely affect the quality of life andprognosis of patients. Esophageal stenting is an effective treatment that has been widely used in clinical practice.However, esophageal stents of different types and parameters have varying adaptability and effectiveness forpatients, and they need to be individually selected according to the patient’s specific situation. The purposeof this study was to provide a reference for clinical doctors to choose suitable esophageal stents. We used 3Dprinting technology to fabricate esophageal stents with different ratios of thermoplastic polyurethane (TPU)/(Poly-ε-caprolactone) PCL polymer, and established an artificial neural network model that could predict the radial forceof esophageal stents based on the content of TPU, PCL and print parameter. We selected three optimal ratios formechanical performance tests and evaluated the biomechanical effects of different ratios of stents on esophagealimplantation, swallowing, and stent migration processes through finite element numerical simulation and in vitrosimulation tests. The results showed that different ratios of polymer stents had different mechanical properties,affecting the effectiveness of stent expansion treatment and the possibility of postoperative complications of stentimplantation.

Effect of surface retaining elements on rock stability:laboratory investigation with sand powder 3D printing

This study aims to investigate the benefcial efects of surface retaining elements (SREs) on the mechanical behaviors of bolted rock and roadway stability. 3D printing (3DP) technology is utilized to create rock analogue prismatic specimens for conducting this investigation. Uniaxial compression tests with acoustic emission (AE) and digital image correlation techniques have been conducted on 3DP specimens bolted with diferent SREs. The results demonstrate that the strength and modulus of elasticity of the bolted specimens show a positive correlation with the area of the SRE;the AE characteristics of the bolted specimens are higher than those of the unbolted specimen, but they decrease with an increase in SRE area, thus further improving the integrity of the bolted specimens. The reinforcement efect of SREs on the surrounding rock of roadways is further analyzed using numerical modelling and feld test. The results provide a better understanding of the role of SREs in rock bolting and the optimization of rock bolting design. Furthermore, they verify the feasibility of 3DP for rock analogues in rock mechanics tests.

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.

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.

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.

A NEW METHOD FOR SOLUTION OF 3D ELASTIC-PLASTIC FRICTIONAL CONTACT PROBLEMS

The solution of 3 D elastic-plastic frictional contact problems belongs to the un specified boundary problems where the interaction between two kinds of nonlinearities should occur. Considering the difficulties for the solution of 3 D frictional contact problems, the key part is the determination of the tangential slip states at the contact points, and a great amount of computing work is needed for a high accuracy result. A new method based on a combination of programming and iteration methods, which are respectively known as two main kinds of methods for contact analysis, was put forward to deal with 3 D elastic-plastic contact problems. Numerical results demonstrate the efficiency of the algorithm illustrated here.

Modeling the 3D Terrain Effect on MT by the Boundary Element Method

A numerical method is put forward in this paper, using the boundary element method （BEM） to model 3D terrain effects on magnetotelluric （MT） surveys, Using vector integral theory and electromagnetic field boundary conditions, the boundary problem of two electromagnetic fields in the upper half space （air） and lower half space （earth medium） was transformed into two vector integral equations just related to the topography ： one magnetic equation for computing the magnetic field and the other electrical equation for computing the electrical field. The topography integral is decomposed into a series of integrals in a triangle element. For the integral in a triangle element, we suppose that the electromagnetic field in it is the stack of the electromagnetic field in the homogeneous earth and the topography response which is a constant; so the computation becomes simple, convenient and highly accurate. By decomposition and computation, each vector integral equation can be calculated by solving three linear equations that are related to the three Cartesian directions. The matrix of these linear equations is diagonally dominant and can be solved using the Symmetric Successive Over-Relaxation （SSOR） method. The apparent resistivity curve of MT on two 3D terrains calculated by BEM is shown in this paper.