CN113962096B - 一种导弹发射筒三向应力确定方法及系统 - Google Patents

一种导弹发射筒三向应力确定方法及系统 Download PDF

Info

Publication number
CN113962096B
CN113962096B CN202111247041.7A CN202111247041A CN113962096B CN 113962096 B CN113962096 B CN 113962096B CN 202111247041 A CN202111247041 A CN 202111247041A CN 113962096 B CN113962096 B CN 113962096B
Authority
CN
China
Prior art keywords
stress
fiber layer
layer
determining
inner liner
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.)
Active
Application number
CN202111247041.7A
Other languages
English (en)
Other versions
CN113962096A (zh
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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202111247041.7A priority Critical patent/CN113962096B/zh
Publication of CN113962096A publication Critical patent/CN113962096A/zh
Application granted granted Critical
Publication of CN113962096B publication Critical patent/CN113962096B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/26—Composites
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Moulding By Coating Moulds (AREA)

Abstract

本发明涉及一种导弹发射筒三向应力确定方法及系统。该方法包括根据导弹发射筒中纤维层及内衬层的材料基本力学性能确定纤维层及内衬层在柱坐标系下的应力与应变关系；根据纤维层及内衬层在柱坐标系下的应力与应变关系确定纤维层和内衬层的应力的分布方程；根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力；根据径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常。本发明能够实现碳纤维增强金属内衬导弹发射筒应力快速计算、校核。

Description

一种导弹发射筒三向应力确定方法及系统
技术领域
本发明涉及导弹发射筒强度分析领域，特别是涉及一种导弹发射筒三向应力确定方法及系统。
背景技术
复合材料具有高强度、高比模量、耐腐蚀和性能可设计强等一系列的优点，在航空航天、船舶、汽车、体育等领域的应用日益广泛。碳纤维增强树脂基复合材料是一种应用广泛的性能优良的先进复合材料，随着复合材料技术的不断革新，这种材料在军工领域的应用越来越广泛，先进复合材料的应用已经成为未来军事技术领域竞争的又一高地。不同于传统的金属材料，碳纤维增强树脂基复合材料是各向异性材料，而金属内衬是各向同性材料，在内压的作用下，内衬的受力情况较易计算，但是复合材料增强层的受力情况较为复杂。在实际制造过程中，增强层的纤维按照设计要求与轴向成一定的角度进行铺放或缠绕，增强层复合材料的材料坐标系和全局坐标系一般不重合，因此复合材料增强层的强度的计算较为复杂。碳纤维复合材料的强度高，但却是典型的脆性材料，因此研究其受力情况判断是否发生失效就成为设计阶段必不可缺的环节。
目前，碳纤维增强金属内衬导弹发射筒应力计算大多采用现有的商用有限元软件，如ABAQUS，ANSYS等，但其在建模时过于复杂以及计算时间过长，大大加长了设计周期。
因此，亟需一种新方法或系统，以实现碳纤维增强金属内衬导弹发射筒应力快速计算、校核，进而对碳纤维增强金属内衬导弹发射筒设计具有指导性的意义。
发明内容
本发明的目的是提供一种导弹发射筒三向应力确定方法及系统，能够实现碳纤维增强金属内衬导弹发射筒应力快速计算、校核。
为实现上述目的，本发明提供了如下方案：
一种导弹发射筒三向应力确定方法，包括：
根据导弹发射筒中纤维层及内衬层的材料基本力学性能确定纤维层及内衬层在柱坐标系下的应力与应变关系；所述材料基本力学性能包括：弹性模量以及泊松比；
根据纤维层及内衬层在柱坐标系下的应力与应变关系确定纤维层的应力的分布方程和内衬层的应力的分布方程；
根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力；
根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常。
可选地，所述根据纤维层及内衬层在柱坐标系下的应力与应变关系确定纤维层的应力的分布方程和内衬层的应力的分布方程，具体包括：
利用公式
Figure BDA0003321415040000021
确定内衬层的应力的分布方程；
利用公式
Figure BDA0003321415040000022
确定纤维层的应力的分布方程；
其中，
Figure BDA0003321415040000023
和
Figure BDA0003321415040000024
分别为内衬层在柱坐标系下的应力，
Figure BDA0003321415040000025
和
Figure BDA0003321415040000026
分别为纤维层在柱坐标系下的应力，w为系数，C3 k和C4 k为系数，μ为内衬层的泊松比，r为沿内衬层径向的坐标，r0≤r≤r1，r0为内衬层的内径，r1是内衬层的外径，C1和C2为系数，f1 k、
Figure BDA0003321415040000027
以及
Figure BDA0003321415040000028
均为系数，rqk为r的(qk)次方，
Figure BDA0003321415040000029
为r的(-qk-1)次方，qk为纤维层径向位移场中的系数。
可选地，所述根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力，之前还包括：
利用公式
Figure BDA0003321415040000031
确定内衬层的径向位移场；
利用公式
Figure BDA0003321415040000032
确定第k层纤维层的径向位移场；
利用公式
Figure BDA0003321415040000033
确定边界约束；
其中，εA＝0，
Figure BDA0003321415040000034
为内衬层的径向位移场，
Figure BDA0003321415040000035
为第k层纤维层的径向位移场，Pw为内压，M为纤维铺层总数，rM为第M层纤维层的内径，M为正数，R为纤维层的总的外径，
Figure BDA0003321415040000036
为第M层纤维层的径向的应力，pk为纤维层径向位移场中的系数。
可选地，所述根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力，具体包括：
利用公式
Figure BDA0003321415040000037
Figure BDA0003321415040000038
确定应力分布方程系数的解的方程；
其中，C1，C2，···，C4 M均为系数。
可选地，所述根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常，具体包括：
利用公式
Figure BDA0003321415040000041
将纤维层和内衬层在柱坐标系下的三向应力转换为笛卡尔直角坐标系下的三向应力；
其中，
Figure BDA0003321415040000042
和
Figure BDA0003321415040000043
分别为笛卡尔直角坐标系下的三向应力，m＝cosα，n＝sinα，α是第k层的缠绕角。
可选地，所述根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常，具体包括：
利用公式
Figure BDA0003321415040000044
验证纤维层；
若ζ>1，则纤维层失效；反之，则纤维层正常；
其中，
Figure BDA0003321415040000045
Figure BDA0003321415040000046
Xt为纤维纵向拉伸强度，Xc为纤维纵向压缩强度，Yt为纤维横向拉伸强度，Yc为纤维横向压缩强度。
一种导弹发射筒三向应力确定系统，包括：
应力与应变关系确定模块，用于根据导弹发射筒中纤维层及内衬层的材料基本力学性能确定纤维层及内衬层在柱坐标系下的应力与应变关系；所述材料基本力学性能包括：弹性模量以及泊松比；
应力的分布方程确定模块，用于根据纤维层及内衬层在柱坐标系下的应力与应变关系确定纤维层的应力的分布方程和内衬层的应力的分布方程；
径向位移和三向应力确定模块，用于根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力；
验证结果确定模块，用于根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常。
可选地，所述应力的分布方程确定模块具体包括：
内衬层的应力的分布方程确定单元，用于利用公式
Figure BDA0003321415040000051
确定内衬层的应力的分布方程；
纤维层的应力的分布方程确定单元，用于利用公式
Figure BDA0003321415040000052
确定纤维层的应力的分布方程；
其中，
Figure BDA0003321415040000053
和
Figure BDA0003321415040000054
分别为内衬层在柱坐标系下的应力，
Figure BDA0003321415040000055
和
Figure BDA0003321415040000056
分别为纤维层在柱坐标系下的应力，w为系数，C3 k和C4 k为系数，μ为内衬层的泊松比，r为沿内衬层径向的坐标，r0≤r≤r1，r0为内衬层的内径，r1是内衬层的外径，C1和C2为系数，f1 k、
Figure BDA0003321415040000057
以及
Figure BDA0003321415040000058
均为系数，rqk为r的(qk)次方，
Figure BDA0003321415040000059
为r的(-qk-1)次方，qk为纤维层径向位移场中的系数。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明所提供的一种导弹发射筒三向应力确定方法及系统，根据建立每一层纤维的柱坐标系下的本构方程，然后根据各层的应力分布方程和径向位移的分布方程、边界约束确定分布方程中的所有系数，即可确定应力-半径和径向位移-半径之间的关系，将各层的半径代入各层的分布方程，即可求出柱坐标系下的三向应力，进而再进行验证。本发明求解过程更加迅速，计算过程更加方便，计算效率有所提高。比现有的有限元计算应力校核强度具有更高的效率，缩短了产品结构设计的周期。进而能够实现碳纤维增强金属内衬导弹发射筒应力快速计算、校核。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明所提供的一种导弹发射筒三向应力确定方法流程示意图；
图2为(内压20MPa，铺层全为环向，S11表示σ1)本发明解析法计算结果与有限元计算结果的对比示意图；
图3为(内压20MPa，铺层全为环向，S22表示σ2)本发明解析法计算结果与有限元计算结果的对比示意图；
图4为(内压20MPa，铺层全为环向，S33表示σ3)本发明解析法计算结果与有限元计算结果的对比示意图；
图5为本发明所提供的一种导弹发射筒三向应力确定系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种导弹发射筒三向应力确定方法及系统，能够实现碳纤维增强金属内衬导弹发射筒应力快速计算、校核。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明所提供的一种导弹发射筒三向应力确定方法流程示意图，如图1所示，本发明所提供的一种导弹发射筒三向应力确定方法，包括：
S101，根据导弹发射筒中纤维层及内衬层的材料基本力学性能确定纤维层及内衬层在柱坐标系下的应力与应变关系；所述材料基本力学性能包括：弹性模量以及泊松比；
作为一个具体的实施例，根据材料基本力学性能(包括弹性模量、泊松比等)建立纤维层及内衬层应力与应变的本构方程，依据铺层设计方案(主要是纤维与中心轴的夹角)，根据偏轴应力-应变关系式，将层合板正轴刚度矩阵转换为偏轴刚度矩阵，建立纤维层柱坐标系下的应力与应变的本构方程。
建立单层复合材料(单层板)应力与应变关系：
Figure BDA0003321415040000071
其中，
Figure BDA0003321415040000072
E1，E2，E3分别表示单层板沿纤维1，2，3方向的弹性模量(1方向是指沿着纤维方向，2、3方向是指垂直于纤维方向)，υ12，υ21，υ13，υ31，υ23，υ32分别表示单层板的六个泊松比，刚度矩阵中的各刚度系数取值如下：
Figure BDA0003321415040000073
金属内衬为各向同性材料，其应力-应变关系只是其中的刚度系数稍有不同：
Figure BDA0003321415040000081
其中，E表示金属内衬的弹性模量，μ是内衬的泊松比；
柱坐标系下第k(k＝1，2，3，···，M，M表示总层数)层纤维的应力-应变关系可表示为：
Figure BDA0003321415040000082
其中，α是缠绕角，m＝cosα，n＝sinα；
Figure BDA0003321415040000083
S102，根据纤维层及内衬层在柱坐标系下的应力与应变关系确定纤维层的应力的分布方程和内衬层的应力的分布方程；
S102具体包括：
利用公式
Figure BDA0003321415040000091
确定内衬层的应力的分布方程；
利用公式
Figure BDA0003321415040000092
确定纤维层的应力的分布方程；
其中，
Figure BDA0003321415040000093
和
Figure BDA0003321415040000094
分别为内衬层在柱坐标系下的应力，
Figure BDA0003321415040000095
和
Figure BDA0003321415040000096
分别为纤维层在柱坐标系下的应力，w为系数，C3 k和C4 k为系数，μ为内衬层的泊松比，r为沿内衬层径向的坐标，r0≤r≤r1，r0为内衬层的内径，r1是内衬层的外径，C1和C2为系数，f1 k、
Figure BDA0003321415040000097
以及
Figure BDA0003321415040000098
均为系数，rqk为r的(qk)次方，
Figure BDA0003321415040000099
为r的(-qk-1)次方，qk为纤维层径向位移场中的系数。。
S103，根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力；
S103之前还包括：
利用公式
Figure BDA00033214150400000910
确定内衬层的径向位移场；
利用公式
Figure BDA00033214150400000911
确定第k层纤维层的径向位移场；
利用公式
Figure BDA00033214150400000912
确定边界约束；
其中，εA＝0，
Figure BDA00033214150400000913
为内衬层的径向位移场，
Figure BDA00033214150400000914
为第k层纤维层的径向位移场，Pw为内压，M为纤维铺层总数，rM为第M层纤维层的内径，M为正数，R为纤维层的总的外径，
Figure BDA00033214150400000915
为第M层纤维层的径向的应力，pk为纤维层径向位移场中的系数。
其中，
Figure BDA0003321415040000101
Figure BDA0003321415040000102
S103具体包括：
利用公式
Figure BDA0003321415040000103
Figure BDA0003321415040000104
确定应力分布方程系数的解的方程；
其中，C1，C2，···，C4 M均为系数。
其中，对于未知数C1，C2，···，C4 M前面的2M×2M阶系数矩阵(这里将其命名为矩阵A，其中的元素为Ai,j，表示第i行第j列的元素)，前面三行的系数是固定的，对于第三行后面的任意行第i行而言，分两种情况：
(1)当i为偶数时
Figure BDA0003321415040000105
A矩阵最后一行的系数只有两个Ai,i-1，Ai,i，第i列的其他系数Ai,j＝0。
(2)当i为奇数时
Figure BDA0003321415040000117
第i列的其他系数Ai,j＝0。
将应力分布方程系数的解的方程简化为A·C＝P，只有矩阵C是未知的，求解矩阵C即可求得系数C1，C2，··，C4 M；
S104，根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常。
S104具体包括：
利用公式
Figure BDA0003321415040000111
将纤维层和内衬层在柱坐标系下的三向应力转换为笛卡尔直角坐标系下的三向应力；
其中，
Figure BDA0003321415040000112
和
Figure BDA0003321415040000113
分别为笛卡尔直角坐标系下的三向应力，m＝cosα，n＝sinα，α是第k层的缠绕角。
S104具体包括：
利用公式
Figure BDA0003321415040000114
验证纤维层；
若ζ>1，则纤维层失效；反之，则纤维层正常；
其中，
Figure BDA0003321415040000115
Figure BDA0003321415040000116
Xt为纤维纵向拉伸强度，Xc为纤维纵向压缩强度，Yt为纤维横向拉伸强度，Yc为纤维横向压缩强度。
如图2至图4所示，本发明专利根据复合材料基本性能参数和铺层方案建立每一层纤维的柱坐标系下的本构方程，然后建立各层的应力分布方程和径向位移的分布方程，分布方程中的系数尚未确定，建立方程的连续边界条件，代入分布方程和本构方程，获得分布方程系数的解的方程，解方程获得分布方程中的所有系数，即可确定应力-半径和径向位移-半径之间的关系，将各层的半径代入各层的分布方程，即可求出柱坐标系下的三向应力σr，σθ，σz，再根据应力转换方程即可将柱坐标系下的三向应力转化为笛卡尔直角坐标系下三向应力σ1，σ2，σ3。
本发明专利提出了应力分布方程系数的解的方程，使得求解过程更加迅速，计算过程更加方便，计算效率有所提高。同时，本发明所提出的这种半解析半数值的求解方法比现有的有限元计算应力校核强度具有更高的效率，缩短了产品结构设计的周期。
图5为本发明所提供的一种导弹发射筒三向应力确定系统结构示意图，如图5所示，本发明所提供的一种导弹发射筒三向应力确定系统，包括：
应力与应变关系确定模块501，用于根据导弹发射筒中纤维层及内衬层的材料基本力学性能确定纤维层及内衬层在柱坐标系下的应力与应变关系；所述材料基本力学性能包括：弹性模量以及泊松比；
应力的分布方程确定模块502，用于根据纤维层及内衬层在柱坐标系下的应力与应变关系确定纤维层的应力的分布方程和内衬层的应力的分布方程；
径向位移和三向应力确定模块503，用于根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力；
验证结果确定模块504，用于根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常。
所述应力的分布方程确定模块502具体包括：
内衬层的应力的分布方程确定单元，用于利用公式
Figure BDA0003321415040000131
确定内衬层的应力的分布方程；
纤维层的应力的分布方程确定单元，用于利用公式
Figure BDA0003321415040000136
确定纤维层的应力的分布方程；
其中，
Figure BDA0003321415040000137
和
Figure BDA0003321415040000138
分别为内衬层在柱坐标系下的应力，
Figure BDA0003321415040000139
和
Figure BDA0003321415040000132
分别为纤维层在柱坐标系下的应力，w为系数，C3 k和C4 k为系数，μ为内衬层的泊松比，r为沿内衬层径向的坐标，r0≤r≤r1，r0为内衬层的内径，r1是内衬层的外径，C1和C2为系数，f1 k、
Figure BDA0003321415040000133
以及
Figure BDA0003321415040000134
均为系数，rqk为r的(qk)次方，
Figure BDA0003321415040000135
为r的(-qk-1)次方，qk为纤维层径向位移场中的系数。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (7)

1.一种导弹发射筒三向应力确定方法，其特征在于，包括：
根据导弹发射筒中纤维层及内衬层的材料基本力学性能确定纤维层及内衬层在柱坐标系下的应力与应变关系；所述材料基本力学性能包括：弹性模量以及泊松比；
根据纤维层及内衬层在柱坐标系下的应力与应变关系确定纤维层的应力的分布方程和内衬层的应力的分布方程；
根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力；
根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常；
所述根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力，具体包括：
利用公式
Figure FDF0000019714260000011
Figure FDF0000019714260000012
确定应力分布方程系数的解的方程；
对于未知数C1，C2，…，C4 M前面的2M×2M阶系数矩阵，将2M×2M阶系数矩阵命名为矩阵A，其中的元素为Ai,j，表示第i行第j列的元素，前面三行的系数是固定的，对于第三行后面的任意行第i行而言，分两种情况：
(1)当i为偶数时
Figure FDF0000019714260000021
A矩阵最后一行的系数只有两个Ai,i-1，Ai,i，第i列的其他系数Ai,j＝0；
(2)当i为奇数时
Figure FDF0000019714260000022
第i列的其他系数Ai,j＝0；
将应力分布方程系数的解的方程简化为A·C＝P，只有矩阵C是未知的，求解矩阵C即可求得系数C1，C2，··，C4 M；
其中，μ为内衬层的泊松比，r为沿内衬层径向的坐标，r0≤r≤r1，r0为内衬层的内径，w为系数，r1是内衬层的外径，Pw为内压，M为纤维铺层总数，rM为第M层纤维层的内径，M为正数，C1，C2，…，C4 M均为系数，f1 k、
Figure FDF0000019714260000023
以及均为系数；
根据复合材料基本性能参数和铺层方案建立每一层纤维的柱坐标系下的本构方程，然后建立各层的应力分布方程和径向位移的分布方程，分布方程中的系数尚未确定，建立方程的连续边界条件，代入分布方程和本构方程，获得分布方程系数的解的方程，解方程获得分布方程中的所有系数，即可确定应力-半径和径向位移-半径之间的关系，将各层的半径代入各层的分布方程，即可求出柱坐标系下的三向应力σr，σθ，σz，再根据应力转换方程即可将柱坐标系下的三向应力转化为笛卡尔直角坐标系下三向应力σ1，σ2，σ3。
2.根据权利要求1所述的一种导弹发射筒三向应力确定方法，其特征在于，所述根据纤维层及内衬层在柱坐标系下的应力与应变关系确定纤维层的应力的分布方程和内衬层的应力的分布方程，具体包括：
利用公式
Figure FDF0000019714260000031
确定内衬层的应力的分布方程；
利用公式
Figure FDF0000019714260000032
确定纤维层的应力的分布方程；
其中，
Figure FDF0000019714260000033
和
Figure FDF0000019714260000034
分别为内衬层在柱坐标系下的应力，
Figure FDF0000019714260000035
和
Figure FDF0000019714260000036
分别为纤维层在柱坐标系下的应力，w为系数，C3 k和C4 k为系数，μ为内衬层的泊松比，r为沿内衬层径向的坐标，r0≤r≤r1，r0为内衬层的内径，r1是内衬层的外径，C1和C2为系数，f1 k、f2 k、f3 k、
Figure FDF0000019714260000037
以及
Figure FDF0000019714260000038
均为系数，rqk为r的(qk)次方，
Figure FDF0000019714260000039
为r的(-qk-1)次方，qk为纤维层径向位移场中的系数。
3.根据权利要求2所述的一种导弹发射筒三向应力确定方法，其特征在于，所述根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力，之前还包括：
利用公式
Figure FDF00000197142600000310
确定内衬层的径向位移场；
利用公式
Figure FDF00000197142600000311
确定第k层纤维层的径向位移场；
利用公式
Figure FDF0000019714260000041
确定边界约束；
其中，εA＝0，
Figure FDF0000019714260000042
为内衬层的径向位移场，
Figure FDF0000019714260000043
为第k层纤维层的径向位移场，Pw为内压，M为纤维铺层总数，rM为第M层纤维层的内径，M为正数，R为纤维层的总的外径，
Figure FDF0000019714260000044
为第M层纤维层的径向的应力，pk为纤维层径向位移场中的系数。
4.根据权利要求3所述的一种导弹发射筒三向应力确定方法，其特征在于，所述根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常，具体包括：
利用公式
Figure FDF0000019714260000045
将纤维层和内衬层在柱坐标系下的三向应力转换为笛卡尔直角坐标系下的三向应力；
其中，
Figure FDF0000019714260000046
和
Figure FDF0000019714260000047
分别为笛卡尔直角坐标系下的三向应力，m＝cosα，n＝sinα，α是第k层的缠绕角。
5.根据权利要求4所述的一种导弹发射筒三向应力确定方法，其特征在于，所述根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常，具体包括：
利用公式
Figure FDF0000019714260000048
验证纤维层；
若ζ>1，则纤维层失效；反之，则纤维层正常；
其中，
Figure FDF0000019714260000051
Figure FDF0000019714260000052
Xt为纤维纵向拉伸强度，Xc为纤维纵向压缩强度，Yt为纤维横向拉伸强度，Yc为纤维横向压缩强度。
6.一种导弹发射筒三向应力确定系统，用于实现权利要求1-5任意一项所述的一种导弹发射筒三向应力确定方法，其特征在于，包括：
应力与应变关系确定模块，用于根据导弹发射筒中纤维层及内衬层的材料基本力学性能确定纤维层及内衬层在柱坐标系下的应力与应变关系；所述材料基本力学性能包括：弹性模量以及泊松比；
应力的分布方程确定模块，用于根据纤维层及内衬层在柱坐标系下的应力与应变关系确定纤维层的应力的分布方程和内衬层的应力的分布方程；
径向位移和三向应力确定模块，用于根据纤维层的径向位移场、内衬层的径向位移场、导弹发射筒中纤维层和内衬层的边界约束以及纤维层的应力的分布方程和内衬层的应力的分布方程确定纤维层和内衬层在柱坐标系下的径向位移和三向应力；
验证结果确定模块，用于根据纤维层和内衬层在柱坐标系下的径向位移和三向应力，采用蔡吴失效准则验证纤维层，确定验证结果；所述验证结果包括：纤维层失效或纤维层正常。
7.根据权利要求6所述的一种导弹发射筒三向应力确定系统，其特征在于，所述应力的分布方程确定模块具体包括：
内衬层的应力的分布方程确定单元，用于利用公式
Figure FDF0000019714260000053
确定内衬层的应力的分布方程；
纤维层的应力的分布方程确定单元，用于利用公式
Figure FDF0000019714260000054
确定纤维层的应力的分布方程；
其中，
Figure FDF0000019714260000061
和
Figure FDF0000019714260000062
分别为内衬层在柱坐标系下的应力，
Figure FDF0000019714260000063
和
Figure FDF0000019714260000064
分别为纤维层在柱坐标系下的应力，w为系数，C3 k和C4 k为系数，μ为内衬层的泊松比，r为沿内衬层径向的坐标，r0≤r≤r1，r0为内衬层的内径，r1是内衬层的外径，C1和C2为系数，f1 k、f2 k、f3 k、
Figure FDF0000019714260000065
以及
Figure FDF0000019714260000066
均为系数，rqk为r的(qk)次方，
Figure FDF0000019714260000067
为r的(-qk-1)次方，qk为纤维层径向位移场中的系数。
CN202111247041.7A 2021-10-26 2021-10-26 一种导弹发射筒三向应力确定方法及系统 Active CN113962096B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111247041.7A CN113962096B (zh) 2021-10-26 2021-10-26 一种导弹发射筒三向应力确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111247041.7A CN113962096B (zh) 2021-10-26 2021-10-26 一种导弹发射筒三向应力确定方法及系统

Publications (2)

Publication Number Publication Date
CN113962096A CN113962096A (zh) 2022-01-21
CN113962096B true CN113962096B (zh) 2022-11-29

Family

ID=79467112

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111247041.7A Active CN113962096B (zh) 2021-10-26 2021-10-26 一种导弹发射筒三向应力确定方法及系统

Country Status (1)

Country Link
CN (1) CN113962096B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112528415A (zh) * 2020-12-18 2021-03-19 沈阳航空航天大学 一种复合材料轴结构宏-细观失效模式分析方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10347890A1 (de) * 2003-10-10 2005-05-04 Abb Patent Gmbh Magnetisch-induktives Messgerät für strömende Stoffe und Verfahren zu dessen Herstellung
FR2909676B1 (fr) * 2006-12-11 2009-03-20 Astrium Sas Soc Par Actions Si Procede pour ameliorer l'adherence de fibres de carbone vis-a-vis d'une matrice organique
CN112036080B (zh) * 2020-08-19 2023-08-22 江苏科技大学 一种复合材料结构参数的确定方法
CN112297464B (zh) * 2020-09-30 2021-08-27 东南大学 一种基于纤维缠绕工艺的深冷高压储氢容器的设计方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112528415A (zh) * 2020-12-18 2021-03-19 沈阳航空航天大学 一种复合材料轴结构宏-细观失效模式分析方法

Also Published As

Publication number Publication date
CN113962096A (zh) 2022-01-21

Similar Documents

Publication Publication Date Title
Onder et al. Burst failure load of composite pressure vessels
US11209112B2 (en) Fiber composite system and method for pipe reinforcement
CN111368389B (zh) 一种预测复合材料层合板失效强度的方法
Cho et al. An empirical formula for predicting the collapse strength of composite cylindrical-shell structures under external pressure loads
Wang et al. Predicting failure pressure of reinforced thermoplastic pipes based on theoretical analysis and experiment
Mistry Theoretical investigation into the effect of the winding angle of the fibres on the strength of filament wound GRP pipes subjected to combined external pressure and axial compression
CN111368439A (zh) 一种基于缠绕成型工艺的压力容器的设计方法
Zacharakis et al. Optimal finite element modeling of filament wound CFRP tubes
CN113962096B (zh) 一种导弹发射筒三向应力确定方法及系统
Zhang et al. Buckling of composite shells with a novel initial imperfection model subjected to hydrostatic pressure
CN112163273B (zh) 复合材料梯形波纹夹芯圆柱壳的多尺度等效建模方法
CN112460464B (zh) 一种基于纤维缠绕工艺的高压储氢容器的设计方法
CN113111517A (zh) 一种非等极孔纤维缠绕压力容器的设计方法
Silva et al. On the effect of the ply stacking sequence on the failure of composite pipes under external pressure
Özaslan et al. Design and validation of a filament wound composite rocket motor case
CN105046076A (zh) 层叠碳纤维复合材料的三层单胞结构微观力学性能计算方法
CN112733341B (zh) 一种海洋纤维增强柔性管的渐进失效分析方法
US11209113B2 (en) Fiber composite system and method for pipe reinforcement
CN112297464B (zh) 一种基于纤维缠绕工艺的深冷高压储氢容器的设计方法
Reddy et al. Optimum design and analysis of filament wound composite tubes in pure and combined loading
Gohari et al. Failure strength of thin-walled cylindrical GFRP composite shell against static internal and external pressure for various volumetric fiber fraction
Tripathy et al. Stiffened composite axisymmetric shells—Optimum lay-up for buckling by ranking
Nazirah et al. Effects of elevated temperatures on glass-reinforced epoxy pipes under multi-axial loadings
Qin et al. Strength and damage analysis of composite tube with metal lining
Petrich et al. Behaviour of FRP-sandwich structures for lightweight composite springs in static and cyclic torsional load cases

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