ansys中函數加載的案例
手把手教你ANSYS的函數加載
最近論壇里很多人問我，如何施加函數載荷。
我今天給大家奉獻的是任意函數加載的操作步驟詳解，手把手教大家操作自定義的函數加載。
如果覺得還不錯，頂下帖子，也算對我的鼓勵了！
大家有什么ANSYS 或 Workbench Mechanical 相關的問題，可以隨時**我 ansys123@qq.com
手把手教你ansys函數加載.doc
Ansys Wrokbench分段復雜函數載荷，加載方式記錄 ￥10
問題：
Ansys Workbench的載荷加載形式有三種，constant/table/function。Constant是在載荷步內給定恒定值；table形式較為便捷，可以在定義每個子步的載荷大小； function形式可以輸入以time/X/Y/Z為變量的簡單方程。
但是仍有某些形式的載荷較難輸入，例如分段復雜函數載荷等。
解決方法：
需要使用Ansys經典界面的function功能編輯分段載荷獲得ADPL載荷命令；再利用Workbench中command的形式施加載荷。
操作方式：
1. Ansys經典中function公式編輯器輸入分段函數。
在function頁卡中選著變量time，在Regime頁卡中逐個定義分段函數；
定義完成后點擊保存，并輸入函數名“TEST3.func”
2. 再次點擊標題欄的Parameters>Functions>Read From files>找到剛才保存的TEST3.func。并在Table Parameter Name中給編輯導入的分段函數命名PForce。此后分段函數即被公式編輯器編譯為表格數組形式，數組的名稱為：PForce。
3. 提取分段函數數值的ADPL命令形式，用于Workbench使用。
完成分段函數導入和命名后，在下拉列表中的File>List>Log file中可以查看經典界面GUI操作對應的ADPL命令。在這里可以將上述function公式編輯器導入的分段函數數組對應ADPL命令顯示出來。（有時log file顯示不及時，再重復一次即可）
4. 在Workbench內創建加載remote point點，并設定加載點的ADPL name為“LoadPoint“，用于加載。
展開 在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）時，即力的方向只沿著theta的反方向。
展開 從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異
如題，《從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異》，形函數對結果的影響大部分人都能聯想到二次單元比線性單元求得的結果更精確，但該文要表達的不僅如此，而是從更一般地討論怎么從單元的形函數來理解節點解與單元解之間的差異。
首先討論單元的階次。作為基礎我們應該明白網格與單元的區別，網格是將幾何體離散化后的結構，即組成幾何體的微元，單元是這些微元的幾何、物理或數學屬性（這里我們并不打算詳細討論單元的這些屬性，但是這些知識會方便對本文的理解）。我們經常在使用ansys或其他CAE軟件時經常會遇到單元的選擇以及單元階次的選擇，一般一種單元包括線性單元和二次單元甚至更高級的單元，比如在ansys中經常被使用的shell181（左）和shell281（右），線性單元使用的形函數是一次的多項式，高次單元使用的形函數是高次的多項式，形函數用于描述相鄰節點之間的位移場，所以高次的單元可以更好的描述形狀復雜的幾何體。
不同于常規材料力學中通過平衡方程求解（首先求得的解是力解），有限元方式求解的特點是首先求解出的結果是節點的位移解，即displacement of nodes，所有的節點位移形成了位移場，在空間上位移場一定是連續的，但是不一定是平滑的。哎哎，是不是特別熟悉的感覺，正是和高數中函數的連續性和可導性兩個性質非常相似，不用奇怪，位移場本來就是用函數描述的，所以自然就存在函數的性質，所以用函數的性質來理解就可以方便解釋一些現象了，下圖分別是用兩種形函數描述的位移場，在有限元求解后得到的首先是節點位移解，即圖中5個節點的位移，假如每個節點的位移用坐標x\y\z的函數來表示，然后通過形函數插值得到相鄰節點之間的位移（也是xyz的函數），上圖是用一次形函數插值，下圖是用二次形函數插值。
展開 
ANSYS 中添加窗函數
比如說我要施加一個周期性的沖擊 希望他作用2個周期后停止
ANSYS中可以使用的數學函數
在ANSYS幫助系統中關于*SET命令的注釋下列出了ANSYS中可以使用的數學函數。所有這些數學函數均可以在ANSYS環境中使用，這些數學函數包括：
ABS(X) 求絕對值
ACOS(X) 反余弦
ASIN(X) 反正弦
ATAN(X) 反正切
ATAN2(X,Y) 反正切, ArcTangent of (Y/X) , 可以考慮變量X,Y 的符號
COS(X) 求余弦
COSH(X) 雙曲余弦
EXP(X) 指數函數
GDIS(X,Y) 求以X為均值,Y為標準差的高斯分布,在使用蒙地卡羅法研究隨機荷載和隨機材料參數時,可以用該函數處理計算結果
LOG(X) 自然對數
LOG10(X) 常用對數(以10為基)
MOD(X,Y) 求 X/Y的余數. 如果 Y=0, 函數值為 0
NINT(X) 求最近的整數
RAND(X,Y) 取隨機數,其中X 是下限, Y是上限
SIGN(X,Y) 取 X的絕對值并賦予Y的符號. Y>=0, 函數值為|X|, Y<0, 函數值為-|X|,.
SIN(X) 正弦
SINH(X) 雙曲正弦
SQRT(X) 平方根
TAN(X) 正切
TANH(X) 雙曲正切
展開 ANSYS中不同形狀的波函數書寫方法
在ANSYS仿真中經常會遇到一些動態的加載方法，加載的載荷（位移、力、電流、溫度等）隨著時間而變化，表示不同的狀態。而相應的在workbench中可以方便的采用表格方法設置不同時間狀態下加載的位移或受力等載荷。但是又實用需要APDL命令的方式書寫不同時刻的載荷，但是函數庫當中又沒有相應的函數，那么如何書寫呢，下面我們選取幾個有代表性的書寫方法
（1）三角波的使用
一個物體在一個平面上移動，從左到右勻速運動，然后再從右到左的勻速運動。如果次數多，則采用表格方式實現比較麻煩，而采用do循環命令的方式就可以方便的加載。
該方法可以采用三角波的形式，給物體一個位移，使它不斷的左右移動，獲取摩擦熱或者應力應變等過程方法如下：
*do,x,1,25,1
cc=ABS(ASIN(ABS(SIN((x+PI()/2)/2)))-PI()/4)!
展開 如何在ANSYS WORKBENCH中施加分段函數激勵
本篇回答一位朋友提出來的問題，說明如何在ANSYS WOKRBENCH中施加分段函數激勵。
假設分段的分布載荷如下
該載荷施加在一長方體的頂面上，作為分布力系施加。
下面說明操作方法。
1. 創建一個瞬態動力學分析系統
2.創建一長方體，尺寸任意。
3.劃分網格
4.分析設置
設置兩個時間步，
第一步終止時間為1秒，打開自動時間步長，通過載荷步來定義載荷子步，初始子步10步，最小5步，最多20步。
再定義第二步如下
其含義是
第2步終止時間為2秒，打開自動時間步長，通過載荷步來定義載荷子步，初始子步10步，最小5步，最多20步.
5.固定左端
6.在上面施加分布載荷1
首先定義第一個載荷步內的函數載荷
接著休眠期第二段（1-2秒內的部分）
得到結果如下
7.在上面施加分布載荷2
接著休眠期第1段（0-1秒內的部分）
得到結果如下
這就可以了。
至于后面的求解就不再贅述了。
來源：宋博士的博客，版權歸作者所有。
展開 『分享』在ANSYS中如何根據函數方程畫曲線
*dim,a,,10
*dim,b,,10
*do,i,1,10
a(i)=i
b(i)=sin(i/5)
*enddo
/prep7
*do,i,1,10
k,i,a(i),b(i),0
*enddo
*do,i,1,9
l,i,i+1
ANSYS的lsdyan中螺栓預緊力Bolt Pretension加載
? 若為同一梁連接同時定義了 Dynamic Relaxation 文件夾中的螺栓預緊力和 LS - DYNA 瞬態分析下的螺栓預緊力，分析時僅使用最后定義的那個。
ANSYS中的循環載荷加載，最易理解的案例來了！
本文的主要目的就是展示在ANSYS中循環加載是如何實現的。
計算結果
橡膠塊循環拉伸變形結果（可以看到有四次循環變形）
本文以一個正方形橡膠塊為例說明，橡膠塊如圖約束（約束XY面節點Z自由度，約束XZ面節點Y自由度，約束YZ面節點X自由度）,在側面施加循環載荷。
計算模型示意圖
循環載荷施加正弦形狀的位移載荷，分為4個正弦周期，四個正弦周期載荷幅值分別為0.1，0.2，0.3，0.4，4個周期加載過后，橡膠內部積累的應力釋放。具體定義分為幾個步驟：
步驟一：首先定義4個周期載荷幅值向量。
*DIM,AMPL,ARRAY,4 ! Amplitude Vector Definition
AMPL(1)=0.01
AMPL(2)=0.02
AMPL(3)=0.03
AMPL(4)=0.04
步驟二：定義離散時間加載點
*DIM,SOLTIME,ARRAY,161 ! Time Vector Definition
SOLTIME(1)=0.0
*DO,I,2,161,1
SOLTIME(I)=SOLTIME(I-1)+0.1
*ENDDO
步驟三：計算每個時間點下的位移激勵大小,也就是正弦曲線上的y值大小。
*DIM,BC_X,ARRAY,161 !
展開 
