CN112800538A - 一种任意薄壁梁结构的层合化分析方法 - Google Patents
一种任意薄壁梁结构的层合化分析方法 Download PDFInfo
- Publication number
- CN112800538A CN112800538A CN202110056368.XA CN202110056368A CN112800538A CN 112800538 A CN112800538 A CN 112800538A CN 202110056368 A CN202110056368 A CN 202110056368A CN 112800538 A CN112800538 A CN 112800538A
- Authority
- CN
- China
- Prior art keywords
- shear
- thin
- section
- formula
- layer
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 35
- 238000005452 bending Methods 0.000 claims abstract description 48
- 238000000034 method Methods 0.000 claims abstract description 18
- 238000003475 lamination Methods 0.000 claims abstract description 13
- 230000009471 action Effects 0.000 claims abstract description 6
- 238000010030 laminating Methods 0.000 claims abstract description 4
- 239000000463 material Substances 0.000 claims description 12
- 230000010354 integration Effects 0.000 claims description 10
- 230000007935 neutral effect Effects 0.000 claims description 8
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 238000005520 cutting process Methods 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 2
- 230000010365 information processing Effects 0.000 claims 1
- 229910000831 Steel Inorganic materials 0.000 abstract description 2
- 239000010959 steel Substances 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 16
- 238000013461 design Methods 0.000 description 7
- 230000007547 defect Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 4
- 238000010008 shearing Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013433 optimization analysis Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/15—Vehicle, aircraft or watercraft design
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种任意薄壁梁结构的层合化分析方法。步骤为:对待分析的薄壁梁的进行建模,并确定待分析的薄壁梁截面的弯曲主轴;在弯曲主轴方向上进行层合化处理,通过等效叠加即可得到总层合结构所对应的在第i层上的宽度;计算每个主轴方向上给定剪切载荷及弯矩时即单向弯剪模型下截面每一层上总的剪流和正应力;计算各个子结构元素在每一层上的剪流和正应力;计算刚心的位置;计算扭转刚度及在扭矩作用下在每个子结构元素上所产生的剪流;计算在刚心处施加两个独立的沿主轴方向的剪切载荷时每个子结构元素上所产生的剪流。方法易于进一步改进,如在基础层合板应用更精确的Timoshenko梁假设,就可容易进一步提高整个方法的精确度。
Description
技术领域
本发明属飞机机体结构设计技术领域,尤其涉及一种适用于具有任意开口的薄壁梁或具有任意多闭室的闭口薄壁梁进行结构分析的任意薄壁梁结构的层合化分析方法。
背景技术
目前,飞机的机翼,机身和尾翼(包括垂直尾翼和水平尾翼)都是各种形状的薄壁结构,而且在飞机的初始剛度和强度的设计阶段很多情况下均可简化为梁结构。因此对于各种薄壁梁结构进行快速有效的分析计算在飞机结构的初始设计阶段起着重要作用。
尽管目前各种通用的商业有限元软件可为薄壁梁结构的分析计算提供有效的工具,但是在初始设计阶段,要考虑各种各样的可能构型并进行各种优化仍然是一个很繁杂的过程。目前参数化有限元的发展,在一定程度上可以代替人工对特定的不同几何参数情况进行自动化分析,但由于其仍然建立在传统有限元方法之上,对于初始的结构设计,在计算时间,结果分析与处理等方面仍然有其不足。
薄壁结构工程梁理论由于模型上的简单长期以来一直被不同程度地广泛应用于飞机设计诸多方面,如刚度计算,强度校核,气动弹性分析,载荷分析等等。然而从国内外目前的文献资料包括各种教科书来看,薄壁梁结构分析目前仍然限于经典的理论分析方法。这种方法尽管模型简单,适合于过去的手工计算,但计算过程较为允长,不同的结构需要不同的处理方法,不太容易编制适用各种结构的统一软件,不仅如此,经典的理论分析方法还存在一些基本缺陷。这些不足具体主要表现在以下几个方面1)需要区分开口薄壁梁和闭口薄壁梁;2)在分析多闭室闭口薄壁梁时,需要先在每个闭室切割开口以变成开口薄壁梁;3)通过引入每个闭室的扭转角相等的关系联立方程,以求解未知量;4)难以对具有大开口的复杂区域进行分析,目前对这些区域的分析需要借助传统有限元方法并伴随复杂的数据结果处理,而且具有很多影响因素;5)目前国内外一直广泛应用的经典的理论分析方法,对扭转刚度的计算存在缺陷,不能反映剛心对扭转刚度的影响。
发明内容
本发明公开了一种薄壁梁结构的优化分析方法,以解决现有技术的上述问题以及其他潜在问题中的任意问题。
为了解决上述技术问题,本发明的技术方案是:一种薄壁梁结构的优化分析方法,所述分析方法具体包括以下步骤:
S1)对待分析的薄壁梁的进行建模,并确定待分析的薄壁梁截面的弯曲主轴;
S2)对所分析的薄壁梁截面在弯曲主轴方向上进行层合化处理,通过等效叠加即可得到总层合结构所对应的在第i层上的宽度;
S3)计算所述薄壁梁截面在每个主轴方向上给定剪切载荷及弯矩时即单向弯剪模型下截面每一层上总的剪流和正应力;
S4)根据S3)得到的薄壁梁截面在单向弯剪下每一层总的剪流和正应力计算所述薄壁梁截面的各个子结构元素在每一层上的剪流和正应力;
S5)根据S4)得到的薄壁梁截面在单向弯剪下各个子结构元素在每一层上的剪流和正应力计算薄壁梁截面的刚心的位置;
S6)根据S4)及S5)计算扭转刚度及在扭矩作用下薄壁梁截面在每个子结构元素上所产生的剪流;
S7)根据S3)、S4)及S5计算在刚心处施加两个独立的沿主轴方向的剪切载荷时每个子结构元素上所产生的剪流。
进一步,所述S2)的具体步骤为:
S2.1)子层合结构的形成;
沿薄壁梁截面的主轴方向在不同的高度将所分析的薄壁结构的截面平行切割成一个层合结构,以此从上到下形成不同种类的元素组,每一个元素组都是一个沿主轴方向的子层合结构;
S2.2)将S2.1)得到每一个子层合结构在第i层上的宽度bik,按照等效原则进行叠加,即可得到总层合结构所对应的在第i层上的宽度bi,公式如下:
式中:Ei为第i层的杨氏模量,Eik为第i层的第k个子结构的杨氏模量,i和k均为大于0的正整数。
进一步,所述S3)的具体步骤为:
S3.1)采用Euler梁假设求解,即平面假设和直法线假设,公式如下:
v=v(x)
σx=Eεx
式中,u、w、εx、σx及σxy分别为x和y方向位移、正应力及剪应力,y0为中性轴的y坐标,
S3.2)根据总体平衡关系及上下边界条件,关系式通过以下公式表达:
σxy=0 at y=ytop
σxy=0 at y=ytop (2a-e)
式中,Mzz为薄壁板截面上的弯矩,Fxy为薄壁板截面上的剪力,ytop和ybot是总层合结构的上表面和下表面的y坐标,b是所考虑层的宽度;
S3.3)求解方程(1a-e)并使用关系式(2a-e)则可求得每一层子层合结构的应力σxi和σxi及中性轴的位置及抗弯刚度EI,公式如下所示:
式中,E为杨氏模量,d为积分点处截面宽度。
进一步,所述S4)的具体步骤为:
S4.1)根据平面假设可知,在同一个层i内不同子结构元素变形一致,所以每个子结构元素的正应变εxik和剪应变γxyik分别都与总的正应变εxi和剪应变γxyi相等,即
εxik=εxi k=1,m
γxyik=γxyi k=1,m (6a-b),
由式6a可得每个子结构的正应力σxi可表示为:
对剪应力σxyik的计算,考虑下面的平衡条件及本构关系式
σxyik=γxyikGik
σxyi=γxyiGi (8a-c),
式中,bik-G为第i层处每个子结构元素在此处承受剪切材料的宽度,
bik-G和Gik为第i层处每个子结构元素在此处承受剪切材料的宽度和剪切模量,Gi为第i层处总的等效剪切模量;
由式(8a-c)及式(6b)最后可得
其中
qyik=σxyikbik-G
qyi=σxyibi (10a-b)
式中,qyik为第i层处每个子结构元素在此处y方向上的剪流,qyi为第i层处总的y方向上的等效剪流。
进一步,所述S5)的具体步骤为:
S5.1)先确定参考点R,再根据剪流对参考点R的扭矩的等值关系即可求得刚心C相对于扭矩参考点R在两个主轴方向上的相对间隔Xp和Yp,由此即可确定刚心的位置,公式如下所示:
VyyXp-VyxYp=TyR
VxyXp-VxxYp=TxR (11a-b)
式中,Vyy为在y方向作用的剪切载荷,Vyx和TyR为Vyy所产生的剪流在x方向的合力及对扭矩参考点R的扭矩,Vxx为在x方向作用的剪切载荷,Vxy和TxR为Vxx所产生的剪流在y方向的合力及扭矩参考点R的扭矩Xp与Yp是刚心在两个平行于弯曲主轴方向上与所给扭矩参考点的间距,
求解(11a-b)即可得到
式中,
αx=Vxy/Vxx
αy=Vyx/Vyy
αxT=TxR/Vxx
αyT=TyR/Vyy (13a-d)
进一步,所述S6)的具体步骤为:
S6.1)假设一个任意薄壁梁截面,相对于其刚心扭转一个角度θ,根据以下公式求出扭矩T,公式如下:
τ=Gγ
式中,ρ为刚心到积分段ds的垂直距离,γ为剪应变,τ为剪应力,G为剪切模量,t为积分点处壁厚,sk为截面所有薄壁中心线段;
由式(14a-c)可得,
式(15b)薄壁结构扭转刚度的计算公式则为本文件所提出的。
由(14a-b)及(15a)即可得剪应力的计算公式
进一步,所述S7)的具体步骤为:
S7.1)给定刚心处沿两个主轴x方向和y方向独立的剪切载荷Vx0和Vy0,单向弯剪模型的输入载荷分别为Vxx和Vyy,
S7.2)在Vxx和Vyy作用下,各自将出现其另外方向上的剪流合力Vxy和Vyx,如式(13a-b)所示。
αx=Vyx/Vxx
αy=Vyx/Vyy (17a-b),
于是根据内力与外载的平衡关系就有
Vyy+αxVxx=Vy0
Vxx+αyVyy=Vx0 (18a-b),
求解式(18a-b),于是就得到下面剪切载荷的转换公式
利用上面的公式就可将给定的剪切载荷转换成单向弯剪模型的输入载荷,从而根据4)就可求得薄壁梁截面的各个子结构元素在每一层上的剪流。考虑到截面非对称时所引起的扭矩耦合,式(19a-b)须作如下修正
其中T为扭矩,系数βTy和βTx为
式中,GJ为截面扭转刚度,G为剪切模量,Sk为截面所有薄壁中心线段,ρ为刚心到积分段ds的垂直距离。
本发明的有益效果是:由于采用上述技术方案,本发明的分析方法具有算统一,不再需要区分开口结构和闭口结构。不再需要对闭口结构尤其是多闭室结构先开口转化成开口结构再进行分析。为开口结构的分析设计提供了有效途径。提出了经典分析方法在计算结构扭转刚度上的理论缺陷;提出了弯剪扭耦合新的处理方法,避免了经典理论在非对称剖面情况下假设纯扭载荷状态所可能带来的不合理性。方法易于进一步改进,如在基础层合板应用更精确的Timoshenko梁假设,就可容易进一步提高整个方法的精确度。
附图说明
图1为本发明一种薄壁梁结构的优化分析方法的流程框图。
图2为一个任意薄壁梁结构示意图。
图3a和图3b为层合化后的子结构示意图。
图4为等效的总层合结构示意图。
图5为第i层上的不同子结构元素。
图6为刚心和扭矩参考点及剪力方向。
图7为一个所分析薄壁截面。
下面结合附图和具体实施例对本发明的技术方案做进一步说明。
如图1所示,本发明一种任意薄壁梁结构的层合化分析方法,所述分析方法具体包括以下步骤:
S1)对待分析的薄壁梁的进行建模,并确定待分析的薄壁梁截面的弯曲主轴;
S2)对所分析的薄壁梁截面在弯曲主轴方向上进行层合化处理,通过等效叠加即可得到总层合结构所对应的在第i层上的宽度;
S3)计算所述薄壁梁截面在每个主轴方向上给定剪切载荷及弯矩时即单向弯剪模型下截面每一层上总的剪流和正应力;
S4)根据S3)得到的薄壁梁截面在单向弯剪下每一层总的剪流和正应力计算所述薄壁梁截面的各个子结构元素在每一层上的剪流和正应力;
S5)根据S4)得到的薄壁梁截面在单向弯剪下各个子结构元素在每一层上的剪流和正应力计算薄壁梁截面的刚心的位置;
S6)根据S4)及S5)计算扭转刚度及在扭矩作用下薄壁梁截面在每个子结构元素上所产生的剪流;
S7)根据S3)、S4)及S5计算在刚心处施加两个独立的沿主轴方向的剪切载荷时每个子结构元素上所产生的剪流,。
所述S2)的具体步骤为:
S2.1)子层合结构的形成;
沿薄壁梁截面的主轴方向在不同的高度将所分析的薄壁结构的截面平行切割成一个层合结构,以此从上到下形成不同种类的元素组,每一个元素组都是一个沿主轴方向的子层合结构;
S2.2)将S2.1)得到每一个子层合结构在第i层上的宽度bik,按照等效原则进行叠加,即可得到总层合结构所对应的在第i层上的宽度bi,公式如下:
式中:Ei为第i层的杨氏模量,Eik为第i层的第k个子结构的杨氏模量,i和k均为大于0的正整数。
所述S3)的具体步骤为:
S3.1)采用Euler梁假设求解,即平面假设和直法线假设,公式如下:
v=v(x)
σx=Eεx
式中,u、w、εx、σx及σxy分别为x和y方向位移、正应力及剪应力,y0为中性轴的y坐标,
S3.2)根据总体平衡关系及上下边界条件,关系式通过以下公式表达:
σxy=0 at y=ytop
σxy=0 at y=ytop (2a-e),
式中,Mzz为薄壁板截面上的弯矩,Fxy为薄壁板截面上的剪力,ybop和ybot是总层合结构的上表面和下表面的y坐标,b是所考虑层的宽度;
S3.3)求解方程(1a-e)并使用关系式(2a-e)则可求得每一层子层合结构的应力σxt和σxi及中性轴的位置及抗弯刚度EI,公式如下所示:
式中,E为杨氏模量,d为积分点处截面宽度。
所述S4)的具体步骤为:
S4.1)根据平面假设可知,在同一个层i内不同子结构元素变形一致,所以每个子结构元素的正应变εxik和剪应变γxyik分别都与总的正应变εxi和剪应变γxyi相等,即
εxik=εxi k=1,m
γxyik=γxyi k=1,m (6a-b),
由式6a可得每个子结构的正应力σxi可表示为:
对剪应力σxyik的计算,考虑下面的平衡条件及本构关系式
σxyik=γxyikGik
σxyi=γxyiGi (8a-c),
式中,bik-G为第i层处每个子结构元素在此处承受剪切材料的宽度,
bik-G和Gik为第i层处每个子结构元素在此处承受剪切材料的宽度和剪切模量,Gi为第i层处总的等效剪切模量;
由式(8a-c)及式(6b)最后可得
其中:
qyik=σxyikbik-G
qyi=σxyibi (10a-b),
式中,qyik为第i层处每个子结构元素在此处y方向上的剪流,qyi为第i层处总的y方向上的等效剪流。
所述S5)的具体步骤为:
S5.1)先确定参考点R,再根据剪流对参考点R的扭矩的等值关系即可求得刚心C相对于扭矩参考点R在两个主轴方向上的相对间隔Xp和Yp,由此即可确定刚心的位置,公式如下所示:
VyyXp-VyxYp=TyR
VxyXp-VxxYp=TxR (11a-b),
式中,Vyy为在y方向作用的剪切载荷,Vyx和TyR为Vyy所产生的剪流在x方向的合力及对扭矩参考点R的扭矩,Vxx为在x方向作用的剪切载荷,Vxy和TxR为Vxx所产生的剪流在y方向的合力及扭矩参考点R的扭矩,
Xp与Yp是刚心在两个平行于弯曲主轴方向上与所给扭矩参考点的间距,
求解(11a-b)即可得到
式中,
αx=Vxy/Vxx
αy=Vyx/Vyy
αxT=TxR/Vxx
αyT=TyR/Vyy (13a-d)。
式中,αx、αy、αxT和αyT均为系数。
所述S6)的具体步骤为:
S6.1)假设一个任意薄壁梁截面,相对于其刚心扭转一个角度θ,根据以下公式求出扭矩T,公式如下:
τ=Gγ
式中,ρ为刚心到积分段ds的垂直距离,γ为剪应变,τ为剪应力,G为剪切模量,t为积分点处壁厚,Sk为截面所有薄壁中心线段;
由式(14a-c)可得,
由(14a-b)及(15a)即可得剪应力的计算公式
所述S7)的具体步骤为:
S7.1)给定刚心处沿两个主轴x方向和y方向独立的剪切载荷Vx0和Vy0,
S7.2)利用S7.1)给定剪切载荷Vx0和Vy0求出单向弯剪模型的输入载荷Vxx和Vyy,公式如下:
S7.3)将S7.2)求出的输入载荷Vxx和Vyy输入单向弯剪模型根据4)就可求得薄壁梁截面的各个子结构元素在每一层上的剪流。
所述S7.2)还可为如果所述薄壁梁截面为非对称时所引起的扭矩耦合,将式(19a-b)须作如下修正
其中T为扭矩,系数βTy和βTx为
式中,GJ为截面扭转刚度,G为剪切模量,Sk为截面所有薄壁中心线段,ρ为刚心到积分段ds的垂直距离。
一种计算机可读存储介质,包括指令,当其在计算机上运行时,使得计算机执行上述的一种任意薄壁梁结构的层合化分析方法。
实施例:
通过“扫描”按照从上到下或从左到右的模式将其分解成一个层合结构(如图2所示)。由于本质上这个层合结构应与普通的层合结构一样,不论在整体上还是在每一层都要遵守相应的平衡条件,变形协调,材料本构关系,及边界条件。所以,结构经过层合化后就可按照层合板的分析方法来求解。通过这样的处理,一个具有任意形状和材料分布的复杂工程梁问题就可转化成为一个具有规则形状和材料分布的层合板问题,从而使问题的求解大大简化,而且便于统一编程计算。
基础层合板结构的分析
如图4所示,一个层合板由n层不同材料不同厚度的材料组成。在此,我们采用Euler梁假设求解,即平面假设和直法线假设。于是就有
v=v(x)
σx=Eεx
式中u、w、εx、σx及σxy分别为x和y方向位移、正应力及剪应力,y0为中性轴的y坐标。
另外还有下面的总体平衡关系及上下边界条件
σxy=0 at y=ytop
σxy=0 at y=ytop (2a-e)
式中,Mzz和Fxy为截面上的弯矩和剪力,ytop和ybot是层合般的上表面和下表面的y坐标,b是所考虑层的宽度。这样,求解方程(1a-e)并使用关系式(2a-e)则可求得每一层的应力σxi和σxi及中性轴或形心位置及抗弯剛度
3任意薄壁结构截面的层合化
3.1子层合结构的形成
对于一个给定的任意薄壁结构截面,沿其一个弯曲中心主轴的方向,如图1所示的厚度方向将其平行切割成一个层合结构。这样,从上到下就形成不同种类的元素组,如左蒙皮(skin L),右蒙皮(skin R),翼梁腹板1(web 1),翼梁腹板2(web 2)以及翼梁腹板3(web3)等。在此,我们每一个元素组都看成是一个沿厚度方向的层合结构,如图3a和3b所示。在物理上它们各自对应一个构件或者构件的一部分,而在理论上它们则各自对应于一个按照如图2所示基础层合结构布局的子层合结构.
这里需要指出的是,所有特殊的节点包括上下和左右的极端点,以及构件之间的交点都需要作为层合化的特殊界面点,否则会引起计算精度的下降。此处,截面中心主轴方向的计算同传统计算方式一样.
图5显示了经过切割后第i层上5个子结构所对应的材料,它们各自的模量为Eik,高度都为hi,宽度为
bik=Liktik/hi k=1,m (4)
这里Lik和tik是元素k的在第i层上的长度和厚度
3.2总层合结构的组装
将每一个子层合结构在第i层上的宽度bik按照E(杨氏模量)*A(截面积)等效原则进行叠加,即可得到总层合板结构所对应的在第i层上的宽度。
式中m为子层合结构的数目,Ei为总层合板结构在第i层所对应的上的模量。
4构件应力应变的计算
得到了总层合板结构每层上的模量、高度及宽度,即可直接应用4.2节基础层合板结构的公式得到此处的层合结构在任意第i层上总的正应力和剪应力σxi及σxyi。同时也可得到弯曲刚度。
根据平面假设可知,在同一个层i内不同子结构元素变形一致,所以每个子结构元素的正应变εxik和剪应变γxyik分别都与总的正应变εxi和剪应变γxyi相等,即
εxik=εxi k=1,m
γxytk=γxyi k=1,m (6a-b)
由式6a可得每个子结构的正应力σxi可表示为
对剪应力σxyik的计算,考虑下面的平衡条件及本构关系式
σxyik=γxyikGik
σxyi=γxyiGi (8a-c)
式中bik-G和Gik表示第i层处每个子结构元素在此处承受剪切材料的宽度和剪切模量,Gi是第i层处总的等效截切模量。
由式(8a-c)及式(6b)最后可得
其中
qyik=σxyikbik-G
qyi=σxyibi (10a-b)
式中qyik表示第i层处每个子结构元素在此处y方向上的剪流,qyi表示第i层处总的y方向上的等效剪流。
式(7)和(9)表明这些子结构元素的正应力在同一层上按照拉伸刚度分配,而剪流按照剪切刚度分配。
5截面刚心的确定
刚心或剪心就是两个方向弯曲剛轴线或剪力线的交点,一般情况下如果截面没有对称轴,则需要分析两个方向的弯曲和剪切特性。按照节1-4所述的步奏,再对应求得薄壁结构在另一个主轴方向受剪时的剪流分布。然后,根矩对剪心处的扭矩平衡关系即可求得剪心的位置。
VyyXp+VyxYp=Typ
VxyXp+VxxYp=Txp (11a-b)
其中Vyy表示在y方向作用的剪切载荷,
Vyx和Typ表示Vyy所产生的剪流在x方向的合力及对剪心的扭矩。Vxx表示在x方向作用的剪切载荷,
Vxy和Txp表示Vxx所产生的剪流在y方向的合力及对剪心的扭矩。表示剪心在弯曲主轴坐标系下的x坐标和y坐标。
6剪切载荷的转换
对于给定的截切载荷,即使相对于弯曲主轴坐标系,也不能直接用于弯剪计算,这里需要有一个转换,因为同弯矩一样,两方向的剪力也有耦合关系。
在此假设在主坐标系下,给定x方向和y方向剪切载荷Vx0和Vy0,对应与弯剪模型的输入载荷分别为Vxx和Vyy。参照上面的分析(11a-b),可知在Vxx和Vyy作用下,各自将出现其另外方向上的剪流合力Vxy和Vyx,我们用下式表示。
αx=Vyx/Vxx
αy=Vyx/Vyy (12a-b)
于是根据内力与外载的平衡关系就有
Vyy+αxVxx=Vy0
Vxx+αyVyy=Vx0 (13a-b)
求解式(13a-b),于是就得到下面剪切载荷的转换公式
利用上面的公式就可将给定的剪切载荷转换成弯剪模型的输入载荷。
7扭转计算
7.1闭剖面扭转计算
经典的薄壁理论分析方法尽管简单,但存在着下面理论上的缺陷:
扭转刚度的计算公式,没有反映剪心的影响。尤其在多闭式情况下,存在明显问题。如下图所示的两种外形尺寸相等但壁厚分布不同的多闭室薄壁梁剖面,由于左右上下对称,可直接确定其剪心位于中心点。根据经典分析方法两种构型具有相同的扭转刚度,但很明显图7所示构型应当具有较大的扭转刚度。
上述说明示出并描述了本申请的若干优选实施例,但如前所述,应当理解本申请并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述申请构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本申请的精神和范围,则都应在本申请所附权利要求书的保护范围内。
Claims (10)
1.一种任意薄壁梁结构的层合化分析方法,其特征在于,所述分析方法具体包括以下步骤:
S1)对待分析的薄壁梁进行建模,并确定待分析的薄壁梁截面的弯曲主轴;
S2)对所分析的薄壁梁截面在弯曲主轴方向上进行层合化处理,通过等效叠加即可得到总层合结构所对应的在第i层上的宽度;
S3)计算所述薄壁梁截面在每个主轴方向上给定剪切载荷及弯矩时即单向弯剪模型下截面每一层上总的剪流和正应力;
S4)根据S3)得到的薄壁梁截面在单向弯剪下每一层总的剪流和正应力计算所述薄壁梁截面的各个子结构元素在每一层上的剪流和正应力;
S5)根据S4)得到的薄壁梁截面在单向弯剪下各个子结构元素在每一层上的剪流和正应力计算薄壁梁截面的刚心的位置;
S6)根据S4)及S5)计算扭转刚度及在扭矩作用下薄壁梁截面在每个子结构元素上所产生的剪流;
S7)根据S3)、S4)及S5计算在刚心处施加两个独立的沿主轴方向的剪切载荷时每个子结构元素上所产生的剪流,完成层合化分析。
3.根据权利要求1所述的分析方法,其特征在于,所述S3)的具体步骤为:
S3.1)采用Euler梁假设求解,即平面假设和直法线假设,公式如下:
v=v(x)
σx=Eεx
式中,u、w、εx、σx及σxy分别为x和y方向位移、正应力及剪应力,y0为中性轴的y坐标,
S3.2)根据总体平衡关系及上下边界条件,关系式通过以下公式表达:
σxy=0 at y=ytop
σxy=0 at y=ytop (2a-e),
式中,Mzz为薄壁板截面上的弯矩,Fxy为薄壁板截面上的剪力,ytop和ybot是总层合结构的上表面和下表面的y坐标,b是所考虑层的宽度;
S3.3)求解方程(1a-e)并使用关系式(2a-e)则可求得每一层子层合结构的应力σxi和σxi及中性轴的位置及抗弯刚度EI,公式如下所示:
式中,E为杨氏模量,d为积分点处截面宽度。
4.根据权利要求3所述的分析方法,其特征在于,所述S4)的具体步骤为:
S4.1)根据平面假设可知,在同一个层i内不同子结构元素变形一致,所以每个子结构元素的正应变εxik和剪应变γxyik分别都与总的正应变εxi和剪应变γxyi相等,即
εxik=εxi k=1,m
γxyik=γxyi k=1,m (6a-b),
由式6a可得每个子结构的正应力σxi可表示为:
对剪应力σxyik的计算,考虑下面的平衡条件及本构关系式
σxyik=γxyikGik
σxyi=γxyiGi (8a-c),
式中,bik-G为第i层处每个子结构元素在此处承受剪切材料的宽度,
bik-G和Gik为第i层处每个子结构元素在此处承受剪切材料的宽度和剪切模量,Gi为第i层处总的等效剪切模量;
由式(8a-c)及式(6b)最后可得
其中:
qyik=σxyikbik-G
qyi=σxyibi (10a-b),
式中,qyik为第i层处每个子结构元素在此处y方向上的剪流,qyi为第i层处总的y方向上的等效剪流。
5.根据权利要求4所述的分析方法,其特征在于,所述S5)的具体步骤为:
S5.1)先确定参考点R,再根据剪流对参考点R的扭矩的等值关系即可求得刚心C相对于扭矩参考点R在两个主轴方向上的相对间隔Xp和Yp,由此即可确定刚心的位置,公式如下所示:
VyyXp-VyxYp=TyR
VxyXp-VxxYp=TxR (11a-b),
式中,Vyy为在y方向作用的剪切载荷,Vyx和TyR为Vyy所产生的剪流在x方向的合力及对扭矩参考点R的扭矩,Vxx为在x方向作用的剪切载荷,Vxy和TxR为Vxx所产生的剪流在y方向的合力及扭矩参考点R的扭矩,
Xp与Yp是刚心在两个平行于弯曲主轴方向上与所给扭矩参考点的间距,
求解(11a-b)即可得到
式中,
αx=Vxy/Vxx
αy=Vyx/Vyy
αxT=TxR/Vxx
αyT=TyR/Vyy (13a-d),
式中,αx、αy、αxT和αyT均为系数。
9.一种实现如权利要求1-8任一项所述的一种任意薄壁梁结构的层合化分析方法的信息处理终端。
10.一种计算机可读存储介质,包括指令,当其在计算机上运行时,使得计算机执行如权利要求1-8任意一项所述的一种任意薄壁梁结构的层合化分析方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110056368.XA CN112800538B (zh) | 2021-01-15 | 2021-01-15 | 一种任意薄壁梁结构的层合化分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110056368.XA CN112800538B (zh) | 2021-01-15 | 2021-01-15 | 一种任意薄壁梁结构的层合化分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112800538A true CN112800538A (zh) | 2021-05-14 |
CN112800538B CN112800538B (zh) | 2023-12-22 |
Family
ID=75809697
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110056368.XA Active CN112800538B (zh) | 2021-01-15 | 2021-01-15 | 一种任意薄壁梁结构的层合化分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112800538B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102844171A (zh) * | 2010-03-26 | 2012-12-26 | 梅西耶-布加蒂-道提公司 | 由拉-压和弯曲下具有改进的机械性能的复合材料制成的机械构件的制造方法 |
CN106547994A (zh) * | 2016-12-01 | 2017-03-29 | 哈尔滨理工大学 | 一种基于层合板理论的复合材料接骨板优化设计方法 |
CN106894328A (zh) * | 2017-02-20 | 2017-06-27 | 重庆大学 | 一种π形结合梁剪力滞的处理方法 |
US20170291328A1 (en) * | 2015-12-02 | 2017-10-12 | Other Lab, Llc | Systems and methods for liner braiding and resin application |
CN107729603A (zh) * | 2017-09-04 | 2018-02-23 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 复合材料壁板螺栓连接的非线性分析方法 |
CN108984841A (zh) * | 2018-06-19 | 2018-12-11 | 威海光威复合材料股份有限公司 | 复合材料层合板等效模量计算及给定载荷下的强度校核 |
CN109141820A (zh) * | 2018-10-08 | 2019-01-04 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 船体模型扭矩及剪力测量方法 |
CN111597638A (zh) * | 2020-05-19 | 2020-08-28 | 上海船舶研究设计院(中国船舶工业集团公司第六0四研究院) | 舷侧大开口船舶总纵剪切强度校核方法 |
-
2021
- 2021-01-15 CN CN202110056368.XA patent/CN112800538B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102844171A (zh) * | 2010-03-26 | 2012-12-26 | 梅西耶-布加蒂-道提公司 | 由拉-压和弯曲下具有改进的机械性能的复合材料制成的机械构件的制造方法 |
US20170291328A1 (en) * | 2015-12-02 | 2017-10-12 | Other Lab, Llc | Systems and methods for liner braiding and resin application |
CN106547994A (zh) * | 2016-12-01 | 2017-03-29 | 哈尔滨理工大学 | 一种基于层合板理论的复合材料接骨板优化设计方法 |
CN106894328A (zh) * | 2017-02-20 | 2017-06-27 | 重庆大学 | 一种π形结合梁剪力滞的处理方法 |
CN107729603A (zh) * | 2017-09-04 | 2018-02-23 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 复合材料壁板螺栓连接的非线性分析方法 |
CN108984841A (zh) * | 2018-06-19 | 2018-12-11 | 威海光威复合材料股份有限公司 | 复合材料层合板等效模量计算及给定载荷下的强度校核 |
CN109141820A (zh) * | 2018-10-08 | 2019-01-04 | 中国船舶科学研究中心(中国船舶重工集团公司第七0二研究所) | 船体模型扭矩及剪力测量方法 |
CN111597638A (zh) * | 2020-05-19 | 2020-08-28 | 上海船舶研究设计院(中国船舶工业集团公司第六0四研究院) | 舷侧大开口船舶总纵剪切强度校核方法 |
Non-Patent Citations (2)
Title |
---|
DIEGO CÁRDENAS等: "Unified theory for curved composite thin-walled beams and its isogeometrical analysis", pages 1 - 17, Retrieved from the Internet <URL:《网页在线公开:https://www.sciencedirect.com/science/article/pii/S0263823117315343》> * |
陈辉等: "硬涂层阻尼结构的随机振动疲劳寿命分析", 《振动、测试与诊断》, vol. 40, no. 6, pages 1212 - 1220 * |
Also Published As
Publication number | Publication date |
---|---|
CN112800538B (zh) | 2023-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108153981B (zh) | 一种基于有限元分析的机身加筋壁板后屈曲分析方法 | |
Chan | Optimal lateral stiffness design of tall buildings of mixed steel and concrete construction | |
CN112417603B (zh) | 一种飞机壁板定位夹紧装配变形预测方法 | |
CN103745066B (zh) | 一种大展弦比机翼结构刚度指标的确定方法 | |
CN112784359B (zh) | 薄壁梁约束扭转极限承载力迭代计算方法 | |
CN110704953B (zh) | 一种大展弦比机翼静气弹性能设计敏度的分析方法 | |
Bechert et al. | Urbach Tower: Integrative structural design of a lightweight structure made of self-shaped curved cross-laminated timber | |
La Magna | Bending-active plates: strategies for the induction of curvature through the means of elastic bending of plate-based structures | |
CN105447269A (zh) | 一种飞机在定常过载状态下的非结构质量配平方法 | |
Ghazlane et al. | Aerostructural adjoint method for flexible wing optimization | |
CN107895086A (zh) | 一种拱坝全坝体的等效应力分析方法 | |
CN108563916B (zh) | 飞行器机翼机身薄壁结构初始尺寸优化设计方法 | |
CN117556568B (zh) | 一种考虑装配偏差的复合材料螺接结构建模方法 | |
CN112800538A (zh) | 一种任意薄壁梁结构的层合化分析方法 | |
Ye | More efficient cold-formed steel elements and bolted connections | |
Groh | Non-classical effects in straight-fibre and tow-steered composite beams and plates | |
Tian et al. | Nonlinear time-varying analysis algorithms for modeling the behavior of complex rigid long-span steel structures during construction processes | |
Hoenlinger et al. | MDO technology needs in aeroelastic structural design | |
CN115221622A (zh) | 一种大尺寸复合材料机身壁板装配定位夹持布局优化方法 | |
Zhao et al. | Reinterpretation of the Dougong joint by the use of parametric tools and robotic fabrication techniques | |
CN110659525B (zh) | 一种复材结构优化设计方法 | |
Dung et al. | Numerical Simulation for the Honeycomb Core Sandwich Panels in Bending by Homogenization Method | |
CN112231864A (zh) | 一种复合材料钉载快速分析系统 | |
Isakson et al. | ASOP-3: A program for optimum structural design to satisfy strength and deflection constraints | |
CN112861069B (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 |