三维全五向编织复合材料弹性性能及热物理性能的有限元分析

在三维全五向(Q5D)编织复合材料细观结构模型的基础上, 建立了其单胞参数化有限元模型。通过施加合理的边界条件, 计算得到了Q5D编织复合材料的弹性常数、 热传导系数和热膨胀系数, 所得结果与现有的实验数据吻合较好。在此基础上, 深入研究了纤维体积分数、 编织角等工艺参数对材料弹性性能和热物理性能的影响规律, 并将计算结果与三维四向(4D)和三维五向(5D)编织复合材料的相应结果进行了对比。结果表明, Q5D编织复合材料具有较好的力学性能和纵向导热性能, 其零膨胀结构的可设计性更强, 为进一步研究此种结构材料的强度问题和热力耦合问题奠定了基础。     

三维全五向编织复合材料弹性性能及热物理性能的有限元分析

在三维全五向(Q5D)编织复合材料细观结构模型的基础上,建立了其单胞参数化有限元模型.通过施加合理的边界条件,计算得到了Q5D编织复合材料的弹性常数、热传导系数和热膨胀系数,所得结果与现有的实验数据吻合较好.在此基础上,深入研究了纤维体积分数、编织角等工艺参数对材料弹性性能和热物理性能的影响规律,并将计算结果与三维四向(4D)和三维五向(5D)编织复合材料的相应结果进行了对比.结果表明,Q5D编织复合材料具有较好的力学性能和纵向导热性能,其零膨胀结构的可设计性更强,为进一步研究此种结构材料的强度问题和热力耦合问题奠定了基础.     

三维五向编织复合材料中纱线截面形状实验分析

采用CCD显微摄像仪获取了树脂基三维五向编织复合材料的一系列截面图像, 分析了编织复合材料内纱线的排列规律及其截面形状的变化。实验结果显示: 沿编织成型方向, 三维五向编织复合材料的不同横截面内纱线的排列方式不同, 且呈周期性变化; 编织纱线的截面形状近似为扁平的凸透镜形; 第5向不动纱的轴线基本保持伸直, 其截面形状沿其轴向近似为扇形和三角形相互过渡变化; 表面编织纱线的排列方式与内部编织纱线的排列有关, 其截面形状近似为椭圆形。     

Effects of braiding angle on modal experimental analysis of three-dimensional and five-directional braided composites

In this paper, the experimental modal analysis of three dimensional (3D) and five (5)-directional braided composite cantilever beams with different braiding angle was conducted systematically with a simple testing system. Preforms were made by four-step 1 × 1 square integral braiding method. The first three order modes of composite specimens were derived, which are different with braiding angles. The experimental result and mechanism were discussed. Experimental observations and analytical predictions show that the natural frequency of specimens decreased and the damping ratio of specimens increased when the braiding angle increased. Furthermore, specimens with smaller braiding angle will be valuable for the better anti-exiting property, and have an opposite effect on dissipation of vibration energy. The comparison between the specimens with different braiding angles reveals that braiding angle is a crucial factor for the vibration performance design of 3D and 5-directional braided composites.     

Progressive failure analysis of bolted single-lap composite joint based on extended finite element method

In this paper, extend finite element method (XFEM) is used to predict the failure of single-lap bolted joints. To simplify calculation of XFEM model, composite laminates of joint are modeled using linear elastic properties. Three-dimensional equivalent material properties are calculated by the MATLAB code written by author. Progressive failure of bolted single-lap composite joint is investigated, and the failure load of joint simulated by XFEM is compared with experiments in literature. Then the influences of geometric parameters on failure load of one bolt single-lap composite joint are studied. Two geometric parameters include plate width-to-hole diameter ratio (W/D) and the edge-to-hole diameter ratio (E/D). At last the failure of single-lap joints with one bolt and two bolts are compared.     

A three-dimensional multilayer composite finite element for stress analysis of composite laminates

