ansys函數(shù)載荷的案例
Ansys Wrokbench分段復(fù)雜函數(shù)載荷，加載方式記錄 ￥10
問題：
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頁卡中選著變量time，在Regime頁卡中逐個(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命令顯示出來。（有時(shí)log file顯示不及時(shí)，再重復(fù)一次即可）
4. 在Workbench內(nèi)創(chuàng)建加載remote point點(diǎn)，并設(shè)定加載點(diǎn)的ADPL name為“LoadPoint“，用于加載。
展開 ansys載荷為位置函數(shù)的一個(gè)例子
設(shè)壓力分布函數(shù)為F=2RSIN(AF)COS(XT)
命令流
*do,i,1,N (共N個(gè)單元）
csys,2 （激活極坐標(biāo)系）
esel,,,i （選擇第i個(gè)單元）
xe=centrx(i) （xe為單元中心X坐標(biāo)）
ye=centry(i) （ye為單元中心y坐標(biāo)）
ze=centrz(i) （ze為單元中心z坐標(biāo)）
*GET,AA,ELEM,I,AREA
f=2*sin(xe)*cos(ye)*cos(ze)
sfe,i,1,pres,,f/AA,,,
*enddo
=====================
先選要加載荷的面，get其最大單元號(hào)和最小單元號(hào)。再做上面提到的循環(huán)，循環(huán)從最小單
元號(hào)到最大單元號(hào)。循環(huán)過程中，判斷所循環(huán)的單元號(hào)在所選的面內(nèi)是否存在。
命令流：
asel,... (要加載荷的面）
esla,,,, (和面相聯(lián)的單元）
*GET,eMAX,ELEM,0,NUM,MAX !獲取最大單元號(hào)
*GET,eMIN,ELEM,0,NUM,MIN !獲取最小單元號(hào)
*do,i,emin,emax
*if,esel(i),eq,1,then (該句語法請(qǐng)查幫助，這里憑記憶寫出）
！！！*do,i,1,N (共N個(gè)單元）
csys,2 （激活極坐標(biāo)系）
esel,,,i （選擇第i個(gè)單元）
xe=centrx(i) （xe為單元中心X坐標(biāo)）
ye=centry(i) （ye為單元中心y坐標(biāo)）
ze=centrz(i) （ze為單元中心z坐標(biāo)）
*GET,AA,ELEM,I,AREA
f=2*sin(xe)*cos(ye)*cos(ze)
sfe,i,1,pres,,f/AA,,,
*endif (判斷結(jié)束）
*enddo
展開 ANSYS知識(shí)普及4——如何施加函數(shù)變化的表面載荷 （ANSYS專家編輯，非原創(chuàng)，歡迎轉(zhuǎn)摘）
本人準(zhǔn)備出一個(gè)ANSYS知識(shí)普及系列，將有用的網(wǎng)上資料歸攏，由于知識(shí)水平有限，不對(duì)之處請(qǐng)諒解。也歡迎各位網(wǎng)友提供好的資料分享，讓我們共同完成這個(gè)ANSYS知識(shí)普及系列。
編輯人：技術(shù)鄰ANSYS專家
業(yè)務(wù)咨詢網(wǎng)址：http://www.yqgqt.org.cn/content/other/402981
（打個(gè)小廣告）
聲 明:1、ANSYS知識(shí)普及系列中所有資料均來自網(wǎng)上；
2、如侵犯知識(shí)產(chǎn)權(quán)，請(qǐng)聯(lián)系ANSYS專家本人或者技術(shù)鄰，我將第一時(shí)間刪除。
小技巧：加本人關(guān)注，可以及時(shí)觀看本人發(fā)布的技術(shù)貼
ANSYS具有函數(shù)加載功能，可以很方便地在模型表面施加函數(shù)變化的各種載荷，在ANSYS中，也可以通過變通的方式來實(shí)現(xiàn)此功能，其思路是：
首先選定所要施加函數(shù)變化表面載荷的表面上的節(jié)點(diǎn)，利用ANSYS的參數(shù)數(shù)組和嵌入函數(shù)知識(shí)寫一簡(jiǎn)單的命令流，定義好相應(yīng)節(jié)點(diǎn)位置的面載荷值，然后通過在節(jié)點(diǎn)上施加面載荷來完成。
下面以在一圓柱表面施加函數(shù)變化載荷為例：
/prep7
et,1,45
cyl4,,,0.5,,,,3
vsweep,all
asel,s,loc,y,0.01,1
nsla
!
*get,nmax,node,,num,max,
*get,nmin,node,,num,min,
*afun,deg
*dim,t1,array,nmax,1,1,
csys,1
*do,k,nmin,nmax
*if,nsel(k),eq,1,then
t1(k)=1000*sin(ny(k))
*else
t1(k)=0
*endif
*enddo
!
sffun,pres,t1(1)
sf,all,pres,0
展開 載荷為振幅的函數(shù)
我想求解一個(gè)自激振動(dòng)時(shí)域響應(yīng)的問題，可是要求其載荷是結(jié)構(gòu)振幅的函數(shù)，有沒有哪位大俠知道此類載荷如何實(shí)現(xiàn)，拜求指點(diǎn)、多謝了
從形函數(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ù)來理解節(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é)中通過平衡方程求解（首先求得的解是力解），有限元方式求解的特點(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)本來就是用函數(shù)描述的，所以自然就存在函數(shù)的性質(zhì)，所以用函數(shù)的性質(zhì)來理解就可以方便解釋一些現(xiàn)象了，下圖分別是用兩種形函數(shù)描述的位移場(chǎng)，在有限元求解后得到的首先是節(jié)點(diǎn)位移解，即圖中5個(gè)節(jié)點(diǎn)的位移，假如每個(gè)節(jié)點(diǎn)的位移用坐標(biāo)x\y\z的函數(shù)來表示，然后通過形函數(shù)插值得到相鄰節(jié)點(diǎn)之間的位移（也是xyz的函數(shù)），上圖是用一次形函數(shù)插值，下圖是用二次形函數(shù)插值。
展開 ansys Workbench螺栓載荷提取時(shí)，如何計(jì)算載荷偏心距離（VDI2230） ￥10
問題：
VDI2230關(guān)于螺栓的計(jì)算中對(duì)于螺栓載荷的提取沒有過多的涉及，本文針對(duì)偏心載荷的提取問題進(jìn)行簡(jiǎn)單說明。
VDI2230中，對(duì)于載荷偏心距a的定義如下，虛擬軸線到截面彎矩為0的點(diǎn)之間的距離。
對(duì)于實(shí)際螺栓連接問題，幾何結(jié)構(gòu)和載荷狀態(tài)復(fù)雜多變，使用經(jīng)驗(yàn)公式估計(jì)并不理想。本文介紹使用有限元仿真的方法確定載荷偏心距離。
示例：
以VDI2230中的案例5為例進(jìn)行對(duì)比計(jì)算，依據(jù)案例5的幾何信息創(chuàng)建仿真模型。
約束筒體底面，在內(nèi)表面施加20Mpa壓力載荷，同時(shí)給螺栓施加約150KN的預(yù)緊力（加不加結(jié)果變化不大），連接面設(shè)定為摩擦面。
將兩個(gè)側(cè)面設(shè)定為，frictionless Support，等效對(duì)稱邊界。（這里沒有使用圓周循環(huán)對(duì)稱邊界，是因?yàn)閳A周對(duì)稱邊界不能支持截面彎矩提取）
注意，在輸出控制中 打開“Nodal Forces”，用于端蓋截面的彎矩提取。
計(jì)算完成后，在結(jié)果提取中，插入Probe——Moment Reaction——使用surface類型進(jìn)行端蓋截面彎矩載荷的提取，這里只需要關(guān)注X軸彎矩。
依次變更截面位置，就可以獲得一條彎矩隨位置變化的曲線，讀取彎矩為0位置的距離值，再進(jìn)一步處理加上螺栓偏心距Ssym，就可以換算到載荷偏心距a。
個(gè)人認(rèn)為仿真結(jié)果17.535，除了在循環(huán)對(duì)稱設(shè)置上與案例給出條件不同外，其余均能反應(yīng)案例邊界。
補(bǔ)充案例：
以機(jī)械設(shè)計(jì)手冊(cè)兩端固支梁，在均布載荷下的反彎點(diǎn)計(jì)算模型為例進(jìn)行驗(yàn)證。
仿真結(jié)果
公式計(jì)算值42.2mm，仿真結(jié)果42.23mm。
展開 ansys的取值函數(shù)
當(dāng)多個(gè)單元共享同一個(gè)節(jié)點(diǎn)時(shí)，在結(jié)果列表中，LOC表示為位置，并在列表的結(jié)束處返回一個(gè)零
NELEM(E,NPOS)
單元E中，NPOS位置處的節(jié)點(diǎn)編號(hào)
與單元面相關(guān)的取值函數(shù)
ELADJ(E,FACE)
與單元E的第FACE面相連接的單元，單元面號(hào)與表面載荷的主常數(shù)是一致的，只有具有相同維數(shù)和形狀的單元才被考慮。如果有多個(gè)單元時(shí)，將返回A-1,如果沒有相鄰近的單元?jiǎng)t返回為A-0
NDFACE(E,FACE,LOC)
單元E中第FACE面中的第LOC位置處的節(jié)點(diǎn)，單元面號(hào)與表面載荷的主常數(shù)是一致的，LOC是面上的節(jié)點(diǎn)位置（對(duì)IJLK面上，LOC=1是在節(jié)點(diǎn)I處，LOC=2是在節(jié)點(diǎn)J處）
NMFACE(E)
是指包含所選擇節(jié)點(diǎn)單元E的面號(hào)，輸出的面號(hào)就是表面載荷的主常數(shù)，如果某個(gè)面上（如線和面單元）存在有多個(gè)載荷常數(shù)，則輸出其最新的主常數(shù)號(hào)
ARFACE(E)
對(duì)于二維平面實(shí)體和三維體積實(shí)體來說，返回容納所選節(jié)點(diǎn)的單元E的面的面積，對(duì)于軸對(duì)稱單元，面積是一個(gè)360的環(huán)單元。
展開 ANSYS 內(nèi)部函數(shù)
VLNEXT(N)
Next higher volume number above N in selected set (or zero if none
found).
30.距離函數(shù)
DISTND(N1,N2) Distance between nodes N1 and N2.
DISTKP(K1,K2) Distance between keypoints K1 and K2.
DISTEN(E,N) Distance between the centroid of element E and node N.
Centroid is determined from the selected nodes on the element.
31.角度函數(shù) (缺省單位為弧度，單位變換用 *AFUN 命令)
ANGLEN(N1,N2,N3) Subtended angle between two lines (defined by
three nodes where N1 is the vertex node). Default is in
radians.
ANGLEK(K1,K2,K3) Subtended angle between two lines (defined by
three keypoints where K1 is the vertex keypoint). Default is in
radians.
32.最近實(shí)體函數(shù)
NNEAR(N) Selected node nearest node N.
KNEAR(K) Selected keypoint nearest keypoint K.
ENEARN(N) Selected element nearest node N.
展開 Ansys Zemax | 什么是點(diǎn)擴(kuò)散函數(shù)（ PSF ）
附件下載
聯(lián)系工作人員獲取附件
本文討論了如何在 OpticStudio 中對(duì)點(diǎn)擴(kuò)散函數(shù)進(jìn)行建模和解釋。使用的分析特征是 Spot Diagram、FFT PSF 和 Huygens PSF。將討論每種工具的優(yōu)點(diǎn)，以及用于最準(zhǔn)確分析的有用特征設(shè)置。
介紹
光學(xué)系統(tǒng)的點(diǎn)擴(kuò)散函數(shù) （PSF） 是單個(gè)點(diǎn)光源產(chǎn)生的輻照度分布。（望遠(yuǎn)鏡拍攝遙遠(yuǎn)恒星的圖像就是一個(gè)很好的例子。盡管源可能是一個(gè)點(diǎn)，但圖像不是。有兩個(gè)主要原因：首先系統(tǒng)中的像差會(huì)將圖像傳播到有限的區(qū)域；其次衍射效果也會(huì)擴(kuò)散圖像，即使在沒有像差的系統(tǒng)中也是如此。
OpticStudio 有三種基本類型的 PSF 計(jì)算：幾何（無衍射）點(diǎn)列圖、基于衍射的 FFT 和 Huygens PSF。本文將討論基本理論，并就正確使用每種類型的 PSF 提供一些指導(dǎo)。
點(diǎn)列圖
OpticStudio 中最基本的分析功能之一是點(diǎn)列圖。此功能從物空間中的單視場(chǎng)點(diǎn)發(fā)射許多光線，通過光學(xué)系統(tǒng)追跡所有光線，并繪制所有光線相對(duì)于某個(gè)公共參考的 （x，y） 坐標(biāo)。因此，點(diǎn)列圖本身就可以看作一個(gè)幾何 PSF。
這里使用的示例光學(xué)系統(tǒng)是一個(gè)焦距為 50 mm 的單拋物面 F/5 反射鏡，物位于無窮遠(yuǎn)處。該系統(tǒng)是一個(gè)簡(jiǎn)化的牛頓望遠(yuǎn)鏡，包含的示例文件為 PSF_Newtonian.ZMX。以下是光學(xué)系統(tǒng)的外觀：
兩個(gè)視場(chǎng)點(diǎn)（一個(gè)在軸上，另一個(gè)呈 2 度角）的點(diǎn)列圖如下所示。
請(qǐng)注意，點(diǎn)列圖是光線落點(diǎn)的集合，每個(gè)點(diǎn)表示一條光線。光線之間沒有相互作用或干擾。點(diǎn)列圖在顯示望遠(yuǎn)鏡的幾何或光線像差的影響方面非常有效。離軸幾何 PSF 清楚地顯示了系統(tǒng)的彗差和像散。然而在軸上，點(diǎn)列圖預(yù)測(cè)了完美的成像。但這是否準(zhǔn)確代表了光學(xué)系統(tǒng)的性能？為了回答點(diǎn)列圖結(jié)果的這個(gè)問題，我們需要將點(diǎn)列分布與衍射極限響應(yīng)進(jìn)行比較。
展開 ANSYS 支持的函數(shù)列表
ANSYS 支持的函數(shù)列表，備用與共享，以后不要老再去找了
SIN(X) Sine
COS(X) Cosine
TAN(X) Tangent
ASIN(X) Arcsine
ACOS(X) Arccosine
ATAN(X) Arctangent
ATAN2(Y,X) Arctangent (Y/X) with the sign of each component considered
SINH(X) Hyperbolic sine
COSH(X) Hyperbolic cosine
TANH(X) Hyperbolic tangent
SQRT(X) Square root
ABS(X) Absolute value
SIGN(X,Y) Absolute value of X with sign of Y.
展開 ANSYS的“get”函數(shù)列表
自己收藏并與大家分享，來自于ＡＮＳＹＳ的ｈｅｌｐ
“get函數(shù)”可用于某些項(xiàng)，并可用于代替*get命令。函數(shù)返回值并在函數(shù)被輸入的地方使用它，繞過了用參數(shù)名存儲(chǔ)值和在要使用值的地方輸入?yún)?shù)名的需要。
例如，假設(shè)要計(jì)算兩個(gè)節(jié)點(diǎn)的平均X位置。使用*GET命令，參數(shù)L1可以指定節(jié)點(diǎn)1的X位置，參數(shù)L2可以指定節(jié)點(diǎn)2的X位置。然后，可以從mid=（L1+L2）/2計(jì)算mid位置：
*GET,L1,NODE,1,LOC,X
*GET,L2,NODE,2,LOC,X
MID=(L1+L2)/2
但是，使用返回節(jié)點(diǎn)N的X位置的節(jié)點(diǎn)位置“get ”函數(shù)NX（N），可以直接計(jì)算MID，而不需要中間參數(shù)L1和L2：
MID=(NX(1)+NX(2))/2
除非另有說明，否則Get函數(shù)返回活動(dòng)坐標(biāo)系中的值。
Get函數(shù)參數(shù)本身可能是參數(shù)或其他Get函數(shù)。get函數(shù)NELEM（E，NPOS）返回元素編號(hào)E的NPOS位置的節(jié)點(diǎn)號(hào)。組合函數(shù)NX（NELEM（E，NPOS））返回該節(jié)點(diǎn)的X位置。
下表列出了按功能分組的可用get函數(shù)。*GET命令還列出GET函數(shù)作為*GET items的替代項(xiàng)（如果適用）
Table 1: *GET - Get Function Summary
"Get Function" Summary
Entity Status Get Function Description
NSEL(N) Status of node N: -1=unselected, 0=undefined, 1=selected.
ESEL(E) Status of element E: -1=unselected, 0=undefined, 1=selected.
展開 
ANSYS 中添加窗函數(shù)
比如說我要施加一個(gè)周期性的沖擊 希望他作用2個(gè)周期后停止
如何利用ANSYS的隨機(jī)分布函數(shù)功能
作者：水哥ANSYS
來源：本文源于ANSYS結(jié)構(gòu)院，上海安世亞太授權(quán)轉(zhuǎn)載
隨機(jī)分布在材料微觀力學(xué)分析中扮演著重要角色，例如混凝土骨料力學(xué)、新型材料纖維力學(xué)分析等內(nèi)容，提及隨機(jī)分布，更多的同學(xué)可能會(huì)聯(lián)想到采用第三方軟件如Matlab來生成，并導(dǎo)入ANSYS計(jì)算，其實(shí)ANSYS本身自帶隨機(jī)分布功能，只是功能略有限制。
ANSYS中產(chǎn)生隨機(jī)分布的一個(gè)重要函數(shù)是 *VFILL，該函數(shù)主要的作用是對(duì)數(shù)組進(jìn)行填充賦值，而在賦值的過程中，用戶既可以選擇自定義數(shù)據(jù)內(nèi)容，也可以選擇利用隨機(jī)函數(shù)產(chǎn)生數(shù)值，ANSYS Help中*VFILL說明如下：
該函數(shù)主要輸入?yún)?shù)為數(shù)組名稱以及輸入數(shù)據(jù)的函數(shù)，當(dāng)選擇為data時(shí)，表示用戶自定義數(shù)據(jù)進(jìn)行填充，當(dāng)選擇其他選項(xiàng)時(shí)，則根據(jù)函數(shù)類型進(jìn)行填充。
*VFILL隨機(jī)數(shù)生成支持均勻分布(Rand)、高斯分布(GDIS)、三角分布(TRIA)、貝塔分布（BETA）、伽馬分布（GRMM），*VFILL用于批量生成，如果需要單獨(dú)生成數(shù)據(jù)，則可以分別使用函數(shù)：
1) num=Rand(con1,con2)
2) num=Gdis(con1,con2)
3) num=Tria(con1,con2,con3)
4) num=Beta(con1,con2,con3,con4)5) num=Gram(con1,con2,con3)
上述con1~con4分別表示函數(shù)參數(shù)，例如針對(duì)均布分布，con1和con2分別表示分布的下限和上限。
下面分別以均布分布、高斯分布、伽馬分布為例進(jìn)行演示。
1、均布分布
APDL代碼：
finish
/clear
/prep7
numA=1000
!
展開 手把手教你ANSYS的函數(shù)加載
最近論壇里很多人問我，如何施加函數(shù)載荷。
我今天給大家奉獻(xiàn)的是任意函數(shù)加載的操作步驟詳解，手把手教大家操作自定義的函數(shù)加載。
如果覺得還不錯(cuò)，頂下帖子，也算對(duì)我的鼓勵(lì)了！
大家有什么ANSYS 或 Workbench Mechanical 相關(guān)的問題，可以隨時(shí)**我 ansys123@qq.com
手把手教你ansys函數(shù)加載.doc
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) 雙曲正切
展開 