CN103678786A - 一种运载火箭蒙皮桁条结构的有限元分析方法 - Google Patents
一种运载火箭蒙皮桁条结构的有限元分析方法 Download PDFInfo
- Publication number
- CN103678786A CN103678786A CN201310629675.8A CN201310629675A CN103678786A CN 103678786 A CN103678786 A CN 103678786A CN 201310629675 A CN201310629675 A CN 201310629675A CN 103678786 A CN103678786 A CN 103678786A
- Authority
- CN
- China
- Prior art keywords
- model
- stringer
- finite element
- ring frame
- covering
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种运载火箭蒙皮桁条结构的有限元分析方法,可用于结构强度快速分析,也可以用于快速建模,属于结构强度分析技术领域。该方法将蒙皮桁条结构的强度分析过程标准化,积累材料数据和结构参数数据,利用分析标准和数据库的效率优势,提高工作效率和分析质量。
Description
技术领域
本发明涉及一种运载火箭蒙皮桁条结构的有限元分析方法,可用于结构强度快速分析,也可以用于快速建模,属于结构强度分析技术领域。
背景技术
蒙皮桁条结构具有成本低、工艺简单、承载能力高等特点,有利于结构轻量化,是国内外运载火箭普遍采用的结构形式。目前蒙皮桁条结构的强度分析有以下两种方法:(1)采用Pro/E等CAD工具建立几何模型,然后通过中间格式文件导入到CAE软件,然后进行赋属性、划网格、加载荷及边界条件、选择求解算法、提交运算等常规有限元分析工作。(2)直接采用CAE软件建立几何模型,然后进行常规有限元分析工作。这两种方法都是CAE工具的初级应用,无数据和经验积累,工作高度重复、极度依赖于分析人员经验、容易出错。
发明内容
本发明的目的是为了提出一种运载火箭蒙皮桁条结构的有限元分析方法,该方法将蒙皮桁条结构的强度分析过程标准化,积累材料数据和结构参数数据,利用分析标准和数据库的效率优势,提高工作效率和分析质量。
本发明的目的是通过以下技术方案实现的。
本发明的一种运载火箭蒙皮桁条结构的有限元分析方法,步骤为:
1)对蒙皮建立有限元模型,有限元模型中所使用的参数为蒙皮的材料属性和蒙皮的几何尺寸;
所述的蒙皮的材料属性是指蒙皮材料的弹性模量、泊松比、屈服应力、极限强度和断后伸长率,并保存在数据库中;
所述的蒙皮的几何尺寸是指蒙皮的中面直径、高度和厚度,并保存在数据库中;
2)对桁条建立有限元模型,有限元模型中所使用的参数为桁条的材料属性和桁条的几何尺寸;
所述的桁条的材料属性是指桁条材料的弹性模量、泊松比、屈服应力、极限强度和断后伸长率,并保存在数据库中;
所述的桁条的几何尺寸是指桁条的截面尺寸,并保存在数据库中;
3)对环框建立有限元模型,有限元模型中所使用的参数为环框的材料属性和环框的几何尺寸;
所述的环框的材料属性是指环框材料的弹性模量、泊松比、屈服应力、极限强度和断后伸长率,并保存在数据库中;
所述的环框的几何尺寸是指环框的截面尺寸,并保存在数据库中;
4)将步骤1)中建立的蒙皮模型、步骤2)中建立的桁条模型和步骤3)中建立的环框模型进行装配,得到装配模型;
5)对步骤4)得到的装配模型进行开孔;
6)采用有限元工具对步骤5)得到的装配模型进行单元划分,得到单元模型,得到的单元模型为应力分析模型或稳定性分析模型;
进行单元划分时要求相邻桁条之间不少于5个单元以及相邻的环框之间不少于5个单元;以达到足以描述失稳波和结构变形模式的要求;
7)如果步骤6)得到的单元模型为应力分析模型,则采用完全积分单元;如果步骤6)得到的单元模型为稳定性分析模型,分两种情况,如果在稳定性分析模型的任意方向的单元数目大于4,则采用减缩积分单元,如果在稳定性分析模型的任意方向的单元数目小于等于4,则采用完全积分单元;
8)采用有限元工具在单元模型上施加边界条件和载荷,然后用有限元工具对单元模型进行应力分析或稳定性分析。
对步骤1)中建立的蒙皮模型、步骤2)中建立的桁条模型和步骤3)中建立的环框模型在数据库中可以进行数据积累,供后续分析调用。
有益效果
对蒙皮桁条结构所用材料和型材数据进行积累,建立数据库供后续分析调用,减少了重复性工作,提高工了作效率。
具体实施方式
下面结合实施例对本发明做进一步说明。
实施例1
一种运载火箭蒙皮桁条结构的有限元分析方法,步骤为:
1)对蒙皮建立有限元模型,蒙皮采用2A12T4状态1mm厚铝板材,弹性模量为70GPa、泊松比为0.3、屈服应力为270MPa、极限强度为400MPa、断后伸长率为13%,保存在数据库中,标记为A1;
蒙皮的中面直径为1700mm、高度为2500mm、厚度为1mm,保存在数据库中,标记为B1;
2)对桁条建立有限元模型,桁条采用2A12T4状态铝型材,弹性模量为70GPa、泊松比为0.3、屈服应力为295MPa、极限强度为390MPa、断后伸长率为10%,保存在数据库中,标记为A2;
桁条的截面为“T型”,尺寸有两个规格:规格1,缘板宽度为30mm、厚度为2mm,腹板高度为25mm、厚度为3mm,保存在数据库中,标记为C1;规格2,缘板宽度为20mm、厚度为1.5mm,腹板高度为22mm、厚度为2mm,保存在数据库中,标记为C2;
3)对环框建立有限元模型,环框采用2A12O状态2mm厚铝板材,弹性模量为70GPa、泊松比为0.3、屈服应力为150MPa、极限强度为250MPa、断后伸长率为30%,保存在数据库中,标记为A3;
环框的截面为方框,长度为30mm、宽度为20mm、厚度为2mm,保存在数据库中,标记为D1;
4)将步骤1)中建立的蒙皮模型、步骤2)中建立的桁条模型和步骤3)中建立的环框模型进行装配,得到装配模型;装配模型共有84根桁条,其中截面为C1的桁条有12根,位于蒙皮内表面,从柱坐标环向0度开始等间距分布,每30度一根;其中截面为C2的桁条有72根,位于蒙皮外表面,从柱坐标环向2.5度开始等间距分布,每5度一根;装配模型共有2个环框,位于蒙皮内表面,截面为D1,1个环框距离蒙皮下端800mm,1个环框距离蒙皮下端2000mm。
5)对步骤4)得到的装配模型进行开孔,共有2个孔:1个为圆孔,半径为300mm,圆心在柱坐标中位于环向310度,轴向距离蒙皮下端400mm;另1个为矩形孔,环向宽度为500mm、轴向高度400mm,矩形中心在柱坐标中位于环向110度,轴向距离蒙皮下端1500mm。
6)采用有限元工具对步骤5)得到的装配模型进行单元划分,得到应力分析模型;相邻桁条之间的最少单元数为6个,相邻的环框之间的最少单元数为100个;
7)步骤6)得到的应力分析模型,采用完全积分单元;
8)采用有限元工具将模型的下端固定约束、上端设为刚性面,在刚性面中心施加轴向压力100kN,然后用有限元工具进行应力分析,得到应力分布,最大Mises应力为159MPa。
实施例2
一种运载火箭蒙皮桁条结构的有限元分析方法,步骤为:
1)对蒙皮建立有限元模型,蒙皮材料采用实施例1中A1,属性从数据库中调取;
蒙皮尺寸采用实施例1中B1,尺寸从数据库中调取;
2)对桁条建立有限元模型,桁条采用7A09T4状态铝型材,弹性模量为70GPa、泊松比为0.3、屈服应力为440MPa、极限强度为500MPa、断后伸长率为8%,保存在数据库中,标记为A4;
桁条的截面为“T型”,采用实施例1中C2,尺寸从数据库中调取;
3)对环框建立有限元模型,桁条材料采用实施例1中A3,属性从数据库中调取;
环框的截面为方框,采用实施例1中D1,尺寸从数据库中调取;
4)将步骤1)中建立的蒙皮模型、步骤2)中建立的桁条模型和步骤3)中建立的环框模型进行装配,得到装配模型;装配模型共有60根桁条,截面为C2,位于蒙皮外表面,从柱坐标环向0度开始等间距分布,每6度一根;装配模型共有1个环框,位于蒙皮内表面,截面为D1,距离蒙皮下端1200mm;
5)对步骤4)得到的装配模型进行开孔,共有1个椭圆孔,长轴沿环向、为300mm,短轴沿轴向、为150mm,椭圆中心在柱坐标中位于环向160度,轴向距离蒙皮下端600mm。
6)采用有限元工具对步骤5)得到的装配模型进行单元划分,得到稳定性分析模型;相邻桁条之间的最少单元数为10个,相邻的环框之间的最少单元数为240个;
7)步骤6)得到的稳定性分析模型,蒙皮采用减缩积分单元,桁条和环框采用完全积分单元;
8)采用有限元工具将模型的下端固定约束、上端设为刚性面,在刚性面中心施加轴向压力200kN,然后用有限元工具进行稳定性分析,得到承载能力为186kN,结构整体失稳。
Claims (1)
1.一种运载火箭蒙皮桁条结构的有限元分析方法,其特征在于步骤为:
1)对蒙皮建立有限元模型,有限元模型中所使用的参数为蒙皮的材料属性和蒙皮的几何尺寸;
所述的蒙皮的材料属性是指蒙皮材料的弹性模量、泊松比、屈服应力、极限强度和断后伸长率,并保存在数据库中;
所述的蒙皮的几何尺寸是指蒙皮的中面直径、高度和厚度,并保存在数据库中;
2)对桁条建立有限元模型,有限元模型中所使用的参数为桁条的材料属性和桁条的几何尺寸;
所述的桁条的材料属性是指桁条材料的弹性模量、泊松比、屈服应力、极限强度和断后伸长率,并保存在数据库中;
所述的桁条的几何尺寸是指桁条的截面尺寸,并保存在数据库中;
3)对环框建立有限元模型,有限元模型中所使用的参数为环框的材料属性和环框的几何尺寸;
所述的环框的材料属性是指环框材料的弹性模量、泊松比、屈服应力、极限强度和断后伸长率,并保存在数据库中;
所述的环框的几何尺寸是指环框的截面尺寸,并保存在数据库中;
4)将步骤1)中建立的蒙皮模型、步骤2)中建立的桁条模型和步骤3)中建立的环框模型进行装配,得到装配模型;
5)对步骤4)得到的装配模型进行开孔;
6)采用有限元工具对步骤5)得到的装配模型进行单元划分,得到单元模型,得到的单元模型为应力分析模型或稳定性分析模型;
进行单元划分时要求相邻桁条之间不少于5个单元以及相邻的环框之间不少于5个单元;以达到足以描述失稳波和结构变形模式的要求;
7)如果步骤6)得到的单元模型为应力分析模型,则采用完全积分单元;如果步骤6)得到的单元模型为稳定性分析模型,分两种情况,如果在稳定性分析模型的任意方向的单元数目大于4,则采用减缩积分单元,如果在稳定性分析模型的任意方向的单元数目小于等于4,则采用完全积分单元;
8)采用有限元工具在单元模型上施加边界条件和载荷,然后用有限元工具对单元模型进行应力分析或稳定性分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310629675.8A CN103678786B (zh) | 2013-11-29 | 2013-11-29 | 一种运载火箭蒙皮桁条结构的有限元分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310629675.8A CN103678786B (zh) | 2013-11-29 | 2013-11-29 | 一种运载火箭蒙皮桁条结构的有限元分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103678786A true CN103678786A (zh) | 2014-03-26 |
CN103678786B CN103678786B (zh) | 2016-06-29 |
Family
ID=50316322
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310629675.8A Active CN103678786B (zh) | 2013-11-29 | 2013-11-29 | 一种运载火箭蒙皮桁条结构的有限元分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103678786B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156498A (zh) * | 2014-05-26 | 2014-11-19 | 北京宇航系统工程研究所 | 一种螺栓的有限元分析方法及强度评估方法 |
CN104408238A (zh) * | 2014-11-11 | 2015-03-11 | 上海宇航系统工程研究所 | 一种桁式半硬壳结构轴压稳定性的数值分析方法 |
CN104732010A (zh) * | 2015-02-06 | 2015-06-24 | 北京宇航系统工程研究所 | 一种多层防热结构快速优化方法 |
CN107451335A (zh) * | 2017-07-06 | 2017-12-08 | 华东建筑设计研究院有限公司 | 一种引入蒙皮效应的铝合金网格结构整体稳定分析方法 |
CN107463716A (zh) * | 2016-06-02 | 2017-12-12 | 宋扬 | 一种地物光谱剖面测量方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103324782A (zh) * | 2013-05-29 | 2013-09-25 | 北京航空航天大学 | 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 |
-
2013
- 2013-11-29 CN CN201310629675.8A patent/CN103678786B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103324782A (zh) * | 2013-05-29 | 2013-09-25 | 北京航空航天大学 | 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 |
Non-Patent Citations (4)
Title |
---|
PAN ZHONGWEN等: "Liquid propellant analogy technique in dynamic modeling of launch vehicle", 《SCIENCE CHINA TECHNOLOGICAL SCIENCES》 * |
张希: "基于显式有限元方法的蒙皮桁条结构稳定性分析方法", 《北京力学会第19届学术年会论文集》 * |
潘忠文等: "运载火箭结构动力学模拟技术研究进展", 《力学进展》 * |
邢誉峰等: "蒙皮加筋圆柱壳弯曲频率的三种计算模型", 《北京航空航天大学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104156498A (zh) * | 2014-05-26 | 2014-11-19 | 北京宇航系统工程研究所 | 一种螺栓的有限元分析方法及强度评估方法 |
CN104408238A (zh) * | 2014-11-11 | 2015-03-11 | 上海宇航系统工程研究所 | 一种桁式半硬壳结构轴压稳定性的数值分析方法 |
CN104408238B (zh) * | 2014-11-11 | 2017-06-06 | 上海宇航系统工程研究所 | 一种桁式半硬壳结构轴压稳定性的数值分析方法 |
CN104732010A (zh) * | 2015-02-06 | 2015-06-24 | 北京宇航系统工程研究所 | 一种多层防热结构快速优化方法 |
CN104732010B (zh) * | 2015-02-06 | 2017-09-26 | 北京宇航系统工程研究所 | 一种多层防热结构快速优化方法 |
CN107463716A (zh) * | 2016-06-02 | 2017-12-12 | 宋扬 | 一种地物光谱剖面测量方法 |
CN107451335A (zh) * | 2017-07-06 | 2017-12-08 | 华东建筑设计研究院有限公司 | 一种引入蒙皮效应的铝合金网格结构整体稳定分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103678786B (zh) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103678786A (zh) | 一种运载火箭蒙皮桁条结构的有限元分析方法 | |
CN204676435U (zh) | 一种高墩混凝土施工可伸缩操作平台 | |
CN204954392U (zh) | 一种建筑管材切割固定装置 | |
CN203141018U (zh) | 工字钢固定工作支架 | |
CN101774188A (zh) | 气压冲床 | |
CN204826714U (zh) | 一种新型钢构支撑胎架 | |
CN204088923U (zh) | 模块化支撑框架 | |
CN104668873A (zh) | 钢筋笼加工模具 | |
CN103382740A (zh) | 带可调节撑杆的弦支门式刚架结构 | |
CN102864746A (zh) | 一种桥梁施工用承重架 | |
CN203381215U (zh) | 双动液压机 | |
CN203901256U (zh) | 一种可调式钢结构加工胎架 | |
CN102151756B (zh) | 刀片成型模具上下冲头的组合结构 | |
CN204564905U (zh) | 组合式冲孔模具 | |
CN203022343U (zh) | 工具式方钢管紧固件 | |
CN203755697U (zh) | 一种钢筋网片样架装置 | |
CN203887537U (zh) | 一种支撑托架 | |
CN202767506U (zh) | 建筑模板施工架体承载钢横杆件 | |
CN203765766U (zh) | 一种核电站生物防护砌块模板 | |
CN204564894U (zh) | 一种铝型材冲孔机 | |
CN205688777U (zh) | 悬挑脚手架的安全防护结构 | |
CN204876088U (zh) | 箱形柱-工形梁预制装配式高延性梁柱节点 | |
CN204156028U (zh) | 用于电池外壳的气缸式限位折边机构 | |
CN103586629A (zh) | 挖掘机驾驶室用异型管材的加工方法 | |
CN202658881U (zh) | 大型简易定型剪力墙竹胶模板 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |