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.

Characteristics of an Axial-flux Permanent Magnet Synchronous Machine with Contra-rotating Rotors under Unbalanced Load Condition from 3-D Finite Element Analysis

During recent years,the axial-flus PMSM with contra-rotating rotors has become a hot topic in academic research due to its high efficiency and simple structure.However,its back-EMF may be distorted under the condition of different angular positions.This paper investigates characteristics of the novel motor used for contra-propeller driving.Considering the torque ripple and current oscillation under unbalanced load condition,this paper analyzes the distorted back-EMF of the machine when its two rotors get different angular positions during rotating.The analysis results are validated by transient-magnetic 3-D FEA method,which the 3-D FEA software is used to model this motor and transient simulations are carried out to obtain its magnetic characteristic and main performances.A main focus is put on the back-EMF characteristic with different angular positions between the two rotors.Furthermore,the characteristic of torque production under unbalanced load is investigated.Finally,a prototype motor is fabricated to validate the analyses of this paper.

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.

2.5-D modeling of cross-hole electromagnetic measurement by nite element method

A finite element method is developed for simulating frequency domain electromagnetic responses due to a dipole source in the 2-D conductive structures. Computing costs are considerably minimized by reducing the full three-dimensional problem to a series of two-dimensional problems. This is accomplished by transforming the problem into y-wave number （Ky） domain using Fourier transform and the y-axis is parallel to the structural strike. In the Ky domain, two coupled partial differential equations for magnetic field Hy and electric field Ey are derived. For a specific value of Ky, the coupled equations are solved by the finite element method with isoparametric elements in the x-z plane. Application of the inverse Fourier transform to the Ky, domain provides the electric and magnetic fields in real space. The equations derived can be applied to general complex two-dimensional structures containing either electric or magnetic dipole source in any direction. In the modeling of the electromagnetic measurement, we adopted a pseudo-delta function to distribute the dipole source current and circumvent the problem of singularity at the source point. Moreover, the suggested method used isoparametric finite elements to accommodate the complex subsurface formation. For the large scale linear system derived from the discretization of the Maxwell＇s equations, several iterative solvers were used and compared to select the optimal one. A quantitative test of accuracy was presented which compared the finite element results with analytic solutions for a dipole source in homogeneous space for different ranges and different wave numbers Ky. to validate the addressed the effects of the distribution range τ of the homogeneous medium. code and check its effectiveness. In addition, we pseudo-delta function on the numerical results in

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.

3-D Modelling of the Confederation Bridge Using Data of Full Scale Tests

Long-span bridges are special structures that require advanced analysis techniques to examine their performance. This paper presents a procedure developed to model the Confederation Bridge using 3-D beam elements. The model was validated using the data collected before the opening of the bridge to the public. The bridge was instrumented to conduct fullscale static and dynamic tests. The static tests were to measure the deflection of the bridge pier while the dynamic tests to measure the free vibrations of the pier due to a sudden release of the static load. Confederation Bridge is one of the longest reinforced concrete bridges in the world. It connects the province of Prince Edward Island and the province of New Brunswick in Canada. Due to its strategic location and vital role as a transportation link between these two provinces, it was designed using higher safety factors than those for typical highway bridges. After validating the present numerical model, a procedure was developed to evaluate the performance of similar bridges subjected to traffic and seismic loads. It is of interest to note that the foundation stiffness and the modulus of elasticity of the concrete have significant effects on the structural responses of the Confederation Bridge.

3D finite element analysis of a two-surface wear model in fretting tests

This article aims at developing a computationally efficient framework to simulate the erosion of two contact surfaces in three-dimensional(3D),depending on the body resistance.The framework involves finite element(FE)resolution of a fretting problem,wear computation via a non-local criterion including a wear distribution parameter(WDP),as well as updating of the geometry and automatic remeshing.Its originality is based on the capability to capture the damage on each surface and obtain local and global results for a quantitative and qualitative analysis.Numerical simulations are carried out for two 3D contact specimens with different values of WDP.The results highlight the importance of correctly modelling wear:One-surface wear model is sufficient from a global point of view(wear volume),or whenever the wear resistance for a body is much higher than that of another one,whereas a 3D two-surface wear model is essential to capturing local effects(contact pressure,wear footprint,etc.)related to the difference in wear resistance of the bodies.