A three-dimensional multilayer composite finite element method has been developed based on a composite variational functional which takes three in-plane strains ε_x, ε_x, ε_xy and three transverse stresses σ_z, σ_yz, σ_xz as the basic variables. The continuity of the transverse stresses σ_z, σ_yz, σ_xz across the laminate thickness is assured a priori by introducing a partial stress field parameter α which is associated with the lower and upper surfaces of a lamina in a laminate. A method has been developed to form the partial stress field based on the assumed displacement field. With this method, a three dimensional (3-D) multilayer composite finite element is formulated for stress analysis of composite laminates. A numerical example is given, which shows some advantages of this composite element.     

BOLJAT: a tool for designing composite bolted joints using three-dimensional finite element analysis

This paper discusses the development of a software tool for design of composite bolted joints, using three-dimensional finite element analysis. The tool allows the user to create the joint geometry through a menu-driven interface and then generate a customised mesh according to the user's needs. Contact parameters are defined automatically, which shields the user from the most difficult part of the process. Boundary conditions, bolt pre-loads, and material properties can also be set. Only a few manual steps are necessary to complete the finite element code generation process. By automating the time-consuming model creation process, the tool facilitates the increased use of three-dimensional finite element analysis in the design of composite bolted joints. A case study is shown to demonstrate the usefulness of the tool.     

Mechanical analysis of 3D braided composites: a finite element model

The analysis of 3D braided composites is more difficult due to its complex microstructure. A new type of finite element method is developed to predict the effective moduli and the local stress within 3D braided composites under the 3D mechanical loading. To verify the present method, the material properties of undamaged 3D braided composites predicted in this paper are compared with the previous work. To demonstrate this method, some examples are analyzed.     

Multi-scale finite element analysis for tension and ballistic penetration damage characterizations of 2D triaxially braided composite

The multi-scale finite element model is presented to analyze tension and ballistic penetration damage characterizations of 2D triaxially braided composite (2DTBC). At the mesoscopic level, the damage of fiber tows is initiated with 3D Hashin criteria, and the damage initiation of pure matrix is predicted by the modified von Mises. The progressive degradation scheme and energy dissipation method are adopted to capture softening behaviors of tow and matrix. The macro-scale damage model is established by maximum-stress criteria and exponential damage evolution. To simulate interface debonding and inter-ply delamination, a triangle traction–separation law is adopted in each scale. Both scale damage models are verified with available experimental results. Based on numerical predictions, the stress–strain responses and damage developments of 2DTBC under axial and transverse tension loading are studied. For ballistic penetration loading, the meso-scale damage mechanisms of 2DTBC are predicted using 1/4 model, 1/2 model, 1-layer model, 2-layer model and 3-layer model. Then, effects of different model and impactor radius on damage modes are analyzed. Additionally, the macro-scale ballistic penetration behaviors of 2DTBC are simulated and compared with experiment. The prediction results for tension and penetration correlate well with experiment results. Both tension and penetration damage characterizations for tow, matrix within tow, pure matrix, interface and inter-ply delamination are revealed. A comparison of penetration damage between meso- and macro-scale presents a similar crack mechanism between two scales.     

三维五向编织复合材料渐进损伤分析的数值方法

基于连续损伤理论, 推导了含损伤裂纹的正交各向异性材料的应力-应变关系。建立了考虑界面脱粘破坏的三维五向编织复合材料细观体胞模型, 并将有限元网格尺寸和单元裂纹尺寸引入损伤演化方程。采用Hashin准则和von-Mises准则分别判断纱线与基体的初始损伤, 结合Eshelby-Mori-Tanaka方法确定材料的刚度退化系数。利用ANSYS有限元软件对材料进行渐进损伤分析, 得到了材料在单向拉伸作用下的应力-应变曲线和极限强度。计算结果表明, 轴纱为材料的主要承力部分, 小编织角材料的破坏模式主要为纱线的拉伸断裂, 界面破坏情况较为严重。模拟计算结果与实验结果吻合较好。     

Comparison between a finite element and a composite method for a three-dimensional potential problem

This paper is concerned with the solution of the three-dimensional potential problem for electromagnetic river gauging. It extends previous ideas of joining finite elements in an interior region to one infinite external element treated by the boundary integral method^1,2 to this case where there are two external infinite elements representing the river and the ground. The boundary continuity conditions on the infinite river–ground interface, as well as the internal–external interfaces, are dealt with by introducing a variational principle with relaxed continuity requirement³.     

