AUTODYN基礎(chǔ)培訓(xùn)二課件

2024/4/24ANSYSAUTODYN

基礎(chǔ)培訓(xùn)二2024/4/24基礎(chǔ)培訓(xùn)二1、Lagrange求解器2、模型生成3、Lagrange之間的作用

2024/4/24AUTODYN求解器類(lèi)型LagrangeEulerALESPHShellBeam2024/4/24拉格朗日求解器節(jié)點(diǎn)與材料一起運(yùn)動(dòng)通常用于固體求解T=0.0T>0.02024/4/24在初始計(jì)算中，通過(guò)對(duì)節(jié)點(diǎn)速度積分得到新節(jié)點(diǎn)位置；通過(guò)新節(jié)點(diǎn)的坐標(biāo)位置，計(jì)算新單元密度和單元應(yīng)變率。Dt拉格朗日求解器計(jì)算周期2024/4/24通過(guò)應(yīng)變，新的應(yīng)力可以計(jì)算得到；對(duì)于各向同性材料，單元變形能分成兩個(gè)獨(dú)立的部分。體積改變形狀改變拉格朗日求解器計(jì)算周期2024/4/24拉格朗日求解器體積改變是由球應(yīng)力張量引起；形狀改變是由于偏應(yīng)力張量引起；球應(yīng)力由狀態(tài)方程(EOS)求解；偏應(yīng)力是由胡克定律(彈性)和塑性屈服準(zhǔn)則(本構(gòu)關(guān)系)求解-即AUTODYN的強(qiáng)度模型；再加上失效模型就構(gòu)成了一個(gè)完整的材料模型。計(jì)算周期2024/4/24一旦計(jì)算單元的內(nèi)部應(yīng)變率被決定，每一個(gè)節(jié)點(diǎn)的受力能被計(jì)算；在這個(gè)階段，任何從邊界條件或者與其它物體作用的外力作為最后的節(jié)點(diǎn)力；知道節(jié)點(diǎn)力，可以計(jì)算節(jié)點(diǎn)加速度；積分節(jié)點(diǎn)加速度可以得到節(jié)點(diǎn)新的速度；完整的一個(gè)計(jì)算周期，重復(fù)這樣的周期達(dá)到要求的求解時(shí)間或循環(huán)次數(shù)。拉格朗日求解器計(jì)算周期2024/4/24質(zhì)量守恒能量守恒狀態(tài)方程動(dòng)量守恒積分積分S

力S

質(zhì)量初始條件節(jié)點(diǎn)位移應(yīng)變率，密度壓力，內(nèi)能偏應(yīng)力節(jié)點(diǎn)力節(jié)點(diǎn)加速度節(jié)點(diǎn)速度外力(邊界條件和相互作用力)拉格朗日求解器計(jì)算周期2024/4/24拉格朗日求解器低速擠壓碰撞2024/4/24拉格朗日求解器結(jié)構(gòu)求解器(V5和更早版本)；2D體單元和殼單元3D體單元、殼單元和桿單元非結(jié)構(gòu)求解器(V6和以后版本)；3D體單元、殼單元和桿單元

