ansys單元函數設置的案例
從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異
如題，《從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異》，形函數對結果的影響大部分人都能聯想到二次單元比線性單元求得的結果更精確，但該文要表達的不僅如此，而是從更一般地討論怎么從單元的形函數來理解節點解與單元解之間的差異。
首先討論單元的階次。作為基礎我們應該明白網格與單元的區別，網格是將幾何體離散化后的結構，即組成幾何體的微元，單元是這些微元的幾何、物理或數學屬性（這里我們并不打算詳細討論單元的這些屬性，但是這些知識會方便對本文的理解）。我們經常在使用ansys或其他CAE軟件時經常會遇到單元的選擇以及單元階次的選擇，一般一種單元包括線性單元和二次單元甚至更高級的單元，比如在ansys中經常被使用的shell181（左）和shell281（右），線性單元使用的形函數是一次的多項式，高次單元使用的形函數是高次的多項式，形函數用于描述相鄰節點之間的位移場，所以高次的單元可以更好的描述形狀復雜的幾何體。
不同于常規材料力學中通過平衡方程求解（首先求得的解是力解），有限元方式求解的特點是首先求解出的結果是節點的位移解，即displacement of nodes，所有的節點位移形成了位移場，在空間上位移場一定是連續的，但是不一定是平滑的。哎哎，是不是特別熟悉的感覺，正是和高數中函數的連續性和可導性兩個性質非常相似，不用奇怪，位移場本來就是用函數描述的，所以自然就存在函數的性質，所以用函數的性質來理解就可以方便解釋一些現象了，下圖分別是用兩種形函數描述的位移場，在有限元求解后得到的首先是節點位移解，即圖中5個節點的位移，假如每個節點的位移用坐標x\y\z的函數來表示，然后通過形函數插值得到相鄰節點之間的位移（也是xyz的函數），上圖是用一次形函數插值，下圖是用二次形函數插值。
展開 ANSYS APDL BEAM 單元的截面設置
選擇梁單元的軸線
latt,1,,1,,7,8,1 !將材料號、截面參考號、實常數（如果有的話）、方向關鍵點等信息分配給
!上面已經選擇好的還沒有劃分單元的梁軸線/
lesize,all,,,10 !指定梁縱向劃分網格的尺寸。由于前面已經用LSEL命令選擇好了的線就是梁的中軸線
!所以不需要再次選擇（ANSYS里，選擇好的實體會有個標志，除非你用命令改變了它們）
lmesh,all !劃分網格，好了，你可以再改變參數，增加荷載項并求解啦。
【附注】
把在ansys中使用梁單元的主意事項列于下：
1. beam188、beam189在section中設定參數；而beam3、beam4則必須在實常數中設置，其中橫截面積、彎曲慣性矩以及扭轉慣性矩是必須填入的，截面厚度（TKY、TKZ）只在圖形顯示中有用，計算的時候并不用到它，看一下梁單元剛度矩陣的推導就可明白，ansys的理論手冊也有梁單元剛度陣元素的詳細介紹。beam188、beam189 是從ansys5.5版本開始出現的兩種新的梁單元，基于Timoshenko梁理論，適于細長梁的計算分析，考慮了剪切變形的影響。
2. 梁單元以實體的形式顯示。大家知道，在ansys中，梁單元默認都是顯示的線條。但是我們可以將賦予了section屬性的梁顯示成實體，這樣做的好處是，更加形象，直觀，可以對梁的布置正確與否作出準確的判斷。方法是：在Utility Menu->PlotCtrls->Style->Size and Shape菜單下，將Display of element后的單選打開，即使其為on的狀態。
3. 單元彎矩圖的繪制。
展開 ansys workbench中設置變厚度殼單元
對于厚度尺寸相對于其他幾何尺寸較小的結構，我們常常采用殼單元來代替三維實體單元進行分析。殼單元模型雖然不像三維實體模型那樣更接近真實模型，但其單元及節點數量少，計算量小，在工程中對復雜模型進行簡化時，采用殼單元能大大降低工作量和計算難度。
在建立殼單元模型時，我們需要輸入殼的厚度值，該厚度值可以在DM中設置，也可以在Mechanical中設置。DM中僅允許輸入常量厚度值（即等厚度），在Mechanical中可以設置隨某一坐標變量變化的厚度值。
等厚度模型
厚度隨坐標變化的模型
大多數情況下，以上厚度設置是能夠滿足工程分析需要的。但是，有一天突發奇想，我想建一個厚度值隨多個坐標值變化的模型，現有的方法以函數進行輸入厚度隨坐標變化時，只允許輸入一個變量，怎么辦？
workbench提供了一個很好的工具—External Data。用它，可以將任意位置的厚度值進行任意編輯，然后導入到Mechanical中。
展開 