CN112052541B - 一种任意形状的钣金型材截面参数化方法 - Google Patents
一种任意形状的钣金型材截面参数化方法 Download PDFInfo
- Publication number
- CN112052541B CN112052541B CN202011028772.8A CN202011028772A CN112052541B CN 112052541 B CN112052541 B CN 112052541B CN 202011028772 A CN202011028772 A CN 202011028772A CN 112052541 B CN112052541 B CN 112052541B
- Authority
- CN
- China
- Prior art keywords
- section
- theta
- angle
- coordinate
- sheet metal
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
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)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Image Generation (AREA)
- Image Analysis (AREA)
Abstract
本发明属于结构动力学技术领域,公开了一种任意形状的钣金型材截面参数化方法。采用矢量追踪的方法和模块化思想,从原点出发对截面轮廓的相邻点依次求解,实现了对带R角的任意形状钣金型材截面的参数化定义,可用于模拟真实非规则截面凸缘的有限元建模以及几何参数的计算,具有很高的普适性。
Description
技术领域
本发明属于结构动力学技术领域,特别涉及一种任意形状的钣金型材截面参数化方法。
背景技术
作为直升机结构中最常用的构型形式,钣金型材常被用于框梁凸缘、平台长绗、加强立筋等结构,与蒙皮组成的壁板加筋结构也是传统航空器高结构效率应用最典型的构型。在强度分析时,型材的截面参数是个关键的设计输入,有限元分析时不可忽略,特别是动力学分析,截面形状的差异对模型的刚度和重量分布具有重要影响,因此,为了模拟真实结构,有限元建模时需要对型材的截面形状进行精确定义。
目前主流有限元商用软件中,截面的定义多采用参数化的方法,通过给定几个常用规则截面,用户可以输入长度、厚度等参数,用于计算截面的所有轮廓点坐标和属性参数。但该方法存在一定缺陷,比如无法考虑R角影响,另外,模块化的参数化方法局限于规则截面,对于航空结构上大量的非规则截面,尚无普适的参数化方法。
发明内容
本发明的目的:提供一种考虑R角的任意形状钣金型材截面参数化方法。
本发明的技术方案:
一种任意形状的钣金型材截面参数化方法,包括以下步骤:
步骤一:确定和输入截面参数;
步骤二,计算R角圆心及角度;
步骤三,计算所有轮廓点坐标;
步骤四,计算截面属性参数。
进一步,所述步骤一中,需要确定和输入的截面参数包括:截面的边数N、钣金的厚度t、钣金的R角半径r、边的长度L1~LN、相邻边之间的角度θ1~θN(|θN|<180°)。
进一步,所述需要确定和输入的截面参数还包括旋转角度θt。
进一步,所述步骤二中R角圆心及角度计算过程如下:
步骤a),选择截面的任意一个起始边的自由端点为原点,起始边为X横轴,建立局部坐标系,起始点坐标为(0,0),
判断起始边与第一个相邻边之间的角度值θ1,如果θ1>0,则以起始边和第一个相邻边的交点为圆心,从起始边位置逆时针旋转θ1得到第一个相邻边的位置;反之,则顺时针旋转θ1得到第一个相邻边的位置;
第一个R角的圆心坐标为C1(L1,-r)(θ1>0)或C1(L1,t+r)(θ1<0);
并构建第一个R角坐标点参数集[C1,r,θstart,θend],θstart=90°(θ1>0)或者θstart=270°(θ1<0),θend=θstart–(180°-θ1)(θ1>0)或者θend=θstart+(180°+θ1)(θ1<0),其中θstart和θend分别为R角轮廓弧线首尾端点与R角圆心连线相对于X横轴的偏转角度;
建立一个幅值为L2-2*t,方向为(180°+θ1)的方向矢量V1,用于求解下一个圆心坐标。
步骤b),将第一个R角圆心坐标按矢量V1平移,得到第二个R角的初始圆心坐标;
判断第二个相邻角θ2,如果θ1*θ2>0时,则以第二个R角的初始圆心坐标作为最终的第二个R角圆心坐标。
如果θ1*θ2<0时,则构建一个幅值为t+2r的转置矢量V12,若θ2<0,则方向为V1+90°,反之,方向为V1-90°。将初始圆心按V12平移,得到第二个R角圆心坐标C2;
构建第二个圆角轮廓点坐标参数集[C2,r,θstart,θend],θstart=180°+θ1+90°(θ2>0)或者θstart=180°+θ1-90°(θ2<0),θend=θstart–(180°-θ2)(θ2>0)或者θend=θstart+(180°+θ2)(θ2<0);
构建一个幅值为L3-2*t,方向为[(180°+θ1)-180°+θ2]的方向矢量V2,用于求解下一个圆心坐标;
步骤c),重复步骤b)依次计算剩余所有圆心坐标并构建对应的圆角轮廓点坐标参数集[CN,r,θstart,θend]。
进一步,在每段圆角轮廓弧线上等分选取5个弧线坐标点,根据每个圆角的轮廓点坐标集的输入参数集,计算每个弧线坐标点的坐标p1~p5,得到弧线坐标点集。
进一步,对于同一个R角圆心,存在内外两个轮廓弧线;
对于外轮廓弧线,CNout参数集为[CN,t+r,θstart,θend];对于内轮廓弧线,CNin参数集为[CN,r,θend,θstart];
建立轮廓点顺序,从原点出发,对截面的所有轮廓点进行排序,按[S1,C1in/C1out,C2in/C2out,……,CNin/CNout,E1,E2,CNout/CNin,……,C2out/C2in,C1out/C1in,S2,S1]的顺序进行排列,其中,S1、S2、E1、E2为截面的四个端点,
组成一个个数为(10*N+5)的坐标集,形成最终的(10*N+5)×2的轮廓点坐标总矩阵。
进一步,若截面参数包括旋转角度θt,则轮廓点坐标总矩阵进行转置,转置矩阵为:[cosθt sinθt;-sinθt cosθt]。
进一步,根据轮廓点坐标总矩阵,采用数学积分的方法,计算截面属性,所述截面属性包括:截面的面积、形心、惯性矩。
本发明的有益技术效果:本发明通过采用矢量追踪的方法和模块化思想,从原点出发对截面轮廓的相邻点依次求解,实现了对带R角的任意形状钣金型材截面的参数化定义,可用于模拟真实非规则截面凸缘的有限元建模以及几何参数的计算,具有很高的普适性。
附图说明
图1是参数化方法流程图;
图2是圆角坐标计算流程图;
图3是截面轮廓点坐标计算示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,一种任意形状的钣金型材截面参数化方法,包括以下步骤:
步骤一:确定和输入截面参数。
对于任意形状的钣金型材截面,需要确定和输入该截面的以下参数,1)截面的边数N;2)钣金的厚度t;3)钣金的R角半径r;4)边的长度L1~LN;5)相邻边之间的角度θ1~θN(|θN|<180°);6)旋转角度θt(必要时)。
步骤二:计算R角圆心及角度。
如图2所示,选择截面的任意一个起始边下端点为原点,起始边为X横轴,建立局部坐标系,起始点坐标为(0,0),然后计算第一个R角的圆心坐标,首先判断输入参数第一个相邻边之间的角度值θ1,如果θ1>0,则表示第二条边是以起始边的另一个端点为圆心,从第一条边位置逆时针旋转θ1后的位置,反之,则绕顺时针旋转θ1。然后根据θ1的值,可以计算第一个R角的圆心坐标C1(L1,-r)(θ1>0)或C1(L1,t+r)(θ1<0),同时构建第一个圆角轮廓点坐标参数集[C1,r,θstart,θend],θstart=90°(θ1>0)或者θstart=270°(θ1<0),θend=θstart–(180°-θ1)(θ1>0)或者θend=θstart+(180°+θ1)(θ1<0),其中θstart和θend分别为R角轮廓弧线首尾端点与R角圆心连线相对于X横轴的偏转角度。最后,建立一个幅值为L2-2*t,方向为(180°+θ1)的方向矢量V1,用于求解下一个圆心坐标。
同理,将第1个R角圆心坐标沿矢量V1方向平移L2-2*t长度,得到第二个R角的初始圆心坐标,然后判断第二个相邻角θ2,如果
θ1*θ2>0时,表示圆心位于边的同侧,则第二个R角圆心即为初始圆心。当θ1*θ2<0时,则表示圆心位于边的不同侧,需再构建一个幅值为t+2r的转置矢量V12,若θ2<0,则方向为V1+90°,反之,方向为V1-90°。将初始圆心沿V12平移,即可得到第二个R角圆心C2,同时构建第二个圆角轮廓点坐标参数集[C2,r,θstart,θend],θstart=180°+θ1+90°(θ2>0)或者θstart=180°+θ1-90°(θ2<0),θend=θstart–(180°-θ2)(θ2>0)或者θend=θstart+(180°+θ2)(θ2<0),最后同样构建一个幅值为L3-2*t,方向为[(180°+θ1)-180°+θ2]的方向矢量V2,用于求解下一个圆心坐标。
同理依次计算剩余所有圆心坐标并构建对应的圆角轮廓点坐标参数集[CN,r,θstart,θend]。
步骤三:计算所有轮廓点坐标。
如图3所示,首先求解轮廓弧的离散点坐标,将每个圆角轮廓弧线进行等分离散,默认为5个点,根据每个圆角的轮廓点坐标集的输入参数,计算每个点的坐标p1~p5,得到弧线坐标点集。对于外轮廓弧线,参数集为[CN,t+r,θstart,θend],对于内轮廓弧线,CNin参数集为[CN,r,θend,θstart]。
然后,建立轮廓点顺序,从原点出发,对截面的所有轮廓点进行排序,按[S1,C1in/C1out,C2in/C2out,……,CNin/CNout,E1,E2,CNout/CNin,……,C2out/C2in,C1out/C1in,S2,S1]的顺序进行排列,其中,S1、S2、E1、E2为截面的四个端点,C N为CN对应的内外轮廓线坐标点集,组成一个个数为(10*N+5)的坐标集,形成最终的(10*N+5)×2的轮廓点坐标总矩阵。必要时,通过参数旋转角度θt,对上述的轮廓点坐标总矩阵进行转置,转置矩阵为:[cosθt sinθt;-sinθtcosθt]。
步骤四:计算截面属性参数。根据上一步的所有轮廓点坐标结果,采用数学积分的方法,可以计算出截面的面积、形心、惯性矩等截面属性,或者直接应用商用软件的内置计算模块,本算例中采用了MSC.Patran中软件的内置函数进行求解和轮廓显示。
以上所述,仅为本发明的具体实施例,对本发明进行详细描述,未详尽部分为常规技术。但本发明的保护范围不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。本发明的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种任意形状的钣金型材截面参数化方法,其特征在于:所述方法包括以下步骤:
步骤一:确定和输入截面参数;包括:截面的边数N、钣金的厚度t、钣金的R角半径r、边的长度L1~LN、相邻边之间的角度θ1~θN,|θN|<180°;
步骤二,计算R角圆心及角度;计算过程如下:
步骤a),选择截面的任意一个起始边的自由端点为原点,起始边为X横轴,建立局部坐标系,起始点坐标为(0,0),
判断起始边与第一个相邻边之间的角度值θ1,如果θ1>0,则以起始边和第一个相邻边的交点为圆心,从起始边位置逆时针旋转θ1得到第一个相邻边的位置;反之,则顺时针旋转θ1得到第一个相邻边的位置;
第一个R角的圆心坐标为C1(L1,-r)(θ1>0)或C1(L1,t+r)(θ1<0);
并构建第一个圆角轮廓点坐标参数集[C1,r,θstart,θend],θstart=90°(θ1>0)或者θstart=270°(θ1<0),θend=θstart–(180°-θ1)(θ1>0)或者θend=θstart+(180°+θ1)(θ1<0),其中θstart和θend分别为R角轮廓弧线首尾端点与R角圆心连线相对于X横轴的偏转角度;
建立一个幅值为L2-2*t,方向为(180°+θ1)的方向矢量V1,用于求解下一个圆心坐标;
步骤b),将第一个R角圆心坐标按矢量V1平移,得到第二个R角的初始圆心坐标;
判断第二个相邻角θ2,如果θ1*θ2>0时,则以第二个R角的初始圆心坐标作为最终的第二个R角圆心坐标;
如果θ1*θ2<0时,则构建一个幅值为t+2r的转置矢量V12,若θ2<0,则转置矢量V12方向为V1+90°,反之,方向为V1-90°;将初始圆心坐标按V12平移,得到第二个R角圆心坐标C2;
构建第二个圆角轮廓点坐标参数集[C2,r,θstart,θend],θstart=180°+θ1+90°(θ2>0)或者θstart=180°+θ1-90°(θ2<0),θend=θstart–(180°-θ2)(θ2>0)或者θend=θstart+(180°+θ2)(θ2<0);
构建一个幅值为L3-2*t,方向为[(180°+θ1)-180°+θ2]的方向矢量V2,用于求解下一个圆心坐标;
步骤c),重复步骤b)依次计算剩余所有圆心坐标并构建对应的圆角轮廓点坐标参数集[CN,r,θstart,θend];
步骤三,计算所有轮廓点坐标;
步骤四,计算截面属性参数。
2.根据权利要求1所述的一种任意形状的钣金型材截面参数化方法,其特征在于:所述需要确定和输入的截面参数还包括旋转角度θt。
3.根据权利要求2所述的一种任意形状的钣金型材截面参数化方法,其特征在于:所述步骤三中,包括:
在每段圆角轮廓弧线上等分选取5个弧线坐标点,根据每个圆角的轮廓点坐标参数集,计算每个弧线坐标点的坐标p1~p5,得到弧线坐标点集。
4.根据权利要求3所述的一种任意形状的钣金型材截面参数化方法,其特征在于:所述步骤三中,还包括
同一个R角圆心,对应内外两个轮廓弧线;
对于外轮廓弧线,CNout参数集为[CN,t+r,θstart,θend];对于内轮廓弧线,CNin参数集为[CN,r,θend,θstart];
从原点出发建立轮廓点顺序,对截面的所有轮廓点按[S1,C1in/C1out,C2in/C2out,……,CNin/CNout,E1,E2,CNout/CNin,……,C2out/C2in,C1out/C1in,S2,S1]的顺序进行排列,其中,S1、S2、E1、E2为截面的四个端点,
组成一个个数为(10*N+5)的坐标集,形成最终的(10*N+5)×2的轮廓点坐标总矩阵。
5.根据权利要求4所述的一种任意形状的钣金型材截面参数化方法,其特征在于:
若截面参数包括旋转角度θt,则轮廓点坐标总矩阵进行转置,转置矩阵为:[cosθt sinθt;-sinθt cosθt]。
6.根据权利要求5所述的一种任意形状的钣金型材截面参数化方法,其特征在于:根据轮廓点坐标总矩阵,采用数学积分的方法,计算截面属性,所述截面属性包括:截面的面积、形心、惯性矩。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011028772.8A CN112052541B (zh) | 2020-09-25 | 2020-09-25 | 一种任意形状的钣金型材截面参数化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011028772.8A CN112052541B (zh) | 2020-09-25 | 2020-09-25 | 一种任意形状的钣金型材截面参数化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112052541A CN112052541A (zh) | 2020-12-08 |
CN112052541B true CN112052541B (zh) | 2022-12-06 |
Family
ID=73605038
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011028772.8A Active CN112052541B (zh) | 2020-09-25 | 2020-09-25 | 一种任意形状的钣金型材截面参数化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112052541B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108592819A (zh) * | 2018-05-09 | 2018-09-28 | 燕山大学 | 一种平面弯曲钣金件截面弯曲轮廓检测装置及方法 |
CN110905852A (zh) * | 2019-11-26 | 2020-03-24 | 北京石油化工学院 | 一种动调轴流风机动叶三维造型方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08174335A (ja) * | 1994-12-26 | 1996-07-09 | Honda Motor Co Ltd | 歯車状ワークの面取り加工方法 |
JP5163331B2 (ja) * | 2008-07-14 | 2013-03-13 | 横浜ゴム株式会社 | タイヤ製造工程の管理方法 |
CN106055794A (zh) * | 2016-05-31 | 2016-10-26 | 浙江科澜信息技术有限公司 | 一种参数化建立三维管线模型的方法 |
CN108332681B (zh) * | 2018-01-03 | 2019-07-16 | 东北大学 | 一种薄壁管材大塑性弯曲截面轮廓曲线的确定方法 |
CN110990951A (zh) * | 2019-12-04 | 2020-04-10 | 中国直升机设计研究所 | 一种直升机外形设计方法 |
CN111274712A (zh) * | 2020-02-22 | 2020-06-12 | 哈尔滨理工大学 | 基于砂轮磨削轨迹的整体式立铣刀容屑槽轴截面建模方法 |
-
2020
- 2020-09-25 CN CN202011028772.8A patent/CN112052541B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108592819A (zh) * | 2018-05-09 | 2018-09-28 | 燕山大学 | 一种平面弯曲钣金件截面弯曲轮廓检测装置及方法 |
CN110905852A (zh) * | 2019-11-26 | 2020-03-24 | 北京石油化工学院 | 一种动调轴流风机动叶三维造型方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112052541A (zh) | 2020-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Leopoldseder et al. | Approximation of developable surfaces with cone spline surfaces | |
CN108074277A (zh) | 一种面向数字化制造的弯管测量和余量定位方法 | |
US6338000B1 (en) | Method of generating shape data method of verifying shape data | |
JP6056016B2 (ja) | 三次元モデル生成方法、システム及びプログラム | |
Pérez et al. | Quasi-developable B-spline surfaces in ship hull design | |
JP2006231345A (ja) | 曲げ加工方法、金属板、加熱位置決定プログラム、及び3次元形状処理装置 | |
CN112052541B (zh) | 一种任意形状的钣金型材截面参数化方法 | |
Marcon et al. | High-order curvilinear hybrid mesh generation for CFD simulations | |
CN113643423B (zh) | 一种具有轴对称特性的Waterbomb衍生折纸结构造型方法及折纸结构 | |
JP2007156604A (ja) | 曲面作成方法及び曲面作成プログラム並びに3次元形状処理装置 | |
JPH10269260A (ja) | 形状データ検証方法 | |
Docampo-Sánchez et al. | A regularization approach for automatic quad mesh generation | |
JP3349917B2 (ja) | 形状データ検証方法 | |
Foufou et al. | Conversion of Dupin cyclide patches into rational biquadratic Bézier form | |
JP5173971B2 (ja) | 曲げ加工方法、金属板、加熱位置決定プログラム及び装置 | |
JP3550949B2 (ja) | プレス成形解析用モデル作成方法 | |
CN117349970B (zh) | 一种船舶装配布置自动合理性验证方法 | |
Delgado-Gutiérrez et al. | An efficient and automated method to generate complex blade geometries for numerical analysis | |
CN116778089A (zh) | 含开口的圆柱形蒙皮桁条结构通用参数化几何建模方法 | |
CN115048825B (zh) | 一种薄壳体曲面仿真方法、装置、设备及介质 | |
JP4269678B2 (ja) | 計測データ比較装置および形状品質評価システム | |
JP3349915B2 (ja) | 形状データ作成方法 | |
CN107203657B (zh) | 一种基于船舶数字总图的舱室传热面积计算方法 | |
Stavric et al. | Ornamental plate shell structures | |
JPH11144093A (ja) | 解析メッシュ生成方法及び装置 |
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 |