LS-dyna能量介紹

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

LS-DYNA FAQ 7.Energy balance 能量平衡

time........................... 4.99735E-03

time step...................... 4.45000E-06

kinetic energy................. 3.80904E+09

internal energy................ 5.15581E+09

spring and damper energy....... 1.00000E-20

hourglass energy .............. 1.34343E+08

system damping energy.......... 0.00000E+00

sliding interface energy....... 1.72983E+07

external work.................. 4.54865E+09

eroded kinetic energy.......... 0.00000E+00

eroded internal energy......... 0.00000E+00

total energy................... 9.11649E+09

total energy / initial energy.. 1.09716E+00

energy ratio w/o eroded energy. 1.09716E+00

global x velocity.............. -6.63878E+01

global y velocity.............. 3.44465E+02

global z velocity.............. -1.86129E+04

time per zone cycle.(nanosec).. 11286

GLSTAT(參見(jiàn)*database_glstat)文件中報(bào)告的總能量是下面幾種能量的和:

內(nèi)能 internal energy

動(dòng)能 kinetic energy

接觸(滑移)能 contact(sliding) energy

沙漏能 houglass energy

系統(tǒng)阻尼能 system damping energy

剛性墻能量 rigidwall energy

GLSTAT中報(bào)告的彈簧阻尼能”Spring and damper energy”是離散單元(discrete elements)、安全帶單元(seatbelt elements)內(nèi)能及和鉸鏈剛度相關(guān)的內(nèi)能(*constrained_joint_stiffness…)之和。而內(nèi)能”Internal Energy”包含彈簧阻尼能”Spring and damper energy”和所有其它單元的內(nèi)能。因此彈簧阻

尼能”Spring and damper energy”是內(nèi)能”Internal energy”的子集。

由SMP 5434a版輸出到glstat文件中的鉸鏈內(nèi)能”joint internal energy”跟*constrained_joing_stiffness不相關(guān)。它似乎與*constrained_joint_revolute(_spherical,etc)的罰值剛度相關(guān)連。這是SMP 5434a之前版本都存在的缺失的能量項(xiàng),對(duì)MPP 5434a也一樣。這種現(xiàn)象在用拉格朗日乘子(Lagrange Multiplier)方程時(shí)不會(huì)出現(xiàn)。

與*constrained_joint_stiffness相關(guān)的能量出現(xiàn)在jntforc文件中,也包含在glstat文件中的彈簧和阻尼能和內(nèi)能中?；叵霃椈勺枘崮堋眘pring and damper energy”,不管是從鉸鏈剛度還是從離散單元而來(lái),總是包含在內(nèi)能里面。

在MATSUM文件中能量值是按一個(gè)part一個(gè)part的輸出的(參見(jiàn)*database_matsum)。沙漏能Hourglass energy僅當(dāng)在卡片*control_energy中設(shè)置HGEN項(xiàng)為2時(shí)才計(jì)算和輸出。同樣,剛性墻能和阻尼能僅當(dāng)上面的卡片中RWEN和RYLEN分別設(shè)置為2時(shí)才會(huì)計(jì)算和輸出。剛性阻尼能集中到內(nèi)能里面。質(zhì)量阻尼能以單獨(dú)的行”system damping energy”出現(xiàn)。由于殼的體積粘性(bulk viscosity)而產(chǎn)生的能量耗散(energy dissipated)在版本970.4748之前是不計(jì)算的。在后續(xù)子版本中,設(shè)置TYPE=-2來(lái)在能量平衡中包含它。

最理想的情況下能量平衡:

總能量total energy = 初始總能量 + 外力功external work 。換句話說(shuō),如果能量比率energy ratio(指的是glstat中的total energy/initial energy,實(shí)際上是total energy/(initial energy + external work)) 等于1.0。注意,質(zhì)量縮放而增加質(zhì)量可能會(huì)導(dǎo)致能量比率增加。

注意在LSprepost的History>Global energies中不包含刪掉的單元(eroded elements)的能量貢獻(xiàn),然而GLSTAT文件中的能量包含了它們。注意它們的貢獻(xiàn)可以通過(guò)ASCII>glstat中的”Eroded Kinetic Energy”& “Eroded Internal Energy”來(lái)繪制。侵蝕能量(Eroded energy)是與刪掉的單元相關(guān)的內(nèi)能和刪掉的節(jié)點(diǎn)相關(guān)的動(dòng)能。典型來(lái)說(shuō),如果沒(méi)有單元?jiǎng)h掉”energy ratio w/o eroded energy”等于1,如果有單元被刪掉則小于1。刪掉的單元與”total energy/initial energy”比率沒(méi)有關(guān)系?？偰芰勘嚷试黾右?dú)w于其它原因,比如增加質(zhì)量。

重述一下,將一個(gè)單元?jiǎng)h掉時(shí),文件glstat中的內(nèi)能和動(dòng)能不會(huì)反映能量的丟失。取而代之的是能量的丟失記錄在glstat文件的”eroded internal energy” & “eroded kinetic energy”中。如果用內(nèi)能減去”eroded internal energy”將得到分析中還存在的單元的內(nèi)能。對(duì)動(dòng)能也一樣。

