Galerkin-based quasi-smooth manifold element(QSME)method for anisotropic heat conduction problems in composites with complex geometry

The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.

In silico optimization of actuation performance in dielectric elastomercomposites via integrated finite element modeling and deep learning

Dielectric elastomers(DEs)require balanced electric actuation performance and mechanical integrity under applied voltages.Incorporating high dielectric particles as fillers provides extensive design space to optimize concentration,morphology,and distribution for improved actuation performance and material modulus.This study presents an integrated framework combining finite element modeling(FEM)and deep learning to optimize the microstructure of DE composites.FEM first calculates actuation performance and the effective modulus across varied filler combinations,with these data used to train a convolutional neural network(CNN).Integrating the CNN into a multi-objective genetic algorithm generates designs with enhanced actuation performance and material modulus compared to the conventional optimization approach based on FEM approach within the same time.This framework harnesses artificial intelligence to navigate vast design possibilities,enabling optimized microstructures for high-performance DE composites.

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.

Low-velocity Impact Damage Analysis of Composite Laminates Using Self-adapting Delamination Element Method

On the basis ofa 2D 4-node Mindlin shell element method, a novel self-adapting delamination finite element method is presented, which is developed to model the delamination damage of composite laminates. In the method, the sublaminate elements are generated automatically when the delamination damage occurs or extends. Thus, the complex process and state of delamination damage can be simulated practically with high efficiency for both analysis and modeling. Based on the self-adapting delamination method, linear dynamic finite element damage analysis is performed to simulate the low-velocity impact damage process of three types of mixed woven composite laminates. Taking the frictional force among sublaminations during delaminating and the transverse normal stress into account, the analytical results are consistent with those of the experimental data.

FINITE ELEMENT MODELING OF MODERATELY THICK COMPOSITE PLATE WITH PIEZOELECTRIC SENSORS AND ACTUATORS

A rectangular finite element for laminated plate with bonded and/or embedded piezoelectric sensors and actuators is developed based on the variational principle and the first order shear deformation theory. The element has four-node, 20-degrees-of-freedom with one potential degree of freedom for each piezoelectric layer to represent the piezoelectric behavior. The higher order derivation of deflection is obtained by using the normal rotation expressions to take the effects of transverse shear deformation into considerations. The finite element can accurately simulate the deformation of both thin and moderately thick plates. A Fortran program is written and a number of benchmark tests are exercised to verify its effectiveness. Results are compared well with the existing data. The unbalanced composite with piezoelectric layers is then analyzed by using the model. Results show that the changes of the ratio between the thickness of positive angle layers and the negative angle layers have an effect on the deformation of the structure under the same electric loading.

Seismic performance of double-skin steel-concrete composite box piers: Part Ⅱ—Nonlinear finite element analysis

An accurate finite element （ FE） model was constructed to examine the hysteretic behavior of double-skin steel-concrete composite box （ DSCB） piers for further understanding the seismic performance of DSCB piers; where the local buckling behavior of steel tubes, the confinement of the in-filled concrete and the interface action between steel tube and in-filled concrete were considered. The accuracy of the proposed FE model was verified by the bidirectional cyclic loading test results. Based on the validated FE model, the effects of some key parameters, such as section width to steel thickness ratio, slenderness ratio, aspect ratio and axial load ratio on the hysteretic behavior of DSCB piers were investigated. Finally, the skeleton curve model of DSCB piers was proposed. The numerical simulation results reveal that the peak strength and elastic stiffness decrease with the increase of the section width to steel thickness ratio. Moreover, the increase of the slenderness ratio may result in a significant reduction in the peak strength and elastic stiffness while the ultimate displacement increases. The proposed skeleton curve model can be taken as a reference for seismic performance analyses of the DSCB piers.

Finite element analysis of stress distribution and burst failure of SiC_f/Ti-6Al-4V composite ring

A three-dimensional cyclic symmetry finite element model of titanium-matrix composites（TMCs） ring was developed to investigate the stress distribution and burst failure. The effects of fiber volume fractions, reinforced areas, thermal residual stresses and two different temperatures on stress distribution were studied. The burst speed was obtained through analyzing the hoop tensile stresses under a series of rotating speeds. The results indicate that at the two different temperatures, the influences of fiber volume fractions and reinforced areas on stress level and distribution are different. Some proposals are provided for the structure design of the TMCs ring. With regard to thermal residual stresses, a larger reinforced area is an advisable choice for design of the ring at higher temperature.

Method to measure tree-ring width,density,elemental composition,and stable carbon and oxygen isotopes using one sample

