功能:一维弹性波动态传输求解

格式:DynamicWave1D(x,t,Model,f,WaveSpeed,DirectionIsLeftToRight,f_x,PeriodTime,StartTime,EndTime)

x:数据变量,要求解点的坐标
t:矩阵变量,表示要求解的对应时刻
Model:1×5的矩阵变量,表示模型参数,其里面参数为 [模型左端点坐标,模型右端点坐标,模型左端点为自由面否,模型右端点为自由面否,模型类型]。具体意思看【注意】事项
f:符号变量存储的符号表达式,表示应力加载函数.此函数里的变量名称必须为x
WaveSpeed:矩阵变量,每个元素对应f里的表达式,表示传输速度
DirectionIsLeftToRight:矩阵变量,每个元素对应f里的表达式，其表示传输方向。为0表示由右到左传输，否则左到右传输
f_x:矩阵变量,每个元素对应f里的表达式，其表示加载函数作用点
PeriodTime:矩阵变量,每个元素对应f里的表达式，其表示加载函数周期
StartTime:矩阵变量,每个元素对应f里的表达式，其表示加载函数开始作用时间
EndTime:矩阵变量,每个元素对应f里的表达式，其表示加载函数结束作用时间

函数执行成功返回x点在t时刻对应的应力值.

原理:
(1)波在自由面反射后波的性质改变，在刚体固定端反射后波的性质不改变，因为是一维问题，则可以直接进行求解，而无需采用近视。
(2)因为应力波的可叠加性，所以本函数支持任意波形、任意加载点，任意个波的作用

注意:
1、本函数求解的是有限、半无限、无限长杆一维波的波动问题。在计算时没有考虑材料破坏，其实一维问题的材料破坏程序也是可以写的，但是目前没有什么时间了。
2、参数Model里的5个参数解释如下
(1)如果你的模型为无限长杆问题,则【模型类型】值为2，其它参数设置任意值即可
(2)如果你的模型为半无限长杆问题，如果是向右无限延伸，则【模型类型】值为1，需要设置【模型左端点坐标】的具体数值，如果左端面为自由面，【模型左端点为自由面否】此值设置为0，否则为固壁边界，则【模型左端点为自由面否】设置为1，其它值可随便设置
(3)如果你的模型为半无限长杆问题，如果是向左无限延伸，则【模型类型】值为-1，需要设置【模型右端点坐标】的具体数值，如果右端面为自由面，【模型右端点为自由面否】此值设置为0，否则为固壁边界，则【模型右端点为自由面否】设置为1，其它值可随便设置
(4)如果你的模型为有限长杆,,则【模型类型】值为0，【模型左端点坐标】和【模型右端点坐标】根据模型设置具体数值，对于【模型右端点为自由面否】和【模型左端点为自由面否】根据如果是自由面则设置为0，否则固壁面设置为1的原则进行设置。

例子:
//1、如下一个长40m的柱形杆，从0时刻开始从左端向右加载一个应力函数y(t)=10*sin(t/1*pi),持续加载2ms,其中y(t)周期为2ms,此波在杆中传播速度为10m/ms,求在距离左端20m处质点在0到15ms时刻内的应力变化曲线图

m=[0 40 0 0 0];//第1个0表示左端坐标,40表示右端面坐标,第2个0表示左端面为自由面,第3个0表示右端面为自由面,第5个0表示此类型为有限长杆模型
t=0:0.075:15;//初始化时刻,表示需要求解0-15ms内的值
f="10*sin(x/1*pi)";//注意函数里的变量必须为x
y=DynamicWave1D(20,t,m,f,10,1,0,2,0,2);
Plot(t,y)//依次回车执行上面几条命令后得到如下结果

{<http://img1.ph.126.net/-MN0msfJKG9sWV_I6m463g==/6619272204956346625.png>}
