如何创建类似于编程语言中的if-elseif-else这种结构的field

hylovegj

如题，请问如何创建类似于编程语言中的if-elseif-else这种结构的标量型空间场（比如压力），需要用到PCL吗？具体如何实现？

ygs2008

你说的我也正在研究，需要用到PCL。如果你的单元和载荷都是自己文件中提取的，用PCL建立FEM场就比较方便。手动建场就比较麻烦。

hylovegj

经过高人指点，终于解决了这个问题，下面以一个实例进行说明。
假设要施加如下压力
f(d)=( (540.1786*(1/d)**k0)**2 + （72.8036*(1/d)**kc)**2 )**0.5/1100
其中，k0和kc的取值为
当k0<=2时，k0=1.8+0.4*d；当k0>时,k0=2.8
当kc<1时，kc=1.7-0.7*d，但是当kc<0.2时，*(1/d)**kc=5；当kc>=1时，kc=1
方法为先创建一个pcl文件，例如pressure_fluctuaction_P.pcl，里面的pressure_fluctuaction_P函数的代码如下

1. 
   
2. $用途：创建左舷浆引起的脉动压力场时需要的PCL函数
   
3. $版本： V1.0
   
4. $开发时间： 20110910
   
5. function pressure_fluctuation_P(x,y,z)
   
6.    real x,y,z,d
   
7.    real k0,kc
   
8.    real pressure
   
9.    
   
10.    d = mth_sqrt((x+300)**2 + (y-6300)**2 + (z-2590)**2)/1100 /*对应霍尔顿法中的d/R*/
    
11.    
    
12.    $计算系数k0
    
13.    if(d <= 2)then
    
14.       k0 = 1.8 + 0.4*d
    
15.    else
    
16.       k0 = 2.8
    
17.    end if
    
18.    
    
19.    $计算系数kc
    
20.    if(d < 1)then
    
21.        kc = 1.7 - 0.7*d
    
22.    else
    
23.       kc = 1.0
    
24.    end if
    
25.    
    
26.    $用简易型霍尔顿法计算脉动压力叶频分量单幅值
    
27.    if(d >= 0.2)then
    
28.       pressure = mth_sqrt((540.1786*(1/d)**k0)**2 + (72.8036*(1/d)**kc)**2)
    
29.    else
    
30.       pressure = mth_sqrt((540.1786*(1/d)**k0)**2 + 132509.1942)
    
31.    end if
    
32.    return pressure
    
33. end function

复制代码

编译该函数，在Field面板的Scalar Function文本框中输入pressure_fluctuaction_P('x,'y,'z)，这样一个标量型空间场就创建好了，如下图所示

这个例子复杂了一些，但仔细看一下不难理解，希望可以帮到有类似需要的朋友，呵呵！

lenovosk3

hylovegj 发表于 2011-9-10 17:05
经过高人指点，终于解决了这个问题，下面以一个实例进行说明。
假设要施加如下压力
f(d)=( (540.1786*(1/d) ...

问下，你这个回复中一行一行的代码效果是怎么弄的？

hylovegj

lenovosk3 发表于 2014-1-8 11:18
问下，你这个回复中一行一行的代码效果是怎么弄的？

看到回复文本框上方的一对尖括号了没？单击后可以在弹出的文本框中粘贴代码。

lenovosk3

hylovegj 发表于 2014-1-8 18:23
看到回复文本框上方的一对尖括号了没？单击后可以在弹出的文本框中粘贴代码。 ...

OK, thanks！

言西早幻城

楼主，*.pcl文件，建立一个文本文件，改后缀就可以了吗？ 建好的*.pcl放在和db同一目录下就可以成功读进去？ 本人pcl小白，望指点！ 祝您一切顺利！

hylovegj

言西早幻城 发表于 2014-1-9 22:42
楼主，*.pcl文件，建立一个文本文件，改后缀就可以了吗？ 建好的*.pcl放在和db同一目录下就可以成功读进去 ...

是的，扩展名改为.pcl，放在db文件所在目录即可。

言西早幻城

hylovegj 发表于 2014-1-12 09:01
是的，扩展名改为.pcl，放在db文件所在目录即可。

谢谢lz。我严格按照您描述的尝试了一下，Field建立成功，不过给模型施加压力的时候报错了“$# (PCL) Missing right parenthesis”，因为没接触过pcl。lz可否把您的*.pcl文件发给我，我直接用，可能复制再粘贴出现了问题。
祝一切顺利！

hylovegj

言西早幻城 发表于 2014-1-13 08:34
谢谢lz。我严格按照您描述的尝试了一下，Field建立成功，不过给模型施加压力的时候报错了“$# (PCL) Missi ...

不好意思，代码文件一时半会找不到。
你先编译一下，看下是不是代码有问题。方法为：
在Patran命令行执行!!INPUT *.pcl命令，看下History Window里是否有错误提示。
如果有错误，它会告诉你错在哪里，你根据提示修改就可以了。
正确的代码文件成功编译后会显示...compiled这样的信息。
根据你贴上来的信息，出错原因是圆括号不匹配，右括号)丢失了，你仔细检查一下。
祝你好运！

言西早幻城

hylovegj 发表于 2014-1-13 18:30
不好意思，代码文件一时半会找不到。
你先编译一下，看下是不是代码有问题。方法为：
在Patran命令行执行 ...

谢谢lz的及时回复。
我又尝试了一下，而且根据您的代码简化了一个test.pcl进行测试，依然没成功。请您在方便的时候帮忙查看一下附件里的代码。祝您一切顺利！

hylovegj

我用你上传的test.pcl试了一下，一切ok呀！脉动压力施加成功。
详细说下你的操作步骤，应该是你哪一步出问题了。