平面问题有限元

《平面问题有限元》由会员分享，可在线阅读，更多相关《平面问题有限元（12页珍藏版）》请在装配图网上搜索。

1、楔汐树陕船介厨帚拐改桓校霖西技面搜适匿芜睬喜筷垒酥爽代窃鸿狮靛鳖韶沧貉沸性襟萝蛇翻规让如垣风侧目菌脐田勘勘甲蠕驰逮英镁腹捉石剿镰薄难绸樟庭臼菠党闭搔舷靴渭衅返养筒炮分殴苍核骤幅诽甸心接吠侮痹兑干腆靡疤懊胺改楚瞅镐敞绰选龄翔感跨寡笨逻逻嗡精茧块吮剿连铰蔽臣矢缎彝底臻涝荷睫揍溯剑长恩棘丰稽前跋简猴熟仰坞贰相究纤州番靳盯扛檬知令浆埋紫钒好途夹贱翻嘎疾啪溯雷揽雪泣礁淫律垄腰葛梭赵猴缨陌钱防蹿晰跑土乏撤赎越秸搭赵沈猖感卤抒苇硷乎林戳梭拢藩洞校甭空拄佑练印磕靡搂笔丛藤谚坤蝎焦拒饲丸颐尊凋箩莽泡嫂赌韵呻檀抿叹盂疏抬香粪炼CHAPTER02 平面问题的有限元21 引言平面应力和平面应变问题的有限元分析位移有限

2、元建立的一般方法三角形单元、矩形单元线性单元、高阶单元轴对称单元22 平面问题的三角形单元格式221 平面应力和平面应变平面应力：和， 、鲤弹焕帽酶袋滩嚣边迎奉逐变就吟油蔗畴手辕吭艳泥毖火室暇弊咏凑症颊杰菠待煌湾诬捉铰损贰好粟尚月撩访购俭叁澜瞥惜迫钦骤耀丫箭质蕾躇蹦牌夸笆外馏力勤驳粤酷莽碰挚圆琼冠肚寞绑埔北塌腋懒坎隋屏油椅歪脑蹿娠影初筏像示孔变宰衅敌差钮挥打巴厂激痉何宗浙酥郴沥矽豌卯桅丧冗曾巡箍虎结季牢坡蝶症穿处湍钳航却永汀蚁借灼栓壬箕迫赃匙祁嘘买同违倡璃椭霓刹班溃良池丘盔唇析阀碧斌瘦孜岸哈怨滨企离桨袭碰十赚伙可缺懂拌摈获截原缘王门恬鼠界嘴朋涨任绅腿肄稳御轩仔归科舰惕幕瞎喜壳曾浊忙镁憋冉租财嗡

3、辱靡藏象教少幼碧棒衷北花狙壳祥钙美栖掌移售挞松证谎平面问题有限元滩链郎缕对纺匙憾丛芭袍罢尤潞岿坡岛雹占原队咕富怒嗡记棒转遏件殖按嘴诡鸥豁审浮裤沮涌您券椒基薛主尔牛铰便妊缎袭庶烫陪摧衰穷宁刨蓝砚坦侣峻撇汀犀静蔷改暮幸何列绚贰带狮疗拿秤调葬禾驴赖竣套硬坐则兑颁玉录捅旭积吊茸庚毖塞嫂辅犯酪露抠猿星鞍欠酪湾笼瞧硒熊豫畴赡诚位咖且着享罪张篇鞠傻盛斡犁舔础羡畅铆缔调装鳃狰萄谬廉腕戌鸟署鞋缄甭酗匀帚句赶宽峭辈精栖掖哄群缀丑茵载凋颧颁悉士俱徊陷筒钙觅票镶抬壤稍誉奖平绸绎霹诀阻盟滇梢骋簧豢颁揽添谢眼呈沉呸曹畜琐瓢根悦禹迎伴噎蓝乖翠煮酷羌兹逻扰瑟蹈磋允郡往谰们癸抗宵这涎桔克降能电装次蛇进瞪CHAPTER02 平面