TEMPORAL AND SPATIAL DISCRETIZATION ON QUASI-3-D GROUNDWATER FINITE ELEMENT MODELLING TO AVOID SPURIOUS OSCILLATION

In this article, the fmite element solution of quasi-three-dimensional （quasi-3-D） groundwater flow was mathematically analyzed. The research shows that the spurious oscillation solution to the Finite Element Model （FEM） is the results choosing the small time step △t or the large element size L and using the non-diagonal storage matrix. The mechanism for this phenomenon is explained by the negative weighting factor of implicit part in the discretized equations. To avoid spurious oscillation solution, the criteria on the selection of △t and L for quasi-3-D groundwater flow simulations were identified. An application example of quasi-3-D groundwater flow simulation was presented to verify the criteria. The results indicate that temporal discretization scale has significant impact on the spurious oscillations in the finite-element solutions, and the spurious oscillations can be avoided in solving practical quasi-3-D groundwater flow problems if the criteria are satisfied.

Benchmarking of two three-dimensional numerical models in time/space domain to predict railway-induced ground vibrations

In the last 30 years,the scientific community has developed and proposed different models and numerical approaches for the study of vibrations induced by railway traffic.Most of them are formulated in the frequency/wave number domain and with a 2.5D approach.Three-dimensional numerical models formulated in the time/space domain are less frequently used,mainly due to their high computational cost.Notwithstanding,these models present very attractive characteristics,such as the possibility of considering nonlinear behaviors or the modelling of excess pore pressure and non-homogeneous and non-periodic geometries in the longitudinal direction of the track.In this study,two 3D numerical approaches formulated in the time/space domain are compared and experimentally validated.The first one consists of a finite element approach and the second one of a finite difference approach.The experimental validation in an actual case situated in Carregado(Portugal)shows an acceptable fitting between the numerical results and the actual measurements for both models.However,there are some differences among them.This study therefore includes some recommendations for their use in practical soil dynamics and geotechnical engineering.

A Three-Dimensional Numerical Model for Predicting the Weld Bead Geometry Characteristics in Laser Overlap Welding of Low Carbon Galvanized Steel

Laser welding (LW) becomes one of the most economical high quality joining processes. LW offers the advantage of very controlled heat input resulting in low distortion and the ability to weld heat sensitive components. To exploit efficiently the benefits presented by LW, it is necessary to develop an integrated approach to identify and control the welding process variables in order to produce the desired weld characteristics without being forced to use the traditional and fastidious trial and error procedures. The paper presents a study of weld bead geometry characteristics prediction for laser overlap welding of low carbon galvanized steel using 3D numerical modelling and experimental validation. The temperature dependent material properties, metallurgical transformations and enthalpy method constitute the foundation of the proposed modelling approach. An adaptive 3D heat source is adopted to simulate both keyhole and conduction mode of the LW process. The simulations are performed using 3D finite element model on commercial software. The model is used to estimate the weld bead geometry characteristics for various LW parameters, such as laser power, welding speed and laser beam diameter. The calibration and validation of the 3D numerical model are based on experimental data achieved using a 3 kW Nd:Yag laser system, a structured experimental design and confirmed statistical analysis tools. The results reveal that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects on the weld quality. The results show great concordance between predicted and measured values for weld bead geometry characteristics, such as depth of penetration, bead width at the top surface and bead width at the interface between sheets, with an average accuracy greater than 95%.

Prediction of Weld Joint Shape and Dimensions in Laser Welding Using a 3D Modeling and Experimental Validation

