以文本方式查看主题

-  声学楼论坛  (http://874085.11480.vipsjym.com.my3w.com/bbs/index.asp)
--  基础理论室  (http://874085.11480.vipsjym.com.my3w.com/bbs/list.asp?boardid=17)
----  拉普拉斯元件实现声阻运算（pspice)  (http://874085.11480.vipsjym.com.my3w.com/bbs/dispbbs.asp?boardid=17&id=24434)
--  作者：liheng_rd
--  发布时间：2013-5-7 0:02:17
--  拉普拉斯元件实现声阻运算（pspice)

如题，

拉普拉斯元件实现声阻运算，参数如何设置
--  作者：nisa
--  发布时间：2013-5-7 11:29:59
--  
Modeling a frequency dependent inductor for AC sweep analysis

 

PROBLEM:
How do I model a frequency dependent inductor for AC sweep analysis?

SOLUTION:
The impedance of an inductor is:

Z = L * s

where

s = j * 2 * PI * Frequency

In a frequency dependent inductor, the "L" in the above formula changes with frequency.

If you know L as a math function of frequency, then:

  1. Use a GLAPLACE device.
  2. Connect the positive and negative inputs to the positive and negative outputs, respectively.
  3. Write the transform expression as 1/(<L>*s), where <L> is replaced by the desired math formula, with frequency replaced by abs(s)/6.283185.
For example, if L = 1 / sqrt(Frequency), then write the transform expression as:

1/(1 / sqrt(abs(s)/6.283185)*s) or sqrt(abs(s)/6.283185)/s



http://www.orcad.com/documents/community.faqs/pspice/020020.aspx

 

更详细的介绍查阅实例(要翻越和谐之墙): Google搜索"Modeling a frequency dependent inductor for AC sweep analysis"的第二个结果:

 http://acousticians.blogspot.com/2012/10/modeling-frequency-dependent-inductor.html
[此贴子已经被作者于2013-05-07 20:03:13编辑过]
--  作者：jma
--  发布时间：2013-5-7 20:09:03
--  
 

1) 将GLAPLACE的输入端正负极与输出端正负极分别对接起来;

2) 改写转换函数为1/<R>形式, 其中<R>为计算后得到的频率f的函数;

3) 考虑到拉普拉斯变量s = jω = j·2πf, 故将上述转换函数中的频率改写成f = abs(s)/6.28 .
--  作者：liheng_rd
--  发布时间：2013-5-7 22:30:39
--  

非常感谢楼上二位的帮助，谢谢
--  作者：liheng_rd
--  发布时间：2013-5-8 0:30:01
--  

举个实例吧：

 

单孔的声阻：R=sqrt(2ηρw)*(l/a+2)/(π*a^2)

 

假设:η=1/2;ρ=1/2;l=1；a=1;

 

则计算的R=[sqrt(π*frequency)*3]/π

 

那么XFORM的参数是1/(([sqrt(π*abs(s)/6.28)*3]/π)*s)这个表达式吗？
--  作者：jeno
--  发布时间：2013-5-9 9:37:09
--  
学习啦
--  作者：nisa
--  发布时间：2013-5-9 11:32:50
--  
以下是引用liheng_rd在2013-05-08 00:30:01的发言：

举个实例吧：

 

单孔的声阻：R=sqrt(2ηρw)*(l/a+2)/(π*a^2)

 

假设:η=1/2;ρ=1/2;l=1；a=1;

 

则计算的R=[sqrt(π*frequency)*3]/π

 

那么XFORM的参数是1/(([sqrt(π*abs(s)/6.28)*3]/π)*s)这个表达式吗？

要把分母最后的s去掉, 即XFORM= 1/Z =1/R= 1/([sqrt(pi*abs(s)/6.28)*3]/pi)

 

假设单孔的声质量表达式子L=sqrt(π*frequency)*3]/π ,

那么XFORM的参数才是XFORM= 1/Z =1/(L*s)=1/(([sqrt(π*abs(s)/6.28)*3]/π)*s)

 

1/Z中的Z表示阻抗
--  作者：lhyjs_63
--  发布时间：2013-6-28 15:33:50
--  [讨论]
小马哥，nisa您们好：
                     第3点中，3) 考虑到拉普拉斯变量s = jω = j·2πf, 故将上述转换函数中的频率改写成f = abs(s)/6.28 

             有不明白请指教，s = jω = j·2πf，那么 f = s /  j·2π, 又因为 j = √-1  所以 f = s / (√-1 ·2π) = -s·√-1 / 2π, 如何得到 f = abs(s)/6.28 的呢，谢谢。
             

                                                 
--  作者：nisa
--  发布时间：2013-6-29 7:54:44
--  
以下是引用lhyjs_63在2013-06-28 15:33:50的发言：
小马哥，nisa您们好：
                     第3点中，3) 考虑到拉普拉斯变量s = jω = j·2πf, 故将上述转换函数中的频率改写成f = abs(s)/6.28 

             有不明白请指教，s = jω = j·2πf，那么 f = s /  j·2π, 又因为 j = √-1  所以 f = s / (√-1 ·2π) = -s·√-1 / 2π, 如何得到 f = abs(s)/6.28 的呢，谢谢。
                                                             