Experimental and three-dimensional global-local finite element analysis of a composite component including degradation process at the interfaces

A specimen is studied in order to analyse the stringer runout effect in a specific design. An experimental test was carried out in several stages until specimen collapse, allowing the damaged regions and their characteristic geometric dimensions to be identified by ultrasound inspections. Several numerical analyses based on finite element method of the component were performed using the commercial software ABAQUS. These numerical analyses aimed to study in detail the onset and propagation of damage in the skin-stringer joint during the load procedure. Due to the high computational cost that this kind of procedure usually requires, a global–local approach was performed, comparing the two techniques offered by the software ABAQUS: Submodeling and Shell-to-Solid Coupling. Shell element typology was used in both cases for the global model. In order to estimate the extension of the critical zone in the specimen, the adhesive layer was modeled by means of solid elements. With reference to the local analysis, laminate continuum elements were selected to discretize the skin and the stringer, whereas the adhesive was represented by interface cohesive elements to take the damage evolution into account. A final correlation was carried out between numerical predictions and experimental results of the damaged regions in the component at several load levels.     

Mesh design in finite element analysis of post-buckled delamination in composite laminates

The role of mesh design in the post-buckling analysis of delamination in composite laminates is addressed in this paper. The determination of the strain energy release rate (SERR) along the crack front is central to the analysis. Frequently, theoretical analysis is limited to treatment of the problem in two dimensions, since considerable complexity is encountered in extending the analysis to three dimensions. However, many practical problems of embedded delamination in composite laminates are inherently three-dimensional in nature. Although in such cases, the finite element (FE) method can be employed, there are some issues that must be examined more closely to ensure physically realistic models. One of these issues is the effect of mesh design on the determination of the local SERR along the delamination front. There are few studies that deal with this aspect systematically. In this paper, the effect of mesh design in the calculation of SERR in two-dimensional (2D) and three-dimensional (3D) FE analyses of the post-buckling behavior of embedded delaminations is studied and some guidelines on mesh design are suggested. Two methods of calculation of the SERR are considered: the virtual crack closure technique (VCCT) and crack closure technique (CCT). The 2D analyses confirm that if the near-tip mesh is symmetric and consists of square elements, then the evaluation of the SERR is not sensitive to mesh refinement, and a reasonably coarse mesh is adequate. Despite agreement in the global post-buckling response of the delaminated part, the SERR calculated using different unsymmetrical near-tip meshes could be different. Therefore, unsymmetrical near-tip meshes should be avoided, as convergence of the SERR with mesh refinement could not be assured. While the results using VCCT and CCT for 2D analyses agree well with each other, these techniques yield different quantitative results when applied to 3D analyses. The reason may be due to the way in which the delamination growth is modeled. The CCT allows simultaneous delamination advance over finite circumferential lengths, but it is very difficult to implement and the results exhibit mesh dependency. Qualitatively, however, the two sets of results show similar distributions of Mode I and Mode II components of the SERR. This is fortunate, since the VCCT is relatively easy to implement.     

基于接触的水下夹桩器卡环连接结构设计

水下夹桩器是一种广泛应用于海洋石油平台建设的液压夹具。卡环式连接结构被应用于水下夹桩器。在水下夹桩器卡环连接结构设计中引入基于多体接触的有限元优化设计过程。阐述了基于扩张的拉格朗日算法有限元分析的接触优化模型建立,利用ANSYS软件参数化设计语言完成了优化计算。优化结果显示,采用非线性有限元法和优化方法对卡环连接结构进行设计能够取得满意的效果,设计结果满足工程要求。     

基于有限元分析的高精密柔性夹具优化方案设计

针对高精密加工过程中柔性夹具的变形问题进行理论计算和实验分析.首先应用CAE软件对现有柔性夹具进行有限元分析,分析夹具变形问题.然后通过CAD软件改变该夹具模型的结构参数和尺寸参数,优化夹具构造,以减少定位过程的夹具变形,并再次通过CAE软件进行分析,验证优化效果.通过增加夹具体主平面刚度、增加夹具体肋板刚度和改变夹具定位点位置等方法,实现夹具体的最大变形量满足加工精度的要求.实验分析以定位点变形为研究对象,测量得到4个定位点在夹具装夹过程中的变形量.实验结果基本符合理论计算结果,得到最终优化方案的4个定位点变形量均小于001 mm,满足了该高精密加工中对夹具体的精度要求.优化方案在实际生产中得到实施,效果良好.     