This paper presents an experimentally validated weld joint shape and dimensions predictive 3D modeling for low carbon galvanized steel in butt-joint configurations. The proposed modelling approach is based on metallurgical transformations using temperature dependent material properties and the enthalpy method. Conduction and keyhole modes welding are investigated using surface and volumetric heat sources, respectively. Transition between the heat sources is carried out according to the power density and interaction time. Simulations are carried out using 3D finite element model on commercial software. The simulation results of the weld shape and dimensions are validated using a structured experimental investigation based on Taguchi method. Experimental validation conducted on a 3 kW Nd: YAG laser source reveals that the modelling approach can provide not only a consistent and accurate prediction of the weld characteristics under variable welding parameters and conditions but also a comprehensive and quantitative analysis of process parameters effects. The results show great concordance between predicted and measured values for the weld joint shape and dimensions.

Improved treatments for evaluating horizontal magnetic components through the 3-D FDM in E-polarization induction problems

To improve the accuracy of the numerical evaluation through the 3-D finite difference method, the surface boundary conditions are added to modify the old program. The author has tested the new program by making calculations for the model constructed by Wanamaker, et al (1984). The comparison between the numerical results obtained from this paper and those by Wannamaker, et al (1984) indicates that a pronounced improvement is realized in the evaluation of the horizontal magnetic components. Moreover, better calculations for the vertical magnetic components are also obtainable by using the new program.

改性1-3型压电复合材料理论模型与仿真

针对高机电耦合系数和曲面换能器的应用需求,该文提出一种由压电陶瓷、环氧树脂和硅橡胶三相复合的改进1-3型压电复合材料,基于Newnhams串、并联分析对复合材料的理论模型进行了推导,并利用Matlab进行了数值计算,分析了改性1-3型复合材料的结构参数对厚度机电耦合系数、等效密度和厚度频率常数的影响。结果表明,陶瓷-环氧复合相的体积分数为0.2~0.5,陶瓷-环氧复合相中的压电陶瓷体积分数为0.6~0.7时,改性1-3型压电复合材料的性能参数较好;利用ANSYS软件对复合材料进行了有限元仿真,将仿真结果与理论计算进行了对比,二者具有较好的一致性,验证了理论模型评估材料性能的准确性。

Mechanical Validation of Perfect Tibia 3D Model Using Computed Tomography Scan

In this paper, the Von Mises stresses and stiffnesses measured by experiments on a human cadaveric tibia and composite ones compared to those predicted by a FE model based on the same bone. Modeling of exact geometrical tibia including cortical and spongy bone using human bone CT scan images and mechanical validating of obtained model, is the aim of this study .The model produced by the current study supplies a tool for simulating mechanical test conditions on human tibia.

Numerical Investigation of Laser Surface Hardening of AISI 4340 Using 3D FEM Model for Thermal Analysis of Different Laser Scanning Patterns

Laser surface hardening is becoming one of the most successful heat treatment processes for improving wear and fatigue properties of steel parts. In this process, the heating system parameters and the material properties have important effects on the achieved hardened surface characteristics. The control of these variables using predictive modeling strategies leads to the desired surface properties without following the fastidious trial and error method. However, when the dimensions of the surface to be treated are larger than the cross section of the laser beam, various laser scanning patterns can be used. Due to their effects on the hardened surface properties, the attributes of the selected scanning patterns become significant variables in the process. This paper presents numerical and experimental investigations of four scanning patterns for laser surface hardening of AISI 4340 steel. The investigations are based on exhaustive modelling and simulation efforts carried out using a 3D finite element thermal analysis and structured experimental study according to Taguchi method. The temperature distribution and the hardness profile attributes are used to evaluate the effects of heating parameters and patterns design parameters on the hardened surface characteristics. This is very useful for integrating the scanning patterns’ features in an efficient predictive modeling approach. A structured experimental design combined to improved statistical analysis tools is used to assess the 3D model performance. The experiments are performed on a 3 kW Nd:Yag laser system. The modeling results exhibit a great agreement between the predicted and measured values for the hardened surface characteristics. The model evaluation reveals also its ability to provide not only accurate and robust predictions of the temperature distribution and the hardness profile as well an in-depth analysis of the effects of the process parameters.

Finite element modelling and experimental investigation of the impact of filling different materials in copper tubes during 3D free bending process