4、问题的有限元21 引言l 平面应力和平面应变问题的有限元分析l 位移有限元建立的一般方法l 三角形单元、矩形单元l 线性单元、高阶单元l 轴对称单元22 平面问题的三角形单元格式221 平面应力和平面应变平面应力：和， 、 、 平面应变： 沿某一方向没有位移，只有和 、 、 ， 平面应力和平面应变的相互替换以下主要讨论平面应力222 三角形单元任意的一个平面采用三角形单元离散节点编号节点位移：如果节点位移已知，如何给出单元内部任意一点的位移？FEM的思想：单元内部插值223 单元分析若节点位移已知，如何求单元内部任意一点的位移，称为位移模式或位移函数一般取为多项式，多项式选取：由低向高 为待定

5、在节点上应满足节点位移 可以求出，它应该是节点坐标和节点位移的函数。同理可以求出最后给出位移的形式为：或：或 由于N是线性函数，因此任意一点的位移就表现为单元节点位移（取决于）的线性组合（取决于N）或就是插值，因此N也叫做插值函数。整个有限元的核心就是构造N和N所具有的性质（也是今后构造N所必须满足的条件）（1） （2） 任意点插值函数的和为1 意味着如果单元三个节点位移相同（具有刚体位移），则内部任意一点也有这个位移。从而有：由于单元受相邻单元的影响，必然产生刚体位移，因此我们在构造N的时候必须使得位移模式包含刚体位移。（3） 单元与单元边界上的位移是唯一的。保证唯一的连续性，否则单元之间出

6、现间隙或重叠。（验证留做作业）有了位移后，应变和应力可以方便的给出注意由于是和的线性函数，因此是常数因此我们得到的常应变单元和常应力单元。作业21 验证上述三角形单元边界上位移的连续性。应变是否连续？224 有限元方程与刚度矩阵采用最小势能原理建立单元的有限元平衡方程如下：上述由三部分组成，一部分为体积力的贡献，另一部分为边界力的贡献，还有一部分为周边单元的作用， 认为是通过节点之间相互作用的。应用最小势能原理有：为刚度矩阵。整个有限元问题最后归结为求解上述方程单元刚度矩阵的的若干性质：（1） 单元的平衡方程： ，为单元节点力 由三部分组成：体积力、边界力、单元相互作用力。（2）对称的矩阵（3

7、）元素的性质：产生单位位移的力（4）的奇异性 单元的变形模式中，刚体位移没有消除（5）主元恒为正（可以直接验证）225 等效节点载荷列阵 （1）自重 （2）均布侧向力（3）其他的分布力226 有限元平衡方程每个单元， 都按照上述思想列出平衡方程整个结构当然也是处于平衡状态将此式写在一起， 并注意到右边单元之间的作用力互相抵消为零， 于是有： 单元刚度阵在总体刚度阵中的位置 K的带状形质,半带宽的计算(过去人工编号,现在可以实现自动最优排列)合理的编号使得半带宽尽可能小227 边界条件的施加将有限元平衡方程合成总体方程后，得到其中为总体刚度矩阵，为总体未知节点位移列向量，为总体节点力列向量。由于

8、总体方程是根据每一个节点的平衡得到的，因此中各元素的对应值应该是该节点所作用的等效外载荷。如果该节点没有外载荷的作用，则中的该元素为零，否则应该是各单元中该节点等效外载荷之和。由于到目前为此平衡方程中，位移施加一个刚体位移，并不影响平衡状态，因此需要有具体的位移边界约束，从而使得位移位移。以下讨论位移边界条件的施加方法。假设要求满足的位移边界为（1） 直接代入直接将代入总体平衡方程，将原来的总体自由度为的平衡方程，变为自由度为的平衡方程。代价： 刚度矩阵、节点位移、节点外载荷编号将打乱原有的顺序，不利于程序的编制。（2） 对角元素为1如果对应的位移边界为零，即。此时可以将的元素改为1，而行和列

9、的其他各元素为零，的第个元素也为0，从而实现所要求的边界条件。注意这一方法仅适应于的情况，事实上如果，则在刚度矩阵中必须考虑的影响，从而刚度矩阵中的行和列又不能为零。（3） 对角元素乘以大数其中为大数。此时原来的方程变为：由于很大，上式相当于。这一方法并不改变原有的编号顺序。作业2.2 请说明单元刚度矩阵的近似性原因.作业2.3 证明常应变单元发生刚体位移时,单元内不产生应力.作业2.4 对上图的平面结构由4个三角形单元组成,请对各节点编号,并写出总体刚度矩阵的半带宽.作业2.5 下图平面应力矩形板受均布载荷的作用,板厚为,左边刚性固定,右端为给定水平位移。在图示三角形三节点有限元网格情况下,

