CN109684694B - 简支条件下壁板承受均布压力中心等效应力评估方法 - Google Patents

简支条件下壁板承受均布压力中心等效应力评估方法 Download PDF

Info

Publication number
CN109684694B
CN109684694B CN201811523876.9A CN201811523876A CN109684694B CN 109684694 B CN109684694 B CN 109684694B CN 201811523876 A CN201811523876 A CN 201811523876A CN 109684694 B CN109684694 B CN 109684694B
Authority
CN
China
Prior art keywords
stress
wall plate
wallboard
simply supported
uniformly distributed
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
CN201811523876.9A
Other languages
English (en)
Other versions
CN109684694A (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.)
AVIC First Aircraft Institute
Original Assignee
AVIC First Aircraft Institute
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 AVIC First Aircraft Institute filed Critical AVIC First Aircraft Institute
Priority to CN201811523876.9A priority Critical patent/CN109684694B/zh
Publication of CN109684694A publication Critical patent/CN109684694A/zh
Application granted granted Critical
Publication of CN109684694B publication Critical patent/CN109684694B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明属于航空结构设计技术领域,涉及一种简支条件下壁板承受均布压力中心等效应力评估方法。本发明的方法通过数学推导,得到简支条件下壁板中心点等效应力的解析解,避免主观因素对计算结果的影响。壁板的等效应力为弯曲应力和薄膜应力的联合作用。弯曲应力于薄膜应力存在关于厚度和壁板中心挠度的一次函数关系,薄膜应力通过壁板材料常数、支持系数以及壁板的中心挠度即可求得,继而得到弯曲应力,最后得到壁板中心处的等效应力。本发明的方法能够精确得到壁板中心点复合应力,避免了目前通过查手册曲线的方法中主观因素引入精度无法保证的技术问题。

Description

简支条件下壁板承受均布压力中心等效应力评估方法
技术领域
本发明属于航空结构设计技术领域,具体而言,涉及一种简支条件下壁板承受均布压力中心等效应力评估方法。
背景技术
在航空结构设计中,壁板承受均布压力是很常见的一种受力形式,壁板的失效往往是壁板中心点复合应力过高造成。目前计算壁板中心点复合应力往往通过查手册曲线来获得,在查曲线过程中对参数的选取主观因素很大,往往不同的人得出的结果不唯一,计算的精度无法保证。因此,亟需解决的技术问题是如何通过数学推导,得到在简支边界下壁板中心点等效应力的解析解,避免主观因素对计算结果的影响。
发明内容
本发明的目的是:设计一种简支条件下壁板承受均布压力中心等效应力评估方法,以解决目前查曲线过程中对参数的选取主观因素很大,计算的精度无法保证的技术问题。
为解决此技术问题,本发明的技术方案是:
一种简支条件下壁板承受均布压力中心等效应力评估方法,所述的简支条件下壁板承受均布压力中心等效应力评估方法主要包括以下步骤:
1、计算壁板的弯曲刚度D、为弯曲应力产生的壁板中心挠度f0
2、计算薄膜应力σm、弯曲应力σbend
Figure BDA0001902689100000011
式中,π代表圆周率,为常量;E代表壁板弹性模型;μ代表壁板材料泊松比;k代表支持系数;b代表壁板宽度;f代表壁板中心挠度;
Figure BDA0001902689100000021
式中,δ代表壁板厚度;
3、计算壁板的复合应力σall
Figure BDA0001902689100000022
所述的步骤1中壁板中心挠度f0根据如下公式计算:
Figure BDA0001902689100000023
式中,P代表均布压力。
所述的步骤2中所述的壁板中心挠度f根据如下公式计算为:
Figure BDA0001902689100000024
式中,a代表无量纲参数。
所述的步骤2中所述的壁板中心挠度f根据如下公式计算为:
Figure BDA0001902689100000025
式中,σm代表薄膜应力。
本发明的有益效果是:本发明的简支条件下壁板承受均布压力中心等效应力评估方法能够精确得到壁板中心点复合应力,避免了目前通过查手册曲线的方法中主观因素引入精度无法保证的技术问题。
附图说明
图1为本发明的方法流程示意图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明:
整体思路:
通过数学推导,得到简支条件下壁板中心点等效应力的解析解,避免主观因素对计算结果的影响。壁板的等效应力为弯曲应力和薄膜应力的联合作用。弯曲应力于薄膜应力存在关于厚度和壁板中心挠度的一次函数关系,薄膜应力通过壁板材料常数、支持系数以及壁板的中心挠度即可求得,继而得到弯曲应力,最后得到壁板中心处的等效应力。
附图1所示为本发明方法的流程示意图,具体计算方法如下:
壁板的强度取决于弯曲应力σbend和薄膜应力σm的联合作用,对承受均布载荷下的壁板中心复合应力σall为:
σall=σmbend
记相邻纵梁的间距为b,壁板厚度为δ,壁板的弯曲刚度为
Figure BDA0001902689100000031
式中:E为壁板材料的弹性模量,μ为壁板材料的泊松比。
壁板中心挠度为
Figure BDA0001902689100000032
式中:f0为弯曲应力产生的中点挠度
Figure BDA0001902689100000033
a为无量纲参数
Figure BDA0001902689100000034
薄膜应力σm按下式计算:
Figure BDA0001902689100000035
式中k为支持系数。
Figure BDA0001902689100000036
其中F为壁板围框元件承压剖面面积。
弯曲应力σband按下式计算:
Figure BDA0001902689100000037
壁板的复合应力为:
Figure BDA0001902689100000041
将(1)代入(3)得:
Figure BDA0001902689100000042
将(2)代入上式得:
Figure BDA0001902689100000043
化简得
Figure BDA0001902689100000044
记为
Figure BDA0001902689100000045
Figure BDA0001902689100000046
对(6)进行求解:
Figure BDA0001902689100000047
其中:
M=B2-3AC
N-BC-9AD
Q=C2-3BD
Figure BDA0001902689100000048
Figure BDA0001902689100000051
具体算例:一宽为200mm,长为500mm,厚度为2mm的壁板,承受50KPa的均布压力,壁板周圈承压纵梁面积为100mm2,壁板材料为7050铝合金,弹性模量为71000MPa,泊松比为0.3,承压纵梁与壁板采用单排螺栓连接,求该壁板中心点最大等效应力。
该壁板与周圈纵梁连接为单排螺栓,壁板支持方式为简支,采用本专利方法进行计算。计算如下:
对应参数为
b=200mm,δ=2mm,P=50KPa,F=100mm2,E=71000MPa,μ=0.3
将以上参数代入公式:
A=12100,B=156200,C=504100,D=1.96e8
M=6.09e9,N=-2.12e13,Q=-9.15e13,
a=3.22,f=4.77
σm=20.80MPa,σbend=17.45MPa
求得最大等效应力为σall=38.25MPa。

