4.弹性力学中的有限元法2-轴对称与空间问题
- 格式:pdf
- 大小:1.51 MB
- 文档页数:66
弹性力学教学中的几个问题李跃宇(淮海工学院机械学院连云港222005)摘要:弹性力学与材料力学是不同的两个力学学科,教师上课时把弹性力学的重要性、与材料力学的区别、在力学中与有限元的关系、弹性力学的几个层次讲清楚,学生才能听得进去。
弹性力学的重要性在于比材料力学有更广的实用范围并为以后有限元的学习打下了基础。
弹性力学放弃了平面假设,方程比较复杂。
弹性力学平面问题有三个层次:应力、应变和位移所满足的一般方程;按应力求解的力法要满足平衡方程、相容方程和应力边界条件;力法中用逆解法或半逆解法求解的应力函数要满足相容条件、用应力函数表示的应力要满足应力边界条件。
关键词:弹性力学;材料力学;有限元法;力法1 弹性力学的重要性弹性力学实用范围比材料力学更广。
材料力学只能解决杆系问题,而且必须是长杆问题。
弹性力学能解决许多材料力学不能解决的问题,比如孔口问题、深梁问题。
有的材料力学书里列出孔口问题的解,但实际上是从弹性力学中移植过来的。
浅梁问题,材料力学可以有较精确的解;但是深梁问题,材料力学只能得到近似解。
梁越深,材料力学解的误差就越大。
弹性力学是断裂力学的基础,裂缝问题就是孔口问题的一个特例。
对高强度的钢材、焊接结构、大型锻件,就需要进行断裂力学的分析,没有弹性力学的基础,就不可能很好地应用断裂力学。
弹性力学也是有限元的基础。
许多重要结构(飞机、船舶、汽车等)的结构设计都少不了有限元的分析。
有限元的分析可以有效地指导设计工作,大大节约了工作量和试验经费。
弹性力学尽管有比材料力学应用更广的理论框架,但是用弹性力学分析的解析解不是很多,这就要用到有限元来分析。
如果不学弹性力学,有限元根本无法正确地理解和应用。
有限元软件的学习是用考题来取得进步的,考题就要用弹性力学解析解与有限元法的结果进行比较。
2 弹性力学与材料力学的区别弹性力学和材料力学所用到的平衡方程都是理论力学中得到的,材料力学是用内力表示的,但是弹性力学平衡方程是用应力表示的。
4. 彈性力學軸對稱問題的有限元法本章包括以下內容:4.1用虛功方程建立有限元方程 4.2三結點單元位移函數 4.3三結點單元剛度矩陣 4.4載荷移置4.5軸對稱分析舉例4.1用虛功方程建立有限元方程物體的幾何形狀、約束情況及所受的外力都對稱於空間的某一根軸,因此在物體中通過該軸的任何平面都是對稱面,所有應力、應變和位移也對稱於該軸,這類問題稱為軸對稱問題。
研究軸對稱問題時通常採用圓柱坐標系(r ,θ,z ),以z 軸為對稱軸。
圖4.1受均布內壓作用的長圓筒如圖4.1所示的受均布內壓作用的長圓筒,通過Z 軸的一個縱截面就是對稱面。
由於對稱性,軸對問題共有4個應力分量:⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧=zrz r τσσσσθ}{ (4-1)其中r σ表示沿半徑方向的正應力,稱為徑向應力;θσ表示沿θ方向的正應力,稱為環向應力或切向應力;z σ表示沿z 方向的正應力,稱為軸向應力;zr τ表示在圓柱面上沿z 方向作用的剪應力。
同樣,軸對稱問題共有4個應變分量:⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧=zrz r γεεεεθ}{ (4-2)其中r ε表示沿半徑方向的正應變,稱為徑向正應變;θε表示沿θ方向的正應變,稱為環向正應變或切向正應變;z ε表示沿z 方向的正應變,稱為軸向正應變;zrγ表示沿r 和z方向的剪應變。
在軸對稱問題中,彈性體內任意一點上,不存在切向位移,只存在徑向位移u 和軸向位移w ,兩個位移分量表示為,⎭⎬⎫⎩⎨⎧=w u f }{(4-3)在討論彈性力學平面問題的有限元法時,我們先由將彈性體劃分為有限個單元的組合體,由虛功方程得到單元剛度矩陣,集成後得到整體剛度矩陣。
在這裏,我們用虛功方程直接得到軸對稱問題的有限元列式。
由虛功方程可得,外力虛功等於內力虛功或虛應變能, ds p f dxdydz F f dxdydz TsTT}{}{}{}{}{}{***⎰⎰⎰⎰⎰⎰⎰⎰+=σε(4-4)其中{F}為體力,{p}為面力。
第四章 轴对称问题有限元法在工程问题中经常会遇到一些实际结构,它们的几何形状、约束条件和外载荷均对称某一固定轴,我们把该固定轴称为对称轴。
则在载荷作用下产生的应力、应变和位移也都对称此轴。
这种问题就称为轴对称问题。
在离心机械、压力容器、矿山机械、飞行器中经常遇到轴对称问题。
第一节 轴对称问题弹性力学基本方程对于轴对称问题,宜采用圆柱坐标系(,,r z θ)。
如果将y弹性体的对称轴作为Z 轴,则所有应力、应变和位移分量都只是r 和Z 轴的函数,而与θ无关,即不随θ变化。
弹性体内任意一点只有两个位移:即沿r 方向的径向位移u 和沿Z 方向的轴向位移w 。
由于轴对称,沿θ方向的环向(周向)位移v 等于零。
因此轴对称问题是二维问题。
在轴对称弹性体内用相距dr 的两个圆柱面和过轴线互成d θ角的两个铅垂面切割出一个高为dz 的微元体,如图2所示。
(a )σ(b)沿r 方向作用的正应力r σ称为径向应力 沿θ方向作用的正应力θσ称为环向应力 沿z 方向作用的正应力z σ称为轴向应力 rz 面内的剪应力 zr τ=rz τ 故轴对称弹性体内任意一点的应力分量{}[]Tr z rz θσσσστ=对应的轴对称弹性体内任意一点的应变分量{}[]Tr z rz θεεεεγ=其中r ε ———-—- 沿r 方向径向线应变θε —----— 沿θ方向环向线应变 z ε ---—-- 沿z 方向轴向线应变rz γ-—---— rz 面内的剪应变与平面问题相比,轴对称问题多了一个环向应变θε。
弹性体受载时,点(,,r z θ)产生径向位移u ,使过点(,,r z θ)的周长增加了2()2r u r ππ+-,因而产生相对伸长,即环向应变:2()22r u r u r rθππεπ+-==轴对称问题的几何方程(应变与位移之间的关系)为,,,r z zr u u w w ur r z r zθεεεγ∂∂∂∂====+∂∂∂∂写成矩阵形式{}r z rz u r u rw z u w z r θεεεεγ⎧⎫⎪⎪⎪⎪⎧⎫⎪⎪⎪⎪⎪⎪⎪⎪⎨⎬⎨⎬⎪⎪⎪⎪⎪⎪⎪⎪⎩⎭⎪⎪⎪⎪⎩⎭==∂∂∂∂∂∂+∂∂根据虎克定律,应力与应变的关系为1()r r z Eθεσμσσ⎡⎤⎣⎦=-+ 1()z r Eθθεσμσσ⎡⎤⎣⎦=-+ 1()z z r Eθεσμσσ⎡⎤⎣⎦=-+ 12(1)rz rz rz r G Eμττ+==由上式得[]10111011(1)(1)(12)101112002(1)r z zr r z rz E θθσσσστμμμμεμμεμμμμμεμμμμγμμ⎧⎫⎪⎪⎪⎪=⎨⎬⎪⎪⎪⎪⎩⎭⎡⎤⎢⎥--⎢⎥⎧⎫⎢⎥⎪⎪⎢⎥---⎪⎪⎢⎥⎨⎬+-⎢⎥⎪⎪⎢⎥⎪⎪--⎩⎭⎢⎥-⎢⎥⎢⎥-⎣⎦=(4—2)这里弹性矩阵[D]为[D ]=10111011(1)(1)(12)101112002(1)E μμμμμμμμμμμμμμμμμ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦-----+-----第二节 三角形截面环单元一、 结构离散化离散化轴对称体时,采用的单元是一些圆环。