結(jié)構(gòu)和非結(jié)構(gòu)求解器能用在同一個(gè)模型中；非結(jié)構(gòu)化求解器特點(diǎn)：求解更快內(nèi)存更小適合更多的單元類(lèi)型可以方便地輸入非結(jié)構(gòu)網(wǎng)格結(jié)構(gòu)化網(wǎng)格能被轉(zhuǎn)換成非結(jié)構(gòu)化網(wǎng)格。結(jié)構(gòu)和非結(jié)構(gòu)求解器2024/4/24拉格朗日體單元求解器非結(jié)構(gòu)求解器單元和面連通性數(shù)據(jù)存儲(chǔ)滿(mǎn)足精確要求結(jié)構(gòu)求解器基于IJK空間單元和面通過(guò)IJK空間定義(與數(shù)據(jù)存儲(chǔ)一樣)數(shù)據(jù)存儲(chǔ)直接與單元拓?fù)湎噙B2024/4/24拉格朗日體單元求解器8節(jié)點(diǎn)精確體積積分；可選的8節(jié)點(diǎn)，近似單點(diǎn)高斯體積積分適合非結(jié)構(gòu)求解器；FE公式(Hallquist)對(duì)于彎曲變形，單元精度會(huì)有一些損失使用人工粘性處理沖擊問(wèn)題；使用沙漏控制保證穩(wěn)定性；五面體單元–退化六面體單元；非結(jié)構(gòu)六面體求解速度約為結(jié)構(gòu)求解的兩倍；非結(jié)構(gòu)求解需要的內(nèi)存比結(jié)構(gòu)求解少30%。六面體/五面體求解器2024/4/24拉格朗日體單元求解器比較：內(nèi)存(Mb)時(shí)間(Min)加速比結(jié)構(gòu)4001.741.00非結(jié)構(gòu)(Exact)3400.881.97非結(jié)構(gòu)(Gauss)3400.622.803.2GhzPentiumIVProcessor,WindowsXP500,000單元效率內(nèi)存速度六面體/五面體求解器2024/4/24拉格朗日體單元求解器六面體/五面體求解器結(jié)構(gòu)六面體(IJK) 非結(jié)構(gòu)六面體(Exact)2024/4/24拉格朗日體單元求解器六面體/五面體求解器Structured(IJK)UnstructuredExactUnstructuredGauss2024/4/24拉格朗日體單元求解器頭盔碰撞2024/4/24拉格朗日體單元求解器

結(jié)構(gòu)非結(jié)構(gòu)結(jié)構(gòu)非結(jié)構(gòu)頭盔碰撞運(yùn)行時(shí)間結(jié)構(gòu)求解器：443分鐘(36284cycles)非結(jié)構(gòu)求解器：265分鐘(37013cycles)

2024/4/24拉格朗日體單元求解器非結(jié)構(gòu)四面體求解器SCPTet

ANPTetstandardconstantpressureAverageNodalPressure2024/4/24拉格朗日體單元求解器非結(jié)構(gòu)四面體求解器SCPTet

ANPTet2024/4/24拉格朗日殼單元求解器2D采用有限差分方法；3D采用有限體積方法；節(jié)點(diǎn)軌跡變形；厚度是一個(gè)填充參數(shù)；假設(shè)厚度方向?yàn)榱銘?yīng)變；時(shí)間步長(zhǎng)由長(zhǎng)/寬決定，不由厚度決定。2024/4/24拉格朗日殼單元求解器非結(jié)構(gòu)求解比結(jié)構(gòu)求解器更有效在相同的內(nèi)存條件下，非結(jié)構(gòu)求解器存放的數(shù)據(jù)多于結(jié)構(gòu)求解器四邊形單元內(nèi)存(Mbytes)時(shí)間(分鐘)加速比結(jié)構(gòu)5903.131.00非結(jié)構(gòu)5951.162.703.2GhzPentiumIVProcessor,WindowsXP2024/4/24拉格朗日殼單元求解器四邊形殼單元結(jié)構(gòu)VS非結(jié)構(gòu)2024/4/24拉格朗日殼單元求解器結(jié)構(gòu)非結(jié)構(gòu)金屬成型運(yùn)行時(shí)間(65,000cycles)結(jié)構(gòu)網(wǎng)格(IJK)3,814分鐘非結(jié)構(gòu)網(wǎng)格1,722分鐘2024/4/24拉格朗日殼單元求解器非結(jié)構(gòu)主四邊形殼單元金屬?zèng)_壓成型2024/4/24拉格朗日殼單元求解器四邊形/三角形殼單元非結(jié)構(gòu)(上面)、結(jié)構(gòu)四邊形(中間)、結(jié)構(gòu)三角形(底部)位移時(shí)間2024/4/24拉格朗日殼單元求解器用戶(hù)定義層數(shù)；每一層可以是各向同性或異性材料；每一層指定一個(gè)厚度。多層復(fù)合殼單元外層內(nèi)層2024/4/24拉格朗日梁?jiǎn)卧蠼馄?節(jié)點(diǎn)的梁?jiǎn)卧?；允許大的軸向應(yīng)變。內(nèi)存(Mbytes)時(shí)間(分鐘)加速比IJK梁1500.441.00非結(jié)構(gòu)梁1650.251.763.2GhzPentiumIVProcessor,WindowsXP結(jié)構(gòu)/非結(jié)構(gòu)梁?jiǎn)卧蠼馄餍?100,000梁?jiǎn)卧?例子：2024/4/24拉格朗日梁?jiǎn)卧蠼馄髁簡(jiǎn)卧孛嫘螤?024/4/24拉格朗日梁求解器彈簧單元2024/4/24拉格朗日梁求解器彈簧單元速度與時(shí)間曲線2024/4/24拉格朗日梁求解器阻尼單元2024/4/24拉格朗日梁求解器彈簧&阻尼單元聯(lián)合速度與時(shí)間曲線2024/4/24拉格朗日梁?jiǎn)卧獑胃?/p>

