复合材料SHELL181单元完全攻略Word文档下载推荐.docx
- 文档编号:18564750
- 上传时间:2022-12-28
- 格式:DOCX
- 页数:20
- 大小:86.86KB
复合材料SHELL181单元完全攻略Word文档下载推荐.docx
《复合材料SHELL181单元完全攻略Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《复合材料SHELL181单元完全攻略Word文档下载推荐.docx(20页珍藏版)》请在冰豆网上搜索。
181壳单元非常适用于线性、大转角和/或大应变非线性的应用。
计算变厚度壳单元应用非线性分析。
在单元范围内支持完全和减缩的积分方法。
181壳单元计及压力分布引起的(载荷刚度)影响。
181壳单元适用于模拟分层的复合壳或夹层结构。
模拟壳的精度取决于第一剪切变形理论(通常称为Mindling-Reisser壳理论)。
对于用43壳单元存在收敛困难问题,可用181壳单元取代43壳单元。
关于这个单元的更详细的内容见ANSYS,Inc.理论参考。
图181.1181壳单元几何形状
xo= Elementx-axisifESYSisnotsupplied.
为单元坐标系的X轴不提供
x =Elementx-axisifESYSissupplied.
为单元坐标系的X轴提供
2.181壳单元数据输入
这个单元的几何形状、节点位置和坐标系示于图181.1:
“181壳单元几何形状”。
此单元由4点I、J、K和L定义。
单元方程式基于对数应变和真实应力方法。
单元动力学考虑到有限膜应变(拉伸)。
然而,在一个时间增量步内的曲率变化假设很小。
你可以用常数或截面定义它的厚度或其他数据。
仅对单层壳选用实常数。
如果一个181壳单元既有实常数设置,又有一个正确有效的截面类型,则实常数将被忽略。
181壳单元也采用预积分壳截面类型(SECTYPE,GENS)。
当此单元使用GENS截面形式,不需要定义厚度或材料。
更详细内容见“UsingPreintgratedGeneralShellSectins”。
用实常数定义厚度
壳单元的厚度可以在单元各节点定义。
假设整个单元厚度平缓变化。
如果单元厚度不变,只输入TK
(1)。
如果厚度变化,则必须输入4个节点的厚度。
层截面定义
另一种选择,可以用截面命令定义壳单元厚度和更一般的特性。
181壳单元可与壳截面联合(见SECTYPE命令说明)。
与选择实常数相比,壳截面是定义壳结构的更通用的方法。
壳截面命令可用于定义分层复合材料壳的定义,而且提供了输入确定的厚度、材料、方位和通过层厚度的积分点的操作。
注意单层壳不排除用壳截面定义,而且提供更灵活的操作,如使用ANSYS函数编码器定义作为整体坐标和采用积分点的函数的厚度。
当采用截面输入时,你可指定经由每层厚度的积分点数(1,3,5,7或9)。
仅当积分点数为1时,积分点总是位于顶面和低面之间。
如果积分点数为3或更多时,其中两个积分点分别位于顶面和底面,其余积分点在上述两点之间等距分布。
当指定积分点数为5时例外,那里四等份的积分点位置向最近的层面移动5%,使与选用实常数输入选定的位置一致。
每层积分点数的默认值为3。
注意,当采用实常数时,ANSYS采用5个积分点。
然而当使用截面定义等效层时,积分点默认值是3。
为了对解法进行比较,用SECDATA命令设置截面积分点数为5。
此单元的默认方位为S1(壳单元面坐标)轴,它与单元中心的单元第一参数方向一致,它连接侧面中线LI和JK。
在最通过的情况下,该轴可确定为:
对于不扭曲单元,默认方位与在CoordinateSystem中描述的一致(第一表面方向与IJ边一致)。
对于空间翘曲或其他扭曲单元,默认方位代表更好的应力状态,因为在单元定义域内,在默认情况下,单元采用单点求积分。
第一方位S1可以THETA角(度)旋转,作为单元的实常数或者使用SECDATA的命令。
对于一个单元,你可以在单元的平面内确定方位的单一值。
当使用截面定义时,可以使用层向方位。
你也能用ESYS确定单元方位。
见CoordinateSystem。
该单元支持退化的三角形形式,然而,除用作网格填充单元或选用薄膜(KEYOPT1=1)外,不推荐使用三角形形式。
当使用具有大变形而选用薄膜时,三角形单元具有更多的优势。
181壳单元用补偿的方法建立独立转动自由度与(对壳表面)平面位移分量之间的联系。
ANSYS程序以默认值选用一个适当的补偿刚度。
然而如果必要,可以通过改变默认值采用第十个实常数(训练刚度因子;
见表181.1“SHELL181RealConstants”)。
这个实常数值是对默认补偿刚度的比例参数。
采用较高值可能在模型中引起较大的非物理能的成分。
因此改变默认值要慎重。
当采用截面定义181单壳元,训练刚度因子可以用SECCONTROLS命令确定。
单元载荷在NodeandElementLoads中说明。
压力可作为单元表面上的面载荷输入,如图181.1;
“181壳单元几何形状”中带圆圈的数字所表示的各表面。
在单元外表面各角和各层(1—1024最大)之间的相交面各角,温度作为体载荷输入。
第一个温度T1默认为TUNIF。
如果其它角的温度都不作规定,它们默认为T1。
如果KEYOPT
(1)=0并被假设精确地输入NL+1温度,每层底面的四角采用一个温度,而最后的温度用于顶层上表面的四角。
假如KEYOPT
(1)=1并假设精确地输入NL个温度值,每层的四角采用一个温度。
这就是:
T1用于T1,T2,T3,和T4;
T2(当输入时)用于T5,T6,T7,和T8等。
对于其它输入形式,未规定的温度默认值为TUNIF。
采用KEYOPT(3),181壳单元支持一致减缩积分和不相容模型的全积分。
通过默认值,这个单元为了效率,在非线性应用中采用一致减缩积分。
带有计时控制的减缩积分的使用,产生一些限制(虽然很小)。
例如,为了获得悬壁梁或刚性构件的平面弯曲(见图181.2“181壳单元典型弯曲应用”),需要若干厚度方向的单元。
采用一致减缩积分所获得的效率足以补偿采用更多单元的需要。
在比较好地改善网格方面,大都与计时的效果无关。
当选用减缩积分时,你可通过对总能(在ETABLE中的SENE标识)和计时控制引起的人工能(在ETABLE中的AENE标识)进行对比,检查解的精度。
如果人工能与总能之比小于5%,一般说来这个解是可接受的。
总能和人工能也可通过使用在解相位中的OUTPR,VENG进行控制。
双线性单元,在全积分时,平面弯曲过硬。
181壳单元在弯曲为主的问题中,使用不相容模型的方法提高精度。
这个方法也叫“额外状态”或“发泡”型方法。
181壳元采用确保补片试验令人满意的方程式(J.C.Simo和F.Amero,“Geometricallynonlinearenhanced strainmixedsandthemethodofincompatiblemethod”IJNME,VOL,33,PP.1413-1449,1992)。
当分析中包括不相容模型时,你必采用全积分。
KEYOPT(3)=2意味着包括不相容模型和采用(2×
2)全积分。
181壳单元采用KEYOPT(3)=2说明没有任何假设的能机理。
这种181壳元的特定形式即使带有粗糙的网格也是高度精确的。
假如你在选用默认值时遇到任何计时有关的困难,我们向你推荐采用KEYOPT(3)=2。
假如网格粗糙和单元平面弯曲支配响应,KEYOPT(3)=2也是必要的。
在所有的分层的应用中我们推荐这种选择。
KEYOPT(3)=2具有最小的使用限制。
你始终可以选择这种方案。
然而对你的问题可以选择最好的方案,提高单元的效率。
对图181.2“181壳元典型弯曲应用”中说明的问题进行研究。
图181.2181壳元典型弯曲应用
KEYOPT(3)=2对平面弯曲只要求用一个贯穿整个厚度的单元
KEYOPT(3)=0对平面弯曲可能至少要用四个贯穿整个厚度的单元
对加强板(平面)模型采用KEYOPT(3)=2
悬壁梁和用壳模拟横截面的梁是平面弯曲为主的典型例子。
在这种情况下,采用KEYOPT(3)=2是最有效的选择。
减缩积分将要求改善网格。
例如悬壁梁问题采用减缩积分要求四个贯穿整个厚度的单元,而带不相容模型的全积分只要一个贯穿整个厚度方向的单元。
对于加强壳,最有效的选择是:
对壳采用KEYOPT(3)=0,对加强板采用KEYOPT(3)=2。
当规定KEYOPT(3)=0,181壳单元对膜和弯曲模型用一个计时控制方法。
默认值,181壳单元对金属和超弹性应用都用计时参数进行计算。
你可以用实常数11和12取代默认值(见表181.1:
“181壳元实常数”)。
取代改变计时的刚度参数,你应该增加网格密度或者选择全积分(KEYOPT(3)=2)。
当采用截面定义,你可以用SECCONTROLS命令规定计时刚度比例因子。
181壳单元包括横截面剪切变形的影响。
采用Bathe-Drorkin的假设剪切应变公式缓解剪切自锁效应。
单元的横截面剪切刚度是下面所示的一个2×
2的矩阵。
在上述矩阵中,R7,R8和R9是实常数7,8和9(见表181.1:
“181壳单元实常数”)你可以用指定不同的实常数取代横截面剪切刚度的默认值。
这个选择对分析分层壳是有效的。
另一种办法是:
用SECCONTROLS命令定义横截面剪切刚度值。
对于各向同性的单层壳,默认的横截面剪切刚度是:
在上述矩阵中,K=5.6,G=剪切膜量,h=壳的厚度。
181壳单元可以与线弹性、弹塑性、蠕变或高弹性材料特性联系。
只有各向同性、各向异性和正交各向异性的线弹性特性可以作为弹性输入。
vonMises各向同性硬化塑性模型可以同BISO(双线性各向同性硬化),MISO(多线性各向同性硬化)和NLISO(非线性各向同性硬化)方案一起引用。
动态硬化塑性模型同BKIN(双线性动态硬化),MKIN(多线性动态硬化),和CHABOCHE(非线性动态硬化)方案一起引用。
引用塑性假设弹性性质是各向同性(那就是,假设正交各向异性的弹性和塑性一起使用,ANSYS假设各向同性弹性膜量=EX和泊松比=NUXY)。
超弹性材料特性(2,3,5或9参数Mooney-Rivlin材料模型,Neo-Hookean模型,Polynomial型式模型,Aradd-Boyce模型和用户定义模型)可以与此单元同时使用。
泊松比用以规定材料的可压缩性。
如果小于0,泊松比设为0;
如果大于或等于0.5,泊松比设为0.5(完全不可压缩)。
各向同性和正交各向异性的热膨胀系数都可用MP,ALPX输入。
当与超弹性一起使用时,假设各向同性膨胀。
用BETAD命令提供总的阻尼值。
如果MP,DAMP规定单元的材料号(与MAT命令一起指定),它是用于单元取代来自BETAD命令的值。
同样,用TREF命令提供参考温度的总值,假如MP,AEFT是规定单元的材料号,它是用于单元取代来自TREF命令的值。
但是,如果MP,REFT是规定层的材料号,它是用以取代总体的或单元的温度值。
在采用减缩积分和计时控制(KEYOPT(3)=0)时,如果使用的质量矩阵不符合求积分规则,可能出现错误的低频模式。
181壳元使用一种设计方案有效地过滤对单元计时模式的惯性影响。
为了有效,必须使用集中质量矩阵。
对于用这种单元进行模态分析,我们推荐设置LUMPM,OFF。
然而,选用集中质量,能用于带全积分方案(KEYOPT(3)=2)。
在单层或多层的壳单元中,KEYOPT(8)=2用来将中面结果储存到成果文件中。
如果你选用SHELL,MID,你将看到这些计算值,而不是顶面或底面的平均值。
你应当使用这个方案储存这些正确的中面结果(膜结果)。
对于这些分析,平均顶面和底面的结果是不确当的;
例子包括用非线性材料性能求得的中面应力和应变以及涉及诸如谱分析的平方运算的膜态综合的中面结果。
KEYOPT(9)=1用于以用户子程序读入初始厚度数据。
你可以通过ISTRESS或ISFILE命令,对这个单元施加初始应力状态。
更多的资料见在“ANSYSBasisAnalysisGuide”中的“InitialStressLoading”。
此外,你可设置KEYOPT(10)=1以用户子程序USTRESS读取初始应力。
关于用户子程序的详情见“GuidetoANSYSUserProgrammableFeaturess”。
这个单元自动计入压力载荷刚度影响。
如果压力载荷刚度的影响需要一个不对称矩阵,用NROPT,UNSYS。
在“SHELL181InputSummary”中给出这个单元输入的摘要。
单元输入总的说明在“ElementInput”中给出。
SHELL181InputSummary:
2.1181壳单元输入摘要
(1)Nodes—节点
I,J,K,L
(2)DegreesofFreedom—自由度
UX,UY,UZ, ROTX,ROTY,ROTZifKEYOPT
(1)=0
UX,UY,UZifKEYOPT
(1)=1
(3)RealConstants—实常数
TK(I), TK(J),TK(K),TK(L),THETA,ADMSUA
E11, E22,E12,DRILL,MEMBRANE,BENDING
SeeTable181.1:
"
SHELL181RealConstants"
formoreinformation.
更多情况见Table181.1:
If aSHELL181elementreferencesa validshellsectiontype,anyrealconstant dataspecifiedwillbeignored.
如果181壳单元输入确定的壳截面类型,任何实常数定义将被忽略。
(4)MaterialProperties—材料特性
EX,EY,EZ,(PRXY, PRYZ,PRXZ,orNUXY,NUYZ,NUXZ),
ALPX,ALPY,ALPZ (orCTEX,CTEY,CTEZorTHSX,THSY,THSZ),
DENS,GXY,GYZ,GXZ
仅为单元提供一次DAMP(阻尼)(使用MAT命令分配材料性能设置),可以为单元提供一次REFT(温度),或可在每层底分配REFT(温度)。
更多的信息见"
SHELL181InputSummary"
中讨论。
(5)SurfaceLoads—面载荷
Pressures--压力
face1 (I-J-K-L)(bottom,in+Ndirection)(底,在+N方向),
face2 (I-J-K-L)(top,in-Ndirection)(顶,在-N方向),
face3(J-I),face4 (K-J),face5(L-K),face6(I-L)
(6)BodyLoads——体载荷
Temperatures--温度
ForKEYOPT
(1)=0(弯曲和膜刚度):
T1,T2,T3,T4(在1层底),T5,T6,T7,T8(在1-2层之间);
以下各层类似,直到最后顶层NL(4*(NL+1))的温度。
因此,每一层单元施加8个温度值。
ForKEYOPT
(1)=1(仅膜刚度):
T1,T2,T3,T4对1层,T5,T6,T7,T8对2层,类似对所有层(4*NL最大层).因此,每一层单元施加4个温度值。
SpecialFeatures—专用名词
Plasticity塑性
Hyperelasticity超弹性
Viscoelasticity粘弹性
Viscoplasticity粘塑性
Creep蠕变
Stressstiffening应力刚化
Largedeflection大位移
Largestrain大应变
Initialstressimport初始应力输入
Birthanddeath活和死
Automaticselectionofelementtechnology单元技术自动选择
Sectiondefinition forlayeredshellsandpreintegratedshellsectionsfor inputofhomogenoussectionstiffnesses
对层壳的截面定义和对相似截面刚度输入的待积分壳截面
Supportsthefollowing typesofdatatablesassociatedwiththeTBcommand:
使用TB命令支持下列数据表类型:
ANEL,BISO,MISO,NLISO,BKIN, MKIN,KINH,CHABOCHE,HILL,RATE,CREEP,HYPER, PRONY,SHIFT,PLASTIC,andUSER.
Note-注:
材料模型的详细说明见ANSYS,Inc.TheoryReference。
关于单元技术选择的更多情况见AutomaticSelectionofElementTechnologiesandETCONTROL。
2.2开关—各种开关的用法
KEYOPT
(1)
Elementstiffness—单元刚度:
0--Bendingandmembranestiffness(default)—弯曲和薄膜刚度(默认)
1--Membranestiffnessonly—只薄膜刚度
KEYOPT(3)
Integrationoption—积分点:
0--Reducedintegrationwithhourglasscontrol(default)—用计时控制的减缩积分(默认)
2--Fullintegrationwithincompatiblemodes—具有不相容模型的全积分
KEYOPT(8)
Specifylayerdatastorage—指定层数储存:
0--Store dataforbottomofbottomlayerand topoftoplayer(multi-layerelements)(default)—储存底层的底面数据和顶层的顶面数据(多层单元)(默认)
1--StoredataforTOP andBOTTOM,foralllayers(multi-layerelements)—储存所有层的顶面数据和底面数据(多层单元)
Volumeofdatamaybeexcessive—体积数据除外
2--Storedata forTOP,BOTTOM,andMIDforall layers;
appliestosingle-andmulti-layerelements—储存顶面、底面数据以及所有层的MID;
适用于单层和多层单元
KEYOPT(9)
Userthicknessoption—用户厚度选项:
0--No usersubroutinetoprovideinitialthickness(default)—不用用户子程序提供初始厚度(默认)
1--Readinitial thicknessdatafromusersubroutineUTHICK—使用用户子程序UTHICK输入初始厚度
SeetheGuidetoANSYS UserProgrammableFeaturesforuserwrittensubroutines—关于用户子程序输出见ANSYS指南
KEYOPT(10)
User-definedinitialstress—用户定义初始应力:
0--No usersubroutinetoprovideinitialstress(default)—没有用户子程序提供初始应力(默认)
1--Readinitial stressdatafromusersubroutineUSTRESS—用户子程序读入初始应力
SeetheGuidetoANSYS UserProgrammableFeaturesforuserwrittensubroutines—关于用户子程序输出见ANSYS指南。
Table181.1SHELL181RealConstants—SHELL181实常数
No.
Name
Description
1
TK(I)
ThicknessatnodeI节点厚度I
2
TK(J)
ThicknessatnodeJ节点厚度J
3
TK(K)
ThicknessatnodeK节点厚度K
4
TK(L)
ThicknessatnodeL节点厚度L
5
THETA
Angleof firstsurfacedirection,indegrees
第一表面方位角(用度)
6
ADMSUA
Addedmassperunitarea
附加单位面质量
7
E11
Transverseshearstiffness[2]横向剪切刚度[2]
8
E22
9
E12
10
DrillStiffnessFactor
In-planerotationstiffness[1,2]平面转
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 复合材料 SHELL181 单元 完全 攻略