三维四向编织复合材料宏观有限元模型冲击损伤仿真及试验验证

三维编织复合材料作为整体编织材料,能够克服层合复合材料层间强度低、易分层的缺陷,相对于金属材料,还具有质量轻、高比刚、高比强度以及良好的抗冲击性能、较高的损伤容限,在汽车、高铁、航海、航空及航天领域中具有广泛应用前景。使用空气炮发射系统开展了钢珠以约210 m/s速度冲击三维四向编织复合材料平板的不同位置试验,基于宏观观察和细观观测,分析了三维编织复合材料在受到钢珠高速冲击下的破坏模式和破坏机制。此外,本文建立了三维四向编织复合材料宏观连续介质损伤（CDM）有限元模型,其中模型计算剩余速度和试验测得剩余速度误差在5%以内,试验和数值模拟的破坏形貌也高度一致,验证了所建立的宏观CDM有限元模型的有效性。     

自适应降噪麦克风壳体声学有限元分析与设计

基于DSP技术的自适应主动降噪麦克风是近年来适应超高噪声下通讯应用而发展起来的麦克风技术。该技术依靠两个独立的无指向型传声元件实现降噪。这两个传声元件通常背靠背地封装在一起。根据自适应滤波理论,系统的降噪性能与这两个元件间对应背景噪声的空间相干函数及对应近场语音的传递函数有关。这两个函数又与麦克风封装壳体的几何设计直接相关。为优化降噪性能,建立了一个以有限元声学仿真为基础的,考虑声场特性及麦克风壳体几何形状的主动降噪性能的模型。依据此模型,提出了新型附加反射板的壳体设计方案。仿真和实验结果显示,相对传统壳体设计,有限元优化后的新方案可以提高麦克风降噪性能大于15dB,提高输出信噪比大于10dB。     

Hyperelastic analysis based on a polygonal finite element method

In this contribution, we present a novel polygonal finite element method applied to hyperelastic analysis. For generating polygonal meshes in a bounded period of time, we use the adaptive Delaunay tessellation (ADT) proposed by Constantinu et al. [1 A. Constantinu, P. Steinmann, T. Bobach, G. Farin, and G. Umlauf, The adaptive Delaunay tessellation: A neighborhood covering meshing technique., Comput. Mech., vol. 42, no. 5, pp. 655–669, 2008.[Crossref], [Web of Science ®] , [Google Scholar]]. ADT is an unstructured hybrid tessellation of a scattered point set that minimally covers the proximal space around each point. In this work, we have extended the ADT to nonconvex domains using concepts from constrained Delaunay triangulation (CDT). The proposed method is thus based on a constrained adaptive Delaunay tessellation (CADT) for the discretization of domains into polygonal regions. We involve the metric coordinate (Malsch) method for obtaining the interpolation over convex and nonconvex domains. For the numerical integration of the Galerkin weak form, we resort to classical Gaussian quadrature based on triangles. Numerical examples of two-dimensional hyperelasticity are considered to demonstrate the advantages of the polygonal finite element method.     

Elasto-plastic finite element analysis with fuzzy parameters

In this paper, an elastoplastic analysis based on fuzzy mathematics and using finite element method will be described. The Drucker–Prager yield criterion and the elasto-plastic matrix are fuzzified for the non-linear analyses which adopt the initial stress method for the solution. A numerical example is given to illustrate the possible variations in the displacements and the extent of plastic zones at the discrete membership functions. A reliability index of a membership grade based on the concept of non-probabilistic entropy is also proposed.     

An isoparametric joint/interface element for finite element analysis

A generally applicable and simple joint/interface element for three- and two-dimensional finite element analysis is presented. The proposed element can model joints/interfaces between solid finite elements and shell finite elements. The derivation of the joint element stiffness is presented and algorithms for the treatment of nonlinear joint behaviour discussed. The performance of the element is tested on typical problems involving shell-to-shell and shell-to-solid interfaces.