函数中的嵌套向量应用于R中的数据帧

问题描述

我有以下基本数据（问题末尾提供了dput）：

P <- 72
E_c <- 80000
head (deflections)
    L   I_comp
1  60 17299.21
2  70 25760.94
3  80 32734.69
4  90 51343.59
5 100 60167.30
6 110 64887.87

我想使用以下公式找到最大挠度：

其中b被定义为b = L - a，而a是可以定义为向量的矢量

a <- seq(10,L,10) 

和x是a这样的嵌套序列，

x <- seq(1,a,1)

我已经编写了下面的函数，以计算所有x值的所有a值的最大偏转量。

delta_less <- function(x){
  a_seq <- seq(10,10)
  L <- x[1]
  I_comp <- x[2]
  X_seq <- seq(a_seq)
  delta <- ((P*(L-a_seq)*X_seq)/(6*E_c*I_comp*L))*(L^2-(L-a_seq)^2-X_seq^2)
  c(val = max(delta),pos = which.max(delta))
}

test <- cbind(deflections,vars = t(apply(deflections,1,delta_less)))

问题是我的代码无法识别嵌套序列。仅当我给a_seq一个奇异的数值时，它才起作用。

提供的数据：

dput(deflections)
structure(list(L = c(60,70,80,90,100,110,120,130,140,60,140),I_comp = c(17299.2063000114,25760.9420686703,32734.6949444858,51343.5889683982,60167.3001695375,64887.8729936444,83451.9473750139,103000.852590734,112631.744898846,23539.3443398425,23544.8558575609,32808.739059987,40880.4964809276,60171.1172692296,64894.6680546358,83469.4777437875,103114.132264558,120223.534252571,25772.7904955165,32810.8231970971,51345.5936366056,74673.7079705815,80752.3694583712,147386.853127916,25767.8092758621,40881.3376639665,55608.9342154021,80694.6568665257,119068.355471205,119402.817753462,147386.853127916)),class = "data.frame",row.names = c(NA,-36L))

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）