CN101691012A - 箱形支撑结构内部加筋板分布优化设计方法 - Google Patents
箱形支撑结构内部加筋板分布优化设计方法 Download PDFInfo
- Publication number
- CN101691012A CN101691012A CN200910197135A CN200910197135A CN101691012A CN 101691012 A CN101691012 A CN 101691012A CN 200910197135 A CN200910197135 A CN 200910197135A CN 200910197135 A CN200910197135 A CN 200910197135A CN 101691012 A CN101691012 A CN 101691012A
- Authority
- CN
- China
- Prior art keywords
- model
- beam element
- design
- box
- loading end
- 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
Landscapes
- Rod-Shaped Construction Members (AREA)
Abstract
本发明涉及一种箱形支撑结构内部加筋板分布优化设计方法,包括以下具体步骤:(1)简化为中空结构的初始模型,(2)中空模型有限元网格划分,(3)形成纤维模型,(4)施加边界条件,(5)承载面材料分布优化设计,(6)确定应保留的梁单元,(7)形成加筋板。本发明提出一种称为纤维模型的简化模型,将复杂的三维结构拓扑优化设计问题转化为简单的二维结构拓扑优化设计问题,利用密度法,自动快速设计箱形支撑结构内部的加筋板分布,使结构达到最好的技术经济综合性能。
Description
技术领域
本发明涉及一种机床的床身及立柱、起重机的大梁与支腿等箱形支撑结构的内部加筋板分布优化设计方法。
背景技术
针对箱形支撑结构内部加筋板的优化设计问题,目前大多数研究都局限于加筋板间隔与截面的尺寸优化设计,而由于加筋板的形式及布局优化设计自由度大,实现困难,现有的设计方法一般根据材料力学理论,结合设计者的经验提出若干种加筋板的分布形式后,在对这些加筋结构进行力学分析比较的基础上,进行选型和改进设计。显然这样的设计方法无法得到真正最优的结构,其最终设计结果很大程度上决定于设计者的经验,而且其设计过程无法实现自动化。
随着近年来结构拓扑优化设计理论的发展,加强筋的布局优化问题,尤其是薄壁板壳结构加强筋的分布优化设计得到了深入的研究。对于箱形支撑结构来说,将箱形结构内部填充材料,然后对其内部的实体采用结构拓扑优化设计方法进行设计,理论上可实现内部加筋板分布优化。但对于工程中常见的大体积箱形结构,对其内部填充的实体进行拓扑优化设计计算量相当大,且其设计结果形态复杂,实际应用中还需要进行进一步的后处理,因此虽然从理论上说该方法行之有效,但实际应用上仍存在问题。
发明内容
本发明是要提供一种箱形支撑结构内部加筋板分布优化设计方法,采用一种称为纤维模型的简化模型,将复杂的三维结构拓扑优化设计问题转化为简单的二维结构设计问题,并采用密度法对结构进行优化设计,自动快速设计箱形支撑结构内部的加筋板分布,使结构达到最好的技术经济综合性能。
为实现上述目的,本发明的技术方案是:一种箱形支撑结构内部加筋板分布优化设计方法,包括以下具体步骤:
1)简化为中空结构的初始设计模型
将箱形支撑结构简化为中空的箱体,形成初始设计模型;
2)有限元网格划分
对初始设计模型进行有限元网格划分,所有外壁均划分为壳单元,且箱体的约束面和设计面的网格需一一对应;
3)形成纤维模型
在约束面和承载面对应的节点上建立梁单元,形成纤维模型;
4)对设计模型施加边界条件
对结构施加载荷和约束,进行静力分析;
5)承载面材料分布优化设计
将承载设置为设计区域,采用式(1)的优化模型进行承载面材料分布的优化设计,
min vol(ρ)=vol(ρ1,ρ2,...,ρn)
subject to:
δmax≤[δ] (1)
σmax≤[σ]
0≤ρi≤1 i=1,2,...n
式中,vol为结构的总体积,ρi为单元i密度,n为单元个数,δmax和σmax分别为结构的最大位移和最大应力,[δ]和[σ]分别为结构的许用位移和许用应力;
6)确定应保留的梁单元
根据第5步骤的设计结果得到承载面的材料最优布置,确定应保留的梁单元,即当某一梁单元的节点在承载面的高密度区时,保留该梁单元;反之,若某一梁单元的节点在承载面的低密度区时,去除该梁单元;
7)形成加筋板
连接各保留的梁单元形成加筋板。
本发明提出一种称为纤维模型的简化模型,将复杂的三维结构拓扑优化设计问题转化为简单的二维结构优化设计问题,利用密度法,自动快速设计箱形支撑结构内部的加筋板分布,使结构达到最好的技术经济综合性能。
附图说明
图1是对角筋板的纤维模型图;
图2是具有内部平行梁单元的箱形结构初始设计模型图;
图3是本发明的流程图;
图4是受扭床身加筋板设计力学模型及设计结果图;
其中,(a)床身承受扭转载荷图,(b)床身力学模型图,(c)有限元设计模型图,(d)床身加筋板最优布置形态图;
图5是床身模型图;其中,(a)优化模型图,(b)横隔板模型图。
具体实施方式
下面结合附图与实施例对本发明作进一步的说明。
一.设计模型的简化
箱形结构内部的加筋板可看成由平行的纤维组成,这些纤维一端连接着承载的外壁,另一端连接着结构的支撑面。为了更清楚地说明,以图1所示的对角筋板箱体为例说明纤维模型的表示方法,可将对角筋板看成由图1所示的平行纤维组成。将平行纤维表示为有限元模型中的梁单元2,则可得图2所示的具有内部平行梁单元的箱形结构初始设计模型。图2中1表示承载面,3表示支撑面。
二.设计流程
设计流程如图3所示。
(1)将设计对象简化为中空的箱体,对其进行有限元网格划分,所有外壁均划分为壳单元,注意箱体的约束面和设计面的网格需一一对应;
(2)在约束面和承载面1对应的节点上建立梁单元2,形成纤维模型;
(3)对结构施加载荷和约束,进行静力分析;
(4)将承载设置为设计区域,采用式(1)的优化模型进行承载面材料分布的优化设计。
min vol(ρ)=vol(ρ1,ρ2,...,ρn)
subject to:
δmax≤[δ] (1)
σmax≤[σ]
0≤ρi≤1 i=1,2,...n
式中,vol为结构的总体积,ρi为单元i密度,n为单元个数,δmax和σmax分别为结构的最大位移和最大应力,[δ]和[σ]分别为结构的许用位移和许用应力。式(1)表示以单元密度为设计变量,在满足结构的强度和刚度要求的条件,寻求使结构体积达到最小的单元密度分布。单元密度在0和1之间连续变化,当单元密度为1时,表示该单元应为实体,当单元密度为0时,表示该单元应为空隙。当单元密度在0和1之间时,需采用下式对中间密度进行惩罚,以使设计结果为仅含有0和1的密度,即单元或为实体或为空隙。
K(ρ)=ρpK (2)
式中,K为惩罚后的单元刚度矩阵,ρ为单元的材料密度,K为实际的单元刚度矩阵,而p称为惩罚因子,通常p的取值在2-3之间。
(5)根据第四步骤的设计结果得到承载面的材料最优布置,确定应保留的梁单元,即当某一梁单元的节点在承载面的高密度区时,保留该梁单元;反之,若某一梁单元的节点在承载面的低密度区时,去除该梁单元。
(6)连接各保留的梁单元形成加筋板。
三.应用例
本发明提出的设计方法可对任意形状箱形支撑结构内部加筋板分布进行优化设计,下面以机床床身结构为例作为说明。
机床床身承受的载荷比较复杂,包括其上支撑的工作台、工件等部件重量,切削力及其本身自重等,这些载荷不仅使床身受压,同时使床身受扭及斜弯曲。试验表明,在各种变形中,扭转变形可占总变形的60%以上,因此本例主要研究如何在床身结构内部布置合理的加筋板,以提高床身的抗扭刚度。
受扭床身的力学模型如图4所示,图4(a)为当床身头尾架上承受绕x轴的扭矩静力当量载荷时的情况,图4(b)为相应的力学模型。根据本发明提出的设计方法,建立床身的有限元模型,并在箱形梁内部建立Z方向的梁单元,如图4(c)所示。将床身底面固定,以上表面为设计面,根据式(1)的优化模型,以体积最小为设计目标,考虑结构的变形和应力约束,得到上表面的材料分布如图4(d)所示。
根据图4(d)所得的设计结果,建立图5(a)所示的带有加筋板的床身结构,该结构为左右两侧对称布置2对斜向加筋板,并与中间布置的一纵向加筋板连接。为了说明设计结果的合理性,将图5(a)的优化设计模型与图5(b)所示的目前常见的带横隔板的床身结构进行比较,两种结构均具有上盖板,并承受图4(b)所示的扭转载荷。为了说明结构的技术经济性能,定义技术经济性能指标μ,
μ=EV (3)
式中,E为结构的弹性应变能,其值越大,说明结构的柔度越大,刚度越小;V为结构的总体积。考虑刚度最优的理想结构是结构的体积最小而应变能也最小,即技术经济性能参数μ的值越小越好。
图5所示2种床身结构的应变能和技术经济性能指标列于表1,表中将带有横隔板的床身结构的各项指标表示为1。由表可见,与常见的横隔板的床身结构相比,优化模型在体积下降1.7%的情况下,结构的总应变能下降25.5%,技术经济性能参数下降26.7%。可见采用本发明提出的设计方法可有效布置床身内部的加筋板,使床身结构的技术经济性能大大提高。
表1各种加筋板的床身性能对比
Claims (1)
1.一种箱形支撑结构内部加筋板分布优化设计方法,其特征在于:包括以下具体步骤:
1)简化为中空结构的初始设计模型
将箱形支撑结构简化为中空的箱体,形成初始设计模型;
2)有限元网格划分
对初始设计模型进行有限元网格划分,所有外壁均划分为壳单元,且箱体的约束面和设计面的网格需一一对应;
3)形成纤维模型
在约束面和承载面对应的节点上建立梁单元,形成纤维模型;
4)对设计模型施加边界条件
对结构施加载荷和约束,进行静力分析;
5)承载面材料分布优化设计
将承载设置为设计区域,采用式(1)的优化模型进行承载面材料分布的优化设计,
min vol(ρ)=vol(ρ1,ρ2,...,ρn)
subject to:
δmax≤[δ](1)
σmax≤[σ]
0≤ρi≤1i=1,2,...n
式中,vol为结构的总体积,ρi为单元i密度,n为单元个数,δmax和σmax分别为结构的最大位移和最大应力,[δ]和[σ]分别为结构的许用位移和许用应力;
6)确定应保留的梁单元
根据第5步骤的设计结果得到承载面的材料最优布置,确定应保留的梁单元,即当某一梁单元的节点在承载面的高密度区时,保留该梁单元;反之,若某一梁单元的节点在承载面的低密度区时,去除该梁单元;
7)形成加筋板
连接各保留的梁单元形成加筋板。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101971350A CN101691012B (zh) | 2009-10-14 | 2009-10-14 | 箱形支撑结构内部加筋板分布优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101971350A CN101691012B (zh) | 2009-10-14 | 2009-10-14 | 箱形支撑结构内部加筋板分布优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101691012A true CN101691012A (zh) | 2010-04-07 |
CN101691012B CN101691012B (zh) | 2011-05-04 |
Family
ID=42079755
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101971350A Expired - Fee Related CN101691012B (zh) | 2009-10-14 | 2009-10-14 | 箱形支撑结构内部加筋板分布优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101691012B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101811257A (zh) * | 2010-04-23 | 2010-08-25 | 上海理工大学 | 机床床身结构优化设计方法 |
CN102063540A (zh) * | 2010-12-30 | 2011-05-18 | 西安交通大学 | 一种机床床身结构优化设计方法 |
CN102779211A (zh) * | 2012-07-31 | 2012-11-14 | 上海理工大学 | 大型电除尘器烟箱结构优化设计方法 |
CN102880761A (zh) * | 2012-10-12 | 2013-01-16 | 上海理工大学 | 复杂装配体的构件拓扑优化设计方法 |
CN103235851A (zh) * | 2013-01-04 | 2013-08-07 | 西安交通大学 | 一种面向高比刚度的机床支撑件筋板布局智能化设计方法 |
CN105202157A (zh) * | 2015-10-21 | 2015-12-30 | 大连理工大学 | 一种功率分流型风电齿轮箱箱体 |
CN105243240A (zh) * | 2015-11-10 | 2016-01-13 | 上海理工大学 | 高刚轻质机床大件结构优化设计方法 |
CN105260567A (zh) * | 2015-11-06 | 2016-01-20 | 上海理工大学 | 三维箱型结构筋板分布仿生优化设计方法 |
CN105631084A (zh) * | 2015-06-04 | 2016-06-01 | 重庆大学 | 行星齿轮减速齿轮箱箱体轻量化结构方法 |
CN106845013A (zh) * | 2017-02-21 | 2017-06-13 | 北京航天发射技术研究所 | 一种齿轮传动箱体内、外筋板结构拓扑优化方法 |
CN108482263A (zh) * | 2018-01-29 | 2018-09-04 | 中国科学院长春光学精密机械与物理研究所 | 车载经纬仪支撑平台 |
CN108549772A (zh) * | 2018-04-14 | 2018-09-18 | 中国科学院云南天文台 | 一种结构优化的偏振片轮及其设计方法 |
CN109145405A (zh) * | 2018-08-01 | 2019-01-04 | 中船第九设计研究院工程有限公司 | 一种钢结构纤维模型损伤评价方法 |
CN109344524A (zh) * | 2018-10-18 | 2019-02-15 | 燕山大学 | 一种薄板结构加强筋分布优化方法 |
CN109408939A (zh) * | 2018-10-18 | 2019-03-01 | 燕山大学 | 一种兼顾应力和位移约束的薄板结构加强筋分布优化的改进方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6188150B1 (en) * | 1999-06-16 | 2001-02-13 | Euv, Llc | Light weight high-stiffness stage platen |
CN101504678B (zh) * | 2009-03-09 | 2011-12-14 | 西北工业大学 | 薄壁加筋结构加筋布局优化设计方法 |
CN101510231A (zh) * | 2009-03-27 | 2009-08-19 | 上海理工大学 | 基于植物根系形成机理的板壳结构加强筋分布设计方法 |
-
2009
- 2009-10-14 CN CN2009101971350A patent/CN101691012B/zh not_active Expired - Fee Related
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101811257A (zh) * | 2010-04-23 | 2010-08-25 | 上海理工大学 | 机床床身结构优化设计方法 |
CN102063540A (zh) * | 2010-12-30 | 2011-05-18 | 西安交通大学 | 一种机床床身结构优化设计方法 |
CN102063540B (zh) * | 2010-12-30 | 2013-03-13 | 西安交通大学 | 一种机床床身结构优化设计方法 |
CN102779211B (zh) * | 2012-07-31 | 2014-09-24 | 上海理工大学 | 大型电除尘器烟箱结构优化设计方法 |
CN102779211A (zh) * | 2012-07-31 | 2012-11-14 | 上海理工大学 | 大型电除尘器烟箱结构优化设计方法 |
CN102880761A (zh) * | 2012-10-12 | 2013-01-16 | 上海理工大学 | 复杂装配体的构件拓扑优化设计方法 |
CN103235851A (zh) * | 2013-01-04 | 2013-08-07 | 西安交通大学 | 一种面向高比刚度的机床支撑件筋板布局智能化设计方法 |
CN103235851B (zh) * | 2013-01-04 | 2016-02-24 | 西安交通大学 | 一种面向高比刚度的机床支撑件筋板布局智能化设计方法 |
CN105631084A (zh) * | 2015-06-04 | 2016-06-01 | 重庆大学 | 行星齿轮减速齿轮箱箱体轻量化结构方法 |
CN105631084B (zh) * | 2015-06-04 | 2019-03-22 | 重庆大学 | 行星齿轮减速齿轮箱箱体轻量化结构方法 |
CN105202157A (zh) * | 2015-10-21 | 2015-12-30 | 大连理工大学 | 一种功率分流型风电齿轮箱箱体 |
CN105260567A (zh) * | 2015-11-06 | 2016-01-20 | 上海理工大学 | 三维箱型结构筋板分布仿生优化设计方法 |
CN105260567B (zh) * | 2015-11-06 | 2018-12-14 | 上海理工大学 | 三维箱型结构筋板分布仿生优化设计方法 |
CN105243240A (zh) * | 2015-11-10 | 2016-01-13 | 上海理工大学 | 高刚轻质机床大件结构优化设计方法 |
CN105243240B (zh) * | 2015-11-10 | 2018-08-28 | 上海理工大学 | 高刚轻质机床大件结构优化设计方法 |
CN106845013A (zh) * | 2017-02-21 | 2017-06-13 | 北京航天发射技术研究所 | 一种齿轮传动箱体内、外筋板结构拓扑优化方法 |
CN106845013B (zh) * | 2017-02-21 | 2021-01-29 | 北京航天发射技术研究所 | 一种齿轮传动箱体内、外筋板结构拓扑优化方法 |
CN108482263A (zh) * | 2018-01-29 | 2018-09-04 | 中国科学院长春光学精密机械与物理研究所 | 车载经纬仪支撑平台 |
CN108482263B (zh) * | 2018-01-29 | 2020-05-19 | 中国科学院长春光学精密机械与物理研究所 | 车载经纬仪支撑平台 |
CN108549772A (zh) * | 2018-04-14 | 2018-09-18 | 中国科学院云南天文台 | 一种结构优化的偏振片轮及其设计方法 |
CN108549772B (zh) * | 2018-04-14 | 2021-11-09 | 中国科学院云南天文台 | 一种结构优化的偏振片轮及其设计方法 |
CN109145405A (zh) * | 2018-08-01 | 2019-01-04 | 中船第九设计研究院工程有限公司 | 一种钢结构纤维模型损伤评价方法 |
CN109145405B (zh) * | 2018-08-01 | 2023-04-14 | 中船第九设计研究院工程有限公司 | 一种钢结构纤维模型损伤评价方法 |
CN109344524A (zh) * | 2018-10-18 | 2019-02-15 | 燕山大学 | 一种薄板结构加强筋分布优化方法 |
CN109408939A (zh) * | 2018-10-18 | 2019-03-01 | 燕山大学 | 一种兼顾应力和位移约束的薄板结构加强筋分布优化的改进方法 |
CN109408939B (zh) * | 2018-10-18 | 2022-11-29 | 燕山大学 | 一种兼顾应力和位移约束的薄板结构加强筋分布优化的改进方法 |
CN109344524B (zh) * | 2018-10-18 | 2022-12-09 | 燕山大学 | 一种薄板结构加强筋分布优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101691012B (zh) | 2011-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101691012B (zh) | 箱形支撑结构内部加筋板分布优化设计方法 | |
CN101950319B (zh) | 一种高速立式加工中心大跨度横梁拓扑优化设计方法 | |
CN102063540B (zh) | 一种机床床身结构优化设计方法 | |
CN103336868B (zh) | 一种立式加工中心机床箱体类支承件结构的建立方法 | |
CN101811257A (zh) | 机床床身结构优化设计方法 | |
CN114756934B (zh) | 一种三维多尺度超材料结构优化设计方法 | |
Chen et al. | Vibration and buckling of truss core sandwich plates on an elastic foundation subjected to biaxial in-plane loads | |
CN103061800A (zh) | 一种矿用救生舱舱体结构优化设计方法 | |
CN111737908B (zh) | 一种基于动载荷静力等效的蒙皮桁条结构快速动态优化设计方法 | |
CN112131770B (zh) | 一种考虑可靠性的功能梯度连续体结构轻量化设计方法 | |
Liu et al. | Bionic design for column of gantry machining center to improve the static and dynamic performance | |
CN115906572A (zh) | 一种基于apdl的钢管混凝土轨枕结构拓扑优化设计方法 | |
Xu et al. | The design and optimization of large-scales heavy gantry NC machining center based on finite element method | |
CN112069715A (zh) | 一种基于多材料结构的拓扑优化方法 | |
CN205057521U (zh) | 一种复合型机床床身 | |
Xiong et al. | The static and modal analysis of 12000KN fine blanking press frame | |
CN217108973U (zh) | 一种gis支架 | |
Zhang et al. | Application and practical validation of topology optimization technology for the frame of biaxial tensile testing machine | |
CN116029043A (zh) | 一种变压器气垫车托盘及其设计方法 | |
Huang et al. | Finite element analysis and topology optimization for the gantry milling machine column structure | |
Qu | Structure Optimal Design Method and Engineering Application | |
Xuanming | STUDY ON OPTIMIZATION DESIGN FOR MECHANICAL SYSTEM OF CRANE BASED ON FINITE ELEMENT ANALYSIS. | |
Fan et al. | Research on structural improvement of worktable of B2-K1018 double grinding wheel frames follow-up CNC crankshaft grinder | |
Tong et al. | Research on the Application of Structural Topology Optimisation in the High-Precision Design of a Press Machine Frame | |
Kumar et al. | Three-dimensional finite element analysis on the flexural behavior of composite beams under linear displacement |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110504 Termination date: 20121014 |