你想多了;   abs（）是求模值，即abs(s)=abs(jω)=ω, ω=2πf → f=ω/2π;

--  作者：mustang
--  发布时间：2013-7-5 19:18:37
--  
学习
--  作者：jqr01
--  发布时间：2013-7-8 2:41:17
--  
阻抗Z= 单孔的声质量 L *S   这里的S是代表什么呢,为什么 要把分母最后的s去掉 , 即XFORM= 1/Z =1/R= 1/([sqrt(pi*abs(s)/6.28)*3]/pi)
--  作者：jqr01
--  发布时间：2013-7-8 2:42:29
--  
为什么这里又会出现S:那么XFORM的参数才是XFORM= 1/Z =1/(L*s)=1/(([sqrt(π*abs(s)/6.28)*3]/π)*s)

 

1/Z中的Z表示阻抗
--  作者：nisa
--  发布时间：2013-7-8 20:29:24
--  
 二、三楼的内容说得够详细了啊；自己多动脑子思考才好理解透。
--  作者：katlin
--  发布时间：2015-5-6 10:07:44
--  
为什么我将参数XFORM设置成 1/([sqrt(pi*abs(s)/6.28)*3]/pi)，提示出错呢，说是无效的传输函数？？
**** INCLUDING ceshi1-SCHEMATIC1.net ****
* source CESHI1
C_Cmd         N16341 N11893  4.5m  
E_pSd         N03344 0 V2 N06338 3.8969e-5
L_Ma1f         OUT 0  163  
X_vSd    N11893 N03344 N06338 V2 SCHEMATIC1_vSd 
R_Rmd         V1 N16144  4.0m  
C_C2         N06242 0  0.094p  
R_Re         N02561 N02595  31.2  
C_C1         N03778 0  0.093p  
X_BLi    N16272 N03050 V1 N03236 SCHEMATIC1_BLi 
R_Ra11         OUT N04192  10  
X_BLv    0 N03236 N03050 0 SCHEMATIC1_BLv 
L_Le         N02595 N16272  80uH  
R_Ra1f         0 OUT  190M  
L_Mmd         N16144 N16341  19uH  
R_Ra1r         N06242 N06138  246M  
G_G1         N03778 N03852 LAPLACE { V(N03778, N03852) } {
+  1/[sqrt(abs(s)/6.28)*0.09] }
---------------------------------------------------------$
ERROR -- Invalid transfer function
L_Maf         V2 N03778  38.6  
R_Ra31         N06138 N05374  10  
L_Mar         N06338 N06242  38.6  
G_G2         0 N05072 LAPLACE { V(0, N05072) } { 1/[sqrt(abs(s)/6.28)*0.09] }
-----------------------------------------------$
ERROR -- Invalid transfer function
R_R20         N03778 0  100G  
L_Ma1         N03852 N04192  210  
R_R21         N06242 0  100G  
L_Ma1r         N06138 N06242  186  
V_Eg         N02561 0 DC 0Vdc AC 1.79 
L_Ma3         N05072 N05374  455  
--  作者：Acoustics
--  发布时间：2015-5-6 15:25:57
--  
* Example Schematics Netlist *
R_R2         $N_0002 $N_0001 0.001
V_V3         $N_0002 0 DC 0V AC 1v
R_R3         $N_0002 0 10meg
R_Rsim     0 $N_0003 10meg
G_Rvar2    $N_0001 0 LAPLACE { V($N_0001, 0) } { 
+ 1/(k2*sqrt(S/(2*3.14159265))) }
R_R1          $N_0004 0 16000
G_Lvar       $N_0003 0 LAPLACE { V($N_0003, 0) } { 
+ 1/(S*1e-6*(k3-(k4*LOG(k5*(S/(2*3.14159265)))))) }
G_Rvar1     $N_0001 $N_0003 LAPLACE { V($N_0001, $N_0003) } { 
+ 1/(k1*sqrt(S/(2*3.14159265))) }
C_C            $N_0001 $N_0004 .64pF
To view the specific effective simulation results
Include the following Macros in your Probe trace analysis to see frequency vs. inductance, impedance, phase angle (in degrees), and Q factor:
PI = 3.14159265
L=(IMG(V(R3:1)/I(R2))/(2*pi*FREQUENCY)
Z = V(R3:1)/I(R2)
ANG = (180/PI)*ARCTAN((IMG(V(R3:1)/I(R2)))/(R(V(R3:1)/I(R2))))
QFACT = ABS((IMG(V(R3:1)/I(R2)))/(R(V(R3:1)/I(R2))))

由于 发布
--  作者：vickyzhang1981
--  发布时间：2017-1-20 16:48:52
--  
我的权限也不够，也来狂顶楼主一通图片点击可在新窗口打开查看
--  作者：XYO104249
--  发布时间：2017-12-25 14:56:47
--  
太深奥了~~收藏学习
--  作者：eie05031
--  发布时间：2023-5-17 13:12:40
--  
么的深奥啊，闷啦