abaqus單元共節點的案例
ANSYS APDL實體單元和殼單元（不共節點）之間的連接 ￥100
實體單元和殼單元之間的連接是ANSYS中常見的問題。即使兩種單元之間共節點，但單元之間不連續（實體單元每個節點有3個平動自由度，而殼單元每個節點有3個平動自由度和3個轉動自由度），對于兩種單元之間面面接觸，可直接定義剛域，本文主要采用MPC法對實體-殼單元的連接方法進行說明。
1 單元類型
算例模型中，實體單元采用SOLID45，殼單元采用SHELL63，接觸位置不共節點。對于兩種單元之間的連接，通過目標單元TARGE170和接觸單元CONTA175實現，定義約束為實體-殼約束，接觸單元為MPC算法，接觸類型為綁定接觸。
2 有限元模型和綁定接觸
圖1 底部固定約束，殼單元施加均布荷載
圖2 目標單元和接觸單元
3 計算結果
圖3 von Mises stress
圖4 X-Component of displacement
付費內容為相關命令流。
展開 LS-DYNA有厚度/零厚度內聚力單元建模案例（共節點法） ￥9.9
本案例為有厚度/零厚度內聚力單元建模示范，操作軟件為LS-DYNA自帶的前處理器Ls-Prepost。
效果如下圖所示：
無厚度
有厚度
付費內容為：
操作過程PDF文件
步驟詳細，演示形式為靜態圖片與動態圖片；給出了詳細的內聚力單元層的材料參數及出處。
2. 對應的K文件及計算結果
展開 Abaqus多單元耦合分析Step by Step (共57頁) ￥3
Abaqus多單元耦合分析-1-5.pdf
ABAQUS UEL - 損傷材料本構簡單應用于4節點平面單元 ￥300
利用ABAQUS自定義單元子程序，既可以開發新的單元，同時也可以定義新的材料本構模型。本文以損傷模型簡單應用于4節點平面單元為案例，介紹ABAQUS UEL的開發和使用。
如上圖所示，該單元包含4個節點，每個節點有兩個自由度，分別在水平（X）和垂直（Y）方向運動。節點1的兩個自由度被固定，節點4的水平自由度被固定，節點2的垂直自由度被固定。節點3和節點4在垂直方向上向上運動，位移為0.1mm。該正方形單元的邊長為100mm。在input文件里，坐標表示為，
定義節點組合與邊界條件為，
為了讓模型收斂性更好，采用quasi-newton 求解器。時間步設置為，
在文件夾中通過Powershell提交job和子程序，
單個單元的變形為，
采用不同的 ??
，在后處理中得到損傷因子的變化，
相對應的力-時間關系為，
對于多個單元的情況，比如9單元組成的模型，
具體介紹見知乎：ABAQUS UEL - 損傷材料本構簡單應用于4節點平面單元 - 知乎 (zhihu.com)
相對應的UEL代碼和input文件在付費內容中，
展開 
ABAQUS中巧用單元節點偏移
前處理過程并不是在ABAQUS中進行，而是采用ansa或者hypermesh網格劃分，再導入ABAQUS做求解和后處理。有時候可能不小心會遇到一個問題，就是設計人員設計的幾何模型提供給我們，它的幾何之間存在微小間隙，如圖1所示。
圖1 幾何模型示例
我們在處理的時候可能或遇到的問題是劃分網格時候因為間隙太小沒有看到就劃分完網格，而我們在做接觸等設置時候發現這個間隙其實不必要，合在一起并不影響分析結果但是能夠大大提升計算和收斂效率，如果修改幾何在返工重畫網格很麻煩。
或者你會說直接在assembly中通過移動幾何體重新裝配，但是都挺麻煩的。尤其是你導入的是孤立網格的話。
這個時候有個很好的方法就是移動單元節點來填充這個空隙。
a.孤立網格已經導入（或者在ABAQUS中已經劃分好網格）如圖2，
圖2 網格模型
b.通過edit mesh選擇某一個面上的所有單元節點（by angle），如圖3所示。
圖3 節點選擇
c. 偏移單元節點，如圖4所示。偏移的距離需要你事先測量出來，測量方法很多，在3D設計軟件中測量或者在ABAQUS中測量都可以。其中CSYS是你選擇的偏移參考系，默認全局坐標系，如果你的模型使用全局坐標系不方便你可以自己在節點面上建立一個垂直的坐標系，那樣偏移就很方便。圖中1,2,3就是對應xyz坐標軸方向，在對應方向輸入對應的距離就行。
圖4 節點偏移
d.完成之后如圖5所示。我這里是偏移了0.5mm。可以看到就自動將節點偏移過去，免去很多返工操作，非常方便。接下來該怎么設置就怎么設置。（這個比設置接觸容差應該好多了，不過具體對比結果我還沒有去試）
圖5 結果
總結：對于微小間隙如果需要處理，又為了避免返工提高分析效率，可以考慮直接節點偏移操作，可理解為將面上網格沿某一方向擴大一定倍數。
展開 ABAQUS 單元-節點排布順序解析（重點講解分析方法）
進行ABAQUS UEL二次開發、或者研究界面問題的時候，比如cohsive單元界面問題，會涉及到單元-節點的排布順序。ABAQUS inp文件中的單元-節點順序排布很有講究，不能搞錯，這是因為節點的排布順序與內部程序有關聯，內部的程序我們不得而知，但是節點順序的規律可以從inp文件中看到，再對比cae節點中的節點編號可以總結規律。以下內容介紹一些經典的實體單元-節點排布順序。
這部分內容其實挺好玩兒的！！！
（一）三角形單元
以一個10x10（無單位）的挖孔板為例，劃分的網格為：（注意到被單元為30、138和23包圍的三角形沒有編號，這是因為那個三角形不是單元，我在那里挖了個孔，因為背景顯示的問題，網格顏色和背景一樣，故沒有區分。作此說明）
這里我們首先查看節點的排布方式。分別以單元1、12、123為例。
編號為“1”的單元，在inp文件中，單元-節點的編號為：
1, 61, 1, 2
在cae界面中，我們只顯示編號為“1”的單元為：
發現節點排布的方式是“逆時針”順序。
編號為“12”的單元，在inp文件中，單元-節點的編號為：
12, 30, 9, 10
在cae界面中，我們只顯示編號為“12”的單元為：
發現節點排布的方式是“逆時針”順序。
編號為“123”的單元，在inp文件中，單元-節點的編號為：
123, 83, 34, 33
在cae界面中，我們只顯示編號為“12”的單元為：
發現節點排布的方式是“逆時針”順序。
至此我們可以得出結論，在一階三角形單元中，單元-節點的索引方式遵循“逆時針”順序，即按照節點順序彎曲右手四個手指，大拇指指向屏幕外。
展開 ABAQUS UEL-梯度損傷模型應用于4節點平面單元 ￥600
傳統損傷模型對于單元的尺寸十分敏感，不同單元尺寸會導致有限元模型精度出現明顯偏差。針對該問題，梯度損傷（Gradient-damage）模型的概念被提了出來。
本文詳細介紹了如何將梯度損傷模型應用于4節點平面單元，并在有限元模型中進行模擬。
ABAQUS提供了UEL（user defined element）給使用者進行開發。筆者利用UEL開發4節點平面單元，其邊界條件如下圖所示。其中，節點1的X、Y方向被限制住，節點2的Y方向被限制，節點4的X方向被限制，節點3、4的Y方向有豎向位移0.1mm。單元為100*100mm的二維正方形。
每個節點除了X和Y方向的位移，還帶有非局部應變（nonlocal strain）。
單個單元模型，
多個單元模型，
具體內容可參見知乎文章：
ABAQUS UEL-梯度損傷模型應用于4節點平面單元 - 知乎 (zhihu.com)
相應的input文件和uel代碼付費可見，
展開 ABAQUS-如何求結構的節點位移單元應力分量和支反力
一.問題
如圖1所示，顯示了四根桁架結構的尺寸與約束情況，材料為鋼，彈性模量設置為2.96Gpa，橫截面積為100，求該模型的
節點位移、單元應力分量、支反力。
圖1
二.部件與材料
首先按照圖1創建部件，選擇二維平面，特征為
線，繪制相應的草圖，并生成實體，命名為link。
圖2
如圖3所示，在屬性模塊創建材料屬性，選擇力學-彈性，在彈框中填寫楊氏模量2960，泊松比0。
圖3
如圖4、5、6所示，創建截面，選擇類型梁-桁架，并賦予材料屬性，填寫截面面積為100。
圖4
圖5
圖6
三.裝配與分析步
如圖7、8所示，將部件進行裝配，創建靜力通用分析步。
圖7
圖8
四.邊界條件與載荷
如圖9所示，創建邊界條件，選擇轉角/位移，約束點1與點2U1和U2的位移；約束點3U2的位移。
圖9
如圖10 所示，創建載荷，選擇點4，給定集中力CF2=-150N。
圖10
五.網格與作業提交
選擇單元族為桁架，單元類型為T2D2，提交作業。
圖11
圖12
六.結果展示與后處理
圖13是放大之后的變形圖，可以看到點4向右下方位移。
圖13
對結果進行處理，提取模型節點位移，單元應力分量和支反力。首先點擊進入可視化模塊，依次點擊選項-通用，在彈框中選擇標簽，勾選顯示單元編號和顯示節點編號，在右側可以自行選擇顏色。
圖14
依次點擊報告-場輸出，在彈框中位置下拉框中選擇唯一節點的，在新窗口選擇RF下拉菜單中的RF1與RF2；選擇U下拉菜單中的U1與U2。
展開 輸出abaqus inp如何去除instance信息，而且單元節點編號也不會有重復
解決的方法：在command line輸入 mdb.models['modelName'].setValues(noPartsInputFile=ON) 重新生成inp文件，這樣在寫出的inp文件里面就沒有instance信息，而且單元節點編號也不會有重復了。
ABAQUS中隱式和顯式的節點和單元的輸出變量解析
FSLIP
Field: yes History: no .fil: no
Length of contact slip path at secondary nodes during contact (FSLIPEQ) and in some
cases (see About contact pairs in Abaqus/Explicit) components of net contact slip in local
tangent directions (FSLIP1 and FSLIP2). These variables remain constant while a
secondary node is not in contact.
FSLIPR
Field: yes History: no .fil: no
Magnitude of contact slip rate at secondary nodes during contact (FSLIPR) and in some
cases (see About contact pairs in Abaqus/Explicit) components of contact slip rate in local
tangent directions (FSLIPR1 and FSLIPR2). These variables are set to zero while a
secondary node is not in contact.
BONDSTAT
Field: no History: yes .fil: no
Spot weld bond status.
展開 【腳本】ABAQUS中判斷某個節點是否在某個單元集 ￥1.99
〇、問題來源
本題目來源于技術鄰ABAQUS主題下的問題
Abaqus Python語言判斷node是否在單元里？
if node1 在集合set1中，請問有這種語句嗎，小白請假不勝感激。
http://www.yqgqt.org.cn/answer/1197962
一、問題分析
這個屬于簡單的腳本判斷，該問題也得到了技術鄰大佬“藍牙”和“君莫”的回答，但大佬的回答并不能夠直接幫助小白解決問題，所以斗膽一試，還請指正。
本問題應該遍歷set-1中所有單元的節點，然后與特定節點編號進行對比，若等于則說明該節點在set-1中，否則則不在。
根據以上思路編寫腳本。
二、腳本測試
在ABAQUS中隨便創建一個模型，劃分網格后形成裝配體，在裝配體中創建一個名為“Set-1”的集合，集合內包含一部分單元，如下
運行編寫的腳本
會彈出如下窗口，要求輸入要查詢的節點ID
然后，程序開始運行，若該節點在Set-1中，則輸出如下信息
否則，輸出如下信息
以上。
最后，感謝關注與支持，如果覺得有用，不如點個關注和收藏？
最后，歡迎關注微信公眾號：ABAQUS在線，讓我們一起愉快的成長吧。
展開 
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題
摘要：
采用基于ABAQUS的UEL子程序開發4節點平面應變等參單元，采用雙線性形函數，4點高斯積分，本構關系為線彈性各向同性材料，得到的單元剛度矩陣和ABABUS自帶的CPE4單元的單元剛度矩陣（剛度矩陣輸出方式為*element matrix output, elset= ALLE, stiffness=yes, OUTPUT FILE=USER DEFINED）不同；對比ANSYS的單元剛度矩陣，結果顯示兩者也不相同。問題出在哪里呢？本文檔將對此問題進行回答。
本文可以作為ABAQUS高級子程序UEL的入門級教程，做UEL的應該關注下！
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題（技術鄰 藍牙）.pdf
展開 ABAQUS Fortran基于提供的代碼 P53，完善 8 節點單元程序，并增加對應力結果的處理。
1.改寫輸入數據格式，使之能適應任意幾何（可利用節點坐標輸 入節點，利用單元-節點關系輸入單元）； 2. 計算節點應力，給出并實現至少一種應力處理方案，提供處理 前后的應力結果（可用表格和云圖表示），可與其它軟件對比； 3.提交總結報告（包括方法/方案描述、帶詳細注釋的代碼、程序框圖、算例描述、結果比較分析等）、可編譯源代碼、可執行文件、 數據文件、結果文件
program p53
!-----------------------------------------------------------------------------
! program 5.3 plane strain of an elastic solid using uniform
! 8-node quadrilateral elements numbered in the x direction
!-----------------------------------------------------------------------------
use new_library ; use geometry_lib ; implicit none
integer::nels,nxe,neq,nband,nn,nr,nip,nodof=2,nod=8,nst=3,ndof,loaded_nodes,&
i,k,iel,ndim=2
real::aa,bb,e,v,det ; character(len=15) :: element = 'quadrilateral'
!
展開 有限元中單元積分點與節點應力相互轉換(CPE4為例)（ABAQUS）
在ABAQUS中,當需要獲取節點上的應力時,可以在后處理中建立路徑或者用查詢功能等獲取.
但是當需要大量的節點上應力數據時,很多人會用Python編程進行大批量的提取應力.但是提取出來的應力為單元積分點上的應力.無法獲取節點上的應力.同時在ABAQUS中的子程序中,也是對積分點上的數據進行操作.
本文基于個人興趣同時想要更加了解有限元背后原理和公式的想法.近日進行了一些初步的探索.希望大家批評指正. 本文基本不涉及原理公式，只在轉換積分點和節點的應力時列出公式。盡可能簡介易懂。
一: 單元類型及節點數目與位移,應變,應力階次的關系
本節內容基于有限元教材及一些網上資料.
(1)有限元求解的思路是:
一: 建立單元節點力與節點位移關系式.
二: 將彈性體上的外載荷等效移置到節點上.
三: 在節點上建立力的平衡方程,求得節點位移.
四: 通過彈性力學基本方程,可求得單元的應力和應變.
(2) 四節點矩形單元
以四節點矩形單元為例，在此只表達有限元教材中的結論，具體公式可參考有限元教材。
（3）ABAQUS中的CPE4單元
CPE4: A 4-node bilinear plane strain quadrilateral.
該單元有四個節點，同時有四個積分點。
對于每個應力分量（注意：在此只看一個應力分量），單元內任一點（x,y）的應力表達式為：
stress=a*x*y+b*x+c*y+d （1）
該表達式有四個未知量：a,b,c,d。
若知道四個積分點的應力分量。將每個積分點帶入上式，則會形成包含四個方程的線性方程組。
展開 