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