matsum文件中的內(nèi)能和動(dòng)能只包含余下(noneroded)的單元的貢獻(xiàn)。

注意,如果在*control_contact卡中將ENMASS設(shè)置為2,則與刪掉的單元的相關(guān)的節(jié)點(diǎn)不會(huì)刪掉,”eroded kinetic energy”是0。

在LSprepost中History>Global 只是動(dòng)能和內(nèi)能的簡(jiǎn)單相加,因此不包含接觸能和沙漏能等的貢獻(xiàn)。

----------

殼的負(fù)內(nèi)能:

為了克服這種不真實(shí)的效應(yīng)

--關(guān)掉考慮殼的減薄(ISTUPD in *control_shell)

--調(diào)用殼的體積粘性(set TYPE=-2 在*control_bulk_viscosity卡中)

--對(duì)在matsum文件中顯示為負(fù)的內(nèi)能的parts使用*damping_part_stiffness;

先試著用一個(gè)小的值,比如0.01

如果在*control_energy中設(shè)置RYLEN=2,因?yàn)閯傂宰枘岫軙?huì)計(jì)算且包含在內(nèi)能中。

----------

正的接觸能:

當(dāng)在接觸定義中考慮了摩擦?xí)r將得到正的接觸能。摩擦將導(dǎo)致正的接觸能。如果沒(méi)有設(shè)置接觸阻尼和接觸摩擦系數(shù),你將會(huì)看到凈接觸能為零或者一個(gè)很小的值(凈接觸能=從邊和主邊能量和)。所說(shuō)的小是根據(jù)判斷-在沒(méi)有接觸摩擦系數(shù)時(shí),接觸能為峰值內(nèi)能的10%內(nèi)可以被認(rèn)為是可接受的。

----------

負(fù)的接觸能:

突然增加的負(fù)接觸能可能是由于未檢測(cè)到的初始穿透造成的。在定義初始幾何時(shí)考慮殼的厚度偏置通常是最有效的減小負(fù)接觸能的步驟。查閱LS-DYNA理論手冊(cè)的23.8.3&23.8.4節(jié)可得到更多接觸能的信息。負(fù)接觸能有時(shí)候因?yàn)閜arts之間的相對(duì)滑動(dòng)而產(chǎn)生。這跟摩擦沒(méi)有關(guān)系,這里說(shuō)的負(fù)接觸能從法向接觸力和法向穿透產(chǎn)生。當(dāng)一個(gè)穿透的節(jié)點(diǎn)從它原來(lái)的主面滑動(dòng)到臨近的沒(méi)有連接的主面時(shí),如果穿透突然檢測(cè)到,則產(chǎn)生負(fù)的接觸能。

如果內(nèi)能為負(fù)接觸能的鏡像,例如glstat文件中內(nèi)能曲線梯度與負(fù)接觸能曲線梯度值相等,問(wèn)題可能是非常局部化的,對(duì)整體求解正確性沖擊較小。你可以在LS-prepost中分離出有問(wèn)題的區(qū)域,通繪制殼單元部件內(nèi)能云圖(Fcomp > Misc > Internal energy)。實(shí)際上,顯示的是內(nèi)能密度,比如內(nèi)能/體積。內(nèi)能密度云圖中的熱點(diǎn)通常表示著負(fù)的接觸能集中于那里。

如果有多于一個(gè)的接觸定義,sleout文件(*database_sleout)將報(bào)告每一個(gè)接觸對(duì)的接觸能量,因此縮小了研究負(fù)接觸能集中處的范圍。

克服負(fù)接觸能的一般的建議如下:

-消除初始穿透(initial penetration)。(在message文件中查找”warning”)

-檢查和排除冗余的接觸條件。不應(yīng)該在相同的兩個(gè)parts之間定義多于一個(gè)的接觸。

-減小時(shí)間步縮放系數(shù)

-設(shè)置接觸控制參數(shù)到缺省值,SOFT=1 & IGNORE=1除外(接觸定義選項(xiàng)卡C)

-對(duì)帶有尖的邊的接觸面,設(shè)置SOFT=2(僅用于segment-to-segment接觸)。而且,在版本970中推薦設(shè)置SBOPT(之前的EDGE)為4對(duì)于部件之間有相對(duì)滑移的SOFT=2的接觸。為了改進(jìn)edge-to-edge SOFT=2接觸行為,設(shè)置DEPTH=5。請(qǐng)注意SOFT=2接觸增加了額外的計(jì)算開(kāi)消,尤其是當(dāng)SBOPT或者DEPTH不是缺省值時(shí),因此應(yīng)該僅在其它接觸選項(xiàng)(SOFT=0或者SOFT=1)不能解決問(wèn)題時(shí)。

模型的細(xì)節(jié)可能會(huì)指示可用其它的一些方法

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