Tree-ring width(RW),density,elemental com-position,and stable carbon and oxygen isotope(δ^(13)C,δ^(18)O)are widely used as proxies to assess climate change,ecology,and environmental pollution;however,a specific pretreat-ment has been needed for each proxy.Here,we developed a method by which each proxy can be measured in the same sample.First,the sample is polished for ring width meas-urement.After obtaining the ring width data,the sample is cut to form a 1-mm-thick wood plate.The sample is then mounted in a vertical sample holder,and gradually scanned by an X-ray beam.Simultaneously,the count rates of the fluorescent photons of elements(for chemical characteriza-tion)and a radiographic grayscale image(for wood density)are obtained,i.e.the density and the element content are obtained.Then,cellulose is isolated from the 1-mm wood plate by removal of lignin,and hemicellulose.After producing this cellulose plate,cellulose subsamples are separated by knife under the microscope for inter-annual and intra-annual stable carbon and oxygen isotope(δ^(13)C,δ^(18)O)analysis.Based on this method,RW,density,elemental composition,δ^(13)C,and δ^(18)O can be measured from the same sample,which reduces sample amount and treatment time,and is helpful for multi-proxy comparison and combination research.

Bending stress of rolling element in elastic composite cylindrical roller bearing

A new structure design method of elastic composite cylindrical roller bearing is proposed, in which PTFE is embedded into a hollow cylindrical rolling element, according to the principle of creative combinations and through innovation research on cylindrical roller bearing structure. In order to systematically investigate the inner wall bending stress of the rolling element in elastic composite cylindrical roller bearing, finite element analysis on different elastic composite cylindrical rolling elements was conducted. The results show that, the bending stress of the elastic composite cylindrical rolling increases along with the increase of hollowness with the same filling material. The bending stress of the elastic composite cylindrical rolling element decreases along with the increase of the elasticity modulus of the material under the same physical dimension. Under the same load, on hollow cylindrical rolling element, the maximum bending tensile stress values of the elastic composite cylindrical rolling element after material filling at 0° and 180° are 8.2% and 9.5%, respectively, lower than those of the deep cavity hollow cylindrical rolling element. In addition, the maximum bending-compressive stress value at 90° is decreased by 6.1%.

Application of the quadrilateral area coordinate method:a new element for laminated composite plate bending problems

Recently, some new quadrilateral finite elements were successfully developed by the Quadrilateral Area Coordinate （QAC） method. Compared with those traditional models using isoparametric coordinates, these new models are less sensitive to mesh distortion. In this paper, a new displacement-based, 4-node 20-DOF （5-DOF per node） quadrilateral bending element based on the first-order shear deformation theory for analysis of arbitrary laminated composite plates is presented. Its bending part is based on the element AC-MQ4, a recent-developed high-performance Mindlin-Reissner plate element formulated by QAC method and the generalized conforming condition method; and its in-plane displacement fields are interpolated by bilinear shape functions in isoparametric coordinates. Furthermore, the hybrid post-rocessing procedure, which was firstly proposed by the authors, is employed again to improve the stress solutions, especially for the transverse shear stresses. The resulting element, denoted as AC-MQ4-LC, exhibits excellent performance in all linear static and dynamic numerical examples. It demonstrates again that the QAC method, the generalized conforming condition method, and the hybrid post-processing procedure are efficient tools for developing simple, effective and reliable finite element models.

Evaluation of interfacial properties in SiC composites using an improved cohesive element method

A two-dimensional axisymmetric finite element model based on an improved cohesive element method was developed to simulate interfacial debonding, sliding friction, and residual thermal stresses in SiC composites during single-fiber push-out tests to extract the interfacial bond strength and frictional stress. The numerical load–displacement curves agree well with experimental curves,indicating that this cohesive element method can be used for calculating the interfacial properties of SiC composites.The simulation results show that cracks are most likely to occur at the ends of the experimental sample, where the maximum shear stress is observed and that the interfacial shear strength and constant sliding friction stress decrease with an increase in temperature. Moreover, the load required to cause complete interfacial failure increases with the increase in critical shear strength, and the composite materials with higher fiber volume fractions have higher bearing capacities. In addition, the initial failure load increases with an increase in interphase thickness.

Three-Dimensional Reconstructed Finite Element Model for C/C Composites by Micro-CT