3D free bending process,an advanced metal forming technology,has attracted much attention due to its unique geometrical flexibility and efficiency.Filling tubes with materials may effectively reduce the distortion in the cross-section area and restrain wrinkling and collapse of tube wall.In the present study,the impacts of filling different materials on the copper tubes during 3D free bending process were investigated.We have found that copper tubes filled with low melting point alloys could induce the most uniform stress distribution,which result in better formability and moderate thickness changing rates when compared to SS304 steel balls and PU rubbers.The lowest ellipticity reached to 1.467%.

Accurate 3D Reconstruction of Subject-Specific Knee Finite Element Model to Simulate the Articular Cartilage Defects

The biomechanical relationship between the articular cartilage defect and knee osteoarthritis (OA) has not been clearly defined. This study presents a 3D knee finite element model (FEM) to determine the effect of cartilage defects on the stress distribution around the defect rim. The complete knee FEM, which includes bones, articular cartilages, menisci and ligaments, is developed from computed tomography and magnetic resonance images. This FEM then is validated and used to simulate femoral cartilage defects. Based on the obtained results, it is confirmed that the 3D knee FEM is reconstructed with high-fidelity level and can faithfully predict the knee contact behavior. Cartilage defects drastically affect the stress distribution on articular cartilages. When the defect size was smaller than 1.00cm2, the stress elevation and redistribution were found undistinguishable. However, significant stress elevation and redistribution were detected due to the large defect sizes ( 1.00cm2). This alteration of stress distribution has important implications relating to the progression of cartilage defect to OA in the human knee joint.

基于ABAQUS的平纹织物同面对向弯曲有限元模拟

为明晰织物结构对其弯曲性能的影响机制,对织物弯曲过程进行了有限元模拟和实验验证。以涤纶织物为例,通过VHX-5000型超景深数码显微镜观测织物,得到纱线几何参数;根据纱线实际尺寸,利用SolidWorks专业建模软件构建涤纶平纹织物的三维几何模型;基于有限元分析软件ABAQUS研究织物同面对向弯曲性能,分析纱线弹性模量、摩擦因数、泊松比等参数对弯曲实验的影响,并将有限元仿真结果与实验结果进行对比。结果表明:模拟所得抗弯力-位移曲线与实验曲线在0.01水平上呈显著相关,模拟结果与实验结果一致,证明用有限元模拟弯曲织物模型的有效性。

磁矢量势的全球电磁感应时间域三维有限元正演计算方法

全球电磁感应测深方法能获得地球深部电导率结构分布,目前被广泛应用于地球内部结构与热状态研究.地磁台站与地磁卫星观测的电磁感应数据为时间序列信号,在时间域分析全球电磁感应数据,特别是分析频谱极宽的磁暴脉冲响应,具有天然的优势.然而,当前的全球电磁感应数据解释技术一般在频率域进行,缺少时间域中的研究成果.为了弥补缺少时间域全球电磁感应数据解释方法的问题,本文开发了一种基于磁场矢量势的全球电磁感应时间域三维有限元并行正演求解方法,具备高精度快速计算源于地球外部时变电流源的地球感应电磁场时间序列的能力,特别适合于计算与分析频谱极宽的磁暴脉冲电磁感应时变响应.首先,建立基于磁场矢量势的全球电磁感应时间域微分控制方程,结合磁层外部电流源的物理属性建立边界条件和初始条件,从而构建出全球电磁感应时间域初始边界值问题.然后,利用四面体矢量有限元技术和无条件稳定的隐式后退欧拉公式,分别实现磁场矢量势的空间域和时间域离散,进而获得不同时刻的实系数大型有限元线性方程组,借助于高性能并行直接求解器,快速高精度地求解不同时刻的磁场矢量势与感应磁场.最后,利用理论模型验证本文算法的正确性.利用Dst磁暴环电流指数建立的时间域电流源与真实地球三维电性模型,研究"澳科一号"后续卫星200 km轨道对中国和澳大利亚下方地幔转换带高导体的探测能力,结果表明这些地幔转换带异常体在200 km卫星轨道高度能产生明显异常.综上所述,本文开发的时间域全球电磁感应方法不仅具备精确快速计算全球感应电磁场时间序列响应的能力,还能为反演与解释"澳科一号"等我国地磁卫星观测数据提供技术支持.