CN114818372A - A Calculation Method of Tsai Modulus of Composite Materials - Google Patents
A Calculation Method of Tsai Modulus of Composite Materials Download PDFInfo
- Publication number
- CN114818372A CN114818372A CN202210553267.8A CN202210553267A CN114818372A CN 114818372 A CN114818372 A CN 114818372A CN 202210553267 A CN202210553267 A CN 202210553267A CN 114818372 A CN114818372 A CN 114818372A
- Authority
- CN
- China
- Prior art keywords
- composite material
- modulus
- sample
- tsai
- stiffness coefficient
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Operations Research (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
Description
技术领域technical field
本公开涉及飞机复合材料结构设计与分析技术领域,尤其涉及一种复合材料Tsai模量的计算方法。The present disclosure relates to the technical field of aircraft composite material structure design and analysis, and in particular, to a method for calculating the Tsai modulus of composite materials.
背景技术Background technique
Tsai模量是2014年由美国斯坦福大学教授、美国国家工程院院士、“复合材料之父”Stephen W.Tsai提出的表征复合材料整体刚度性能的新的弹性常数,即复合材料弹性矩阵的迹(Trace=Q11+Q22+2Q66),2020年,为了纪念该弹性常数的发现,同时为了表彰Stephen W.Tsai在复合材料领域的历史贡献,国际上正式将该弹性常数命名为“Tsai’sModulus”,即Tsai模量。The Tsai modulus is a new elastic constant proposed by Stephen W. Tsai, a professor at Stanford University, an academician of the National Academy of Engineering, and the "father of composite materials" in 2014, to characterize the overall stiffness performance of composite materials, that is, the trace of the elastic matrix of composite materials (Trace). =Q 11 +Q 22 +2Q 66 ), in 2020, in order to commemorate the discovery of the elastic constant, and in recognition of Stephen W. Tsai's historical contribution in the field of composite materials, the elastic constant was officially named "Tsai's Modulus" internationally. , which is the Tsai modulus.
依据Stephen W.Tsai的原始理论,仅需要通过0°铺层的单向板纵向拉伸试验即可近似得出某种材料的Tsai模量,具体方法为,通过某材料0°铺层的单向板纵向拉伸试验获得该材料纵向拉伸模量E1,然后即可以通过以下方法近似计算出该材料的Tsai模量:According to Stephen W.Tsai's original theory, the Tsai modulus of a certain material can be approximated by only the longitudinal tensile test of a 0° layered unidirectional plate. The longitudinal tensile modulus E 1 of the material is obtained by the longitudinal tensile test of the plate, and then the Tsai modulus of the material can be approximately calculated by the following method:
上述方法中的0.88由多种碳纤维复合材料力学性能统计得到。但当铺层角度φ=0°时,不同材料的归一化刚度系数仍有较大分散,因此,采用0°层压板试验计算碳纤维复合材料的Tsai模量时,误差会比较大。0.88 in the above method is obtained from the statistics of the mechanical properties of various carbon fiber composite materials. However, when the layup angle φ=0°, the normalized stiffness coefficients of different materials are still widely dispersed. Therefore, when the Tsai modulus of carbon fiber composites is calculated by the 0° laminate test, the error will be relatively large.
因此,有必要提供一种新的技术方案改善上述方案中存在的一个或者多个问题。Therefore, it is necessary to provide a new technical solution to improve one or more problems existing in the above solutions.
需要说明的是,在上述背景技术部分公开的信息仅用于加强对本公开的背景的理解,因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。It should be noted that the information disclosed in the above Background section is only for enhancement of understanding of the background of the present disclosure, and therefore may contain information that does not form the prior art that is already known to a person of ordinary skill in the art.
发明内容SUMMARY OF THE INVENTION
本公开的目的在于提供一种复合材料Tsai模量的计算方法,进而至少在一定程度上克服由于相关技术的限制和缺陷而导致的一个或者多个问题。The purpose of the present disclosure is to provide a method for calculating the Tsai modulus of a composite material, so as to at least to a certain extent overcome one or more problems caused by limitations and defects of the related art.
基于复合材料的纵向拉伸模量与铺层角、归一化刚度系数之间的预设对应关系,根据所述纵向拉伸模量范围来确定复合材料试样的铺层角和归一化刚度系数,其中,所述复合材料为碳纤维复合材料;Based on the preset correspondence between the longitudinal tensile modulus of the composite material, the layup angle, and the normalized stiffness coefficient, the layup angle and normalization of the composite material sample are determined according to the longitudinal tensile modulus range stiffness coefficient, wherein the composite material is a carbon fiber composite material;
根据确定出的所述铺层角按照正负铺层角进行复合材料的试样制备;According to the determined layup angle, the sample preparation of the composite material is carried out according to the positive and negative layup angles;
对制备好的所述试样进行拉伸试验测定,得到所述试样在拉伸载荷作用下的纵向应变,并根据所述纵向应变和所述试样的宽度确定试样的面内纵向刚度;Perform a tensile test on the prepared sample to obtain the longitudinal strain of the sample under tensile load, and determine the in-plane longitudinal stiffness of the sample according to the longitudinal strain and the width of the sample ;
根据所述试样的面内纵向刚度和所述试样的厚度确定复合材料的归一化刚度;determining a normalized stiffness of the composite material from the in-plane longitudinal stiffness of the specimen and the thickness of the specimen;
基于所述复合材料的归一化刚度和归一化刚度系数确定所述复合材料的Tsai模量。The Tsai modulus of the composite material is determined based on the normalized stiffness and normalized stiffness coefficient of the composite material.
本公开实施例中,所述试样为矩形平板拉伸试样或螺旋缠绕圆管试样。In the embodiment of the present disclosure, the sample is a rectangular flat plate tensile sample or a helically wound round tube sample.
本公开实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到。In the embodiment of the present disclosure, the layup angle and the normalized stiffness coefficient in the preset corresponding relationship are statistically obtained from a plurality of known elastic constants of composite materials.
本公开实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:In the embodiment of the present disclosure, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
统计每种已知复合材料的纵向拉伸模量、横向拉伸模量、面内泊松比、面内剪切模量,并根据经典层压板理论计算每种复合材料的弹性矩阵系数。The longitudinal tensile modulus, transverse tensile modulus, in-plane Poisson's ratio, and in-plane shear modulus of each known composite were counted, and the elastic matrix coefficients of each composite were calculated according to classical laminate theory.
本公开实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:In the embodiment of the present disclosure, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
根据每种所述复合材料的弹性矩阵系数计算出每种已知复合材料的Tsai模量。The Tsai modulus of each known composite was calculated from the elastic matrix coefficients of each of the composites.
本公开实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:In the embodiment of the present disclosure, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
设定具有所述正负铺层角两种铺层组成的虚拟层压板,将所述铺层角从0°逐渐增加到90°,根据经典层压板理论计算出任意所述正负铺层角的虚拟层压板的面内纵向刚度。Set a virtual laminate composed of two layers of the positive and negative layup angles, gradually increase the layup angle from 0° to 90°, and calculate any of the positive and negative layup angles according to the classical laminate theory The in-plane longitudinal stiffness of the virtual laminate.
本公开实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:In the embodiment of the present disclosure, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
根据所述虚拟层压板的厚度、面内纵向刚度和所述复合材料的Tsai模量得到所述归一化刚度系数。The normalized stiffness coefficient is derived from the thickness of the virtual laminate, the in-plane longitudinal stiffness, and the Tsai modulus of the composite.
本公开实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:In the embodiment of the present disclosure, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
根据所述虚拟层压板的铺层角绘制所述归一化刚度系数随所述铺层角的变化曲线。The normalized stiffness coefficient is plotted according to the layup angle of the virtual laminate as a function of the layup angle.
本公开实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:In the embodiment of the present disclosure, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
根据所述归一化刚度系数随量所述铺层角的变化曲线中的密集交汇点,得到该密集交汇点的横坐标和纵坐标,所述横坐标对应的数值为所述铺层角,所述纵坐标的数值为所述归一化刚度系数,并对得到的所述铺层角和所述归一化刚度系数进行修正。According to the dense intersection point in the variation curve of the normalized stiffness coefficient with the ply angle, the abscissa and the ordinate of the dense intersection point are obtained, and the value corresponding to the abscissa is the ply angle, The value of the ordinate is the normalized stiffness coefficient, and the obtained layup angle and the normalized stiffness coefficient are corrected.
本公开实施例中,所述复合材料为热塑性复合材料或热固性复合材料。In the embodiments of the present disclosure, the composite material is a thermoplastic composite material or a thermosetting composite material.
本公开的实施例提供的技术方案可以包括以下有益效果:The technical solutions provided by the embodiments of the present disclosure may include the following beneficial effects:
本公开的一种实施例中,通过上述方法,根据纵向拉伸模量范围确定复合材料的铺层角和归一化刚度系数,根据铺层角进行制备试样,对制备好的试样进行拉伸试验测定,得到拉伸载荷作用下的纵向应变,基于纵向应变和试样宽度确定试样的面内纵向刚度,通过面内纵向刚度和试样厚度确定归一化刚度,并基于归一化刚度和归一化刚度系数得到复合材料的Tsai模量。该方法能够解决原始方法误差大的问题,可简化传统的复合材料层压结构刚度分析方法。In an embodiment of the present disclosure, by the above method, the layup angle and the normalized stiffness coefficient of the composite material are determined according to the range of the longitudinal tensile modulus, the sample is prepared according to the layup angle, and the prepared sample is subjected to The tensile test is determined to obtain the longitudinal strain under the tensile load, the in-plane longitudinal stiffness of the specimen is determined based on the longitudinal strain and the specimen width, the normalized stiffness is determined by the in-plane longitudinal stiffness and the specimen thickness, and the normalized stiffness is determined based on the in-plane longitudinal stiffness and specimen thickness. The normalized stiffness and the normalized stiffness coefficient were used to obtain the Tsai modulus of the composite. This method can solve the problem of large error of the original method, and can simplify the traditional stiffness analysis method of composite laminated structure.
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本公开。It is to be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present disclosure.
附图说明Description of drawings
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本公开的实施例,并与说明书一起用于解释本公开的原理。显而易见地,下面描述中的附图仅仅是本公开的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the disclosure and together with the description serve to explain the principles of the disclosure. Obviously, the drawings in the following description are only some embodiments of the present disclosure, and for those of ordinary skill in the art, other drawings can also be obtained from these drawings without creative effort.
图1示意性示出本公开示例性实施例中复合材料Tsai模量的计算方法的步骤流程图;FIG. 1 schematically shows a flow chart of steps of a method for calculating the Tsai modulus of a composite material in an exemplary embodiment of the present disclosure;
图2示意性示出本公开示例性实施例中复合材料矩形平板拉伸试样结构及加载示意图;FIG. 2 schematically shows the structure and loading schematic diagram of the composite rectangular flat plate tensile specimen in the exemplary embodiment of the present disclosure;
图3示意性示出本公开示例性实施例中复合材料螺旋缠绕圆管试样结构及加载示意图;FIG. 3 schematically shows the structure and loading schematic diagram of the composite material helically wound round tube sample in the exemplary embodiment of the present disclosure;
图4示意性示出本公开示例性实施例中虚拟层压板的归一化刚度系数随铺层角φ的变化曲线。FIG. 4 schematically illustrates the normalized stiffness coefficients of virtual laminates in an exemplary embodiment of the present disclosure Variation curve with ply angle φ.
具体实施方式Detailed ways
现在将参考附图更全面地描述示例实施方式。然而,示例实施方式能够以多种形式实施,且不应被理解为限于在此阐述的范例;相反,提供这些实施方式使得本公开将更加全面和完整,并将示例实施方式的构思全面地传达给本领域的技术人员。所描述的特征、结构或特性可以以任何合适的方式结合在一个或更多实施方式中。Example embodiments will now be described more fully with reference to the accompanying drawings. Example embodiments, however, can be embodied in various forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of example embodiments to those skilled in the art. The described features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.
此外,附图仅为本公开的示意性图解,并非一定是按比例绘制。图中相同的附图标记表示相同或类似的部分,因而将省略对它们的重复描述。附图中所示的一些方框图是功能实体,不一定必须与物理或逻辑上独立的实体相对应。可以采用软件形式来实现这些功能实体,或在一个或多个硬件模块或集成电路中实现这些功能实体,或在不同网络和/或处理器装置和/或微控制器装置中实现这些功能实体。Furthermore, the drawings are merely schematic illustrations of the present disclosure and are not necessarily drawn to scale. The same reference numerals in the drawings denote the same or similar parts, and thus their repeated descriptions will be omitted. Some of the block diagrams shown in the figures are functional entities that do not necessarily necessarily correspond to physically or logically separate entities. These functional entities may be implemented in software, or in one or more hardware modules or integrated circuits, or in different networks and/or processor devices and/or microcontroller devices.
但由图4可以看出,当铺层角度φ=0°时,不同材料的归一化刚度系数仍有较大分散,因此,采用0°层压板试验计算碳纤维复合材料的Tsai模量时,误差会比较大。However, it can be seen from Figure 4 that when the layup angle φ=0°, the normalized stiffness coefficients of different materials are still widely dispersed. Therefore, when the Tsai modulus of carbon fiber composites is calculated by the 0° laminate test, the error will be larger.
本示例实施方式中提供了一种复合材料Tsai模量的计算方法,参考图1中所示,该方法可以包括:This example embodiment provides a method for calculating the Tsai modulus of a composite material. Referring to Figure 1, the method may include:
步骤S101:基于复合材料的纵向拉伸模量与铺层角、归一化刚度系数之间的预设对应关系,根据所述纵向拉伸模量范围来确定复合材料试样的铺层角和归一化刚度系数,其中,所述复合材料为碳纤维复合材料。Step S101: Based on the preset correspondence between the longitudinal tensile modulus of the composite material, the layup angle, and the normalized stiffness coefficient, determine the layup angle and Normalized stiffness coefficient, wherein the composite material is a carbon fiber composite material.
步骤S102:根据确定出的所述铺层角按照正负铺层角进行复合材料的试样制备。Step S102: According to the determined layup angle, the sample preparation of the composite material is performed according to the positive and negative layup angles.
步骤S103:对制备好的所述试样进行拉伸试验测定,得到所述试样在拉伸载荷作用下的纵向应变,并根据所述纵向应变和所述试样的宽度确定试样的面内纵向刚度。Step S103 : perform a tensile test on the prepared sample, obtain the longitudinal strain of the sample under the action of the tensile load, and determine the surface of the sample according to the longitudinal strain and the width of the sample Inner longitudinal stiffness.
步骤S104:根据所述试样的面内纵向刚度和所述试样的厚度确定复合材料的归一化刚度。Step S104: Determine the normalized stiffness of the composite material according to the in-plane longitudinal stiffness of the sample and the thickness of the sample.
步骤S105:基于所述复合材料的归一化刚度和归一化刚度系数确定所述复合材料的Tsai模量。Step S105: Determine the Tsai modulus of the composite material based on the normalized stiffness and the normalized stiffness coefficient of the composite material.
通过上述方法,根据纵向拉伸模量范围确定复合材料的铺层角和归一化刚度系数,根据铺层角进行制备试样,对制备好的试样进行拉伸试验测定,得到拉伸载荷作用下的纵向应变,基于纵向应变和试样宽度确定试样的面内纵向刚度,通过面内纵向刚度和试样厚度确定归一化刚度,并基于归一化刚度和归一化刚度系数得到复合材料的Tsai模量。该方法能够解决原始方法误差大的问题,可简化传统的复合材料层压结构刚度分析方法。By the above method, the layup angle and normalized stiffness coefficient of the composite material are determined according to the range of longitudinal tensile modulus, the sample is prepared according to the layup angle, and the prepared sample is subjected to tensile test measurement to obtain the tensile load. The longitudinal strain under the action, the in-plane longitudinal stiffness of the specimen is determined based on the longitudinal strain and the specimen width, the normalized stiffness is determined by the in-plane longitudinal stiffness and the specimen thickness, and is obtained based on the normalized stiffness and the normalized stiffness coefficient Tsai modulus of composites. This method can solve the problem of large error of the original method, and can simplify the traditional stiffness analysis method of composite laminated structure.
下面,将参考图1至图4对本示例实施方式中的上述方法的各个步骤进行更详细的说明。Hereinafter, each step of the above method in this example embodiment will be described in more detail with reference to FIGS. 1 to 4 .
在步骤S101中,基于复合材料的纵向拉伸模量与铺层角、归一化刚度系数之间的预设对应关系,根据所述纵向拉伸模量范围来确定复合材料试样的铺层角和归一化刚度系数,其中,所述复合材料为碳纤维复合材料。具体的,铺层角φ、归一化刚度系数与复合材料的纵向拉伸模量E1之间的预设对应关系如表1所示。根据纵向拉伸模量E1范围通过表1查找对应的铺层角φ和归一化刚度系数 In step S101, based on the preset correspondence between the longitudinal tensile modulus of the composite material, the layup angle, and the normalized stiffness coefficient, the layup of the composite material sample is determined according to the longitudinal tensile modulus range angle and normalized stiffness coefficient, wherein the composite material is a carbon fiber composite material. Specifically, the layup angle φ, normalized stiffness coefficient The preset correspondence with the longitudinal tensile modulus E1 of the composites is shown in Table 1 . According to the range of longitudinal tensile modulus E1, look up the corresponding ply angle φ and normalized stiffness coefficient through Table 1
表1复合材料的纵向拉伸模量E1与铺层角φ、归一化刚度系数之间的预设对应关系Table 1 Longitudinal tensile modulus E1 of composites versus ply angle φ, normalized stiffness coefficient preset correspondence between
在步骤S102中,根据确定出的所述铺层角按照正负铺层角进行复合材料的试样制备。具体的,根据表1中预设对应关系确定出的铺层角大小,按照正负铺层角±φ两种铺层角进行复合材料的试样制备,便于得到具有对应该铺层角大小的试样,方便后续进行试样的试验测定。In step S102, the sample preparation of the composite material is performed according to the determined layup angles according to the positive and negative layup angles. Specifically, according to the size of the layup angle determined by the preset corresponding relationship in Table 1, the sample preparation of the composite material is carried out according to the positive and negative layup angles ±φ two kinds of layup angles, so that it is convenient to obtain a sample with a size corresponding to the layup angle. The sample is convenient for subsequent testing and measurement of the sample.
在步骤S103中,对制备好的所述试样进行拉伸试验测定,得到所述试样在拉伸载荷作用下的纵向应变,并根据所述纵向应变和所述试样的宽度确定试样的面内纵向刚度。具体的,按照确定出的铺层角制备好试样后,在对试样进行拉伸试样测定时,对试样一端固定,另一端施加拉伸载荷F,应保证在当前载荷下试样层压板内不发生任何形式的损伤,根据试验测定该试样层压板在拉伸载荷F作用下的纵向应变ε0。然后根据纵向应变和试样宽度计算试样的面内纵向刚度。具体按照如下公式(1)进行计算。In step S103, a tensile test is performed on the prepared sample to obtain the longitudinal strain of the sample under the action of the tensile load, and the sample is determined according to the longitudinal strain and the width of the sample in-plane longitudinal stiffness. Specifically, after the sample is prepared according to the determined layup angle, when the sample is subjected to tensile sample measurement, one end of the sample is fixed, and the other end is subjected to a tensile load F to ensure that the sample is under the current load. No damage of any kind occurred in the laminate, and the longitudinal strain ε 0 of the sample laminate under tensile load F was determined experimentally. The in-plane longitudinal stiffness of the specimen is then calculated from the longitudinal strain and the specimen width. Specifically, it is calculated according to the following formula (1).
在步骤S104中,根据所述试样的面内纵向刚度和所述试样的厚度确定复合材料的归一化刚度。具体的,根据如下公式(2)计算试样的归一化刚度 In step S104, the normalized stiffness of the composite material is determined according to the in-plane longitudinal stiffness of the sample and the thickness of the sample. Specifically, the normalized stiffness of the specimen is calculated according to the following formula (2)
在步骤S105中,基于所述复合材料的归一化刚度和归一化刚度系数确定所述复合材料的Tsai模量。具体的,根据如下公式(3)计算复合材料的Tsai模量。In step S105, the Tsai modulus of the composite material is determined based on the normalized stiffness and the normalized stiffness coefficient of the composite material. Specifically, the Tsai modulus of the composite material was calculated according to the following formula (3).
可选的,在一些实施例中,所述试样为矩形平板拉伸试样或螺旋缠绕圆管试样。具体的,矩形平板拉伸试样为复合材料以正负铺层角±φ两种铺层角进行铺放制备而成,螺旋缠绕圆管试样为复合材料以正负铺层角±φ两种铺层角进行螺旋缠绕制备而成,其中,矩形平板拉伸试样如图2所示,其中,图中L表示试样的长度,螺旋缠绕圆管试样如图3所示。Optionally, in some embodiments, the sample is a rectangular flat plate tensile sample or a helically wound round tube sample. Specifically, the rectangular flat plate tensile specimen is prepared by laying the composite material at the positive and negative layup angles ±φ, and the spirally wound round tube sample is the composite material with the positive and negative layup angles ±φ. The seed ply corners are prepared by spiral winding. The rectangular flat plate tensile sample is shown in Figure 2, where L represents the length of the sample, and the spiral wound round tube sample is shown in Figure 3.
可选的,在一些实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到。具体的,通过大量已知的复合材料弹性常数统计得到如表1所示的纵向拉伸模量与铺层角、归一化刚度系数之间的预设对应关系。Optionally, in some embodiments, the layup angle and the normalized stiffness coefficient in the preset correspondence are obtained by statistics of a plurality of known elastic constants of composite materials. Specifically, the preset corresponding relationship between the longitudinal tensile modulus, the layup angle, and the normalized stiffness coefficient as shown in Table 1 is obtained through statistics of a large number of known elastic constants of composite materials.
可选的,在一些实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:Optionally, in some embodiments, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
统计每种已知复合材料的纵向拉伸模量、横向拉伸模量、面内泊松比、面内剪切模量,并根据经典层压板理论计算每种复合材料的弹性矩阵系数。具体的,统计每种已知复合材料的纵向拉伸模量E1、横向拉伸模量E2、面内泊松比μ12、面内剪切模量G12,根据经典层压板理论计算每种复合材料的弹性矩阵系数Q11、Q22、Q12、Q66等。其中,经典层压板理论基于以下3个假设:(1)层间变形一致假设:层压板各单层之间粘和牢固,层间变形一致,无相对位移;(2)直法线假设:变形前垂直于板中面的直线在变形后仍保持垂直,且长度不变;(3)平面应力状态假设:层压板中各单层都可近似地认为处于平面应力状态。The longitudinal tensile modulus, transverse tensile modulus, in-plane Poisson's ratio, and in-plane shear modulus of each known composite were counted, and the elastic matrix coefficients of each composite were calculated according to classical laminate theory. Specifically, the longitudinal tensile modulus E 1 , the transverse tensile modulus E 2 , the in-plane Poisson’s ratio μ 12 , and the in-plane shear modulus G 12 of each known composite material are calculated according to the classical laminate theory. The elastic matrix coefficients Q 11 , Q 22 , Q 12 , Q 66 , etc. of each composite material. Among them, the classical laminate theory is based on the following 3 assumptions: (1) The interlayer deformation consistency hypothesis: the single layers of the laminate are sticky and firm, the interlayer deformation is consistent, and there is no relative displacement; (2) The straight normal line assumption: deformation The straight line perpendicular to the middle surface of the board remains vertical after deformation, and its length remains unchanged; (3) Assumption of plane stress state: each single layer in the laminate can be approximately considered to be in a plane stress state.
可选的,在一些实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:Optionally, in some embodiments, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
根据每种所述复合材料的弹性矩阵系数计算出每种已知复合材料的Tsai模量。具体的,根据如下公式(4)计算每种已知复合材料的Tsai模量,其中每种复合材料的Tsai模量即为该复合材料弹性矩阵的迹。The Tsai modulus of each known composite was calculated from the elastic matrix coefficients of each of the composites. Specifically, the Tsai modulus of each known composite material is calculated according to the following formula (4), wherein the Tsai modulus of each composite material is the trace of the elastic matrix of the composite material.
Trace=Q11+Q22+2Q66 (4)Trace=Q 11 +Q 22 +2Q 66 (4)
可选的,在一些实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:Optionally, in some embodiments, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
设定具有所述正负铺层角两种铺层组成的虚拟层压板,将所述铺层角从0°逐渐增加到90°,根据经典层压板理论计算出任意所述正负铺层角的虚拟层压板的面内纵向刚度A11。具体的,虚拟层压板由+φ和-φ两种铺层角进行铺层而成,当铺层角从0°逐渐增加到90°,根据经典层压板理论计算出任意±φ的虚拟层压板的面内纵向刚度A11。Set a virtual laminate composed of two layers of the positive and negative layup angles, gradually increase the layup angle from 0° to 90°, and calculate any of the positive and negative layup angles according to the classical laminate theory The in-plane longitudinal stiffness A 11 of the virtual laminate. Specifically, the virtual laminate is formed by laying up two ply angles of +φ and -φ. When the ply angle gradually increases from 0° to 90°, the virtual laminate of any ±φ is calculated according to the classical laminate theory. In-plane longitudinal stiffness A 11 .
可选的,在一些实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:Optionally, in some embodiments, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
根据所述虚拟层压板的厚度、面内纵向刚度和所述复合材料的Tsai模量得到所述归一化刚度系数。具体的,通过面内纵向刚度A11和虚拟层压板的厚度h计算得到归一化纵向刚度然后再根据得到的归一化纵向刚度和复合材料的Tsai模量得到归一化刚度系数 The normalized stiffness coefficient is derived from the thickness of the virtual laminate, the in-plane longitudinal stiffness, and the Tsai modulus of the composite. Specifically, the normalized longitudinal stiffness is obtained by calculating the in-plane longitudinal stiffness A 11 and the thickness h of the virtual laminate The normalized stiffness coefficient is then obtained based on the obtained normalized longitudinal stiffness and the Tsai modulus of the composite
可选的,在一些实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:Optionally, in some embodiments, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
根据所述虚拟层压板的所述铺层角绘制所述归一化刚度系数随所述铺层角的变化曲线。具体的,虚拟层压板的归一化刚度系数与铺层角φ的关系绘制虚拟层压板的归一化刚度系数随铺层角φ的变化曲线,如图4所示。The normalized stiffness coefficient as a function of the layup angle is plotted according to the layup angle of the virtual laminate. Specifically, the normalized stiffness coefficient of the virtual laminate Normalized stiffness coefficients for virtual laminates plotted as a function of ply angle φ The variation curve with the layup angle φ is shown in Figure 4.
可选的,在一些实施例中,所述预设对应关系中的所述铺层角和所述归一化刚度系数通过多个已知的复合材料弹性常数统计得到的步骤包括:Optionally, in some embodiments, the step of statistically obtaining the layup angle and the normalized stiffness coefficient in the preset correspondence relationship through a plurality of known elastic constants of composite materials includes:
根据所述归一化刚度系数随所述铺层角的变化曲线中的密集交汇点,得到该密集交汇点的横坐标和纵坐标,所述横坐标对应的数值为所述铺层角,所述纵坐标的数值为所述归一化刚度系数,并对得到的所述铺层角和所述归一化刚度系数进行修正。具体的,通过图4,找到多种复合材料归一化刚度系数随铺层角φ的变化曲线中的密集交汇点A,该密集交汇点A可称为焦点A,密集交汇点A的横坐标对应的数值即为铺层角φ,密集交汇点A的纵坐标对应的数值即为归一化刚度系数此时得到的铺层角φ和归一化刚度系数分别可以通过大量复合材料纵向拉伸模量、横向拉伸模量、泊松比及剪切模量的样本数据进行修正。修正后的铺层角φ、归一化刚度系数和纵向拉伸模量E1范围之间形成预设对应关系,如表1所示。The abscissa and ordinate of the dense junction are obtained according to the dense junction in the variation curve of the normalized stiffness coefficient with the ply angle, and the value corresponding to the abscissa is the ply angle, so The value of the ordinate is the normalized stiffness coefficient, and the obtained layup angle and the normalized stiffness coefficient are corrected. Specifically, through Figure 4, find the normalized stiffness coefficients of various composite materials The dense junction A in the curve with the change of the ply angle φ, the dense junction A can be called the focus A, the value corresponding to the abscissa of the dense junction A is the ply angle φ, and the ordinate of the dense junction A The corresponding value is the normalized stiffness coefficient The ply angle φ and normalized stiffness coefficient obtained at this time It can be corrected by a large number of sample data of longitudinal tensile modulus, transverse tensile modulus, Poisson's ratio and shear modulus of composite materials, respectively. Corrected ply angle φ, normalized stiffness coefficient A preset correspondence is formed between the range of the longitudinal tensile modulus E1, as shown in Table 1 .
可选的,在一些实施例中,所述复合材料为热塑性复合材料或热固性复合材料。具体的,不论是热塑性复合材料,或是热固性复合材料,都可以计算出其Tsai模量。Optionally, in some embodiments, the composite material is a thermoplastic composite material or a thermosetting composite material. Specifically, whether it is a thermoplastic composite material or a thermosetting composite material, its Tsai modulus can be calculated.
通过上述方法,根据纵向拉伸模量范围确定复合材料的铺层角和归一化刚度系数,根据铺层角进行制备试样,对制备好的试样进行试验测定,得到等效纵向拉伸模量,基于等效纵向拉伸模量、试样宽度和试样厚度得到复合材料的归一化刚度,再通过归一化刚度和归一化刚度系数得到复合材料的Tsai模量。该方法能够解决原始方法误差大的问题,可简化传统的复合材料层压结构刚度分析方法。Through the above method, the layup angle and normalized stiffness coefficient of the composite material are determined according to the range of longitudinal tensile modulus, the samples are prepared according to the layup angle, and the prepared samples are tested and measured to obtain the equivalent longitudinal tensile force. Modulus, the normalized stiffness of the composite material is obtained based on the equivalent longitudinal tensile modulus, specimen width and specimen thickness, and then the Tsai modulus of the composite material is obtained by the normalized stiffness and normalized stiffness coefficient. This method can solve the problem of large error of the original method, and can simplify the traditional stiffness analysis method of composite laminated structure.
本领域技术人员在考虑说明书及实践这里公开的发明后,将容易想到本公开的其它实施方案。本申请旨在涵盖本公开的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本公开的真正范围和精神由所附的权利要求指出。Other embodiments of the present disclosure will readily occur to those skilled in the art upon consideration of the specification and practice of the invention disclosed herein. This application is intended to cover any variations, uses, or adaptations of the present disclosure that follow the general principles of the present disclosure and include common knowledge or techniques in the technical field not disclosed by the present disclosure . The specification and examples are to be regarded as exemplary only, with the true scope and spirit of the disclosure being indicated by the appended claims.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210553267.8A CN114818372B (en) | 2022-05-20 | 2022-05-20 | A method for calculating Tsai modulus of composite materials |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210553267.8A CN114818372B (en) | 2022-05-20 | 2022-05-20 | A method for calculating Tsai modulus of composite materials |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114818372A true CN114818372A (en) | 2022-07-29 |
CN114818372B CN114818372B (en) | 2024-10-15 |
Family
ID=82517682
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210553267.8A Active CN114818372B (en) | 2022-05-20 | 2022-05-20 | A method for calculating Tsai modulus of composite materials |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114818372B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114818374A (en) * | 2022-05-20 | 2022-07-29 | 中国航空工业集团公司西安飞机设计研究所 | A Calculation Method of Tsai Modulus of Composite Materials |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20050112788A (en) * | 2004-05-28 | 2005-12-01 | 학교법인 대양학원 | Method for estimating effective stiffness and optimizing layer angle of composite material |
CN106126773A (en) * | 2016-06-12 | 2016-11-16 | 北京航空航天大学 | A kind of intensity prediction method containing uncertain parameter composite laminated plate lost efficacy based on whole layer |
CN108984841A (en) * | 2018-06-19 | 2018-12-11 | 威海光威复合材料股份有限公司 | Strength check under the calculating of composite laminated plate concrete moduli and specified load |
CN109063280A (en) * | 2018-07-16 | 2018-12-21 | 北京航空航天大学 | A kind of needle thorn C/C Modulus of Composites calculation method |
CN114818374A (en) * | 2022-05-20 | 2022-07-29 | 中国航空工业集团公司西安飞机设计研究所 | A Calculation Method of Tsai Modulus of Composite Materials |
-
2022
- 2022-05-20 CN CN202210553267.8A patent/CN114818372B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20050112788A (en) * | 2004-05-28 | 2005-12-01 | 학교법인 대양학원 | Method for estimating effective stiffness and optimizing layer angle of composite material |
CN106126773A (en) * | 2016-06-12 | 2016-11-16 | 北京航空航天大学 | A kind of intensity prediction method containing uncertain parameter composite laminated plate lost efficacy based on whole layer |
CN108984841A (en) * | 2018-06-19 | 2018-12-11 | 威海光威复合材料股份有限公司 | Strength check under the calculating of composite laminated plate concrete moduli and specified load |
CN109063280A (en) * | 2018-07-16 | 2018-12-21 | 北京航空航天大学 | A kind of needle thorn C/C Modulus of Composites calculation method |
CN114818374A (en) * | 2022-05-20 | 2022-07-29 | 中国航空工业集团公司西安飞机设计研究所 | A Calculation Method of Tsai Modulus of Composite Materials |
Non-Patent Citations (3)
Title |
---|
LIYONG JIA等: "Validation and development of trace-based approach fro composite laminates", COMPOSITES SCIENCE AND TECHNOLOGY, vol. 221, 12 April 2022 (2022-04-12) * |
伍春波;陈秀华;余音;汪海;: "复合材料厚板双轴非线性刚度特性分析", 航空学报, no. 07, 6 December 2011 (2011-12-06) * |
马子广;陈庆童;王卫卫;: "复合材料层压板开孔拉伸力学性能探究", 直升机技术, no. 01, 15 March 2015 (2015-03-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114818374A (en) * | 2022-05-20 | 2022-07-29 | 中国航空工业集团公司西安飞机设计研究所 | A Calculation Method of Tsai Modulus of Composite Materials |
CN114818374B (en) * | 2022-05-20 | 2024-10-18 | 中国航空工业集团公司西安飞机设计研究所 | A calculation method for Tsai modulus of composite materials |
Also Published As
Publication number | Publication date |
---|---|
CN114818372B (en) | 2024-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Al Abadi et al. | Elastic properties of 3D printed fibre-reinforced structures | |
CN110672417B (en) | A method for obtaining elastic-plastic properties of ultra-thin materials by small punch test | |
CN109991077B (en) | A Prediction Method for Hybrid Delamination Resistance Curves of Composite Materials | |
CN104794299B (en) | A kind of composite interference coordinates joint stress distribution computational methods | |
CN107966354B (en) | Fatigue life prediction method and device for composite material and electronic equipment | |
Groo et al. | In situ damage detection for fiber‐reinforced composites using integrated zinc oxide nanowires | |
CN114818372A (en) | A Calculation Method of Tsai Modulus of Composite Materials | |
CN112163273B (en) | Multi-scale equivalent modeling method for trapezoidal corrugated sandwich cylindrical shell made of composite material | |
CN114839082A (en) | Method for measuring and calculating interlayer shear strength of heterogeneous multilayer composite board | |
CN112733341B (en) | A Progressive Failure Analysis Method for Marine Fiber Reinforced Flexible Pipes | |
CN114818374A (en) | A Calculation Method of Tsai Modulus of Composite Materials | |
Yang et al. | Improved tensile and single-lap-shear mechanical-electrical response of epoxy composites reinforced with gridded nano-carbons | |
CN110161132A (en) | The construction method of the defect model of carbon fibre composite | |
CN117433990B (en) | Method, system and equipment for determining interlayer shear strength of composite material laminated plate | |
CN111553102A (en) | A Finite Element Prediction Method for Static Strength of Composite Stiffened Panels Containing Moisture Heat Effect | |
CN111261239B (en) | Composite material structure damage identification method based on damage mechanics | |
CN114297842B (en) | An efficient method for estimating torsional stiffness of carbon fiber composite wound pipe fittings | |
CN114297839B (en) | Method for designing and checking torsional rigidity of carbon fiber composite material wound pipe fitting | |
Guo et al. | Failure mode analysis and validation of hierarchical sandwich structures under compression | |
CN111959095B (en) | An online health monitoring method for fiber-reinforced metal laminate materials | |
CN115688523A (en) | Method for establishing single-layer equivalent structure segmented mechanical model of multilayer structure of aramid fiber honeycomb wall | |
Zhang et al. | Experiments and analysis on stability of the sandwich structures with soft core | |
Chen et al. | Elastic constants identification of symmetric angle-ply laminates via a two-level optimization approach | |
CN105043728A (en) | Method of calculating tensile bearing capacity of joint to be measured | |
CN114184474B (en) | Tensile-shear coupling test method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |