2 弹性力学问题的有限单元法思考题2.1 有限元法离散结构时为什么要在应力变化复杂的地方采用较密网格,而在其他地方采用较稀疏网格?答:在应力变化复杂的地方每一结点与相邻结点的应力都变化较大,若网格划分较稀疏,则在应力突变处没有设置结点,而使得所求解的误差很大,若网格划分较密时,则应力变化复杂的地方可以设置更多的结点,从而使得所求解的精度更高一些。
2.2 因为应力边界条件就是边界上的平衡方程,所以引用虚功原理必然满足应力边界条件,对吗?答:对。
2.3 为什么有限元只能求解位移边值问题和混合边值问题?弹性力学中受内压和外压作用的圆环能用有限元方法求解吗?为什么?答:有限元法是一种位移解法,故只能求解位移边值问题和混合边值问题。
而应力边值问题没有确定的位移约束,不能用位移法求解,所以也不能用有限元法求解。
2.4 矩形单元旋转一个角度后还能够保持在单元边界上的位移协调吗?答:能。
矩形单元的插值函数满足单元内部和单元边界上的连续性要求,是一个协调元。
矩形的插值函数只与坐标差有关,旋转一个角度后各个结点的坐标差保持不变,所以插值函数保持不变。
因此矩形单元旋转一个角度后还能够保持在单元边界上的位移协调。
2.5 总体刚度矩阵呈带状分布,与哪些因素有关?如何计算半带宽? 答:因素:总体刚度矩阵呈带状分布与单元内最大结点号与最小结点号的差有关。
计算:设半带宽为B ,每个结点的自由度为n ,各单元中结点整体码的最大差值为D ,则B=n(D+1),在平面问题中n=2。
2.6 为什么单元尺寸不要相差太大,如果这样,会导致什么结果? 答:由于实际工程是一个二维或三维的连续体,将其分为具有简单而规则的几何单元,这样便于网格计算,还可以通过增加结点数提高单元精度。
在几何形状上等于或近似与原来形状,减小由于形状差异过大带来的误差。
若形状相差过大,使结构应力分析困难加大,误差同时也加大。
2.7 剖分网格时,在边界出现突变和有集中力作用的地方要设置结点或单元边界,试说明理由。
答:有限元处于弹性力学问题的方法是离散法。
它将一个受外力作用的连续弹性体离散成一定数量的有限小的单元集合体,单元之间只在结点上相互联系,即只有结点才能传递力。
所以在边界出现突变和有集中力作用的地方要设置结点和单元边界。
2.8 为什么说三角形三结点单元是常应变单元,如果在每边中点增加一个结点,那么单元内应力如何分布?答:(1)应变矩阵[B]中的参数m j i m j i c c c b b b 、、、、、由坐标变量x 、y 之差确定。
当单元的坐标差确定之后,这些参数与坐标变量x 、y 无关,因此[B]为常量阵。
当单元的结点位移{a}确定后,由[B]转换求得的单元应变都是常量,也就是说在荷载作用下单元中各点具有统一的xy y x γεε、、值。
因此三结点三角形单元称为常应变单元。
(2)如果在每边中点增加一个结点,单元内的应力为线性分布。
习 题2.1试证明x 、y 与面积坐标的关系 证明:设P 点坐标为(x,y )j jpijy x yxy x A ii11121=()()()()[]()y c x b a y x x x y y y x y x xy y x y x y x x y y x m m m i j j i i j j i j i j i j i j i ++=-+-+-=---++=212121同理可求得:由面积坐标定义得:()y c x b a AA A L i i i ijmPjm i ++==21()()y c x b a y x y x y x A y c x b a y x y x yxA j j j iim m pmii i i mm j jpjm++==++==21111212111121()y c x b a AA A L j j j ijm Pmi j ++==21()y c x b a AA A L m m m ijmPij m ++==21由此推出坐标y x 、与面积坐标的函数关系:()()22j i i j j i i j i j j i m j j m m j m j m j j m A c L c L a c a c x b c b c A b L b L a b b a y b c b c ⎧-+-=⎪-⎪⎨-+-⎪=⎪-⎩式(2.1)面积:m i i m j m m j i j j i m j i c b c b c b c b c b c b a a a A -=-=-=++=2代入式(2.1)有:mj j m m j j m j m m j i j j i j i i j j i i j c b c b b a b a L b L b y c b c b c a c a L c L c x --+-=--+-=其中形状参数由下式确定:mj mji m j mji j m m j mm jji x x x x c y y y y b y x y x y x y x a +-==-=-=-==1111代入上式(2.1)可转化为:m m j j i i m m j j i i L y L y L y y L x L x L x x ++=++=再加上 m j i L L L ++=1 所以用面积坐标表示直角坐标矩阵形式如下:⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧m j i m jim j i L L L y y y x x x y x 1111 2.2 试证明两相似三角形的单元刚度矩阵相同。
证明:由于两个三角形相似,故设h A A =21, h 为一常数。
三角形:()111121i j ji i c b c b A -=111111i m j m j i y y b y y b -=-=111111i m j m j i x x c x x c +-=+-=参数Λj i j i c c b b 、、、,只与坐标差有关,所以hb b i i 121=hc c i i 121=单元刚度矩阵通式为:hb b b b s r s r 12211=h A A 121= 故[][]21rs rs K K =所以[][][][][][][][][][][][][]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==mm mjmijm jjjiimij ii Te K K K K K K K K K tA B D B K [][]eeK K 21=因此两相似三角形的单元刚度矩阵相同。
2.3 直角三角形固定在刚性基础上,受齐顶的水压力和自重作用,如图2.14所示。
若按一个单元计算,水的容重g γ,三角形平面构件容重g ρ,取泊松比v =1/6,试求顶点位移和固定面上的反力。
解:按逆时针编码,局部编码与整体编码相同:1-2-3建立坐标())0,0(3)3,0(20,21:a a xoy(1) 求形函数矩阵:[]()⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-+-+-+-+=-=s r s r s r s r s r s r s r s r rsb bc c c b b c b c c b c c b b A Et K 21212121142νννννννa a a a 600321===a b b a b 303321-===a c a c c 220321-===图(2.14) 形函数:)(21y c x b a A N i i i i ++=233221a a a A =⨯⨯=所以:ay a x N a yN a xN 32132321--===形函数的矩阵为:[][][][]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----==a y a x a y a x a y a x a y a x N N N N m ji3210302003210302(2) 刚度矩阵[][][][][][][][][][]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333231232221131211K K K K K K K K K K e[]()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+-+-+-+=-=s r s r s r sr s r s r s r s r rsb bc c c b b c b c c b c c b b A Et K 21212121142ννννννν()125213531416122=-=-==νννa EA Et t可得:[][]⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡=400353534150093532211EK E K [][]⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0251035343127273323531233E K E K[]⎥⎥⎦⎤⎢⎢⎣⎡----=215251935313E K []⎥⎥⎦⎤⎢⎢⎣⎡----=41253535323E K(3)位移列向量和右端项由边界条件可确定:{}{}Te u a 000022υ=水压力和构件厚分别为:10==t ghp γ[]⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=43127425215127332135259414001253503525021525025415019109353E K e{}T Te t l q h q h q R ⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧=03203100203600001自重为W 与支座反力:{}Ty x y x eW R R W WR R R ⎭⎬⎫⎩⎨⎧---=330333112所以:{}Ty x y x eW R h q R W h q W R R R ⎭⎬⎫⎩⎨⎧-+--=33363303011由[]{}{}eeeR a K =得到下列矩阵方程组:⎪⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎨⎧-+--=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧33363000030301122W R h q R W h q W R R u y x y x υ 化简得:⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧-=⎭⎬⎫⎩⎨⎧⎥⎥⎦⎤⎢⎢⎣⎡3640035353022W h q u E υ]⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=431274252151273321352594140012535035250215250254150191009353E K e可得:⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧-=⎭⎬⎫⎩⎨⎧E W E h q u 363567022υ将⎭⎬⎫⎩⎨⎧22υu 代入下式: ⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧-+-=⎭⎬⎫⎩⎨⎧⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡----333425135025103533031122W R h q R W R R u E y x y x υ 固定面上的反力:a h ga gh q 330===γγ从而可得支座反力为:4322123412033011h q W R h q W R W h q R WR y x y x -=-=+=-=2.4 试从式(2.69)说明对角线元素改1法只能用于给定零位移的情形,而对角元素充大数法可以适用任意指定位移的情形。