CN102622472A - Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins - Google Patents

Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins Download PDF

Info

Publication number
CN102622472A
CN102622472A CN2012100460222A CN201210046022A CN102622472A CN 102622472 A CN102622472 A CN 102622472A CN 2012100460222 A CN2012100460222 A CN 2012100460222A CN 201210046022 A CN201210046022 A CN 201210046022A CN 102622472 A CN102622472 A CN 102622472A
Authority
CN
China
Prior art keywords
nail
model
stress
load
geometric
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012100460222A
Other languages
Chinese (zh)
Other versions
CN102622472B (en
Inventor
李彪
李亚智
姚振华
刘向东
李玺
刘兴科
舒怀
束一秀
谢强
王启
王海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN 201210046022 priority Critical patent/CN102622472B/en
Publication of CN102622472A publication Critical patent/CN102622472A/en
Application granted granted Critical
Publication of CN102622472B publication Critical patent/CN102622472B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及一种分析复合材料机械连接钉载和应力分布及应力水平的方法,技术特征在于:建立二维壳梁模型计算连接件钉载分配,并从二维模型中取出一个目标钉附近矩形区域的旁路载荷;在二维壳梁模型的基础上,根据矩形区域的大小及螺栓尺寸等几何信息,新建仅包含目标钉的三维实体模型,将二维模型取出的旁路载荷作为三维模型的力边界条件,并可施加螺栓预紧力,计算该三维模型即可得到该区域应力分布及应力水平。本方法发明充分发挥有限元模型中二维建模方法与三维建模方法的优点,为工程上计算多钉连接钉载分配、及应力分析提供快速而准确的计算方法,降低设计成本。

Figure 201210046022

The invention relates to a method for analyzing the nail load, stress distribution and stress level of the mechanical connection of composite materials. The technical feature is: establish a two-dimensional shell beam model to calculate the nail load distribution of the connector, and take out a rectangle near the target nail from the two-dimensional model Area bypass load; on the basis of the 2D shell beam model, according to the geometric information such as the size of the rectangular area and the bolt size, a new 3D solid model containing only target nails is created, and the bypass load extracted from the 2D model is used as the 3D model The force boundary condition can be applied, and the bolt preload can be applied, and the stress distribution and stress level of the area can be obtained by calculating the three-dimensional model. The invention of the method gives full play to the advantages of the two-dimensional modeling method and the three-dimensional modeling method in the finite element model, provides a fast and accurate calculation method for engineering calculation of multi-nail connection nail load distribution and stress analysis, and reduces design cost.

Figure 201210046022

Description

分析复合材料机械连接钉载和应力分布及应力水平的方法Method for analyzing nail loads and stress distributions and stress levels in mechanical connections of composite materials

技术领域 technical field

本发明涉及一种分析复合材料机械连接钉载和应力分布及应力水平的方法,可用于航空航天飞行器结构设计,机械设计与制造领域。  The invention relates to a method for analyzing the mechanical connection nail load, stress distribution and stress level of composite materials, which can be used in the fields of aerospace vehicle structure design, mechanical design and manufacturing. the

背景技术 Background technique

复合材料结构由于具有比强度大、比刚度高、抗腐蚀性能好以及优异的可设计性,近些年来在航空和航天领域得到越来越广泛的应用,已成为飞机结构的重要材料之一。复合材料用量通常用其所占飞机机体结构重量的百分比表示,并已成为飞机结构设计先进水平的重要指标之一。  Due to its high specific strength, high specific stiffness, good corrosion resistance and excellent designability, composite material structures have been more and more widely used in the field of aviation and aerospace in recent years, and have become one of the important materials for aircraft structures. The amount of composite materials is usually expressed as a percentage of the weight of the aircraft body structure, and has become one of the important indicators of the advanced level of aircraft structure design. the

复合材料主要承力结构都采用大量螺栓进行连接,该类机械连接结构是飞机结构中的主要连接构造和载荷传递途径。因此,复合材料机械连接接头成为飞机结构的关键部位,其连接力学性能、强度和寿命情况值得关注。复合材料强度、刚度的各向异性和缺乏延性,在多钉机械连接中连接板上各孔的载荷重新分配能力差,各个连接孔上的承载比例相差很大。因此,复合材料层合板多钉机械连接钉载分配问题在工程设计与应用中一直受到高度重视,目前有两大类方法可供应用:试验方法和有限元数值分析方法。试验方法耗资大,试件加工要求严格,难以对影响钉载分配和强度的众多因素进行系统的研究,而且试验结果带有一定的分散性;有限元数值分析方法耗资小,并可考虑不同参数的影响。  The main load-bearing structures of composite materials are connected by a large number of bolts. This type of mechanical connection structure is the main connection structure and load transmission path in aircraft structures. Therefore, the mechanical connection joint of composite material has become a key part of the aircraft structure, and its connection mechanical properties, strength and life are worthy of attention. Due to the anisotropy of strength and stiffness of composite materials and lack of ductility, the load redistribution ability of each hole in the connection plate is poor in the multi-nail mechanical connection, and the load ratio of each connection hole varies greatly. Therefore, the problem of nail load distribution in multi-nail mechanical connections of composite laminates has always been highly valued in engineering design and application. At present, there are two types of methods available for application: experimental methods and finite element numerical analysis methods. The cost of the test method is high, the processing requirements of the specimen are strict, and it is difficult to conduct systematic research on many factors that affect the nail load distribution and strength, and the test results have a certain degree of dispersion; the finite element numerical analysis method is low in cost and can consider different parameters. Impact. the

确定多钉连接结构中各个钉的承载比例、应力分布及应力水平是计算多钉连接强度的基础。国内外关于复合材料机械连接问题广泛采用的有限元方法包括以下两种:一种是建立二维板壳模型的方法;另一种是建立三维实体模型的方法。单纯的二维板壳模型方法是一种简化的模拟方法,计算周期短、能够得到满足工程精度的钉载分配结果;但得到的层合板中各个铺层的应力分布不能考虑实际情况中的螺栓预紧力、钉孔装配间隙、摩擦等因素的影响,且不能获得层间应力。建立三维实体模型的方法能 够考虑上述实际情况中的多种影响因素,但其缺点是建模过程复杂,计算周期过长,不具有通用性,接触问题非线性迭代计算的收敛性问题常常成为难点,不适于工程实际运用。  Determining the bearing ratio, stress distribution and stress level of each nail in the multi-nail connection structure is the basis for calculating the strength of the multi-nail connection. The finite element methods widely used at home and abroad on the mechanical connection of composite materials include the following two types: one is the method of establishing a two-dimensional plate and shell model; the other is the method of establishing a three-dimensional solid model. The simple two-dimensional plate and shell model method is a simplified simulation method, which has a short calculation period and can obtain nail load distribution results that meet engineering accuracy; Preload, nail hole assembly clearance, friction and other factors, and the interlayer stress cannot be obtained. The method of establishing a three-dimensional solid model can take into account various influencing factors in the above-mentioned actual situation, but its disadvantages are that the modeling process is complicated, the calculation cycle is too long, and it is not universal. The convergence problem of nonlinear iterative calculation of contact problems often becomes a problem. Difficulties, not suitable for practical engineering applications. the

发明内容 Contents of the invention

要解决的技术问题  technical problem to be solved

为了避免现有技术的不足之处,本发明提出一种分析复合材料机械连接钉载和应力分布及应力水平的方法,发挥有限元模型中二维建模方法与三维建模方法的优点,为工程上计算多钉连接钉载分配、及应力分析提供快速而准确的计算方法,优化设计多钉连接结构,降低设计成本。  In order to avoid the deficiencies of the prior art, the present invention proposes a method for analyzing the nail load, stress distribution and stress level of the mechanical connection of composite materials, and utilizes the advantages of the two-dimensional modeling method and the three-dimensional modeling method in the finite element model, for In engineering, the calculation of multi-nail connection nail load distribution and stress analysis provide fast and accurate calculation methods, optimize the design of multi-nail connection structures, and reduce design costs. the

技术方案  Technical solutions

一种分析复合材料机械连接钉载的方法,其特征在于步骤如下:  A method for analyzing composite material mechanical connection nail load, characterized in that the steps are as follows:

步骤1:根据各个被连接板的实际长度和宽度,在各自的中面位置处建立对应的几何面,在钉位置处建立几何线段;  Step 1: According to the actual length and width of each connected plate, establish the corresponding geometric plane at the respective mid-plane position, and establish the geometric line segment at the nail position; 

步骤2:在几何面上,分别划分出以各个钉为中心,排距为边长的正方形区域;  Step 2: On the geometric surface, divide the square area with each nail as the center and the row distance as the side length;

步骤3:将几何面赋予壳单元属性,并划分成四边形网格;将几何线段赋予梁元属性,并划分成线网格,且梁元的两个端点在被连接板的网格节点上,正方形边上将存在节点;完成网格划分,得到二维有限元板壳模型;  Step 3: Assign the geometric surface to the shell element attribute and divide it into a quadrilateral mesh; assign the geometric line segment to the beam element attribute and divide it into a line mesh, and the two end points of the beam element are on the grid nodes of the connected plate, There will be nodes on the side of the square; the grid division is completed to obtain a two-dimensional finite element plate and shell model;

步骤4:在二维模型中,将连接件的固支端进行固支约束,加载端施加所需的载荷,载荷在不同被连接板之间通过梁元来传递;  Step 4: In the two-dimensional model, the fixed end of the connector is restrained by the fixed support, and the required load is applied to the loading end, and the load is transmitted between different connected plates through beam elements;

步骤5:通过有限元软件求解计算模型,当计算完成后,在得到的计算结果中提取出梁元的节点力,该节点力大小为钉载。  Step 5: Solve the calculation model through the finite element software. After the calculation is completed, extract the nodal force of the beam element from the obtained calculation results, and the magnitude of the nodal force is the nail load. the

一种分析复合材料机械连接钉载的周围应力分布及应力水平的方法,其特征在于步骤如下:  A method for analyzing the surrounding stress distribution and stress level of the mechanical connection nail load of composite materials, characterized in that the steps are as follows:

步骤1:根据各个被连接板的实际长度和宽度,在各自的中面位置处建立对应的 几何面,在钉位置处建立几何线段;  Step 1: According to the actual length and width of each connected plate, establish the corresponding geometric plane at the position of the respective middle plane, and establish the geometric line segment at the position of the nail;

步骤2:在几何面上,分别划分出以各个钉为中心,排距为边长的正方形区域;  Step 2: On the geometric surface, divide the square area with each nail as the center and the row distance as the side length;

步骤3:将几何面赋予壳单元属性,并划分成四边形网格;将几何线段赋予梁元属性,并划分成线网格,且梁元的两个端点在被连接板的网格节点上,正方形边上将存在节点;完成网格划分,得到二维有限元板壳模型;  Step 3: Assign the geometric surface to the shell element attribute and divide it into a quadrilateral mesh; assign the geometric line segment to the beam element attribute and divide it into a line mesh, and the two end points of the beam element are on the grid nodes of the connected plate, There will be nodes on the side of the square; the grid division is completed to obtain a two-dimensional finite element plate and shell model;

步骤4:在二维模型中,在任意一个钉为目标钉的被连接板上的正方形区域上,分别提取并保存正方形区域四周的旁路载荷,所述旁路载荷包括力和弯矩;  Step 4: In the two-dimensional model, on any square area on the connected plate where any nail is the target nail, extract and save the bypass loads around the square area, and the bypass load includes force and bending moment;

步骤5:根据步骤2中正方形区域的边长和各个被连接板的真实厚度,建立被连接板的三维几何模型;选择其中的任意一块被连接板,将它在轴线方向的一端延长至模型长度的10%,为弹性假体;  Step 5: According to the side length of the square area in step 2 and the real thickness of each connected plate, establish a 3D geometric model of the connected plate; select any one of the connected plates, and extend one end of it in the axial direction to the length of the model 10% of the elastic prosthesis;

步骤6:根据钉和垫圈的尺寸建立紧固件的三维几何模型,将钉与垫圈当成整体建模,并以圆柱体表示螺纹形状;  Step 6: Establish a three-dimensional geometric model of the fastener according to the size of the nail and washer, model the nail and washer as a whole, and use a cylinder to represent the thread shape;

步骤7:将步骤5和步骤6得到的三维几何模型赋予实体单元属性,扫略生成六面体网格;在钉与钉孔之间、垫圈与被连接板之间、不同的被连接板之间建立接触对,完成三维实体模型的建立;  Step 7: Assign the 3D geometric model obtained in Step 5 and Step 6 to the attribute of the solid unit, and scan to generate a hexahedral mesh; establish between nails and nail holes, between gaskets and connected plates, and between different connected plates Contact pairs to complete the establishment of the 3D solid model;

步骤8:以步骤5中的弹性假体的一端作固支约束;  Step 8: Use one end of the elastic prosthesis in step 5 as a fixed support constraint;

步骤9:以步骤4中取出的旁路载荷作为三维实体模型对应面上的力边界条件,且弹性假体的侧边不加载;力按照均布力的方式加载,弯矩按照弯曲正应力公式折算成应力在面上梯度加载,并在钉上通过预紧力单元施加螺栓预紧力;  Step 9: Use the bypass load taken out in step 4 as the force boundary condition on the corresponding surface of the 3D solid model, and the side of the elastic prosthesis is not loaded; the force is loaded in a uniformly distributed manner, and the bending moment is in accordance with the bending normal stress formula It is converted into stress and loaded gradiently on the surface, and the bolt pretightening force is applied on the nail through the pretightening force unit;

步骤10:利用有限元软件求解计算该三维模型,得到目标钉周围应力分布及应力水平。  Step 10: use finite element software to solve and calculate the three-dimensional model, and obtain the stress distribution and stress level around the target nail. the

依次将所有钉作为目标钉,并反复重复步骤4~步骤10得到所有钉周围的应力分布及应力水平。  Take all the nails as the target nails in turn, and repeat steps 4 to 10 repeatedly to obtain the stress distribution and stress level around all nails. the

有益效果  Beneficial effect

本发明提出的一种分析复合材料机械连接钉载和应力分布及应力水平的方法,建立二维壳梁模型计算连接件钉载分配,并从二维模型中取出一个目标钉附近矩形区域的旁路载荷;在二维壳梁模型的基础上,根据矩形区域的大小及螺栓尺寸等几何信息,新建仅包含目标钉的三维实体模型,将二维模型取出的旁路载荷作为三维模型的力边界条件,并可施加螺栓预紧力,计算该三维模型即可得到该区域应力分布及应力水平。  A method for analyzing the nail load, stress distribution and stress level of the mechanical connection of composite materials proposed by the present invention, establishes a two-dimensional shell beam model to calculate the nail load distribution of the connector, and takes out a side of a rectangular area near the target nail from the two-dimensional model road load; on the basis of the 2D shell-beam model, according to the geometric information such as the size of the rectangular area and the bolt size, a new 3D solid model containing only target nails is created, and the bypass load extracted from the 2D model is used as the force boundary of the 3D model Conditions, and the bolt pre-tightening force can be applied, and the stress distribution and stress level of the area can be obtained by calculating the three-dimensional model. the

本发明与现有技术相比有益效果:充分发挥二维与三维有限元模型的优点,应用性强,计算时间短,并能保证计算精度,提高了计算效率。本发明方法中,有限元分析从二维模型到三维模型结合的方法来分析问题的思想较新颖,对其他大型复杂模型都可以使用。本发明方法对复合材料板与复合材料板连接,金属板与金属板连接,复合材料板与金属板混合连接均适用,通用性强。  Compared with the prior art, the present invention has beneficial effects: the advantages of the two-dimensional and three-dimensional finite element models are fully utilized, the applicability is strong, the calculation time is short, the calculation accuracy can be ensured, and the calculation efficiency is improved. In the method of the invention, the idea of combining the finite element analysis from the two-dimensional model to the three-dimensional model to analyze the problem is relatively novel, and can be used for other large and complex models. The method of the invention is applicable to the connection of composite material boards to composite material boards, metal boards to metal boards, and composite material boards to metal boards, and has strong universality. the

附图说明 Description of drawings

图1:本发明技术方案流程图;  Fig. 1: flow chart of technical scheme of the present invention;

图2:一列三钉单搭接连接件的二维板壳模型;  Figure 2: A two-dimensional plate and shell model of a row of three-nail single-lap joints;

1-复合材料板,2-金属板,I-第一个钉,II-第二个钉,III-第三一个钉;  1-composite plate, 2-metal plate, I-the first nail, II-the second nail, III-the third nail;

图3:二维板壳模型中壳单元与梁元共用节点的示意图;  Figure 3: Schematic diagram of joints shared by shell elements and beam elements in the 2D plate and shell model;

图4:三维实体模型网格形式及其约束方式示意图;  Figure 4: Schematic diagram of the grid form of the 3D solid model and its constraints;

图5:三维实体模型中不同类型板旁路载荷的加载内容及方式示意图;  Figure 5: Schematic diagram of the loading contents and methods of different types of plate bypass loads in the 3D solid model;

a:固支板加载内容及方式;b:非固支板加载内容及方式。  a: Loading content and method of fixed support plate; b: Loading content and method of non-fixed support plate. the

具体实施方式 Detailed ways

现结合实施例、附图对本发明作进一步描述:  Now in conjunction with embodiment, accompanying drawing, the present invention will be further described:

本发明的计算流程如图1所示,为了更好的理解本发明的技术方案,将以上算法应用于常见的复合材料板与金属板机械连接分析,分析内容为:钉载分配计算,及最大钉载孔周围的应力场分析。连接件是一列三钉的单搭接板,如图2所示,以下步骤基于有限元软件ANSYS:  The calculation process of the present invention is shown in Figure 1. In order to better understand the technical solution of the present invention, the above algorithm is applied to the common mechanical connection analysis of composite material plates and metal plates. The analysis content is: nail load distribution calculation, and maximum Stress field analysis around the nail loading hole. The connector is a single lapped plate with three nails, as shown in Figure 2. The following steps are based on the finite element software ANSYS:

步骤1根据不同的被连接板的实际长度和宽度,在各自的中面位置处建立几何面,将工作平面的原点移动到钉所在位置,转动工作平面,分别沿X方向和Y方向在两个几何面上切割出关键点。把不同板上同一个钉的关键点用线段连接起来,记录保存关键点编号。复合材料板是图2a中的板1,金属板是图2a中的板2;  Step 1. According to the actual length and width of different plates to be connected, establish a geometric surface at the position of the respective middle plane, move the origin of the working plane to the position of the nail, rotate the working plane, and move the plane along the X direction and the Y direction at two Keypoints are cut out of the geometry face. Connect the key points of the same nail on different boards with line segments, and record and save the key point numbers. The composite plate is plate 1 in Fig. 2a, and the metal plate is plate 2 in Fig. 2a;

步骤2以钉为正方形中心,排距为边长,用工作平面在各钉周围划分正方形区域;如图2b所示;  Step 2. Take the nail as the center of the square, and the row distance as the side length, and use the working plane to divide the square area around each nail; as shown in Figure 2b;

步骤3定义材料属性,并对复合材料板赋予SHELL99壳单元属性,金属板赋予SHELL63壳单元属性,将几何面划分生成四边形网格。线段赋予BEAM188梁元属性,并划分生成一个线网格。梁元的两个端点在被连接板的网格节点上,即梁元与壳单元是共用节点的(如图3所示);  Step 3 defines the material properties, and assigns SHELL99 shell element attributes to the composite material plate, and SHELL63 shell element attributes to the metal plate, and divides the geometric surface to generate a quadrilateral mesh. The line segment is assigned the beam element attribute of BEAM188, and is divided to generate a line mesh. The two end points of the beam element are on the grid node of the connected plate, that is, the beam element and the shell element share the node (as shown in Figure 3);

步骤4对连接件施加图2a中的边界条件;  Step 4 applies the boundary conditions in Figure 2a to the connector;

步骤5在被连接板的加载端施加所需载荷,在Solution模块中计算该二维板壳模型;  Step 5 Apply the required load on the loading end of the connected plate, and calculate the two-dimensional plate and shell model in the Solution module;

步骤6计算完成后,保存该模型为2D.db格式的数据文件。在后处理模块分别取出各个公共节点的节点力大小,该值就是钉载;  After the calculation in step 6 is completed, save the model as a data file in 2D.db format. The nodal force of each common node is taken out in the post-processing module, and this value is the nail load;

步骤7将钉载分配比例最大的I钉作为目标钉。在二维板壳模型中,选取出每块被连接板上目标钉附近的正方形区域,即选取出图2b中的阴影区域;  In step 7, the I nail with the largest nail load distribution ratio is taken as the target nail. In the two-dimensional plate and shell model, select the square area near the target nail on each connected plate, that is, select the shaded area in Figure 2b;

步骤8提取正方形区域四周的旁路载荷。先选中阴影区域内的单元,再选择正方形边上的节点,在命令框内输入FSUM命令,可以得到选中单元对这些节点的作用力;用该命令分别提取出每个正方形四条边上的作用力。垂直于X轴的正方形边提取的旁路载荷包括NX,QY,QZ,MY,MZ;平行于X轴的正方形边提取的旁路载荷包括NY,QX,QZ,MX,MZ;将提取出的旁路载荷存入数据文件之中;  Step 8 extracts the bypass loads around the square area. First select the elements in the shaded area, then select the nodes on the sides of the square, and enter the FSUM command in the command box to get the force of the selected elements on these nodes; use this command to extract the forces on the four sides of each square . The bypass loads extracted from the square sides perpendicular to the X axis include N X , Q Y , Q Z , M Y , M Z ; the side loads extracted from the square sides parallel to the X axis include N Y , Q X , Q Z , M X , M Z ; store the extracted bypass load into the data file;

步骤9新建作业,并根据二维板壳模型中正方形的尺寸、不同被连接板各自的厚度建立被连接板的几何体。根据钉和垫圈的实际尺寸,建立紧固件的几何体,并忽 略螺纹形状,将钉与垫圈当做整体。选择金属板方块,将其沿X方向的一端延长模型长度的10%作为弹性假体;弹性假体的目的是消除固支这个强约束对应力分布造成的不利影响;  Step 9 Create a new job, and establish the geometry of the connected plates according to the size of the square in the two-dimensional plate and shell model and the respective thicknesses of different connected plates. According to the actual size of the nail and washer, the geometry of the fastener is established, and the thread shape is ignored, and the nail and washer are regarded as a whole. Select a metal plate square, and extend its end along the X direction by 10% of the model length as an elastic prosthesis; the purpose of the elastic prosthesis is to eliminate the adverse effects of the strong constraint of the fixation on the stress distribution;

步骤10对上述三维几何模型赋予SOLID185实体单元属性,并扫略生成六面体网格。在钉与钉孔之间、垫圈与被连接板之间、不同的被连接板之间分别定义接触对,生成的三维实体模型,模型网格如图4所示;  Step 10 assigns SOLID185 solid element attributes to the above-mentioned 3D geometric model, and scans to generate a hexahedral mesh. Define contact pairs between nails and nail holes, between gaskets and connected plates, and between different connected plates, and generate a three-dimensional solid model. The model grid is shown in Figure 4;

步骤11将步骤8中取出的旁路载荷作为应力边界条件分别加到三维实体模型对应的面上。其中,力按均布加载的方式:垂直于面的力(NX或NY)直接加于单元面上,剪力用表面效应单元SURF154加载;弯矩按弯曲正应力公式折算成正应力,以梯度的形式加到对应面上;预紧力使用ANSYS内部提供的预紧力单元施加。三维实体模型中弹性假体部分的侧边不包含应力边界条件,如在图5a中,仅在阴影部分施加应力边界条件;  In step 11, the bypass load taken out in step 8 is added as a stress boundary condition to the corresponding surface of the three-dimensional solid model. Among them, the force is loaded in a uniform manner: the force perpendicular to the surface (N X or N Y ) is directly applied to the element surface, and the shear force is loaded by the surface effect unit SURF154; the bending moment is converted into a normal stress according to the bending normal stress formula, as The form of the gradient is added to the corresponding surface; the preload is applied using the preload element provided inside ANSYS. The side of the elastic prosthesis in the three-dimensional solid model does not contain stress boundary conditions, as shown in Figure 5a, only stress boundary conditions are applied to the shaded part;

步骤12在Solution模块中计算该三维实体模型,当计算完成后,可以在后处理模块中提取出各个铺层在各自单元坐标系下的结果,并绘制应力云图,亦可查看整体变形图。  Step 12 Calculate the 3D solid model in the Solution module. After the calculation is completed, the results of each layer in the respective unit coordinate system can be extracted in the post-processing module, and the stress cloud diagram can be drawn, and the overall deformation diagram can also be viewed. the

步骤13若要改变目标钉重新分析,则需恢复步骤6中保存的2D.db文件,再将步骤7到步骤12重复一遍即可。  Step 13 If you want to change the target nail and re-analyze, you need to restore the 2D.db file saved in step 6, and then repeat steps 7 to 12. the

Claims (3)

1.一种分析复合材料机械连接钉载的方法,其特征在于步骤如下:1. A method for analyzing composite material mechanical connection nail loading, characterized in that the steps are as follows: 步骤1:根据各个被连接板的实际长度和宽度,在各自的中面位置处建立对应的几何面,在钉位置处建立几何线段;Step 1: According to the actual length and width of each connected plate, establish the corresponding geometric surface at the position of the respective middle plane, and establish the geometric line segment at the position of the nail; 步骤2:在几何面上,分别划分出以各个钉为中心,排距为边长的正方形区域;Step 2: On the geometric surface, divide the square area with each nail as the center and the row distance as the side length; 步骤3:将几何面赋予壳单元属性,并划分成四边形网格;将几何线段赋予梁元属性,并划分成线网格,且梁元的两个端点在被连接板的网格节点上,正方形边上将存在节点;完成网格划分,得到二维有限元板壳模型;Step 3: Assign the geometric surface to the shell element attribute and divide it into a quadrilateral mesh; assign the geometric line segment to the beam element attribute and divide it into a line mesh, and the two end points of the beam element are on the grid nodes of the connected plate, There will be nodes on the side of the square; the grid division is completed to obtain a two-dimensional finite element plate and shell model; 步骤4:在二维模型中,将连接件的固支端进行固支约束,加载端施加所需的载荷,载荷在不同被连接板之间通过梁元来传递;Step 4: In the two-dimensional model, the fixed support end of the connector is fixedly supported, and the required load is applied to the loading end, and the load is transmitted between different connected plates through beam elements; 步骤5:通过有限元软件求解计算模型,当计算完成后,在得到的计算结果中提取出梁元的节点力,该节点力大小为钉载。Step 5: Solve the calculation model through the finite element software. After the calculation is completed, extract the nodal force of the beam element from the obtained calculation results, and the magnitude of the nodal force is the nail load. 2.一种利用权利要求1所述方法分析复合材料机械连接钉载的周围应力分布及应力水平的方法,其特征在于步骤如下:2. a method utilizing the described method of claim 1 to analyze the surrounding stress distribution and the stress level of the composite material mechanical connection nail load, it is characterized in that the steps are as follows: 步骤1:根据各个被连接板的实际长度和宽度,在各自的中面位置处建立对应的几何面,在钉位置处建立几何线段;Step 1: According to the actual length and width of each connected plate, establish the corresponding geometric surface at the position of the respective middle plane, and establish the geometric line segment at the position of the nail; 步骤2:在几何面上,分别划分出以各个钉为中心,排距为边长的正方形区域;Step 2: On the geometric surface, divide the square area with each nail as the center and the row distance as the side length; 步骤3:将几何面赋予壳单元属性,并划分成四边形网格;将几何线段赋予梁元属性,并划分成线网格,且梁元的两个端点在被连接板的网格节点上,正方形边上将存在节点;完成网格划分,得到二维有限元板壳模型;Step 3: Assign the geometric surface to the shell element attribute and divide it into a quadrilateral mesh; assign the geometric line segment to the beam element attribute and divide it into a line mesh, and the two end points of the beam element are on the grid nodes of the connected plate, There will be nodes on the side of the square; the grid division is completed to obtain a two-dimensional finite element plate and shell model; 步骤4:在二维模型中,在任意一个钉为目标钉的被连接板上的正方形区域上,分别提取并保存正方形区域四周的旁路载荷,所述旁路载荷包括力和弯矩;Step 4: In the two-dimensional model, on any square area on the connected plate where any nail is the target nail, extract and save the bypass loads around the square area, the bypass loads include force and bending moment; 步骤5:根据步骤2中正方形区域的边长和各个被连接板的真实厚度,建立被连接板的三维几何模型;选择其中的任意一块被连接板,将它在轴线方向的一端延长至模型长度的10%,为弹性假体;Step 5: According to the side length of the square area in step 2 and the real thickness of each connected plate, establish a 3D geometric model of the connected plate; select any one of the connected plates, and extend one end of it in the axial direction to the length of the model 10% of the elastic prosthesis; 步骤6:根据钉和垫圈的尺寸建立紧固件的三维几何模型,将钉与垫圈当成整体建模,并以圆柱体表示螺纹形状;Step 6: Establish a three-dimensional geometric model of the fastener according to the size of the nail and washer, model the nail and washer as a whole, and use a cylinder to represent the thread shape; 步骤7:将步骤5和步骤6得到的三维几何模型赋予实体单元属性,扫略生成六面体网格;在钉与钉孔之间、垫圈与被连接板之间、不同的被连接板之间建立接触对,完成三维实体模型的建立;Step 7: Assign the 3D geometric model obtained in Step 5 and Step 6 to the attribute of the solid unit, and scan to generate a hexahedral mesh; establish between nails and nail holes, between gaskets and connected plates, and between different connected plates Contact pairs to complete the establishment of the 3D solid model; 步骤8:以步骤5中的弹性假体的一端作固支约束;Step 8: Use one end of the elastic prosthesis in step 5 as a fixed support constraint; 步骤9:以步骤4中取出的旁路载荷作为三维实体模型对应面上的力边界条件,且弹性假体的侧边不加载;力按照均布力的方式加载,弯矩按照弯曲正应力公式折算成应力在面上梯度加载,并在钉上通过预紧力单元施加螺栓预紧力;Step 9: Use the bypass load taken out in step 4 as the force boundary condition on the corresponding surface of the 3D solid model, and the side of the elastic prosthesis is not loaded; the force is loaded in a uniformly distributed manner, and the bending moment is in accordance with the bending normal stress formula It is converted into stress and loaded on the surface gradient, and the bolt pretightening force is applied on the nail through the pretightening force unit; 步骤10:利用有限元软件求解计算该三维模型,得到目标钉周围应力分布及应力水平。Step 10: use finite element software to solve and calculate the three-dimensional model, and obtain the stress distribution and stress level around the target nail. 3.根据权利要求2所述的方法,其特征在于:在步骤4中改变目标钉,依次将所有钉作为目标钉,并反复重复步骤4~步骤10得到所有钉周围的应力分布及应力水平。3. The method according to claim 2, characterized in that in step 4, the target nails are changed, all nails are used as target nails in turn, and steps 4-10 are repeated repeatedly to obtain the stress distribution and stress level around all nails.
CN 201210046022 2012-02-27 2012-02-27 Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins Expired - Fee Related CN102622472B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210046022 CN102622472B (en) 2012-02-27 2012-02-27 Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210046022 CN102622472B (en) 2012-02-27 2012-02-27 Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins

Publications (2)

Publication Number Publication Date
CN102622472A true CN102622472A (en) 2012-08-01
CN102622472B CN102622472B (en) 2013-12-25

Family

ID=46562388

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210046022 Expired - Fee Related CN102622472B (en) 2012-02-27 2012-02-27 Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins

Country Status (1)

Country Link
CN (1) CN102622472B (en)

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102841969A (en) * 2012-09-19 2012-12-26 中国航天科工集团第二研究院二十三所 Finite element modeling method for shell and beam strengthening structure
CN103353360A (en) * 2012-11-14 2013-10-16 西安交通大学 Method for measuring residual stress of welded rotor
CN103455679A (en) * 2013-09-06 2013-12-18 湘潭电机股份有限公司 Finite element analysis method of composite material
CN103837361A (en) * 2014-03-06 2014-06-04 卿新林 Comprehensive monitor system for composite material mechanical connecting structure
CN103925267A (en) * 2014-03-19 2014-07-16 刘龙权 Composite material and metal connection structure and forming method thereof
CN104699872A (en) * 2013-12-06 2015-06-10 中国飞机强度研究所 Numerical analysis method for connection between three-dimensional and two-dimensional finite element units
CN105277661A (en) * 2015-11-02 2016-01-27 西北工业大学 Delamination damage analyzing method of composite material laminated plate interference bolt mounting process
CN105447237A (en) * 2015-11-16 2016-03-30 西北工业大学 Stress analysis method for composite material interference connection under action of tensile load
CN105588759A (en) * 2016-01-22 2016-05-18 北京航空航天大学 Test method for indirectly measuring bolt load distribution proportion in damage process of composite multi-bolt joint structure
CN106709183A (en) * 2016-12-26 2017-05-24 中国航空工业集团公司西安飞机设计研究所 Detail stress simulation method of composite mixed joint replaced by adopting metal plate
CN106768916A (en) * 2016-11-21 2017-05-31 北京航空航天大学 It is a kind of to use the method uniformly followed closely and be loaded into row compound material bolt attachment structure ultimate failure loading prediction
CN106874545A (en) * 2017-01-05 2017-06-20 武汉船用机械有限责任公司 A kind of method for designing of flat board with holes
CN107220422A (en) * 2017-05-19 2017-09-29 上海宇航系统工程研究所 A kind of extensive connector intensity all standing finite element modeling method
CN107526879A (en) * 2017-08-04 2017-12-29 华中科技大学 A kind of threedimensional model gradient finite element solving method based on ANSYS Parametric Design Languages
CN107885946A (en) * 2017-11-23 2018-04-06 江苏理工学院 A kind of method based on progressive damage model prediction composite Multi-fasteners joint intensity
CN109388868A (en) * 2018-09-25 2019-02-26 陕西飞机工业(集团)有限公司 A kind of contact of consideration nail hole and follow closely curved nail load calculation method
CN109614713A (en) * 2018-12-13 2019-04-12 北京航空航天大学 Prediction method of nail load distribution for composite multi-nail connections based on modified stiffness method
CN109682680A (en) * 2018-11-29 2019-04-26 中国商用飞机有限责任公司北京民用飞机技术研究中心 A kind of complex composite material mechanical connected nail load measurement test method
CN109726435A (en) * 2018-12-04 2019-05-07 中国航空工业集团公司西安飞机设计研究所 A Method for Calculating the Loads of Joint Nail Groups
CN109813587A (en) * 2019-03-08 2019-05-28 南昌航空大学 A test piece for the study of nail load distribution of rivet eccentric connection joints
CN111159879A (en) * 2019-12-25 2020-05-15 中国航空工业集团公司西安飞机设计研究所 Nail load processing method of centralized load transmission structure
WO2020107529A1 (en) * 2018-11-29 2020-06-04 中国商用飞机有限责任公司北京民用飞机技术研究中心 Test method for pin load distribution measurement for mechanical connection structure
CN112231864A (en) * 2020-11-16 2021-01-15 西北工业大学 A Rapid Analysis System for Composite Nail Loading
CN114741810A (en) * 2022-04-20 2022-07-12 中国船舶重工集团海装风电股份有限公司 High-precision bolt internal force extraction method and device based on finite element large deformation solution
CN115031942A (en) * 2022-05-27 2022-09-09 南京航空航天大学 Method for determining bolt tightening process parameters for assembling composite material structure
CN115098972A (en) * 2022-07-05 2022-09-23 浙江吉利控股集团有限公司 Multilayer structure bolted connection modeling method, device, equipment and storage medium
CN115828707A (en) * 2022-12-29 2023-03-21 中国航空工业集团公司成都飞机设计研究所 Method for calculating service life of fastening hole under complex loading condition

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1991846A (en) * 2005-12-29 2007-07-04 比亚迪股份有限公司 Method for creating three-dimensional finite element mesh
US20100131256A1 (en) * 2008-11-26 2010-05-27 Livermore Software Technology Corporation Spot Weld Failure Determination Method in a Finite Element Analysis
CN101877021A (en) * 2009-12-28 2010-11-03 中国农业大学 A method for automatic data conversion among dynamic analysis, three-dimensional modeling and finite element analysis software

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1991846A (en) * 2005-12-29 2007-07-04 比亚迪股份有限公司 Method for creating three-dimensional finite element mesh
US20100131256A1 (en) * 2008-11-26 2010-05-27 Livermore Software Technology Corporation Spot Weld Failure Determination Method in a Finite Element Analysis
CN101877021A (en) * 2009-12-28 2010-11-03 中国农业大学 A method for automatic data conversion among dynamic analysis, three-dimensional modeling and finite element analysis software

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘向东 等: "多钉连接钉载分配特性研究方法", 《复合材料学报》 *

Cited By (43)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102841969A (en) * 2012-09-19 2012-12-26 中国航天科工集团第二研究院二十三所 Finite element modeling method for shell and beam strengthening structure
CN103353360A (en) * 2012-11-14 2013-10-16 西安交通大学 Method for measuring residual stress of welded rotor
CN103353360B (en) * 2012-11-14 2015-07-01 西安交通大学 Method for measuring residual stress of welded rotor
CN103455679A (en) * 2013-09-06 2013-12-18 湘潭电机股份有限公司 Finite element analysis method of composite material
CN104699872A (en) * 2013-12-06 2015-06-10 中国飞机强度研究所 Numerical analysis method for connection between three-dimensional and two-dimensional finite element units
CN103837361A (en) * 2014-03-06 2014-06-04 卿新林 Comprehensive monitor system for composite material mechanical connecting structure
CN103837361B (en) * 2014-03-06 2016-03-23 常州飞智传感科技有限公司 Compound substance machinery syndeton comprehensive monitor system
CN103925267A (en) * 2014-03-19 2014-07-16 刘龙权 Composite material and metal connection structure and forming method thereof
CN103925267B (en) * 2014-03-19 2016-09-07 刘龙权 A kind of composite and metal connecting structure and the method forming described attachment structure
CN105277661A (en) * 2015-11-02 2016-01-27 西北工业大学 Delamination damage analyzing method of composite material laminated plate interference bolt mounting process
CN105277661B (en) * 2015-11-02 2017-03-22 西北工业大学 Delamination damage analyzing method of composite material laminated plate interference bolt mounting process
CN105447237A (en) * 2015-11-16 2016-03-30 西北工业大学 Stress analysis method for composite material interference connection under action of tensile load
CN105447237B (en) * 2015-11-16 2019-01-08 西北工业大学 Tensile load acts on lower composite material interference connection stress analysis method
CN105588759A (en) * 2016-01-22 2016-05-18 北京航空航天大学 Test method for indirectly measuring bolt load distribution proportion in damage process of composite multi-bolt joint structure
CN105588759B (en) * 2016-01-22 2018-02-23 北京航空航天大学 The test method of pin load distribution ratio in a kind of indirect determination composite Multi-fasteners joint structure failure history
CN106768916B (en) * 2016-11-21 2018-11-20 北京航空航天大学 A method of row compound material bolt connection structure ultimate failure loading prediction is loaded into using uniformly nail
CN106768916A (en) * 2016-11-21 2017-05-31 北京航空航天大学 It is a kind of to use the method uniformly followed closely and be loaded into row compound material bolt attachment structure ultimate failure loading prediction
CN106709183A (en) * 2016-12-26 2017-05-24 中国航空工业集团公司西安飞机设计研究所 Detail stress simulation method of composite mixed joint replaced by adopting metal plate
CN106874545A (en) * 2017-01-05 2017-06-20 武汉船用机械有限责任公司 A kind of method for designing of flat board with holes
CN106874545B (en) * 2017-01-05 2020-05-19 武汉船用机械有限责任公司 A design method of a flat plate with holes
CN107220422A (en) * 2017-05-19 2017-09-29 上海宇航系统工程研究所 A kind of extensive connector intensity all standing finite element modeling method
CN107526879B (en) * 2017-08-04 2020-05-19 华中科技大学 Three-dimensional model gradient finite element solving method based on ANSYS parametric design language
CN107526879A (en) * 2017-08-04 2017-12-29 华中科技大学 A kind of threedimensional model gradient finite element solving method based on ANSYS Parametric Design Languages
CN107885946A (en) * 2017-11-23 2018-04-06 江苏理工学院 A kind of method based on progressive damage model prediction composite Multi-fasteners joint intensity
CN109388868A (en) * 2018-09-25 2019-02-26 陕西飞机工业(集团)有限公司 A kind of contact of consideration nail hole and follow closely curved nail load calculation method
WO2020107529A1 (en) * 2018-11-29 2020-06-04 中国商用飞机有限责任公司北京民用飞机技术研究中心 Test method for pin load distribution measurement for mechanical connection structure
CN109682680B (en) * 2018-11-29 2021-07-20 中国商用飞机有限责任公司北京民用飞机技术研究中心 Complex composite material mechanical connection nail load measurement test method
CN109682680A (en) * 2018-11-29 2019-04-26 中国商用飞机有限责任公司北京民用飞机技术研究中心 A kind of complex composite material mechanical connected nail load measurement test method
CN109726435B (en) * 2018-12-04 2022-11-22 中国航空工业集团公司西安飞机设计研究所 A Method for Calculating the Load of Joint Nail Group
CN109726435A (en) * 2018-12-04 2019-05-07 中国航空工业集团公司西安飞机设计研究所 A Method for Calculating the Loads of Joint Nail Groups
CN109614713A (en) * 2018-12-13 2019-04-12 北京航空航天大学 Prediction method of nail load distribution for composite multi-nail connections based on modified stiffness method
CN109614713B (en) * 2018-12-13 2022-10-11 北京航空航天大学 Composite material multi-nail connection nail load distribution prediction method based on corrected stiffness method
CN109813587B (en) * 2019-03-08 2021-04-02 南昌航空大学 A test piece that is used for eccentric connection of rivet to connect nail to carry distribution research
CN109813587A (en) * 2019-03-08 2019-05-28 南昌航空大学 A test piece for the study of nail load distribution of rivet eccentric connection joints
CN111159879A (en) * 2019-12-25 2020-05-15 中国航空工业集团公司西安飞机设计研究所 Nail load processing method of centralized load transmission structure
CN111159879B (en) * 2019-12-25 2023-07-21 中国航空工业集团公司西安飞机设计研究所 A nail-loading processing method for a centralized load-carrying structure
CN112231864A (en) * 2020-11-16 2021-01-15 西北工业大学 A Rapid Analysis System for Composite Nail Loading
CN114741810A (en) * 2022-04-20 2022-07-12 中国船舶重工集团海装风电股份有限公司 High-precision bolt internal force extraction method and device based on finite element large deformation solution
CN114741810B (en) * 2022-04-20 2025-01-28 中船海装风电有限公司 High-precision bolt internal force extraction method and device based on finite element large deformation solution
CN115031942A (en) * 2022-05-27 2022-09-09 南京航空航天大学 Method for determining bolt tightening process parameters for assembling composite material structure
CN115098972A (en) * 2022-07-05 2022-09-23 浙江吉利控股集团有限公司 Multilayer structure bolted connection modeling method, device, equipment and storage medium
CN115098972B (en) * 2022-07-05 2025-04-15 浙江吉利控股集团有限公司 Modeling method, device, equipment and storage medium for multi-layer structure bolt connection
CN115828707A (en) * 2022-12-29 2023-03-21 中国航空工业集团公司成都飞机设计研究所 Method for calculating service life of fastening hole under complex loading condition

Also Published As

Publication number Publication date
CN102622472B (en) 2013-12-25

Similar Documents

Publication Publication Date Title
CN102622472A (en) Method for analyzing load and stress distribution and stress levels of composite mechanical connection pins
CN105701296B (en) A Finite Element Modeling Method for Racetrack Bolted Connection Structure
CN104268342A (en) Method for analyzing vibration characteristics of bearing retainer based on finite elements
CN103559390B (en) A kind of composite π shape based on average inefficacy index glueds joint attachment structure hot strength Forecasting Methodology
CN108345742B (en) Dynamic equivalent modeling method suitable for corrugated sandwich structure
CN103714216A (en) Composite material mechanics property evaluation method based on three-dimensional microcosmic crystal whisker composition
CN108344634A (en) A kind of armored concrete carefully sees threedimensional model modeling method
CN111159943B (en) A Buckling Treatment Method for Sealing Structure of Moving Airfoil
CN106650045A (en) A Method for Optimum Design and Strength Analysis of the Structural Section Between Secondary Boxes
CN105302994A (en) Finite element simulation method of beam_plate shell structure
CN115495962B (en) A bolt connection modeling method based on layered virtual materials
CN112528415A (en) Macro-meso failure mode analysis method for composite material shaft structure
CN109726438A (en) A Method for Calculating Temperature Stress of Hybrid Structures
CN111159879A (en) Nail load processing method of centralized load transmission structure
CN107368660A (en) A kind of analysis method of predetermined period lattice material yield surface
CN109726435A (en) A Method for Calculating the Loads of Joint Nail Groups
CN114297877A (en) Multi-working-condition simulation automation system and method for rod structure metamaterial structure
CN105468826A (en) Design method of composite material
CN109388868B (en) Nail load calculation method considering nail hole contact and nail bending
CN103177168A (en) Reverse engineering method of gear tooth surface friction factor
CN117292777A (en) Method, device, medium and equipment for predicting interlayer stress of composite material laminated plate based on sawtooth effect
CN117171864A (en) A linear vibration prediction method for beam structures
CN116090113A (en) A Mesh Finite Element Model of Block Structure with Variable Node Suspension Elements
Bournival et al. A mesh-geometry based method for coupling 1D and 3D elements
CN109902341A (en) An Improved Finite Element Method for Analysis of Nail Load Distribution in Composite Multi-Nail Connections

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20131225

Termination date: 20180227