ansys 焊縫單元的案例
基于ANSYS的某焊接件兩焊縫在順序焊接過程中的分析（生死單元應(yīng)用案例）
關(guān)于生死單元的簡(jiǎn)單介紹
在ansys計(jì)算過程中，如果需要向模型中加入（或刪除）實(shí)體，模型中對(duì)應(yīng)實(shí)體部位的單元就“存在”（或消亡）。單元生死選項(xiàng)就用于在這種情況下殺死或重新激活選擇的單元。例如，在焊接分析過程中，隨著高溫焊料的加入，坡口處的單元需要不斷地被激活；在材料斷料分析中，隨著裂紋的延伸，斷裂處的單元需要不斷的被殺死；在隧道挖掘和橋梁建立分析中，材料也需要不斷的被殺死或激活。因此，單元的生死應(yīng)用技術(shù)廣泛的存在于ansys仿真分析中，是一項(xiàng)應(yīng)用非常廣泛的技術(shù)。
單元的生死并不是ansys程序?qū)⑺?em>單元對(duì)應(yīng)的實(shí)體從模型中刪除，或者激活重新生成材料，而是通過將其剛度矩陣，或者傳導(dǎo)矩陣（對(duì)應(yīng)于不同的分析），乘以很小的因子（ESTIF），默認(rèn)值為1E-6。死單元的單元載荷將為0，從而不對(duì)載荷向量生效，等效于將單元殺死；同樣，當(dāng)一個(gè)單元被重新激活時(shí)，其剛度，單元載荷等恢復(fù)其原始的數(shù)值，重新激活的單元也沒有應(yīng)變記錄，在熱分析里面沒有熱量存儲(chǔ)。需要注意的是，生死單元對(duì)大部分單元可以應(yīng)用，然而對(duì)某些單元卻是不可用的。
在一些情況下，單元生死狀態(tài)可以根據(jù)ansys的計(jì)算結(jié)果決定。如在斷裂分析中，我們需要將應(yīng)力值大于材料屈服強(qiáng)度的單元殺死，可以利用Etable選擇相應(yīng)的單元進(jìn)行殺死，繼而返回到求解器進(jìn)行求解，如果如此循環(huán)，則可觀察到裂紋的生長(zhǎng)過程。
可以在大多數(shù)靜態(tài)和非線性瞬態(tài)分析中使用單元生死，其基本分析與相應(yīng)的分析過程是一致的，主要包括三個(gè)步驟：建模，施加載荷并求解，查看結(jié)果。
今年隨著ANSYS19.0的推出，也帶來了一個(gè)好消息：ANSYS V19.0在Workbench界面下新增了網(wǎng)格生死功能。以往我們只能在經(jīng)典界面下進(jìn)行網(wǎng)格生死操作，或者在Workbench界面下借助APDL來實(shí)現(xiàn)網(wǎng)格生死，這種操作既不方便又容易出錯(cuò)。
展開 ANSYS Mechanical聯(lián)合ANSYS nCode DesignLife 在實(shí)體焊縫疲勞分析
三、ANSYS nCode DesignLife實(shí)體焊縫疲勞分析流程
1、搭建有限元分析基本流程
根據(jù)疲勞載荷工況，搭建靜力學(xué)分析項(xiàng)目，例如圖8中搭建兩個(gè)載荷工況，值得注意的是，nCode DesignLife求解計(jì)算需要網(wǎng)格信息的一致性。
圖8
靜力學(xué)分析結(jié)果如圖9所示，注意這個(gè)等效應(yīng)力結(jié)果不能作為焊縫靜力學(xué)評(píng)價(jià)標(biāo)準(zhǔn)（僅能作為疲勞分析中應(yīng)力分布趨勢(shì)觀測(cè)使用）。
圖9
2、定義焊縫SCL文件
依據(jù)第2節(jié)說明進(jìn)行焊縫文件的定義，注意在焊縫疲勞分析的配置文件中焊縫表達(dá)線的命名需要彼此不同。多條焊縫下“Existing Weldline file selected”需要選擇“Append”。鼠標(biāo)右鍵點(diǎn)選“Evaluate All Results”生成焊縫信息到指定存放路徑，如圖10所示。
圖10
3、實(shí)體焊縫疲勞分析概述
圖11
ANSYS Workbench平臺(tái)并沒有直接的nCode DesignLife焊縫疲勞預(yù)定義分析模塊，但可以間接更改疲勞求解引擎進(jìn)行求解，且采用ANSYS Workbench數(shù)據(jù)管理系統(tǒng)便于管理。
例如按照?qǐng)D11所示引入A6和B6單元格求解信息進(jìn)入C3單元格，修改項(xiàng)目名稱為【nCode SN Solidseam(DesignLife)】，方便識(shí)別工程分析內(nèi)容。
ANSYS nCode DesignLife分析環(huán)境替換預(yù)定義求解引擎為焊縫疲勞引擎“Seam weld CAE Fatigue”，并完成輸入和輸出通道之間連接，如圖11所示。
展開 ANSYS Mechanical聯(lián)合ANSYS nCode DesignLife 在實(shí)體焊縫疲勞分析
作者 | 付穌昇 安世中德結(jié)構(gòu)仿真咨詢專家
首發(fā) | 仿真秀（ID:fangzhenxiu2018）
引言：ANSYS nCode DesigenLife具有強(qiáng)大的焊縫疲勞分析能力，由于分析過程的復(fù)雜性， ANSYS Workbench工作平臺(tái)預(yù)定義7類nCode DesignLife疲勞分析模塊并不包括對(duì)于焊縫疲勞的相關(guān)分析，需要間接完成。
ANSYS nCode DesigenLife焊縫疲勞分析能夠?qū)Ρ”诮Y(jié)構(gòu)進(jìn)行，同時(shí)也能夠基于非薄壁結(jié)構(gòu)進(jìn)行實(shí)體焊縫疲勞模擬，如圖1所示。
實(shí)體焊縫疲勞分析，基于結(jié)構(gòu)應(yīng)力法，對(duì)于實(shí)體網(wǎng)格建立的焊縫分析具有相當(dāng)?shù)钠者m性，相對(duì)于熱點(diǎn)應(yīng)力法，無需對(duì)網(wǎng)格進(jìn)行強(qiáng)制控制。
限于篇幅，本文僅對(duì)實(shí)體焊縫疲勞分析一般流程進(jìn)行概述。
① 基于“DesignLife theory”對(duì)實(shí)體焊縫疲勞分析方法進(jìn)行概述；
② 基于ANSYS Mechanical創(chuàng)建有限元求解；
③ 基于nCode Weldline創(chuàng)建實(shí)體焊縫信息；
④ 基于ANSYS nCode DesignLife進(jìn)行實(shí)體焊縫疲勞求解引擎求解。
圖1
一、實(shí)體焊縫模型創(chuàng)建準(zhǔn)則
1、ANSYS nCode DesignLife實(shí)體焊縫分析方法
ANSYS nCode DesignLife實(shí)體焊縫分析理論中對(duì)于實(shí)體焊縫評(píng)估采用結(jié)構(gòu)應(yīng)力法，與熱點(diǎn)應(yīng)力法（距離焊趾表面一定距離的兩點(diǎn)或三點(diǎn)，進(jìn)行線性或二次插值計(jì)算來確定焊趾處的熱點(diǎn)應(yīng)力值，如圖2所示）相比較，結(jié)構(gòu)應(yīng)力法對(duì)于網(wǎng)格無需特殊考慮，對(duì)網(wǎng)格敏感程度相對(duì)較低。
展開 4/28 Ansys nCode DesignLife焊縫疲勞分析詳解
內(nèi)容簡(jiǎn)介
首先，介紹焊縫疲勞行為特點(diǎn)；進(jìn)而，說明焊縫疲勞分析的名義應(yīng)力法（如：BS7608）和結(jié)構(gòu)應(yīng)力法（如：Volvo (Shell單元)& ASME (Solid單元)基本原理，在Ansys系列軟件中的實(shí)現(xiàn)流程及案例；最后，介紹Ansys Mechanical 近年在處理焊縫建模的功能改進(jìn)以及在Mechanical UI下調(diào)用nCode DesignLife開展焊縫疲勞分析的方法、流程及案例。
面向受眾
重型機(jī)械、風(fēng)電、汽車（零部件）、航空航天、造船、橋梁、電子信息、海洋鉆探及高層建筑等行業(yè)需要對(duì)焊縫結(jié)構(gòu)進(jìn)行強(qiáng)度及疲勞分析的仿真工程師，相關(guān)科研人員及高校師生。
展開 
干貨 | ANSYS Ncode焊縫疲勞壽命評(píng)估方法簡(jiǎn)介
若采用網(wǎng)格點(diǎn)力方法，依據(jù)網(wǎng)格點(diǎn)力數(shù)據(jù)去推導(dǎo)焊趾和焊根單元邊中間點(diǎn)的平均膜應(yīng)力和彎曲應(yīng)力，進(jìn)而可以計(jì)算焊趾和焊根上下表面的法向應(yīng)力，用于做法向方向的疲勞壽命評(píng)估。當(dāng)需要對(duì)焊喉部位進(jìn)行壽命評(píng)估計(jì)算時(shí)，Ncode將基于焊縫單元的兩個(gè)焊縫邊計(jì)算應(yīng)力值，然后平均到中心位置。網(wǎng)格力方法要求采用線性單元。
7. Ncode焊縫疲勞壽命評(píng)估算法評(píng)估了彎曲應(yīng)力對(duì)總應(yīng)力的貢獻(xiàn)度，根據(jù)占比大小取確定，焊縫為剛性或柔性，不同的彎曲力占比，需要采用不同的S-N材料曲線，軟件會(huì)根據(jù)彎曲應(yīng)力比重S-N曲線進(jìn)行自動(dòng)插值處理。
ANSYS Ncode Designlife焊縫疲勞仿真流程
展開 角焊縫（殼體）疲勞在ANSYS nCode DesigenLife的創(chuàng)建與計(jì)算原則淺述
ANSYS nCode DesigenLife焊縫疲勞分析最初用于汽車行業(yè)薄板結(jié)構(gòu)(1-3 mm) 的焊接分析模擬，采用薄殼搭建有限元模型，相關(guān)工業(yè)應(yīng)用也都針對(duì)于此類結(jié)構(gòu)進(jìn)行。ANSYS nCode DesigenLife焊縫疲勞分析采用結(jié)構(gòu)應(yīng)力法進(jìn)行計(jì)算，具有好的網(wǎng)格不敏感性，目前該方法也適用于以實(shí)體建模的焊縫疲勞分析。
限于篇幅本文僅針對(duì)角焊縫（殼體）焊縫單元創(chuàng)建和計(jì)算的準(zhǔn)則基于ANSYS nCode Theory手冊(cè)進(jìn)行編寫，關(guān)于搭接焊縫、激光焊等請(qǐng)參考相關(guān)文獻(xiàn)資料。
兩名筆者水平極為有限，錯(cuò)誤必然較多，另原稿成稿較早且截取原稿部分并非完整，某種程度未能緊跟相關(guān)技術(shù)發(fā)展，因此嚴(yán)禁直接應(yīng)用于企業(yè)項(xiàng)目的產(chǎn)品分析以免造成重大事故和傷害。另外本文建立的焊縫有限元模型不能作為評(píng)估焊縫極限強(qiáng)度的方法進(jìn)行使用。
一、殼體焊縫有限元建模通用原則
不同類型的焊縫形式具有不同的分析方式，需要根據(jù)焊縫種類進(jìn)行分組，每一個(gè)有限元輸入分組應(yīng)對(duì)應(yīng)疲勞引擎中對(duì)應(yīng)的有限元焊縫類型，并設(shè)置一個(gè)合理的參數(shù)數(shù)值。
對(duì)于以薄殼單元建立焊縫有限元建模具有一定的通用準(zhǔn)則：
① 網(wǎng)格應(yīng)以4節(jié)點(diǎn)四邊形單元為主，表達(dá)金屬薄板的中面。
② 以單排或雙排殼單元進(jìn)行焊縫建模表達(dá)。
③ 焊縫網(wǎng)格規(guī)整，尺寸以5mm為最好，規(guī)避三角形網(wǎng)格出現(xiàn)。
④ 疲勞分析焊縫單元需設(shè)置特殊焊接屬性。
⑤ 焊縫單元法向保證設(shè)置法向朝外。
⑥ 毗鄰焊縫的單元的非平均化節(jié)點(diǎn)應(yīng)力被提取作為焊趾和焊根疲勞計(jì)算評(píng)估使用，該應(yīng)力也可以是平均化的或在單元邊長(zhǎng)的中點(diǎn)處進(jìn)行計(jì)算，通過在“ANSYS Group Properties”中設(shè)置“WeldLocation = MidElementEdge”進(jìn)行考慮。
展開 斯姆勒精品案例：基于ANSYS子模型技術(shù)的焊縫結(jié)構(gòu)的精細(xì)化計(jì)算
基于ANSYS子模型技術(shù)的焊縫結(jié)構(gòu)精細(xì)化計(jì)算
掌握ANSYS焊縫子模型分析技巧
●技術(shù)背景
焊縫（welded seam）利用焊接熱源的高溫，將焊條和接縫處的金屬熔化連接而成的縫。焊縫金屬冷卻后，即將兩個(gè)焊件連接成整體。根據(jù)焊縫金屬的形狀和焊件相互位置的不同，分對(duì)接焊縫、角焊縫、塞焊縫和電鉚焊等；
焊接失效就是焊接接頭由于各種因素，在一定條件下斷裂（如：應(yīng)力、溫度、材質(zhì)、焊接質(zhì)量和實(shí)際使用工況條件等）。接頭一旦失效，就會(huì)使相互緊密聯(lián)系成一體的構(gòu)件局部分離、撕裂并擴(kuò)展，造成焊接結(jié)構(gòu)損壞，致使設(shè)備停機(jī)，影響正常生產(chǎn)。；
焊接失效
（1）因設(shè)計(jì)不合理，存在局部剛性過大，應(yīng)力集中的現(xiàn)象。
（2）材料缺陷。鑄鋼件相對(duì)于軋制板材存在著沖擊韌度差，屈服強(qiáng)度低的特點(diǎn)，還有焊接工藝制定不合理、焊接規(guī)范的運(yùn)用不當(dāng)、焊接方法的選擇不正確等。
（3）焊工技術(shù)水平高低與焊接位置的好壞；還有焊接檢驗(yàn)水平，包括對(duì)材質(zhì)的檢驗(yàn)和焊縫檢驗(yàn)等。另外，環(huán)境溫度對(duì)焊接質(zhì)量也是一個(gè)重要的影響因素。
展開 基于ansys的梁單元、實(shí)體單元徐變精細(xì)化分析（含各參數(shù)解釋） ￥25
2、改網(wǎng)格模型，改成自己對(duì)應(yīng)的網(wǎng)格模型，網(wǎng)格用ansys，hypermesh，ansa等前處理軟件都沒問題。
3、改材料參數(shù)，改成你想要的徐變模型，對(duì)著規(guī)范或者是你做出來的試驗(yàn)擬合曲線。
以上即可實(shí)際應(yīng)用。
ANSYS單元類型選擇方法 附ansys結(jié)構(gòu)單元與材料應(yīng)用手冊(cè)下載
六、單元類型選擇方法
7.進(jìn)行完前面的選擇工作，單元類型就基本上已經(jīng)定位在2-3種單元類型上了，接下來打開這幾種單元的幫助手冊(cè)，進(jìn)行以下工作：
仔細(xì)閱讀其單元描述，檢查是否與分析問題的背景吻合、
了解單元所需輸入的參數(shù)、單元關(guān)鍵項(xiàng)和載荷考慮；
了解單元的輸出數(shù)據(jù)；
下載地址：ansys結(jié)構(gòu)單元與材料應(yīng)用手冊(cè)
ANSYS中桿單元和殼單元的單元耦合問題
在比較復(fù)雜的結(jié)構(gòu)的有限元分析中，不同的結(jié)構(gòu)部件通常使用不同類型的單元來模擬。
通常情況下，不同類型的單元的各個(gè)節(jié)點(diǎn)的自由度數(shù)目是不同的，不同類型單元的連接節(jié)點(diǎn)處的自由度的耦合問題，是一個(gè)比較令人頭疼的問題。
在ANSYS中通?？梢杂民詈厦頒P來耦合不同類型單元在連接節(jié)點(diǎn)處的自由度(DOF)。
也可以用CE命令來認(rèn)為添加自由度之間的約束方程來達(dá)到耦合的目的。
下面是一個(gè)簡(jiǎn)單的算例，使用了CE命令來耦合連接節(jié)點(diǎn)處的自由度。
模型是航天器的機(jī)翼的一個(gè)Section的某一個(gè)隔框。上下表皮是薄殼結(jié)構(gòu)，用Shell63單元來模擬，在上下表皮之間有起支撐作用的桿件，用link8單元來模擬。
建模的時(shí)候，link8單元和shell63單元在連接有各自獨(dú)立的節(jié)點(diǎn)。即：link8單元和shell63單元的節(jié)點(diǎn)在連接處是重合的，但是，節(jié)點(diǎn)編號(hào)是各自獨(dú)立的。
link8單元在每個(gè)節(jié)點(diǎn)有 ux,uy,uz3個(gè)平動(dòng)自由度；
shell63在每個(gè)節(jié)點(diǎn)有ux,uy,uz這3個(gè)平動(dòng)自由度和rotx,roty,rotz這3個(gè)轉(zhuǎn)個(gè)自由，共6個(gè)自由度。
在耦合節(jié)點(diǎn)處，兩個(gè)耦合節(jié)點(diǎn)的ux,uy,uz自由度應(yīng)該是相等的。
這個(gè)等式可以用CE命令來描述。
完整的命令流如下：
finish
/clear,start
/prep7
!定義第一種材料屬性；
mp,ex,1,30e6
mp,prxy,1,0.3
!定義shell63單元和實(shí)常數(shù)；
et,1,shell63
r,1,1e-3
!建立幾何模型；
rectng,31.8,33.2,0,0.3556
agen,2,1,1,1,0,0,1
a,1,4,8,5
a,6,7,3,2
KL,7,0.5, ,
KL,3,0.5, ,
在關(guān)鍵點(diǎn)處生成節(jié)點(diǎn)；
nkpt,100,4 !與編號(hào)為117的節(jié)點(diǎn)耦合
nkpt,101,9 !
展開 ANSYS各類型單元連接專題講解(五)之3D梁單元與殼單元剛接
例如采用ANSYS模擬一個(gè)多層混凝土框架結(jié)構(gòu)，一般除計(jì)算整體指標(biāo)外，我們?cè)谟?jì)算具體荷載作用時(shí)（如風(fēng)荷載、地震作用、恒載、活載等），樓板一般采用彈性版，此時(shí)可用殼單元模擬，主梁、次梁采用梁單元模擬，此時(shí)變?yōu)榱?em>單元包含在殼面內(nèi)的情況，當(dāng)然此類情況是否需要考慮截面偏置，可根據(jù)具體工程而定。
對(duì)這中梁單元包含在殼單元面內(nèi)的情況，只需要將梁單元與殼單元共用節(jié)點(diǎn)即可，而無須格外建立約束方程。
三、梁單元在殼單元內(nèi)但不包含
此種情況為梁與殼位于同一面內(nèi)，但其中面不包含梁線，適用于多尺度建模分析（如下圖）。梁單元與殼單元的連接在端部可以通過剛性梁和剛性區(qū)域兩種方式連接。剛性梁采用MPC184單元，剛性區(qū)域采用Cerig命令，具體使用方法下期文章討論。
展開 
ANSYS APDL實(shí)體單元和殼單元（不共節(jié)點(diǎn)）之間的連接 ￥100
實(shí)體單元和殼單元之間的連接是ANSYS中常見的問題。即使兩種單元之間共節(jié)點(diǎn)，但單元之間不連續(xù)（實(shí)體單元每個(gè)節(jié)點(diǎn)有3個(gè)平動(dòng)自由度，而殼單元每個(gè)節(jié)點(diǎn)有3個(gè)平動(dòng)自由度和3個(gè)轉(zhuǎn)動(dòng)自由度），對(duì)于兩種單元之間面面接觸，可直接定義剛域，本文主要采用MPC法對(duì)實(shí)體-殼單元的連接方法進(jìn)行說明。
1 單元類型
算例模型中，實(shí)體單元采用SOLID45，殼單元采用SHELL63，接觸位置不共節(jié)點(diǎn)。對(duì)于兩種單元之間的連接，通過目標(biāo)單元TARGE170和接觸單元CONTA175實(shí)現(xiàn)，定義約束為實(shí)體-殼約束，接觸單元為MPC算法，接觸類型為綁定接觸。
2 有限元模型和綁定接觸
圖1 底部固定約束，殼單元施加均布荷載
圖2 目標(biāo)單元和接觸單元
3 計(jì)算結(jié)果
圖3 von Mises stress
圖4 X-Component of displacement
付費(fèi)內(nèi)容為相關(guān)命令流。
展開 在ANSYS 中3維坐標(biāo)下的 shell structure 使用2D 平面單元劃分，應(yīng)該使用哪個(gè)單元型號(hào)的單元
在ANSYS 中3維坐標(biāo)下的 shell structure 使用2D 平面單元（僅考慮平面內(nèi)的位移）劃分，應(yīng)該使用哪個(gè)單元型號(hào)的單元？
ANSYS梁單元與實(shí)體單元的耦合與約束方程
ANSYS梁單元與實(shí)體單元的耦合與約束方程
By長(zhǎng)安CAE
1 概述
在ANSYS計(jì)算過程中，有時(shí)候會(huì)遇到不同單元之間進(jìn)行連接，由于不同的單元自由度不同，連接時(shí)通常需要通過耦合和約束方程建立節(jié)點(diǎn)自由度的關(guān)系，保證結(jié)果的準(zhǔn)確性。
耦合可以理解成是將耦合的對(duì)象某個(gè)自由度作相等處理，而約束方程則不局限于相等這個(gè)關(guān)系，其可以描述具有某種關(guān)系的自由度。如圖1所示，為梁單元與平面單元的連接。如果不采用約束方程，力矩的傳遞無法完成，因?yàn)槠矫?em>單元沒有轉(zhuǎn)動(dòng)自由度。
圖1 梁單元與平面單元連接
為使節(jié)點(diǎn)2具有力矩傳遞的能力，要求1、2、3節(jié)點(diǎn)之間的自由度滿足以下關(guān)系：
ROTZ2 = (UY3 - UY1)/10
再通過CE命令，即可將此關(guān)系通過約束方程的形式施加給1、2、3節(jié)點(diǎn)。
2 命令
查看ANSYS的幫助文檔，查詢CE命令的解釋，如圖2所示。
圖2 ANSYS的CE命令解釋
CE, NEQN, CONST, NODE1, Lab1, C1, NODE2, Lab2, C2, NODE3, Lab3, C3
其中，NEQT表示常數(shù)，用于區(qū)別約束方程，一般可以用數(shù)字1、2、3表示即可，表示第幾個(gè)約束方程；
CONST表示方程的常數(shù)項(xiàng)，一般為0；
NODE1，表示第一個(gè)節(jié)點(diǎn)；
Lab1，表示自由度標(biāo)簽，對(duì)于結(jié)構(gòu)而言，就是三個(gè)平移和三個(gè)轉(zhuǎn)動(dòng)自由度；
C1，表示該自由度的系數(shù)；
同理，后面的也一樣。
展開 ANSYS特殊單元——Follw201（隨動(dòng)荷載）單元
ANSYS的Follw201單元是ANSYS的幾個(gè)特殊單元（比如mesh200）之一，稱為隨動(dòng)荷載單元。都知道在ansys里面施加壓力載荷pressure時(shí)，其實(shí)載荷是可以隨動(dòng)的，也就是能夠一直保持著面的法線方向，而施加集中了或者力矩時(shí)則不能保證。Follw201單元便能解決這個(gè)問題。
Follw201單元是一個(gè)單節(jié)點(diǎn)的3D單元，具有六個(gè)自由度，只能夠覆蓋在既有單元節(jié)點(diǎn)上，而且節(jié)點(diǎn)必須具有3個(gè)平移自由度和3個(gè)轉(zhuǎn)動(dòng)自由度，也即是只能用在梁單元和殼單元上，實(shí)體單元僅有三個(gè)自由度。
Follw201單元主要用于幾何非線性分析問題中，在這類問題的分析過程中幾何會(huì)發(fā)生比較大的變形，面或線的法線方向可能發(fā)生比較大的變化，施加的載荷的方向是否隨動(dòng)對(duì)結(jié)果的影響非常大。也就是用到此單元時(shí)會(huì)配對(duì)使用Nlgeon,on命令以打開大變形開關(guān)。如下圖所示為單元示意圖。
圖1
每個(gè)單元有兩個(gè)面，面1用于設(shè)定集中力的大小，面2用于設(shè)置力矩的大小，面的方向在應(yīng)用時(shí)是通過單元的實(shí)常數(shù)進(jìn)行定義的。
另外還需要注意，有限元求解的時(shí)候大部分是求解對(duì)稱矩陣，但隨動(dòng)荷載單元的應(yīng)用則包括了隨動(dòng)荷載剛化效應(yīng)，使剛度矩陣為非對(duì)稱的，因此需要采用非對(duì)稱求解器進(jìn)行計(jì)算。
下面是具體應(yīng)用，建一根梁單元，在梁的端部施加隨動(dòng)集中力。
/prep7
!定義參數(shù)
EE=207E3
B=10
LCD=300
AA=B*B
IZ=B**4/12
PHZ=EE*IZ/LCD/LCD
!定義單元和材料
!201單元不需要定義材料
et,1,beam4
et,2,follw201
mp,ex,1,ee
mp,prxy,1,0.3
!定義實(shí)常數(shù)，實(shí)常數(shù)1設(shè)置梁單元的參數(shù)
r,1,aa,iz,iz,b,b
r,2,,1.0
!
展開 