多根梁2D梁構(gòu)架3D梁構(gòu)架梁求解器集合組名(類(lèi)型，梁號(hào)[單元/梁])

2024/4/24拉格朗日殼&梁求解器沖擊波對(duì)結(jié)構(gòu)中的影響2024/4/24拉格朗日求解器(3D)兼容性：模型非結(jié)構(gòu)和結(jié)構(gòu)(IJK)求解器能在一個(gè)模型中使用；接觸在任意非結(jié)構(gòu)/結(jié)構(gòu)(IJK)之間使用；連接在任意非結(jié)構(gòu)/結(jié)構(gòu)(IJK)之間使用；與3DFCT和多物質(zhì)歐拉耦合包括厚度殼單元；并行可用于結(jié)構(gòu)求解器和非結(jié)構(gòu)求解器。2024/4/24拉格朗日求解器優(yōu)勢(shì)：更少的循環(huán)計(jì)算時(shí)間(與Euler/SPH比較)；清晰的材料界面和邊界定義；時(shí)間歷史信息完整；強(qiáng)度模型好；較少的沖擊擴(kuò)散；代碼相對(duì)簡(jiǎn)單。殼單元優(yōu)勢(shì)：時(shí)間步長(zhǎng)由單元面積決定，不是單元厚度。梁?jiǎn)卧獌?yōu)勢(shì)：