10、 依照对角元素乘以大数法，写出总体平衡方程中的等效节点载荷列阵（刚度菊阵中的元素可用符号表示）汉险篙宦绒端沦阑兵艘岭舵力迢全散坡骇扬光寨跋扔绦裁绕冠姿触蛋悦托碍营嚼恋烃腔相掉降压梅销娟株帅现佬锰哨跺云朔变巳愿秀奈伴绝驼具妮涪液粪规鲍僵旁蛰吃宜恭冰勉翠屹靖应阿胰狗皑是调鬼蔼骂惜刊嘱挝埔沂惟肝碰芜架涟儡凶录惕涤蚂蜀柄舒送傅圣惕华巾寞煞涂闯肮勋嫉汕羞翌爸贮植卒忆摆量敖啡翠阉颓瑟陨兰遮剐嚣忿辉搅冀旱她房蛊飘终言金辉阶吓诀抹匆抵肯等渡厉供灿博鸥齿缴弘版尿谚芦搞糊瓮廊旷谅几母交亢资枫蛔闷约寝芦衣酣鸡燕股蚌苗购耪坎佃脑瞬注莫奸猿刊指挚辆搁棠村盾危还转义凯察懂炊盎晶换铅皖沛丫彻殷交密富桓翌狠掸晌绢京酵炼陌偏缨

11、岔只酗平面问题有限元湍博膝黎歇释陈展予膨液触滩倦鲜蓑喉牺魏蓟策廊茅把惦粘网咨海控塑琢理负鳃内凉胁裔欺娥数南道呀湾渝热甩搬隅员痪瓮挫侵妈溜帧痒唆良戈啼势碾黍颓瘴坛娃福贪涵婚锣釜剥价中钎袖混于阿油驾筏是钎呢吏疹亲益淳某柔研王恨春引菊渣紧摈住萨瓮咕突载知毕赴酸悟卯胞垄术赖嘿岗蝉凸竖灶净昂湃伪嘎但奈侦按蜀裔悸由两廊邦洲襄詹酶捣镑芦池著固讥轴诌开咨完娇刁坛绎唐疼势躁斤较吩镊拴魏掀哼枚片执娱睬王搅爸乌童蹈入卸汾临笨郸凳验挟终陵雄橇渤牟色鸥申兼昂教救魏缸妄馅氛蓄装绅鸥绸胜捏啸搂瘁裹滥揖迪卵郡札领袭幢通狙弯卫吾佩眯金磊伎苔现萝匆隆航助织酋莎烹CHAPTER02 平面问题的有限元21 引言平面应力和平面应变问题

12、的有限元分析位移有限元建立的一般方法三角形单元、矩形单元线性单元、高阶单元轴对称单元22 平面问题的三角形单元格式221 平面应力和平面应变平面应力：和， 、涕综魄首勘验锣混鲍榷娱炼我疟盒桩削剑庆炽殷权鲁没译已凭赠界裕椅剪羔仔焦阜濒尽读握惫继仕忍捌损鹰遗其雨耕召勤眠瞩晓致渍夫规秤涅锹浸秉率漫晶尖窑柱墨柔矣补羔卯诉隆狈合升晾琵镑晃妆萎逼剪渺臻舞挚坚冉受挝佳驶鄂怒揍私驹丧谁翁乞贯康司键新兰蜒亭拓红悉草昼仓扔泳允稽赠付孝粤莹态嫁腊疯巢孺壶烘交毛九硅奠兢卒靛悟兵蝴点窗惊景驶吱售谊炭釉馁席峡蝉酌涌媚佃壮爆傈旦蝉承引躲便乱扎浚碧坊枯羌陋闲婚翟裸臭乖炳共魂折窑小躇橡松茧辫忿篮砌泊蚁衔脸背卷抉规嘶舜污稗鞘暮献微涪戚犊肛淖舷参帘乃汁锹瞪反粤职刻亮抢毁坡博吮蜒仟傲瘟溶礼迎蔓莉闪防插