abaqus主表面從表面的案例
ABAQUS粗糙表面模型生成插件
ABAQUS粗糙表面模型生成插件
abaqus根據規范更改表面粗糙度 ￥50
python代碼：依據FFT變換生成不同等級的路面粗糙度
Abaqus前處理插件-裝配體外表面提取 ￥100
在許多分析任務的前處理過程中，我們需要提取結構的外表面，例如施加壓強等分布載荷，添加對流、輻射、溫度邊界條件等等。
對于較為復雜的模型提取外表面的工作十分繁瑣。尤其當模型中存在多個instance實例，其結構表面之間往往相互交叉接觸，選取模型外表面就變得十分繁瑣。
本插件可以實現“一鍵提取”所有外表面。
插件簡易界面通過RSG功能制作，內核函數包含對單元面的遍歷識別和處理，可作為前處理二次開發進階技巧供大家學習。
ABAQUS隨機粗糙度表面地形建模
本案例介紹在ABAQUS內建立三維隨機粗糙度表面或地形圖模型，并通過隨機粗糙度表面進行簡單的動力學模擬。
首先采用CAD隨機粗糙度表面插件建立三維隨機粗糙度實體幾何模型，并將模型導出為iges格式文件。
在ABAQUS內將隨機粗糙度表面文件以部件的形式進行導入。
為了動力學模擬的需要，這里新建一個球體部件，并將其與粗糙度表面進行裝配，球體置于粗糙度表面的任意位置。
設置球體與粗糙度表面間的相互作用，切向行為設置罰，法向行為設置硬接觸，并在載荷中設置重力并將模型下表面固定。
為模型劃分網格，單元形狀設置為四面體。
提交作業并查看球體在隨機粗糙度表面或特定地形中的運動路徑情況。
展開 
基于abaqus的鋼球振動沖擊金屬表面
基于abaqus的鋼球振動沖擊金屬表面。鋼球在X方向做勻速平動，Z軸做順時針轉動，Y方向做簡諧振動，沖擊金屬表面?？稍诮Y果中觀察金屬塊的應力應變、溫度分布、表面形貌，晶粒細化，殘余應力等情況。
abaqus粗糙表面的微動磨損分析
本文基于Abaqus分析了粗糙表面的微動磨損行為。
進行粗糙表面的微動磨損分析，首先需要建立粗糙表面的幾何模型。試驗表明分形理論可以有效表征粗糙面的幾何特征。二維表面的輪廓由W-M分形函數確定
通過python結合式(1)可以得到模型輪廓如下。
圖 1 Python生成的輪廓
圖 2 粗糙面網格
磨損模型如下
通過umeshmotion子程序將式(2)磨損模型引入有限元分析。
壓頭上，法向施加固定載荷，切向施加周期性位移。計算得到的結果如下所示。
圖 3 光滑表面和粗糙表面磨損后的變形對比
展開 Abaqus子程序系列：FRIC（定義接觸表面的摩擦行為）
如果定義的是接觸表面，傳遞進來的是0
CINAME
與摩擦定義關聯的用戶定義的表面的名稱，左對齊。對于接觸單元，它是與摩擦定義相關聯的界面定義所給定的單元集名稱。如果給界面定義分配了一個可選名稱，那么CINAME將作為這個名稱傳入，左對齊。
SLNAME
從表面名稱。如果使用接觸單元，傳入空格。
NSNAME
主表面名稱。如果使用接觸單元，傳入空格。
NPT
接觸單元的積分點數量。如果接觸表面被定義，傳入0。
NODE
與這個接觸點相關的用戶定義的全局從節點號（或根據部件實例的組裝定義的模型的內部節點編號)。如果使用面對面接觸公式，則對應于約束的主要從節點。如果從接觸單元調用，則作為零傳遞。
COORDS(MCRD)
包含這個點，此時坐標的數組
RCOORD(MCRD)
如果主表面被定義為剛性表面，則傳遞此數組，其中包含剛性表面上對應點的當前位置和方向的坐標。
DROT(2,2)
旋轉增量矩陣。對于與三維剛性表面的接觸，這個矩陣表示相對于剛性表面的表面方向的增量旋轉。這樣做是為了在這個子程序中，將向量值或張力值狀態變量進行正確的旋轉。在調用FRIC之前，應力和滑動分量已經旋轉了這個量。該矩陣作為單位矩陣傳入二維和軸對稱接觸問題。
TEMP（2）
從節點和對應的主表面，此時的溫度
PREDEF（2，NPRED）
一個數組，包含當前增量步結束時的所有用戶指定場變量的一對值(分析開始時的初始值和分析期間的當前值)。從接觸對調用FRIC，一對值中的第一個值對應于從節點，第二個值對應于主表面上最近的點。如果從一個大滑動接觸單元調用FRIC，則PREDEF(1,NPRED)對應于該單元積分點處的值，PFREDEF(2,NPRED)對應于對應表面上的最近點。
展開 Abaqus中多相材料不規則表面輪廓提取
1 前言
在Abaqus切削仿真中，目前多采用二維正交模型來轉化和代替各種加工形式。目前對于Abaqus切削仿真的可查資料中，多是模型建立和一些注意事項，對于其后處理過程較少提及。加工表面的粗糙度是表面質量評價的一項重要指標，仿真得到的微觀結構的細觀變化也是切削仿真的一大研究重點。因此，對切削表面的輪廓提取是有必要的。
在Abaqus中，在后處理過程中，輪廓提取可以采用多種方式，例如建立路徑和導出連續節點坐標。但是多相材料建模通常采用不同的Part最后Assembly得到，建立路徑只能在獨立的Part中進行。除此之外，網格劃分與一般的規則形狀得到的網格也不相同，多相材料劃分得到的網格往往并不規則，因此導出連續節點也是不現實的。因此，要想導出多相復合材料的表面輪廓需要尋求一種別的方式。
本篇小節只要針對多相復合材料的切削表面輪廓進行講述，所使用的軟件包括Abaqus、AutoCAD、Excel和Origin。除了Abaqus切削仿真表面輪廓提取，也可對一些其它復雜形狀和結構的輪廓提取做出指導意義。
2 提取過程
如圖1所示為所選擇的一個案例的最終切削結果，其中顆粒和基體是單獨的Part，切削表面并不平整。提取目標既最上面的一條輪廓曲線。要求：其中某一條線的實際長度。
圖1 樣件的切削結果
第一步：在Abaqus Visualization界面下選擇命令Plot Deformed Shape，再點擊命令Render Model: Wireframe，得到如圖2所示的線框圖。
展開 Abaqus 如何處理接觸表面的初始過盈問題 ￥5
目錄
初始過盈產生的原因
解決初始過盈的方法
在 Abaqus 中指定不同選項
壓力工況下的性能驗證
總結
在本文中，我們將以軸對稱 O 型圈為例，闡述并展示 Abaqus 處理接觸表面初始重疊（即初始過盈）的多種方法。
ABAQUS任意單元表面加入膜單元或加入復合材料纖維層
以上內容來自360百科
本期是教大家如何在ABAQUS有限元模型中在任意實體單元表面加入殼單元作為纖維增強材料來模擬復合材料：
孔眼壁上的膜單元來模擬壁面加固材料
內加入纖維增強材料
轉自公眾號——ABAQUS大世界
旨在分享，若侵即刪.
abaqusO形環接觸表面初始過盈的幾種處理方法
O形環接觸表面初始過盈的幾種處理方法
本文以O形環為例，介紹Abaqus接觸表面初始過盈的幾種處理方法
初始過盈原因
接觸定義中，接觸面可能在分析開始時就有部分干涉，即初始過盈。其由不同原因引起，CAD模型較差導致零件之間的重疊，或離散誤差導致，如圖中2兩個圓，同樣半徑的接觸幾何圖形被不同單元尺寸離散，它們以不同的方式偏離原來的圓形邊。
初始過盈還有一個原因是干涉。在實際中，干涉部分在初始位置時就會變形，產生應力和應變。
初始過盈解決方法
Abaqus中，主要有2種方法解決初始過盈問題：無應變調整、設置過盈配合。
無應變調整
無應變調整即移動從面上的節點，以消除初始過盈，不會產生任何應力或應變。這種方法適用于非預期的過盈，模型在第一個分析步0位置時移動節點，改變了網格，如果干涉部分過大，會產生網格扭曲。
設置過盈配合
設置過盈配合后，在從面節點移動到非干涉位置的過程中，有對應的應力和應變產生。這個方法適用于預期過盈。在第1個分析步0位置時，接觸從面的節點不進行調整，在第1個增量步中開始進行過盈求解。
ABAQUS設置
上述兩種方法在不同接觸屬性中的處理方式不同。
通用接觸
一般情況下，通用接觸中使用無應變調整處理小的初始過盈。Abaqus自動判斷主面和從面，確定一個限定值來調整小過盈量，超過此限定值的過盈部分不會被調整。例如，圖3中使用通用接觸的默認設置，只調整了一個節點。
圖3 使用默認設置的通用接觸
可通過定義和指定接觸初始化任務，更改限定值。在“Edit Interaction”對話框中，指定初始化任務，如圖4所示。
圖4 通用接觸初始化任務
點擊“Edit Interaction”對話框的初始化任務后面的鉛筆圖標，將新創建的接觸初始化任務分配給一對表面。
展開 
abaqus計算受重力作用下的土體表面重應力為什么不是零
abaqus計算受重力作用下的土體表面重應力為什么不是零呢，這個是下面的個人一些想法，可以供小白參考。
舉一個最簡單例子，假設土體大小10X10X10米，材料密度2000kg/m3.彈性模量100Mpa，泊松比0.3，摩擦角30度，粘聚力30Kpa，只受重力作用，重力加速度取10。單元尺寸大小分別取0.5、1、2、5m。
計算地表豎向應力分別為0.5X104pa、1X104pa、2X104pa、5X104pa，可以看出，單元尺寸越小，地表單元的應力就越小，結果偏于更準確。這是因為重力是作用在每個單元的重心位置，該模型標準矩形，單元也規整，第一層每個單元的標高是單元網格尺寸的一半，第一層重心位置的應力就是密度X重力加速度X該層單元格重心深度，再通過有限元原理轉化到每個單元的節點上，可想而知，要想地表網格節點尺寸為0，必須是單元網格大小足夠小，接近于0，這就是為什么abaqus模擬巖土工程不準確的地方，不可能做到足夠小，一般巖土工程的模型都是比較大的（幾十米幾百米幾千米），模型越大網格尺寸會劃分的很大，精確度也越低。更多案例可以關注抖音abaquser。
展開 高爾夫球的秘密-Abaqus撞擊試驗仿真；揭秘球體分層和凹凸表面背后的湍流
根據網上搜到的試驗數據與材料參數，使用Abaqus對試驗中的3層球進行撞擊響應建模分析。球體按層切分，并賦予指定材料的截面屬性。
通過Abaqus/Explicit分析的高爾夫球撞擊過程、球體的應力以及速度云圖如下：
對于球體在球桿打擊作用下的響應（速度與旋轉速率），在設計環節也會進行大量的計算分析，通常會計算球桿不同表面特征（U型開槽、V型開槽）下的出球響應。
如下圖所示，設計部門在仿真前期會做一些基于試驗參數的對標工作，以矯正仿真分析時高應變率條件下的材料本構模型參數。
在參數修正的基礎上，再進行仿真計算，以更準確地對高爾夫球的動態響應進行預測，從而指導產品設計，縮短研發周期。對標后的仿真基本上可以做到和高速攝影同步。
先歇會，整點啤酒。
現在再來談談，球面上的凹槽怎么回事。上面提到高爾夫球的出球響應中，有個變量是旋轉速率，原來，球在飛行的過程中，不同旋轉速率下，由于凹凸的氣動外形，導致球體產生氣動阻力、升力是完全不一樣的，這也就決定了高爾夫球的運動軌跡。
對于高速飛行的高爾夫球，凹凸的表面會導致湍流，影響球體受力，下面這個視頻是Youtube上ID為CFD Support的團隊通過OpenFOAM計算的不同旋轉速率條件下高爾夫球的升力和阻力系數。有沒有旋轉，差別還是挺顯著的，所以球桿的擊球面要開槽，這樣在出球時，球才會更容易轉起來。
高爾夫球CFD
~上期inp文件下載~
老鼠夾子鏈接：https://pan.baidu.com/s/1TUgt76E8nxz1g3tpjBy44Q 密碼：d73f
展開 在ABAQUS中基于圓柱坐標系設置關于坐標函數的表面力（keyword 曲面加載,圓柱坐標,面力）
物體受力在單一應力狀態下，內部各點也是單一狀態。在部分科研理論分析中，需要在物體內部挖孔，利用逆推法推導物理量。例如下圖所示，受Y方向某拉力作用，各點應力狀態為：
在圓孔中心位置建立圓柱坐標系，該應力狀態在圓柱坐標系下的公式為：
在這種情況下反推物理量，需要對曲面施加基于圓柱坐標系的面力。
案例如下：在圓弧面基于圓柱坐標系施加等效于單向應力狀態的面力。
加載前先建立圓柱坐標系（注意R軸方向為0度位置，T軸方向為角度增大方向，示意圖見文后的加載圖）
具體設置方法為：Load>Create Load>Mechanical>surface traction
選中中間曲面后，先設置徑向力，按以下參數設置：
Distribution：應力分配，點擊后面的f(x)創建一個基于圓柱坐標系的表達式，Local system 要選擇圓柱坐標系，Th為角度變量。
Traction：選擇General，為一般力。
Vector：點擊選擇圖標后，依次選擇（0，0，0） （-1，0，0） ，坐標選擇建立的圓柱坐標系。
注：面力方向矢量是基于所選坐標系，（-1，0，0）就是沿圓柱坐標系下的R軸反向。
Magnitude：選擇應力大小為1。
然后在創建一個Load,設置切向力，如下圖所示，也是基于圓柱坐標系。
再創建一個Load，在整體坐標系下對兩側的平面施加Y方向的面力，大小為1，同時對后面的面施加全約束。
最后加載形式為下圖所示：
求解結果如下圖：
大部分位置應力在0.99~1.01之間，為單向應力狀態，加載方式正確。
本問題的關鍵是面力的方向問題，在選擇面力的方向矢量時，是基于所選坐標系。對于圓柱坐標系，切向力矢量為（0，-1，0）時，即力的方向只沿著
展開 