ansys 殼單元應力的案例
使用等效結構應力法預測殼單元/實體單元焊趾的疲勞壽命
建立了以殼單元和實體單元建模的焊縫模型，并標記了焊趾點位置。
2.在不同的兩個工況天下對模型施加兩種載荷，并計算焊趾處的節點結構應力。
3.提取兩種模型焊趾處的節點力。
4.使用自己編寫的代碼計算兩種模型的焊趾等效結構應力，并計算損傷。
有意咨詢代碼或算法相關問題的可私聊我。
WB14.0水上滑道結構應力分析（殼單元，梁單元在WB中應用）
復件 水上滑道結構應力分析報告.doc
這是我做的一個游樂設施的 仿真分析的報告，由于涉及企業的隱私，刪除了數據和名稱。
報告的關鍵點：殼單元和梁單元在workbench中的應用，殼單元和梁單元是DM建模的顯著特征，在WB中建模，分析，方便，快捷。
關于dyna殼單元應力問題
當我使用nip為5的殼單元時，我的殼單元輸出的應力是5個積分點的和還是平均值？
ANSYS中桿單元和殼單元的單元耦合問題
在比較復雜的結構的有限元分析中，不同的結構部件通常使用不同類型的單元來模擬。
通常情況下，不同類型的單元的各個節點的自由度數目是不同的，不同類型單元的連接節點處的自由度的耦合問題，是一個比較令人頭疼的問題。
在ANSYS中通常可以用耦合命令CP來耦合不同類型單元在連接節點處的自由度(DOF)。
也可以用CE命令來認為添加自由度之間的約束方程來達到耦合的目的。
下面是一個簡單的算例，使用了CE命令來耦合連接節點處的自由度。
模型是航天器的機翼的一個Section的某一個隔框。上下表皮是薄殼結構，用Shell63單元來模擬，在上下表皮之間有起支撐作用的桿件，用link8單元來模擬。
建模的時候，link8單元和shell63單元在連接有各自獨立的節點。即：link8單元和shell63單元的節點在連接處是重合的，但是，節點編號是各自獨立的。
link8單元在每個節點有 ux,uy,uz3個平動自由度；
shell63在每個節點有ux,uy,uz這3個平動自由度和rotx,roty,rotz這3個轉個自由，共6個自由度。
在耦合節點處，兩個耦合節點的ux,uy,uz自由度應該是相等的。
這個等式可以用CE命令來描述。
完整的命令流如下：
finish
/clear,start
/prep7
!定義第一種材料屬性；
mp,ex,1,30e6
mp,prxy,1,0.3
!定義shell63單元和實常數；
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, ,
在關鍵點處生成節點；
nkpt,100,4 !與編號為117的節點耦合
nkpt,101,9 !
展開 
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
付費內容為相關命令流。
展開 ANSYS各類型單元連接專題講解(五)之3D梁單元與殼單元剛接
例如采用ANSYS模擬一個多層混凝土框架結構，一般除計算整體指標外，我們在計算具體荷載作用時（如風荷載、地震作用、恒載、活載等），樓板一般采用彈性版，此時可用殼單元模擬，主梁、次梁采用梁單元模擬，此時變為梁單元包含在殼面內的情況，當然此類情況是否需要考慮截面偏置，可根據具體工程而定。
對這中梁單元包含在殼單元面內的情況，只需要將梁單元與殼單元共用節點即可，而無須格外建立約束方程。
三、梁單元在殼單元內但不包含
此種情況為梁與殼位于同一面內，但其中面不包含梁線，適用于多尺度建模分析（如下圖）。梁單元與殼單元的連接在端部可以通過剛性梁和剛性區域兩種方式連接。剛性梁采用MPC184單元，剛性區域采用Cerig命令，具體使用方法下期文章討論。
展開 天方地圓結構-梁殼單元建模實例！再次驗證應力奇異的可怕性！
應力集中和應力奇異的驗證
▲▲▲
因此模型存在矩形截面，即有90°夾角部分，從理論上來說，此部分位置必然會產生應力集中，筆者通過不斷細化網格求解驗證發現，網格越細化，應力值則越大，如果僅僅出現應力集中的話，即使網格不斷細化，那么最終的應力值也會趨于一個穩定值，但此模型則是隨著網格的不斷細化，應力值似乎有無限變大的趨勢，所以此時，應力集中的部位已然產生了應力奇異，但此應力奇異是實際模型本身的幾何缺陷造成的，因而是無可避免的。
因而，對于有幾何結構突變部位并非網格越細化越好的，網格細化卻帶來了無法判斷結果準確性的一個可怕后果，建議：
1. 在建模的時候如果存在幾何突變部位，務必要采用倒圓角進行圓滑過渡，以避免應力集中變成了應力奇異（即使真實模型中并未有圓角過渡）。
2. 梁殼單元在突變部位隨網格細化產生應力奇異是無可避免的，所以盡量避免在有幾何突變的部位采用梁殼單元，最好采用實體單元。
3. 實體建模時，四面體網格也會隨網格細化產生無可避免的應力奇異，所以最好采用六面體網格，這一點壓力容器行業基本上都采用六面體網格。
應力分析僅僅有漂亮的計算云圖和計算結果是遠遠不夠的，分析設計人必須具有判斷計算結果準確性的能力，試想如果以并非準確的應力結果作為評定依據，那將是多么可怕的一件事情！
有從事壓力容器分析設計行業的朋友如感興趣可掃描上方二維碼關注一下公眾號！
展開 ANSYS殼單元分析箱梁
選定殼中部截面
shell,mid
! Define a path to plot the stress
! 定義一個路徑來繪制應力
PATH,MidPath1,4,30,20,
! Define path with nodes
! 用節點來定義路徑
PPATH,1,89
PPATH,2,34
PPATH,3,570
PPATH,4,1066
PATH,STAT
!*
AVPRIN,0,0,
!*
! Map stress Z to the path
! 將Stress Z映射到路徑
PDEF, ,S,Z,AVG
/PBC,PATH, ,0
! Plot the stress
! 繪制應力
PLPATH,SZ
FINISH
展開 ANSYS中薄殼厚殼分類及單元特性
②平行于中曲面的面素上的正應力與其它應力相比可忽略不計。
但上述假定同時假定了兩種不相容的變形狀態，即平面應變和平面應力狀態。因此許多學者提出了許多修正理論，但是只要是基于 Kirchhoff-Love 假定為基礎的薄殼理論，其精度都不會超過 Kirchhoff-Love 理論的精度范圍。
為構造協調的薄板殼單元，可采用多種方法，如增加自由度法、再分割法（也稱復合法）、離散克希霍夫（Discrete Kirchhoff Theory）法等，但都適用于薄板殼結構，也不考慮橫向剪切變形的影響。
5. 考慮橫向剪切變形的殼理論
可考慮橫向剪切變形影響的理論，一般稱為 Mindlin-Reissner 理論，是將 Reissner 關于中厚板理論的假定推廣到殼中。
ANSYS殼單元
薄板殼單元基于 Kirchhoff-Love 理論，即不計橫向剪切變形的影響；中厚板殼單元則基于 Mindlin-Reissner 理論，考慮橫向剪切變形的影響。
在 ANSYS中，SHELL 單元采用平面應力單元和板殼彎曲單元的疊加。除SHELL63、SHELL51、SHELL61 不計橫向剪切變形外（可用于薄板殼分析），其余均計入橫向剪切變形的影響（可用于中厚板殼分析）。
對于板殼單元還應注意以下幾個問題：
⑴ 面內行為
由于面內采用平面應力狀態，因此不存在“體積鎖死”問題，但“剪切自鎖”問題依然存在，因此許多單元采用了 ESF 以響應面內行為， 如 SHELL41、SHELL43 和SHELL63 單元等，SHELL181 支持橫向剪切剛度的讀入。
展開 ansys workbench中設置變厚度殼單元
對于厚度尺寸相對于其他幾何尺寸較小的結構，我們常常采用殼單元來代替三維實體單元進行分析。殼單元模型雖然不像三維實體模型那樣更接近真實模型，但其單元及節點數量少，計算量小，在工程中對復雜模型進行簡化時，采用殼單元能大大降低工作量和計算難度。
在建立殼單元模型時，我們需要輸入殼的厚度值，該厚度值可以在DM中設置，也可以在Mechanical中設置。DM中僅允許輸入常量厚度值（即等厚度），在Mechanical中可以設置隨某一坐標變量變化的厚度值。
等厚度模型
厚度隨坐標變化的模型
大多數情況下，以上厚度設置是能夠滿足工程分析需要的。但是，有一天突發奇想，我想建一個厚度值隨多個坐標值變化的模型，現有的方法以函數進行輸入厚度隨坐標變化時，只允許輸入一個變量，怎么辦？
workbench提供了一個很好的工具—External Data。用它，可以將任意位置的厚度值進行任意編輯，然后導入到Mechanical中。
展開 ANSYS巧用殼單元給實體劃分六面體網格
ANSYS巧用殼單元給實體劃分六面體網格
1 概述
眾所周知，ANSYS經典劃分網格的功能比較弱，映射劃分（Map）和掃掠劃分(Sweep)對幾何形狀的要求都十分高。而四面體網格一方面導致單元數目多余六面體，一方面給計算后處理帶來一定的不便。
有些情況下，幾何模型的結構導致即使再怎么對模型進行切分都不可能掃掠出六面體網格，這種情況下，可以巧妙地利用殼單元。
ANSYS經典里對于一個平面，劃分網格非常簡單，而且幾何形狀約束很少，即使是自由劃分的網格，一般情況下都比較規整。利用這個特點，用殼單元對面進行網格劃分，然后再對整個實體模型進行網格劃分。本次以一個例子示意此過程。
2 過程
首先在ANSYS經典界面定義兩個單元類型，分別是shell181和solid185。如圖1所示。
圖1 單元類型
建立幾何模型，采用block命令，建立100x40x10的長方體：
block,-50,50,0,10,-20,20
如圖2所示。
圖2 幾何模型
之后為了演示網格劃分，將模型切分成幾塊，如圖3所示。
圖3 切分模型
再然后選擇殼單元shell181，如圖4所示：
圖4 選擇shell181單元
然后設置模型最上層的面各個線條的分數：
圖5 操作
份數分別如圖6.
圖6 線條份數
之后點擊MeshTool，如圖7所示。
圖7 劃分面網格設置
如圖7設置，點擊Mesh，選中模型的最上一層表面劃分，得到圖8的結果。
圖8 面網格劃分
再重復前面的選擇單元的操作，選擇單元類型為solid185，并且在MeshTool里選擇Volumes 的掃掠（sweep）劃分，如圖9所示。
圖9 操作
點擊Sweep選中所有的體，即得到如圖10所示的網格。
展開 
ANSYS巧用殼單元給實體劃分六面體網格
ANSYS巧用殼單元給實體劃分六面體網格
1 概述
眾所周知，ANSYS經典劃分網格的功能比較弱，映射劃分（Map）和掃掠劃分(Sweep)對幾何形狀的要求都十分高。而四面體網格一方面導致單元數目多余六面體，一方面給計算后處理帶來一定的不便。
有些情況下，幾何模型的結構導致即使再怎么對模型進行切分都不可能掃掠出六面體網格，這種情況下，可以巧妙地利用殼單元。
ANSYS經典里對于一個平面，劃分網格非常簡單，而且幾何形狀約束很少，即使是自由劃分的網格，一般情況下都比較規整。利用這個特點，用殼單元對面進行網格劃分，然后再對整個實體模型進行網格劃分。本次以一個例子示意此過程。
2 過程
首先在ANSYS經典界面定義兩個單元類型，分別是shell181和solid185。如圖1所示。
建立幾何模型，采用block命令，建立100x40x10的長方體：
block,-50,50,0,10,-20,20
如圖2所示。
圖2 幾何模型
之后為了演示網格劃分，將模型切分成幾塊，如圖3所示。
圖3 切分模型
再然后選擇殼單元shell181，如圖4所示：
圖4 選擇shell181單元
然后設置模型最上層的面各個線條的分數：
圖5 操作
份數分別如圖6.
圖6 線條份數
之后點擊MeshTool，如圖7所示。
圖7 劃分面網格設置
如圖7設置，點擊Mesh，選中模型的最上一層表面劃分，得到圖8的結果。
圖8 面網格劃分
再重復前面的選擇單元的操作，選擇單元類型為solid185，并且在MeshTool里選擇Volumes 的掃掠（sweep）劃分，如圖9所示。
圖9 操作
點擊Sweep選中所有的體，即得到如圖10所示的網格。
展開 ANSYS巧用殼單元給實體劃分六面體網格
ANSYS巧用殼單元給實體劃分六面體網格
1 概述
眾所周知，ANSYS經典劃分網格的功能比較弱，映射劃分（Map）和掃掠劃分(Sweep)對幾何形狀的要求都十分高。而四面體網格一方面導致單元數目多余六面體，一方面給計算后處理帶來一定的不便。
有些情況下，幾何模型的結構導致即使再怎么對模型進行切分都不可能掃掠出六面體網格，這種情況下，可以巧妙地利用殼單元。
ANSYS經典里對于一個平面，劃分網格非常簡單，而且幾何形狀約束很少，即使是自由劃分的網格，一般情況下都比較規整。利用這個特點，用殼單元對面進行網格劃分，然后再對整個實體模型進行網格劃分。本次以一個例子示意此過程。
2 過程
首先在ANSYS經典界面定義兩個單元類型，分別是shell181和solid185。如圖1所示。
圖1 單元類型
建立幾何模型，采用block命令，建立100x40x10的長方體：
block,-50,50,0,10,-20,20
如圖2所示。
圖2 幾何模型
之后為了演示網格劃分，將模型切分成幾塊，如圖3所示。
圖3 切分模型
再然后選擇殼單元shell181，如圖4所示：
圖4 選擇shell181單元
然后設置模型最上層的面各個線條的分數：
圖5 操作
份數分別如圖6.
圖6 線條份數
之后點擊MeshTool，如圖7所示。
圖7 劃分面網格設置
如圖7設置，點擊Mesh，選中模型的最上一層表面劃分，得到圖8的結果。
圖8 面網格劃分
再重復前面的選擇單元的操作，選擇單元類型為solid185，并且在MeshTool里選擇Volumes 的掃掠（sweep）劃分，如圖9所示。
圖9 操作
點擊Sweep選中所有的體，即得到如圖10所示的網格。
展開 ANSYS中如何獲取采用殼單元模擬時的截面內力
部分朋友反應在采用殼單元進行仿真計算時不知如何提取殼單元的截面內力，今日水哥就殼單元的截面內力提取方法簡單說明下，供諸君參考一二。
首先講講殼單元的應力和內力輸出。
薄殼單元和中厚板殼單元應力和內力的輸出項目不盡相同，對于薄殼單元如 SHELL63 就不輸出次要應力（τxz、τyz）和內力（Nx、Ny），而中厚板殼單元則輸出這些應力和內力。
注意，殼單元的內力輸出均是相對于單元坐標系，單元各邊內力相同，為該單元單位長度上的內力，如 Mx 的單位為“力×長度/長度”，如需該單元的總彎矩則再乘以單元邊長即可。單元的內力可通過單元表輸出，例如shell181的結果輸出示意圖如圖，單元表選項如下：
上述方法針對的是單個單元，然而實際計算過程中，我們常常需要獲取某個截面的總內力，此時可通過計算獲取。一般而言，有兩種方式，一種是路徑積分法，另外一種是單元節點力求和法。水哥個人建議采用單元節點力求和法，簡單快捷。
單元節點力求和法需要掌握兩個命令：Spoint \ Fsum
Spoint,node,x,y,z
該命令定義力矩求和的位置點，如果求和不位于總體直角坐標系下，可輸入node定義或采用Rsys命令定義。
Fsum,lab,Item
該命令計算所選擇單元集中選擇節點集的所有節點力的合力和合力矩。因而在求具體某截面的內力時，應選擇該截面附件的單元以及節點。
下面以某懸臂板為例，闡述基本思路。
某混凝土懸臂板，板厚100mm，尺寸為900mmX2000mm，混凝土等級為C30，在板的端部100mm范圍內受到均布荷載0.5KN/m^2，求板跨中間截面的剪力以及彎矩。
展開 ansys Workbench 靜應力模塊，利用生死單元技術結合APDL命令，模擬轉軸最大扭力 ￥10
?
ansys Workbench 靜應力模塊，利用生死單元技術結合APDL命令，模擬轉軸最大扭力
示例：要求計算轉軸所能承受的最大扭轉力矩，轉軸抗拉強度1230MPa
模型如下： 中間最細位置R=3
Workbench計算時，左側固定。右側面施加圓轉位移。
效果展示
?
操作過程：
首先，初步計算轉軸旋轉多少會接近許用最大值1000Mpa。確定初始載荷大小。
當加載1° ——0.0174 弧度 ，時 轉軸約945Mpa。
其次，利用APDL命令分載荷步逐步增大轉角載荷，并在每個載荷步中進入后處理中查看是否有單元應力超過許用值1000Mpa。當有單元超過許用值時記錄該單元，在下一步載荷過程中將該單元抑制。繼續加載直到循環結束。
1.創建加載點——remotePoint
在Pilot Node APDL Name 中定義名稱：后期將在插入的APDL命令中使用該名稱，更改載荷大小。
創建單元組——Name Selection
在每個載荷步的后處理中需要篩選單元結果，查看是否超過許用應力。為了縮小查詢范圍可以先根據經驗判斷危險截面位置，將危險截面附近的單元定義為一個組。在后期結果查看時，僅在該組內查找單元應力。從而提高計算效率。
注意：選著的是單元組，可以使用框選功能。
在Analysis setting 中插入Command 命令
插入命令如下所示，同時注意單位制的選著，本例使用mm kg N。 命令見附錄
命令中包含有三種 應力評估方法，一:剪應力失效。二：等效應力失效。三：第一主應力失效。應根據實際工況條，結合零部件失效模式，自主選著。
!!!!!1.使用剪切應力判斷是否失效*********************
!
展開 