CN107451360A - 复合材料直角转弯连续铺层的数值模拟方法 - Google Patents
复合材料直角转弯连续铺层的数值模拟方法 Download PDFInfo
- Publication number
- CN107451360A CN107451360A CN201710633097.3A CN201710633097A CN107451360A CN 107451360 A CN107451360 A CN 107451360A CN 201710633097 A CN201710633097 A CN 201710633097A CN 107451360 A CN107451360 A CN 107451360A
- Authority
- CN
- China
- Prior art keywords
- edge strip
- laying
- composite
- strip area
- area
- 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
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及飞机结构设计技术领域,具体提供了一种复合材料直角转弯连续铺层的数值模拟方法,首先建立复合材料3D壳几何模型,包括水平缘条区、转角区和垂直缘条区,然后建立水平缘条区直角坐标系和垂直缘条区直角坐标系,并建立转角区柱坐标系,再分别定义水平缘条区、转角区和垂直缘条区的复合材料铺层属性,并根据法向分别调试上述三个区的铺层基准参考方向,选取上述三个区铺层对应坐标系,并保证上述三个区域的铺层连续且同一铺层的角度一致,然后对复合材料3D壳几何模型进行网格划分,对所述几何模型施加约束和载荷,得到有限元模型,最后求出所述有限元模型的静强度场变量。
Description
技术领域
本发明涉及飞机结构设计技术领域,特别涉及复合材料直角转弯连续铺层的数值模拟方法。
背景技术
复合材料层压板结构相互垂直缘板交界处的R区结构特征,广泛存在于各种飞机复合材料零部件中,如长桁、肋及梁等。在R区预浸料存在转折铺设,坐标系将发生变化,在一般的总体有限元分析中,各独立面参考各自的直角坐标系,未考虑转角及铺层连续问题,其输出结果不能反应真实结构的受力情况,如某型飞机复合材料翼肋由于缘条和腹板转角区提早出现分层而导致放弃采用复合材料翼肋方案。
发明内容
为克服上述现有技术存在的缺陷,本发明提供了一种复合材料直角转弯连续铺层的数值模拟方法,包括以下步骤:
步骤一,建立复合材料3D壳几何模型,所述壳几何模型包括水平缘条区、转角区和垂直缘条区,水平缘条区对应水平缘条面集合,转角区对应转角面集合,垂直缘条区对应垂直缘条面集合;
步骤二,建立水平缘条区直角坐标系和垂直缘条区直角坐标系,并建立转角区柱坐标系;
步骤三,分别定义水平缘条区、转角区和垂直缘条区的复合材料铺层属性,并根据法向分别调试上述三个区的铺层基准参考方向,选取上述三个区铺层对应坐标系,并保证上述三个区域的铺层连续且同一铺层的角度一致;
步骤四,对复合材料3D壳几何模型采用八节点减缩壳单元类型进行网格划分,对所述几何模型施加约束和载荷,得到有限元模型;
步骤五,求出所述有限元模型的静强度场变量。
优选的,所述复合材料为预浸料,所述复合材料铺层属性包括材料力学性能、铺层角度、铺层顺序。
优选的,所述静强度场变量包括铺层应力、应变、位移。
本发明提供的复合材料直角转弯连续铺层的数值模拟方法,将R区连接区域进行分割,并分别建立各自区域的局部坐标系,实现了复合材料铺层转弯后的连续性,将复合材料转折铺层连接成整体铺层,该建模方法为复合材料空间结构及R区的细节分析提供了可靠实用的分析方法,可算出整个铺层的应力应变等结果。
附图说明
图1是复合材料3D壳几何模型的示意图;
图2是复合材料的坐标系示意图;
图3是复合材料的有限元网格示意图;
图4是复合材料的45°铺层剪应力云图。
附图标记:垂直缘条区1,R转角区2,水平缘条区3。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行更加详细的描述。
需要说明的是:下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例,而不是全部的实施例,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,均仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明保护范围的限制。
本发明提供的复合材料直角转弯连续铺层的数值模拟方法,旨在建立垂直缘板R区铺层连续铺设下的细节模型,为该细节结构特征的力学行为(如静强度、稳定性及损伤容限)评估提供了一种数值模拟方法,该方法的具体实施步骤如下:
步骤一,如图1所示,建立复合材料3D壳几何模型,所述壳几何模型包括水平缘条区、转角区和垂直缘条区,水平缘条区对应水平缘条面集合,转角区对应转角面集合,垂直缘条区对应垂直缘条面集合,本实施例中优选的是,所述复合材料为预浸料。
步骤二,如图2所示,建立水平缘条区直角坐标系和垂直缘条区直角坐标系,并建立转角区柱坐标系。
步骤三,分别定义水平缘条区、转角区和垂直缘条区的复合材料铺层属性,并根据法向分别调试上述三个区的铺层基准参考方向,选取上述三个区铺层的对应坐标系,并保证上述三个区域的铺层连续且同一铺层的角度一致,本实施例中优选的是,所述复合材料铺层属性包括材料力学性能、铺层角度、铺层顺序等属性。
步骤四,如图3所示,对复合材料3D壳几何模型采用八节点减缩壳单元类型进行网格划分,对所述几何模型施加约束和载荷,得到有限元模型,本实施例中优选的是,所述有限元模型采用Abaqus软件实现,本实施例中:
加载及约束:一端顶点加载F=10N,另一边固支;
层板铺层:[45/0/90-45];
单向预浸料材料性能:E1=136GPa,E2=E3=9.8GPa,G12=G13=4.7GPa,G23=4.261GPa,v12=v13=0.28,v23=0.15。
步骤五,如图4所示,求出所述有限元模型的静强度场变量,图中的编码为单元编号,本实施例中优选的是,所述静强度场变量包括铺层应力、应变、位移等变量。
复合材料的45°铺层剪应力数据如表1所示:
表1复合材料的45°铺层剪应力数据
单元 | τ12(MPa) |
264 | -15.8377 |
361 | -10.1633 |
362 | -11.8414 |
363 | -11.4608 |
620 | -19.6173 |
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (3)
1.一种复合材料直角转弯连续铺层的数值模拟方法,其特征在于:包括以下步骤:
步骤一,建立复合材料3D壳几何模型,所述壳几何模型包括水平缘条区、转角区和垂直缘条区,水平缘条区对应水平缘条面集合,转角区对应转角面集合,垂直缘条区对应垂直缘条面集合;
步骤二,建立水平缘条区直角坐标系和垂直缘条区直角坐标系,并建立转角区柱坐标系;
步骤三,分别定义水平缘条区、转角区和垂直缘条区的复合材料铺层属性,并根据法向分别调试上述三个区的铺层基准参考方向,选取上述三个区铺层对应坐标系,并保证上述三个区域的铺层连续且同一铺层的角度一致;
步骤四,对复合材料3D壳几何模型采用八节点减缩壳单元类型进行网格划分,对所述几何模型施加约束和载荷,得到有限元模型;
步骤五,求出所述有限元模型的静强度场变量。
2.根据权利要求1所述的数值模拟方法,其特征在于:所述复合材料为预浸料,所述复合材料铺层属性包括材料力学性能、铺层角度、铺层顺序。
3.根据权利要求1所述的数值模拟方法,其特征在于:所述静强度场变量包括铺层应力、应变、位移。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710633097.3A CN107451360B (zh) | 2017-07-28 | 2017-07-28 | 复合材料直角转弯连续铺层的数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710633097.3A CN107451360B (zh) | 2017-07-28 | 2017-07-28 | 复合材料直角转弯连续铺层的数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107451360A true CN107451360A (zh) | 2017-12-08 |
CN107451360B CN107451360B (zh) | 2021-03-09 |
Family
ID=60490572
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710633097.3A Active CN107451360B (zh) | 2017-07-28 | 2017-07-28 | 复合材料直角转弯连续铺层的数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107451360B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583040A (zh) * | 2018-11-07 | 2019-04-05 | 中国航空工业集团公司西安飞机设计研究所 | 一种考虑复合材料结构参数连续性的优化方法 |
CN111976168A (zh) * | 2020-08-12 | 2020-11-24 | 常州启赋安泰复合材料科技有限公司 | 一种复合材料进气道的下料方法 |
CN113779850A (zh) * | 2021-09-28 | 2021-12-10 | 上海航空工业(集团)有限公司 | 复合材料结构模型的检测方法、检测装置、设备和介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103605865A (zh) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | 一种基于内聚力模型的复合材料多向铺层板分层扩展行为模拟方法 |
CN104268317A (zh) * | 2014-09-12 | 2015-01-07 | 西北工业大学 | 机械零件圆角结构形状优化方法 |
-
2017
- 2017-07-28 CN CN201710633097.3A patent/CN107451360B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103605865A (zh) * | 2013-11-28 | 2014-02-26 | 北京航空航天大学 | 一种基于内聚力模型的复合材料多向铺层板分层扩展行为模拟方法 |
CN104268317A (zh) * | 2014-09-12 | 2015-01-07 | 西北工业大学 | 机械零件圆角结构形状优化方法 |
Non-Patent Citations (1)
Title |
---|
文友谊 等: "L型复合材料结构固化变形有限元模拟及其应用", 《PROCEEDINGS OF THE 11TH CHINA-JAPAN JOINT CONFERENCE ON COMPOSITE MATERIALS》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109583040A (zh) * | 2018-11-07 | 2019-04-05 | 中国航空工业集团公司西安飞机设计研究所 | 一种考虑复合材料结构参数连续性的优化方法 |
CN109583040B (zh) * | 2018-11-07 | 2022-11-18 | 中国航空工业集团公司西安飞机设计研究所 | 一种考虑复合材料结构参数连续性的优化方法 |
CN111976168A (zh) * | 2020-08-12 | 2020-11-24 | 常州启赋安泰复合材料科技有限公司 | 一种复合材料进气道的下料方法 |
CN113779850A (zh) * | 2021-09-28 | 2021-12-10 | 上海航空工业(集团)有限公司 | 复合材料结构模型的检测方法、检测装置、设备和介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107451360B (zh) | 2021-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107609268B (zh) | 一种基于tin的三维地质模型与bim模型耦合的方法 | |
CN107451360A (zh) | 复合材料直角转弯连续铺层的数值模拟方法 | |
US7406403B2 (en) | Systems and methods for automatically generating 3D wireframe CAD models of aircraft | |
Peeters et al. | High-fidelity finite element models of composite wind turbine blades with shell and solid elements | |
US8762110B2 (en) | Interactive truss objects in computer-aided design models | |
CN103412996B (zh) | 在VA-One中快速建立声学部件模型的方法 | |
KR20120082112A (ko) | 복합재 구조 분석 시스템과 분석방법 | |
US20130325417A1 (en) | Numerical Simulation Of A Structure Having A Heat-Affected Zone Using A Finite Element Analysis Model | |
Loseille et al. | Boundary layer mesh generation and adaptivity | |
CN104346828A (zh) | 一种Tribon M3零件级船体三维模型重构方法 | |
Kussmaul et al. | A novel computational framework for structural optimization with patched laminates | |
Längst et al. | The potentials of isogeometric analysis methods in integrated design processes | |
CN117217049A (zh) | 一种基于PloughMesh的有限元分析前处理方法和系统 | |
CN115238419A (zh) | 一种基于复合材料的轴流压气机叶片设计方法 | |
CN104899381B (zh) | 一种焊接箱型截面钢节点多尺度有限元建模方法 | |
KR101696330B1 (ko) | Am 도면 환경에서의 3d 의장 모델링 장치 및 그것을 이용한 3d 의장 모델링 방법 | |
CN105138750A (zh) | 一种可扩展多线交叉机身线框模型构建方法 | |
US20140222395A1 (en) | Contact surface definition creation involving low order and quadratic finite elements in a numerical simulation of an impact event | |
Autry et al. | Automated top level aircraft structural sizing tool (Atlass): A framework for preliminary aircraft design and optimization | |
US20190258762A1 (en) | Methods And Systems For Manufacturing Products/Parts Made Of Carbon Fiber Reinforced Composite Based On Numerical Simulations | |
Vijayachandran et al. | Steered fiber paths for improved in-plane compressive response of aerostructural panels: experimental studies and numerical modeling | |
Li et al. | Automated Generation of Finite-Element Meshes for Aircraft Conceptual Design | |
KR20060027477A (ko) | 선박의 선체 상세 설계 시점의 컴퓨터 에디드 디자인모델을 이용한 선박 선체 상세 해석을 위한 컴퓨터 에디드엔지니어링 모델의 생성 시스템 및 그 방법 | |
Peeters et al. | Development of automated high fidelity finite element models for large wind turbine blades | |
Wang et al. | Automated hex meshing for turbomachinery secondary air system |
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 |