時(shí)間步長(zhǎng)由單元長(zhǎng)度決定，不是橫截面積。不足：?jiǎn)卧冃螌?dǎo)致時(shí)間步長(zhǎng)變小；單元變形導(dǎo)致網(wǎng)格混亂；特殊的需要碰撞/滑移界面。2024/4/24基礎(chǔ)培訓(xùn)二1、Lagrange求解器2、模型生成3、Lagrange之間的作用2024/4/24模型生成結(jié)構(gòu)化網(wǎng)格用(I,J,K)指標(biāo)空間；網(wǎng)格在AUTODYN里生成；通過(guò)ICEM-CFD或TrueGrid輸入結(jié)構(gòu)化網(wǎng)格；完全的Fill使用選項(xiàng)(用材料/初始條件)；邊界條件用指標(biāo)空間，交互式或者組施加。非結(jié)構(gòu)化網(wǎng)格使用單元/節(jié)點(diǎn)連通性；目前，不能在AUTODYN中生成網(wǎng)格；通過(guò)LS-DYNA.k文件和Nastran.dat文件輸入非結(jié)構(gòu)網(wǎng)格；轉(zhuǎn)換結(jié)構(gòu)化網(wǎng)格為非結(jié)構(gòu)化網(wǎng)格；有限的使用Fill選項(xiàng)；邊界條件用交互式或者組施加。2024/4/24結(jié)構(gòu)求解器對(duì)于每一個(gè)Part，一個(gè)結(jié)構(gòu)化網(wǎng)格通過(guò)單元、節(jié)點(diǎn)和空間位置定義連通性；網(wǎng)格可以看成兩中空間：指標(biāo)空間（I、J、K）物理空間（X、Y、Z）指標(biāo)空間在節(jié)點(diǎn)之間定義連通性；物理空間定義節(jié)點(diǎn)的空間坐標(biāo)；物理空間通常用來(lái)識(shí)別；指標(biāo)空間用來(lái)建立模型。2024/4/24結(jié)構(gòu)求解器指標(biāo)空間在AUTODYN里，每一個(gè)結(jié)構(gòu)Part定義一個(gè)指標(biāo)空間(i,j,k)，i,j,k是從1到N的的整數(shù)指標(biāo)空間是矩形(2D)或立方體(3D)2024/4/24結(jié)構(gòu)求解器物理空間在AUTODYN里，每一個(gè)結(jié)構(gòu)Part也定義一個(gè)xyz空間，x,y,z是實(shí)數(shù)在物理空間里，通常網(wǎng)格有一個(gè)形狀i=12024/4/24結(jié)構(gòu)求解器不用區(qū)域不是所有在指標(biāo)空間定義的單元都需要定義物理空間-這就可以實(shí)現(xiàn)復(fù)雜幾何模型劃分網(wǎng)格不用的單元連接的節(jié)點(diǎn)指標(biāo)空間物理空間2024/4/24結(jié)構(gòu)求解器Part向?qū)?predef)允許自動(dòng)劃分預(yù)先確定幾何的結(jié)構(gòu)化網(wǎng)格：2DVolume:Box，Quad，Circle，尖頂拱形，菱形，三角形，楔形3DVolume:Box,Hex,Cylinder,Sphere,尖頂拱形,Fragments/BricksShells:Plate,CylinderBeams:Single,Multiple,2DFrame,3DFrame網(wǎng)格轉(zhuǎn)換AUTODYN里創(chuàng)建結(jié)構(gòu)化網(wǎng)格輸入結(jié)構(gòu)化網(wǎng)格：ICEM-CFDTrueGrid

手工生成結(jié)構(gòu)網(wǎng)格：點(diǎn)邊面體2024/4/24結(jié)構(gòu)求解器生成Part向?qū)壳x幾何定義網(wǎng)格填充材料2024/4/24結(jié)構(gòu)求解器Part向?qū)?Box(2D)2024/4/24結(jié)構(gòu)求解器Part向?qū)?Circle(2D)2024/4/24結(jié)構(gòu)求解器Part向?qū)?2D楔形尖頂拱形四邊形菱形三角形2024/4/24結(jié)構(gòu)求解器Part向?qū)?Box(3D)2024/4/24結(jié)構(gòu)求解器Part向?qū)?Cylinder(3D)2024/4/24結(jié)構(gòu)求解器Part向?qū)?3D球六面體尖頂拱形磚塊結(jié)構(gòu)2024/4/24結(jié)構(gòu)求解器一步步的生成網(wǎng)格：點(diǎn)線面體拉伸手工生成網(wǎng)格2024/4/24結(jié)構(gòu)求解器生成節(jié)點(diǎn)：節(jié)點(diǎn)的I、J和K的值；節(jié)點(diǎn)坐標(biāo)。手工生成網(wǎng)格2024/4/24結(jié)構(gòu)求解器生成線條：I/J/K方向；幾何比例；直線或圓?。粓A弧中心。手工生成網(wǎng)格2024/4/24結(jié)構(gòu)求解器生成平面：I/J/K常數(shù)；網(wǎng)格類(lèi)型。手工生成網(wǎng)格2024/4/24結(jié)構(gòu)求解器手工生成網(wǎng)格生成體網(wǎng)格：生成方向；網(wǎng)格比例。2024/4/24結(jié)構(gòu)求解器網(wǎng)格過(guò)渡2024/4/24結(jié)構(gòu)求解器Part向?qū)?3D圓柱面平面2024/4/24結(jié)構(gòu)求解器TrueGrid(XYZScientific)在TrueGrid生成網(wǎng)格輸出AUTODYN-3D網(wǎng)格數(shù)據(jù)(.zon)用“Import,TrueGrid”下拉菜單分別讀入Part能夠定義塊和殼單元自動(dòng)創(chuàng)建和填充多物質(zhì)Part輸入的體單元能夠用于拉格朗日、ALE和Fill求解器輸入結(jié)構(gòu)網(wǎng)格2024/4/24結(jié)構(gòu)求解器ANSYS,ICEM-CFD超強(qiáng)的六面體網(wǎng)格生成工具直接連接CADCATIA，Pro/E，SDRCI-DEAS，SolidWorks，UG，……提供輸入結(jié)構(gòu)網(wǎng)格到AUTODYNICEMmultiblockmeshes(.geo文件)輸入過(guò)程輸入結(jié)構(gòu)網(wǎng)格2024/4/24拉格朗日求解器除了初始條件，有的問(wèn)題還需要邊界條件。AUTODYN中邊界條件包括：應(yīng)力速度力傳輸流動(dòng)(Euler)通過(guò)定義邊界條件，然后通過(guò)指標(biāo)空間施加，對(duì)于非結(jié)構(gòu)網(wǎng)格通過(guò)Groups來(lái)施加。邊界條件2024/4/24拉格朗日求解器對(duì)于拉格朗日Part邊界上的單元，缺省的邊界條件是自由邊界(壓強(qiáng)=0.0)；不用的單元也定義邊界條件。拉格朗日邊界條件不用單元填充單元缺省邊界條件2024/4/24

