知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較

2016-10-29 by:CAE仿真在線(xiàn) 來(lái)源:互聯(lián)網(wǎng)

Civil釋放梁端約束的含義及輸入方式比較 MIDAS 朱鋒

1具體問(wèn)題

Civil 程序提供釋放梁端約束功能用于模擬顯示結(jié)構(gòu)中可能出現(xiàn)的鉸接或者半鉸接狀態(tài)。程序提供兩種輸入方式(相對(duì)值、數(shù)值)如下圖所示:

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil分析圖片1

圖 1.1 程序輸入窗口

中文版幫助文檔提供如下說(shuō)明:

選擇類(lèi)型釋放比率決定選定節(jié)點(diǎn)在單元局部坐標(biāo)系中各自由度方向的約束條件。

類(lèi)型:當(dāng)選擇"相對(duì)值"時(shí),輸入釋放后殘留的約束能力的百分比;當(dāng)選擇"數(shù)值"時(shí),輸入釋放后殘留的約束能力的絕對(duì)數(shù)值大小。

選擇某個(gè)方向自由度時(shí),表示將釋放該自由度方向上的約束。在后面的輸入框中可以輸入釋放后殘留的約束能力。例如左圖中,i節(jié)點(diǎn)(N1端節(jié)點(diǎn))彎矩My的值為0.3,當(dāng)選擇"相對(duì)值"時(shí),表示My抗彎剛度的30%有效;當(dāng)選擇"數(shù)值"時(shí),表示My的抗彎剛度在0.3kN*m/[rad]內(nèi)有效。右側(cè)j節(jié)點(diǎn)(N2端節(jié)點(diǎn))彎矩My值為0,表示My抗彎剛度為0,即成為鉸支。

關(guān)于幫助文檔中所描述的殘余約束等詞匯很多用戶(hù)表示不甚理解,究竟具體輸入的數(shù)值大小及程序內(nèi)部如何處理,此幫助中沒(méi)有明確說(shuō)明。故本文通過(guò)簡(jiǎn)單測(cè)試模型詳細(xì)說(shuō)明輸入值的含義及注意事項(xiàng)。

解決方法:

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil分析圖片2

圖 1.2 英文版幫助文檔摘錄

從上述英文版幫助文檔不難看出,程序處理釋放梁端約束實(shí)質(zhì)是在釋放節(jié)點(diǎn)位置建立虛擬彈簧單元,彈簧剛度矩陣為上文中k ?的第一個(gè)矩陣(由于虛擬彈簧的間距無(wú)限小,所以扭轉(zhuǎn)和剪切的相關(guān)性不存在,體現(xiàn)在矩陣中的0項(xiàng))。

考慮到測(cè)試模型需要輸出變形和內(nèi)力,從而可以反算彈簧剛度。測(cè)試模型如下圖所示:

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil分析圖片3

圖 1.3 測(cè)試模型

由于最終結(jié)果輸出以i,j端為輸出位置,加載彎矩荷載采用單元荷載中的集中彎矩,加載位置選擇距離i端很小距離:

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil分析圖片4

圖 1.4 輸入荷載表格

顯然通過(guò)查看計(jì)算后的轉(zhuǎn)角結(jié)果(Ry)以及輸入的外力My可以反算虛擬彈簧的內(nèi)部剛度值(K=My/Ry)通過(guò)輸入“數(shù)值”剛度檢查測(cè)試模型是否有效:

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil技術(shù)圖片5

圖 1.5 輸入固定剛度值

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil技術(shù)圖片6

圖 1.6 位移計(jì)算結(jié)果

輸入剛度值和位移結(jié)果達(dá)到預(yù)期效果:10=1000000/10000。至此可得到程序輸入固定剛度的含義:彎矩和轉(zhuǎn)角之間的比值(剛度的基本定義)。

通過(guò)輸入“相對(duì)值”剛度確定程序內(nèi)部實(shí)際采用剛度“數(shù)值”

通過(guò)上述方法,將輸入方式調(diào)整為按“相對(duì)值”方式輸入剛度。通過(guò)程序計(jì)算得到的位移結(jié)果和已知的外荷載(1000000kNm)可進(jìn)一步測(cè)試相對(duì)值方式程序內(nèi)部采用的彈簧剛度。