Claims (1)

1.一种简支条件下壁板承受均布压力中心等效应力评估方法,其特征在于:所述的简支条件下壁板承受均布压力中心等效应力评估方法包括以下步骤:
1.1、计算壁板的弯曲刚度D、为弯曲应力产生的壁板中心挠度f0
1.2、计算薄膜应力σm、弯曲应力σbend
Figure FDA0003882774830000011
式中,π代表圆周率,为常量;E代表壁板弹性模型;μ代表壁板材料泊松比;k代表支持系数;b代表壁板宽度;f代表壁板中心挠度;
Figure FDA0003882774830000012
式中,δ代表壁板厚度;
1.3、计算壁板的复合应力σall
Figure FDA0003882774830000013
所述的步骤1.1中为弯曲应力产生的壁板中心挠度f0根据如下公式计算:
Figure FDA0003882774830000014
式中,P代表均布压力;
所述的步骤1.2中所述的壁板中心挠度f根据如下公式计算为:
Figure FDA0003882774830000015
式中,a代表无量纲参数;
所述的无量纲参数a根据如下公式计算为:
Figure FDA0003882774830000016
式中,σm代表薄膜应力。
CN201811523876.9A 2018-12-12 2018-12-12 简支条件下壁板承受均布压力中心等效应力评估方法 Active CN109684694B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811523876.9A CN109684694B (zh) 2018-12-12 2018-12-12 简支条件下壁板承受均布压力中心等效应力评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811523876.9A CN109684694B (zh) 2018-12-12 2018-12-12 简支条件下壁板承受均布压力中心等效应力评估方法

Publications (2)

Publication Number Publication Date
CN109684694A CN109684694A (zh) 2019-04-26
CN109684694B true CN109684694B (zh) 2023-01-13

Family

ID=66186582

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811523876.9A Active CN109684694B (zh) 2018-12-12 2018-12-12 简支条件下壁板承受均布压力中心等效应力评估方法

Country Status (1)