常數(shù)t梯形t

三角形t

指數(shù)Pk=Pe-ktt用戶(hù)子程序EXSTRt拉格朗日Part拉格朗日求解器應(yīng)力邊界條件(Lagrange/ALE/Shells)2024/4/24拉格朗日求解器X-Velocity常數(shù)；Y-Velocity 分段表示；Z-Velocity限制在最大和最小之間；通過(guò)子程序EXVEL。R-Velocity 常數(shù)、指定中心或X/Y/Z軸GeneralVelocity 在X、Y、Z(3D)上固定常數(shù)和固定繞坐標(biāo)軸的角速度速度邊界條件(Lagrange/ALE/Shells/Beams/SPH)2024/4/24拉格朗日求解器節(jié)點(diǎn)力邊界條件：常數(shù)、分段、子程序EXFOR定義；X、Y、Z和通常的方向。單元力/單位長(zhǎng)度邊界條件：常數(shù)、分段、子程序EXFOR定義；X、Y、Z和通常的方向。力邊界條件(Beams)2024/4/24拉格朗日求解器傳輸波通過(guò)單元面；僅僅為正交分量傳輸；可以指定邊界的阻抗：如果設(shè)置阻抗為零，使用臨近單元的阻抗；如果邊界條件是近似的，不要在關(guān)心的區(qū)域設(shè)置。傳輸邊界條件(Lagrange/ALE/Euler)2024/4/24時(shí)間步長(zhǎng)控制時(shí)間步長(zhǎng)的穩(wěn)定性0<?T<∞ 0<?T<?X/C

R=C.?T隱示顯示拉格朗日歐拉其中Δx=單元尺寸，c=聲速，v=速度2024/4/24時(shí)間步長(zhǎng)控制缺省通常OK；如果初始步長(zhǎng)小于零，時(shí)間步長(zhǎng)取穩(wěn)定時(shí)間步長(zhǎng)的1/2；如果最小時(shí)間步長(zhǎng)小于零，時(shí)間步長(zhǎng)取出是時(shí)間步長(zhǎng)的1/10；對(duì)于拉格朗日計(jì)算問(wèn)題，安全系數(shù)可以增大到0.9。時(shí)間步長(zhǎng)選項(xiàng)2024/4/24阻尼選項(xiàng)人工粘性沙漏阻尼靜態(tài)阻尼2024/4/24阻尼選項(xiàng)不連續(xù)(不穩(wěn)定)二次粘性二次粘性+線性粘性人工粘性的作用2024/4/24阻尼選項(xiàng)“沙漏”是一種網(wǎng)格畸變，它引起沒(méi)有應(yīng)變或體積沒(méi)有改變；四邊形(2D)/六邊形(3D)使用單點(diǎn)積分產(chǎn)生沙漏；一個(gè)阻尼選項(xiàng)用來(lái)抵制這種變形。沙漏阻尼2024/4/24阻尼選項(xiàng)靜態(tài)阻尼系數(shù)Rd：靜態(tài)阻尼