具體操作可以采用如下技巧:

(1) 通過(guò)釋放梁端約束表格修改釋放比值:

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil技術(shù)圖片7

(2) 通過(guò)查看表格結(jié)果并修改輸出位移小數(shù)位數(shù)

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil技術(shù)圖片8

圖 1.8 位移結(jié)果表格(右鍵菜單選擇增加顯示小數(shù)位數(shù))

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil仿真分析圖片9

圖 1.9 修改Ry輸出小數(shù)位數(shù)

最終對(duì)比結(jié)果

顯然從單元?jiǎng)偠染仃囍锌梢钥闯隽簡(jiǎn)卧猧端的扭轉(zhuǎn)剛度取值為:4EI/L。根據(jù)上述方法可整理對(duì)比結(jié)果如

下表所示:

表 1.1 反算剛度與桿件剛度對(duì)比表格

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil仿真分析圖片10

結(jié)論:

再次審視英文幫助文檔提供的最終處理方法表格:

釋放梁端約束的含義及輸入方式的比較

10 / 41

知識(shí)庫(kù)-Civil釋放梁端約束的含義及輸入方式比較Midas Civil學(xué)習(xí)資料圖片11

圖 1.10 英文幫助文檔

從表中公式不難看出,程序處理相對(duì)值方式輸入的剛度值,是以定義釋放梁端約束的桿件剛度為基準(zhǔn)的,

本次測(cè)試中以扭轉(zhuǎn)剛度為基準(zhǔn)(4EI/L),修正系數(shù)不是簡(jiǎn)單的乘以用戶(hù)輸入的相對(duì)值,而是通過(guò)r1/(1-r1)

的方式輸入(此處r1為用戶(hù)輸入的相對(duì)值)。

現(xiàn)將整體測(cè)試結(jié)論整理如下:

(1) 絕對(duì)值輸入的剛度就是簡(jiǎn)單的扭轉(zhuǎn)剛度,輸入單位量綱同EI/L。

(2) 相對(duì)值輸入的剛度與選擇的單元有關(guān),特別注意此處由于用到了單元?jiǎng)偠?故剛度取值與單元長(zhǎng)度有關(guān)。如果縮短和增長(zhǎng)單元的長(zhǎng)度,顯然程序內(nèi)部采用的釋放剛度(虛擬彈簧剛度)將發(fā)生變化,整個(gè)結(jié)構(gòu)的位移結(jié)果也將發(fā)生變化(整個(gè)模型的效應(yīng)結(jié)果由于剛度分配原理發(fā)生改變)。

(3) 相對(duì)值輸入的比值不是簡(jiǎn)單的乘以桿件剛度后作為虛擬彈簧的剛度,而是需要用到上述表格的換算公式r1/(1-r1)。

(4) 最后想特別說(shuō)明,單元的剛度與長(zhǎng)度直接相關(guān)(EI/L),虛擬單元由于長(zhǎng)度非常小,其剛度值將非常大,故沒(méi)有釋放梁端約束的節(jié)點(diǎn)位置,如果考慮節(jié)點(diǎn)本身的剛度將類(lèi)似于表格中輸入比例1的情況,即此處相當(dāng)于有一個(gè)剛度非常大的彈簧,其起到連接左右桿件的作用。正因如此對(duì)比表格結(jié)果會(huì)發(fā)現(xiàn),當(dāng)輸入>0.5的比例值后,實(shí)際采用的節(jié)點(diǎn)剛度較原桿件的剛度是增大的。這可能是很多用戶(hù)不甚理解的地方,究其根本主要是理解節(jié)點(diǎn)本身扭轉(zhuǎn)剛度無(wú)窮大的基本想法(現(xiàn)實(shí)中單獨(dú)一個(gè)節(jié)點(diǎn)的剛度是沒(méi)有意義的,此處僅為便于理解)。

(5) 實(shí)際操作中輸入1和0是比較容易理解的。輸入中間值建議根據(jù)比例公式r1/(1-r1)或工程經(jīng)驗(yàn)輸入絕對(duì)值即可,但程序內(nèi)部采用的連接剛度本文已經(jīng)比較詳實(shí)的闡述完畢。

開(kāi)放分享：優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才