Country Link
CN (1) CN109684694B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007102599A (ja) * 2005-10-06 2007-04-19 Sumitomo Kinzoku Kozan Siporex Kk 許容下地部材判定システム及び許容下地部材判定プログラム
CN103955593A (zh) * 2013-11-29 2014-07-30 大连船舶重工集团有限公司 一种自升式钻井平台悬臂梁载荷图谱获取方法
CN106650168A (zh) * 2017-01-03 2017-05-10 山东理工大学 基于限位挠度的非等偏频一级渐变板簧应力强度的校核法
KR20170064760A (ko) * 2015-12-02 2017-06-12 세메스 주식회사 프로브 검사 장치
CN108254949A (zh) * 2018-01-26 2018-07-06 惠州市华星光电技术有限公司 一种曲面状态下液晶面板的寿命及最大曲率的估算方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3763697A (en) * 1970-04-15 1973-10-09 Sturm Stress Inc Method and apparatus for determining stress
DE3415754A1 (de) * 1984-03-20 1985-09-26 Josef Gartner & Co, 8883 Gundelfingen Verfahren zur erhoehung der tragfaehigkeit duennwandiger verkleidungselemente
JP4198821B2 (ja) * 1999-04-02 2008-12-17 日本フルハーフ株式会社 ウイング車両の荷箱構造
JP5602455B2 (ja) * 2009-12-24 2014-10-08 大成建設株式会社 梁部材および建物構造体
CN104280155B (zh) * 2013-07-14 2019-03-12 付康 一种确定多层膜薄膜应力的系统与方法
CN104034449B (zh) * 2014-06-03 2016-04-13 东南大学 薄膜材料残余应力测试结构及方法
CN105787149B (zh) * 2015-12-27 2019-03-08 北京航空航天大学 一种由弧齿锥齿轮传动系统轴上功率谱向齿面应力谱精确转换的方法
CN105928810B (zh) * 2016-06-21 2018-11-30 巨石集团有限公司 纤维增强材料弯曲疲劳测试方法
CN106407557B (zh) * 2016-09-18 2019-07-09 常州大学 无模液压准静态薄板鼓胀变形极限应变及载荷的图算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007102599A (ja) * 2005-10-06 2007-04-19 Sumitomo Kinzoku Kozan Siporex Kk 許容下地部材判定システム及び許容下地部材判定プログラム
CN103955593A (zh) * 2013-11-29 2014-07-30 大连船舶重工集团有限公司 一种自升式钻井平台悬臂梁载荷图谱获取方法
KR20170064760A (ko) * 2015-12-02 2017-06-12 세메스 주식회사 프로브 검사 장치
CN106650168A (zh) * 2017-01-03 2017-05-10 山东理工大学 基于限位挠度的非等偏频一级渐变板簧应力强度的校核法
CN108254949A (zh) * 2018-01-26 2018-07-06 惠州市华星光电技术有限公司 一种曲面状态下液晶面板的寿命及最大曲率的估算方法

Also Published As

Publication number Publication date
CN109684694A (zh) 2019-04-26

Similar Documents

Publication Publication Date Title
CN1289220C (zh) 用于对带材进行平直度测量的测量辊
Adeoti et al. Stability of 6082-T6 aluminium alloy columns with H-section and rectangular hollow sections
CN105675484B (zh) 均布载荷下中心带刚性板的环形薄膜弹性能的确定方法
CN109726451B (zh) 固支条件下壁板承受均布压力中心等效应力评估方法
CN109323924A (zh) 最大挠度受弹性限制下圆形薄膜最大应力的确定方法
CN109684694B (zh) 简支条件下壁板承受均布压力中心等效应力评估方法
CN113092041B (zh) 横向均布载荷下环形薄膜最大挠度的确定方法
CN111442985A (zh) 一种确定横向均布载荷下圆形薄膜最大挠度的方法
CN105956379B (zh) 均布载荷下中心带刚性板的环形薄膜最大应力的确定方法
CN105138806B (zh) 油气弹簧不等厚环形阀片的强度校核方法
CN108268675B (zh) 疲劳损伤确定方法以及疲劳损伤确定装置
CN110303070B (zh) 一种双金属复合板矫直工艺新方法
CN104237037A (zh) 基于Berkovich压痕的材料弹塑性参数仪器化压入测试方法
CN111442980A (zh) 均布载荷下圆形薄膜最大挠度的确定方法
JP5655394B2 (ja) 絞りビード試験方法及びその試験方法で求めた物性値を用いたプレス成形解析方法
CN114674683B (zh) 获取材料单轴应力应变关系的锥压入式小冲杆试验方法
CN105956311B (zh) 非端部接触式少片端部加强型副簧根部厚度的设计方法
CN111842547B (zh) 一种具有拉压不对称性金属板材矫直工艺新方法
DE102015113147A1 (de) Einströmdüse mit einer Beschichtung in einem Messringraumbereich
CN113551977A (zh) 带有刚性内边缘的环形薄膜的挠度的确定方法
CN107983772A (zh) 一种板材冷轧过程中可靠性有限元模型的建立方法
CN113051780A (zh) 一种平板结构轴压屈曲载荷判断方法
CN114386234B (zh) 一种厚板jco成形工艺压下量的计算新方法
CN105890883B (zh) 非端部接触式少片斜线型主副簧各片强度的校核方法
CN104252964A (zh) 一种变压器辐射型铁心柱制作方法

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