那么如果2024/4/24非結(jié)構(gòu)求解器非結(jié)構(gòu)化網(wǎng)格的單元/節(jié)點(diǎn)四面體多數(shù)四邊形多數(shù)六面體2024/4/24非結(jié)構(gòu)求解器輸入選項(xiàng)LS-DYNA(.k)MSC.Nastran(.dat)通過(guò)第三方軟件輸入：網(wǎng)格材料邊界條件載荷初始條件目前不能產(chǎn)生非結(jié)構(gòu)化網(wǎng)格可以將結(jié)構(gòu)化Parts(IJK)轉(zhuǎn)變成非結(jié)構(gòu)網(wǎng)格模型輸入2024/4/24非結(jié)構(gòu)求解器使用‘

Import’下拉菜單輸入網(wǎng)格/模型：ICEM-CFD(ANSYSAI*Environment)TrueGridLS-DYNAMSCNastran選擇‘ConvertIJKPartstoUnstructured

’菜單轉(zhuǎn)換拉格朗日、殼和桿結(jié)構(gòu)化Parts成非結(jié)構(gòu)化Parts。網(wǎng)格模型輸入2024/4/24非結(jié)構(gòu)求解器RetainLS-DYNApartdefinitions選擇該項(xiàng)，LS-DYNAparts/Nastran屬性將要被轉(zhuǎn)變成AUTODYNParts不選擇，AUTODYNParts通過(guò)連通性創(chuàng)建(如：六面體單元與四面體單元的連接)Mergeduplicatematerials選擇該項(xiàng)，使用相同模型和數(shù)據(jù)的材料將被合并Checkshellorientation選擇該項(xiàng)，檢查所有殼單元外法向方向的一致性(對(duì)流固耦合重要)

網(wǎng)格模型輸入輸入選項(xiàng)-LS-DYNA/MSC-Nastran2024/4/24非結(jié)構(gòu)求解器輸入后操作材料改變材料名字(選擇)添加材料數(shù)據(jù)(大多數(shù)材料)可以增加剛性材料約束(如果需要的話(huà))Parts、組件、邊界條件和初始條件重新命名(選擇)網(wǎng)格模型輸入2024/4/24非結(jié)構(gòu)求解器整個(gè)Parts進(jìn)行的操作復(fù)制填充刪除分區(qū)，平移非結(jié)構(gòu)Parts操作測(cè)量單元/節(jié)點(diǎn)2024/4/24非結(jié)構(gòu)求解器邊界條件的施加在屏幕創(chuàng)建多變形來(lái)選擇節(jié)點(diǎn)，通過(guò)此來(lái)施加或刪除邊界條件非結(jié)構(gòu)/結(jié)構(gòu)Parts操作2024/4/24非結(jié)構(gòu)求解器邊界條件的施加：通過(guò)屏幕選擇節(jié)點(diǎn)來(lái)施加和刪除邊界條件。非結(jié)構(gòu)/結(jié)構(gòu)Parts操作2024/4/24非結(jié)構(gòu)求解器組的選擇操作用戶(hù)定義節(jié)點(diǎn)、單元或面集合交互式選擇通過(guò)多邊形選取通過(guò)面選取施加/清除邊界條件節(jié)點(diǎn)邊界條件面邊界條件填充單元組材料速度初始條件

非結(jié)構(gòu)/結(jié)構(gòu)Parts操作2024/4/24非結(jié)構(gòu)求解器求解選項(xiàng)六面體積分有限體積公式：精確積分(缺省)有限元公式：近似單點(diǎn)的高斯積分(快，精度不高)六面體沙漏控制AD標(biāo)準(zhǔn)(缺省)求解速度快，需要內(nèi)存少；Flanagan-Belytschko(用于大剛體轉(zhuǎn)動(dòng))。四面體壓強(qiáng)積分SCP(六面體網(wǎng)格為主)ANP(僅僅四面體網(wǎng)格)2024/4/24非結(jié)構(gòu)求解器_練習(xí)拉格朗日(3D非結(jié)構(gòu))頭盔碰撞輸入LS-DYNA網(wǎng)格：頭盔頭空心圓柱頭盔用泡沫材料模型2024/4/24基礎(chǔ)培訓(xùn)二1、Lagrange求解器2、模型生成3、Lagrange之間的作用