The precise microscopic feature of carbon-carbon（C/C） composites is essential {or an accurate predic tion of their mechanical behavior. After fabrication, actual microscopic feature differs from simple ideal spatial model. Micro computed lomography（CT） scan can well describe internal microstruetures of composites. Therefore, a reconstructed model is developed based on mireo-CT, by a series of prodcedures including extrac tlng components, generating new binary images and establishing a finite element （FE） model. Compared with the model designed by reconstructed commercial software MIMICS. the presented reconstructed FE model is superior in terms of high mesh quality and eontrollable mesh cluantity. The precision of the model is verified by experiment.

A new higher-order shear deformation theory and refined beam element of composite laminates

A new higher-order shear deformation theory based on global-local superposition technique is developed. The theory satisfies the free surface conditions and the geometric and stress continuity conditions at interfaces. The global displacement components are of the Reddy theory and local components are of the internal first to third-order terms in each layer. A two-node beam element based on this theory is proposed. The solutions are compared with 3D-elasticity solutions. Numerical results show that present beam element has higher computational efficiency and higher accuracy.

Damage prediction for magnesium matrix composites formed by liquid-solid extrusion process based on finite element simulation

A damage prediction method based on FE simulation was proposed to predict the occurrence of hot shortness crocks and surface cracks in liquid-solid extrusion process. This method integrated the critical temperature criterion and Cockcroft ＆ Latham ductile damage model, which were used to predict the initiation of hot shortness cracks and surface cracks of products, respectively. A coupling simulation of deformation with heat transfer as well as ductile damage was carried out to investigate the effect of extrusion temperature and extrusion speed on the damage behavior of Csf/AZ91D composites. It is concluded that the semisolid zone moves gradually toward deformation zone with the punch descending. The amplitude of the temperature rise at the exit of die from the initial billet temperature increases with the increase of extrusion speed during steady-state extrusion at a given punch displacement. In order to prevent the surface temperature of products beyond the incipient melting temperature of composites, the critical extrusion speed is decreased with the increase of extrusion temperature, otherwise the hot shortness cracks will occur. The maximum damage values increase with increasing extrusion speed or extrusion temperature. Theoretical results obtained by the Deform^TM-2D simulation agree well with the experiments.

Regulating element distribution to improve magnetic properties of sintered Nd–Fe–B/Tb–Fe–B composite magnets

Nd content was varied in Nd13.2-xFe80.8+xB6(x = 0, 0.5, 1, and 1.5) to optimize the magnetic properties of sintered Nd–Fe–B/Tb–Fe–B composite magnets, which were prepared by mixing 9 g of Nd–Fe–B with 1 g of Tb17Fe75B8 powder.In conventional magnets, by reducing Nd content, the coercivity of 10.4 kOe in Nd13.2Fe80.8B6 decreases to 7.2 kOe in Nd12.2Fe81.8B6;meanwhile, in Nd–Fe–B/Tb–Fe–B magnets the coercivity does not decrease when reducing Nd content.In the intergranular phase, the Tb content increases owing to the reducing Nd content of the Nd–Fe–B alloy in the sintered composite magnets.Therefore, the excess Tb in Tb17Fe75B8 enters the intergranular phase, and more Tb atoms can substitute for Nd at the grain boundary of the Nd–Fe–B phase, leading to a more significant increase in coercivity.The remanence increases with reducing Nd content, and the energy product of 39.1 MGOe with a high coercivity of 21.0 kOe is obtained in Nd12.2Fe81.8B6/Tb17Fe75B8 magnets.These investigations show that magnetic properties can be further improved by regulating the element distribution in sintered composite magnets.

Theoretical and numerical research on the dynamic launch response of carbon fiber composite cartridges

Understanding the dynamic response of composite material cartridges during the firing process is of great significance for improving their reliability and safety.A theoretical model describing the dynamic response of composite material cartridges is established based on the thick-walled cylinder theory and rate-dependent constitutive model of composite materials.The correctness of the theoretical model is validated through finite element simulations of cartridge deformation.The influence of chamber pressure and cartridge wall thickness on the cartridge's deformation process and stress distribution is analyzed.The results indicate that the primary deformation of composite material cartridges inside the chamber is elastic deformation.Compared to metal cartridges,composite material cartridges require higher pressure for touching-chamber and are more prone to developing gaps after unloading to ensure smooth extraction.During the deformation process,the touching-chamber behavior of the cartridge can improve the stress distribution.Under the same chamber pressure,the touching-chamber behavior can reduce the circumferential stress by approximately 30%.The inner wall surface of the cartridge is a critical area that requires attention.The touching-chamber behavior can be facilitated by appropriately reducing the cartridge wall thickness while ensuring overall strength.This study can provide guidance for the optimization design of composite material cartridges.

Effects of cement-enhanced soil on the ultimate lateral resistance of composite pile in clayey soil

The composite pile consisting of core-pile and surrounding cement-enhanced soil is a promising pile foundation in recent years.However,how and to what extent the cement-enhanced soil influences the ultimate lateral resistance has not been fully investigated.In this paper,the ultimate lateral resistance of the composite pile was studied by finite element limit analysis(FELA)and theoretical upper-bound analysis.The results of FELA and theoretical analysis revealed three failure modes of laterally loaded composite piles.The effects of the enhanced soil thickness,strength,and pile-enhanced soil interface characteristics on the ultimate lateral resistance were studied.The results show that increasing the enhanced soil thickness leads to a significant improvement on ultimate lateral resistance factor(N P),and there is a critical thickness beyond which the thickness no longer affects the N P.Increasing the enhanced soil strength induced 6.2%-232.6%increase of N P.However,no noticeable impact was detected when the enhanced soil strength was eight times higher than that of the natural soil.The maximum increment of N P is only 30.5%caused by the increase of interface adhesion factor(a).An empirical model was developed to calculate the N P of the composite pile,and the results show excellent agreement with the analytical results.

THE ANALYSIS OF THIN WALLED COMPOSITE LAMINATED HELICOPTER ROTOR WITH HIERARCHICAL WARPING FUNCTIONS AND FINITE ELEMENT METHOD

In the present paper, a series of hierarchical warping functions is developed to analyze the static and dynamic problems of thin walled composite laminated helicopter rotors composed of several layers with single closed cell. This method is the development and extension of the traditional constrained warping theory of thin walled metallic beams, which had been proved very successful since 1940s. The warping distribution along the perimeter of each layer is expanded into a series of successively corrective warping functions with the traditional warping function caused by free torsion or free beading as the first term, and is assumed to be piecewise linear along the thickness direction of layers. The governing equations are derived based upon the variational principle of minimum potential energy for static analysis and Rayleigh Quotient for free vibration analysis. Then the hierarchical finite element method. is introduced to form a,. numerical algorithm. Both static and natural vibration problems of sample box beams axe analyzed with the present method to show the main mechanical behavior of the thin walled composite laminated helicopter rotor.

Effect of Silver Element on Microstructure and Properties of W-30Cu/TiC Composites

W-30 wt%Cu and TiC-50 wt%Ag were successfully synthesized by a novel simplified pretreatment followed by electroless plating. The 0 wt% TiC, 0.5 wt% TiC, and 0.5 wt%TiC-0.5 wt%Ag composite powders were added to W-30 wt%Cu composite powders by blending, and then reduced. The reduced W-30 Cu, W-30 Cu/0.5 TiC, and W-30 Cu-0.5 Ag/0.5 TiC composite powders were then compacted and sintered at 1 300 ℃ in protective hydrogen for 60 min. The phase and morphology of the composite powders and materials were analyzed using X-ray diffraction and field emission scanning electron microscopy. The relative density, electrical conductivity, and hardness of the sintered samples were examined. Results showed that W-30 Cu and TiC-Ag composite powders with uniform structure were obtained using simplified pretreatment followed by electroless plating. The addition of TiC particles can significantly increase the compressive strength and hardness of the W-30 Cu composite material but decrease the electrical conductivity. Next, 0.5 wt% Ag was added to prepare W-30 Cu-0.5 Ag/TiC composites with excellent electrical conductivity. The electrical conductivity of these composites(61.2%) is higher than that in the national standard(the imaginary line denotes electrical conductivity of GB IACS 42%) of 45.7%.

3D FINITE ELEMENT ANALYSIS OF THE DAMAGE EFFECTS ON THE DENTAL COMPOSITE SUBJECT TO WATER SORPTION

The damage effects of water sorption on the mechanical properties of the hydroxyapatite particle reinforced Bis-GMA/TEGDMA copolymer （HA/Bis-GMA/TEGDMA） h-ave been predicted using 3D finite cell models. The plasticizer effect on the polymer matrix was considered as a variation of its Young＇s modulus. Three different cell models were used to determine the influence of varying particle contents, interphase strength and moisture concentration on the debonding damage. The stress distribution pattern has been examined and the stress transfer mode clarified. The Young＇s modulus and fracture strength of the Bis-GMA/TEGDMA composite were also predicted using the model with and without consideration of the damage. ine Iormer results with consideration of the debonding damage are in good agreement with existing literature experimental data. The shielding effect of our proposed model and an alternative approach were discussed. The FCC cell model has also been extended to predict the critical load for the damaged and the undamaged composite subject to the 3-point flexural test.