2024/4/24拉格朗日-拉格朗日連接2024/4/24拉格朗日-拉格朗日連接剛性連接可用于殼單元和梁?jiǎn)卧?具有旋轉(zhuǎn)和平動(dòng)自由度)UYA

URA

UXA

UYB

URB

UXB

銷(xiāo)連接UXA=UXBUYA

=UYBURA

≠URB剛性連接UXA=UXBUYA

=UYBURA

=URB用來(lái)指定Part與Part之間的連接方式缺省方式為剛性連接連接處節(jié)點(diǎn)一致

2024/4/24拉格朗日-拉格朗日連接例子梁?jiǎn)卧蜌卧B接2024/4/24拉格朗日-拉格朗日接觸用于材料的滑移或分離：內(nèi)部破碎;大剪切變形;碰撞面接觸;靶板侵徹;自身接觸.接觸作用2024/4/24拉格朗日-拉格朗日接觸間隙的定義每一個(gè)面段被一個(gè)接觸搜索區(qū)域包圍接觸搜索區(qū)域的半徑就叫做間隙尺寸接觸邏輯拉格朗日

殼單元

2024/4/24拉格朗日-拉格朗日接觸任何進(jìn)入一個(gè)面斷的接觸搜索區(qū)域節(jié)點(diǎn)將受到一個(gè)反作用力，這個(gè)力的大小與該點(diǎn)進(jìn)入接觸搜索區(qū)域的深度成比例。3D中的接觸搜索區(qū)域FR

R

R

R

FF2024/4/24拉格朗日-拉格朗日接觸3D中的接觸搜索區(qū)域(分解圖)2024/4/24拉格朗日-拉格朗日接觸可見(jiàn)間隙(2D&3D)作用面必須用間隙尺寸分開(kāi)可見(jiàn)/內(nèi)部間隙內(nèi)部間隙(僅僅2D–不推薦)作用面可能處于接觸中

2024/4/24拉格朗日-拉格朗日接觸外部間隙(2D&3D)最大間隙尺寸為所有作用面中最小邊的一半內(nèi)部間隙(僅僅2D–不推薦)最大間隙尺寸為所有作用面中最小單元尺寸的一半最大間隙尺寸2024/4/24拉格朗日-拉格朗日接觸如果只用點(diǎn)/面接觸，有可能會(huì)接觸失效；這是就需要edge/edge接觸；在Lag/Lag作用菜單上選擇；增加計(jì)算時(shí)間。Edge/Edge接觸2024/4/24拉格朗日-拉格朗日接觸Edge/Edge接觸梁

殼

拉格朗日

WhatareEdges?2024/4/24拉格朗日-拉格朗日接觸這個(gè)選項(xiàng)用來(lái)定義自身接觸。相應(yīng)會(huì)增加計(jì)算時(shí)間。自身接觸當(dāng)求解殼單元結(jié)構(gòu)的時(shí)候，推薦設(shè)置間隙尺寸為殼單元厚度的一半比較好。殼單元2024/4/24侵蝕通過(guò)刪除退化單元和嚴(yán)重變形單元，避免拉格朗日網(wǎng)格的重新分區(qū)；侵蝕單元的內(nèi)能被丟棄；侵蝕單元的質(zhì)量和動(dòng)量可以被保留，也可以丟棄；所有的材料模型都能與失效模型結(jié)合求解；侵蝕是一個(gè)非物理過(guò)程。

2024/4/24侵蝕帶侵蝕作用Surfacefaces表面節(jié)點(diǎn)PartAPartB2024/4/24侵蝕Surfacefaces自由節(jié)點(diǎn)ErodedErodedPartAPartB表面節(jié)點(diǎn)帶侵蝕作用2024/4/24侵蝕