CN111210880B - 基于非线性力电耦合生长模型的软物质生长自组装方法 - Google Patents

基于非线性力电耦合生长模型的软物质生长自组装方法 Download PDF

Info

Publication number
CN111210880B
CN111210880B CN201911400313.5A CN201911400313A CN111210880B CN 111210880 B CN111210880 B CN 111210880B CN 201911400313 A CN201911400313 A CN 201911400313A CN 111210880 B CN111210880 B CN 111210880B
Authority
CN
China
Prior art keywords
growth
incremental
electric field
stress
displacement
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911400313.5A
Other languages
English (en)
Other versions
CN111210880A (zh
Inventor
吕朝锋
杜洋坤
陈伟球
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201911400313.5A priority Critical patent/CN111210880B/zh
Publication of CN111210880A publication Critical patent/CN111210880A/zh
Application granted granted Critical
Publication of CN111210880B publication Critical patent/CN111210880B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/70—Machine learning, data mining or chemometrics

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Medical Informatics (AREA)
  • Evolutionary Computation (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Prostheses (AREA)

Abstract

本发明公开一种基于非线性力电耦合生长模型的软物质生长自组装方法，该方法可以实现通过力电耦合效应来控制生长软组织的自组装方法，所述方法包括将生长前后的软组织的几何特征模型化，进而得到生长前后的残余应力和电场强度；对于生长前后的残余应力和电场强度进行稳定性分析，获得满足增量平衡方程的模态，进而可以实现通过外部电场和机械力场实现指导调控生长自组装。本方法可以被广泛应用在软体机器人设计制造，生物医学研究，以及临床假体移植等科学及应用领域。

Description

基于非线性力电耦合生长模型的软物质生长自组装方法
技术领域
本发明涉及软组织自组装领域，具体涉及一种基于非线性力电耦合生长模型的软物质生长自组装方法。
背景技术
生物体的生长和重塑最基本和必要的生物活动之一，不仅赋予了生物体的多样性而且确保了很多生物体的正常的生物功能和自修复过程。广为人知的是生物遗传信息，化学刺激以及营养物质等因素控制着生物体的生长和修复过程。不过也有很多实验证实了外部分力电等物理场因素也能够影响生物体的生长发育过程。进一步的实验也证实了生物材料例如骨骼，牙齿，皮肤，血管，胶原蛋白，弹性蛋白等都具有力电耦合效应。
现阶段基于乘法分解的体生长理论，已有通过调控生物材料的初始材料的弹性参数和几何尺寸的方式来控制生长的方法提出。但是此方法的实际操作性不高，对于已有的生物组织，很难实现改变随意组织的厚度和材料参数。
发明内容
针对现有技术的不足，本发明提出一种基于非线性力电耦合生长模型的软物质生长自组装方法，具体包括如下步骤：
S1：将生长前后的软组织的几何特征模型化，进而得到生长前后的残余应力和电场强度；
S2：对于生长前后的残余应力和电场强度进行稳定性分析，获得满足增量平衡方程的模态，进而根据模态通过力电耦合效应控制生物组织的生长以及自组装。
进一步地，所述的S1具体为：
S1.1：将生长前后的软组织的几何特征模型化，进而得到生长前后的总体变形梯度F的具体形式；
S1.2：假设软组织为各向同性生长，即软组织的生长变形梯度Fg＝diag[g,g,g]，其中，g3为软组织的体积变化率；进而利用体生长模型F＝F+Fg，得到材料弹性变形梯度F+；
S1.3：对于不可压缩的各向同性材料，假设其应变能函数为Ω(Fe,g-3FgDl)，由
Figure BDA0002347318230000011
得到生长之后的应力表达式
Figure BDA0002347318230000012
其中，Dl是相对于软组织生长之前的构型下的电位移张量，
Figure BDA0002347318230000021
为弹性变形的左柯西格林张量，D为相对于软组织生长之后的构型的电位移张量，Ωi为应变能函数对第i张量不变量的偏导数，p为拉格朗日乘子，I1为右柯西格林张量
Figure BDA0002347318230000022
的第一张量不变量；
S1.4：由
Figure BDA0002347318230000023
得到
Figure BDA0002347318230000024
这里E是电场强度张量，对于理想介电材料，进一步得到对应的电压值为
Figure BDA0002347318230000025
S1.5：建立应力和电场强度的平衡方程divσ＝0,divD＝0,curlE＝0，计算出软组织生长之后的几何尺寸、残余应力和电场强度的值。
进一步地，所述的S2具体如下：
S2.1：通过线性增量理论分析，在生长后的构型上施加一个微小的增量位移场
Figure BDA0002347318230000026
和增量电位移
Figure BDA0002347318230000027
求偏导数得到增量位移梯度
Figure BDA0002347318230000028
进而代入
Figure BDA0002347318230000029
和
Figure BDA00023473182300000210
Figure BDA00023473182300000211
得到相对于生长后构型的增量应力
Figure BDA00023473182300000212
和增量电场强度
Figure BDA00023473182300000213
其中，
Figure BDA00023473182300000214
是对应于生长后构型中的物质点x摄动的位移；
Figure BDA00023473182300000215
为增量的拉格朗日乘子，
Figure BDA00023473182300000216
为增量电位移的相对于变形后构型的分量，
Figure BDA00023473182300000217
ΓI，KI为相对于生长后构型的增量即时模量；
S2.2：将增量应力的具体表达式代入到增量应力和电场的平衡方程
Figure BDA00023473182300000218
和边界条件
Figure BDA00023473182300000219
得到增量场的控制方程η′＝Gη，这里η是增量的位移应力以及电位仪向量，G是取决于生长和电场的Stroh矩阵；
S2.3：通过求解一个数学上存在的增量位移应力向量能够恰好满足控制方程η′＝Gη，即得到软物质生长之后失稳的临界时刻以及相应的模态，从而实现通过力电耦合效应控制生物组织的生长以及自组装。
本发明的有益效果如下：
相对于传统自组装方法，通过改变材料弹性参数以及几何尺寸调控软物质的生长模态，本发明可以通过外电场和机械场影响生长过程中的残余应力的积聚以及外电场强度的变化来改变软组织的生长过程状态的稳定性，进而实现通过外部作用影响调控软组织的自组装模态，有着更好的可操作性和实践性。特别是对于生物组织，几乎无法实现通过改变材料参数和几何尺寸调控软物质的生长模态。同时本发明也可以很好地被利用在软体机器人领域，通过设定外电场的加载模式，产生不同的生长失稳模态，进而实现不同的或者自适应的模态。同时本方法也可以被用于解释生物体力电对生长过程的共同影响作用。
附图说明
图1表示的是外力和外电场作用的体生长模型的原理图，其中，总体生长变形梯度分解为生长变形梯度和弹性变形梯度；
图2表示的是本发明的方法的流程图，首先给定施加的外力和电场值，其次通过生长理论得到生长之后的残余应力，进而对生长之后状态的稳定性，最后找出临界失稳模态也即是自组装的目标模态。
具体实施方式
下面根据附图和优选实施例详细描述本发明，本发明的目的和效果将变得更加明白，应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
如图1-2所示，本发明的基于非线性力电耦合生长模型的软物质生长自组装方法，首先，将生长前后的软组织的几何特征模型化，得到生长前后的总体变形梯度F的具体形式，进而假设软组织为各向同性生长，即软组织的生长变形梯度Fg＝diag[g,g,g]，其中，g3为软组织的体积变化率；进而利用体生长模型F＝F+Fg，得到材料弹性变形梯度F+；
对于不可压缩的各向同性材料，基于目标材料选择其合适的应变能函数为Ω(Fe,g-3FgDl)，其中Dl是相对于软组织生长之前的构型下的电位移张量，进而，由
Figure BDA0002347318230000031
Figure BDA0002347318230000032
可以得到生长之后的应力张量的表达式
Figure BDA0002347318230000033
Figure BDA0002347318230000034
其中
Figure BDA0002347318230000035
为弹性变形的左柯西格林张量，D为相对于软组织生长之后的构型的电位移张量，Ωi为应变能函数对第i张量不变量的偏导数，p为拉格朗日乘子，I1为右柯西格林张量
Figure BDA0002347318230000036
的第一张量不变量，I1＝trCe；
由
Figure BDA0002347318230000037
得到
Figure BDA0002347318230000038
这里E是电场强度张量，对于理想介电材料，进一步得到对应的电压值为
Figure BDA0002347318230000039
最后将应力和电场表达式代入应力和电场强度的平衡方程divσ＝0,divD＝0,curlE＝0，计算出来生长之后的几何尺寸和残余应力以及电场强度分量的值；
基于求解的生长后的几何尺寸，残余应力和电场强度，进而对生长后构型进行线性增量理论分析，在生长后的构型上施加一个微小的增量位移场
Figure BDA00023473182300000310
和增量电位移
Figure BDA00023473182300000311
这里
Figure BDA00023473182300000312
是对应于生长后构型中的物质点x摄动的位移；求偏导数即可以得到增量位移梯度
Figure BDA00023473182300000313
进而代入
Figure BDA00023473182300000314
和
Figure BDA00023473182300000315
可以得到相对于生长后构型的增量应力
Figure BDA00023473182300000316
和增量电场强度
Figure BDA00023473182300000317
其中
Figure BDA00023473182300000318
为增量的拉格朗日乘子，
Figure BDA00023473182300000319
ΓI，KI为相对于生长后构型的增量即时模量，其具体表达式为
Figure BDA0002347318230000041
Figure BDA0002347318230000042
依赖于生长过程中的弹性变形梯度和电位移分量可以得到增量方程中的材料即时模量
Figure BDA0002347318230000043
ΓI，KI；
将增量应力的具体表达式代入到增量应力和电场的平衡方程
Figure BDA0002347318230000044
和边界条件
Figure BDA0002347318230000045
得到增量场的控制方程η′＝Gη，这里η是增量的位移应力以及电位仪向量，G是取决于生长和电场的Stroh矩阵；
通过求解一个数学上存在的增量位移应力向量能够恰好满足控制方程η′＝Gη，即可得到软物质生长之后失稳的临界时刻以及相应的模态，也就可以实现通过力电耦合效应控制生物组织的生长以及自组装。
本领域普通技术人员可以理解，以上所述仅为发明的优选实例而已，并不用于限制发明，尽管参照前述实例对发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实例记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在发明的精神和原则之内，所做的修改、等同替换等均应包含在发明的保护范围之内。

Claims (1)

1.一种基于非线性力电耦合生长模型的软物质生长自组装方法，其特征在于，
S1：将生长前后的软组织的几何特征模型化，进而得到生长前后的残余应力和电场强度；所述的S1具体为：
S1.1：将生长前后的软组织的几何特征模型化，进而得到生长前后的总体变形梯度F的具体形式；
S1.2：软组织为各向同性生长，即软组织的生长变形梯度Fg＝diag[g，g，g]，其中，g3为软组织的体积变化率；进而利用体生长模型F＝FeFg，得到材料弹性变形梯度Fe；
S1.3：对于不可压缩的各向同性材料，假设其应变能函数为Ω(Fe，g-3FgDl)，由
Figure FDA0003536996830000011
得到生长之后的应力表达式
Figure FDA0003536996830000012
其中，Dl是相对于软组织生长之前的构型下的电位移张量，
Figure FDA0003536996830000013
为弹性变形的左柯西格林张量，D为相对于软组织生长之后的构型的电位移张量，Ωi为应变能函数对第i张量不变量的偏导数，p为拉格朗日乘子，I1为右柯西格林张量
Figure FDA0003536996830000014
的第一张量不变量；
S1.4：由
Figure FDA0003536996830000015
得到
Figure FDA0003536996830000016
这里E是电场强度张量，对于理想介电材料，进一步得到对应的电压值为
Figure FDA0003536996830000017
S1.5：建立应力和电场强度的平衡方程div σ＝0，div D＝0，curl E＝0，计算出软组织生长之后的几何尺寸、残余应力和电场强度的值；
S2：对于生长前后的残余应力和电场强度进行稳定性分析，获得满足增量平衡方程的模态，进而根据模态通过力电耦合效应控制生物组织的生长以及自组装；
所述的S2具体如下：
S2.1：通过线性增量理论分析，在生长后的构型上施加一个微小的增量位移场
Figure FDA00035369968300000111
和增量电位移
Figure FDA0003536996830000018
求偏导数得到增量位移梯度
Figure FDA0003536996830000019
进而代入
Figure FDA00035369968300000110
Figure FDA0003536996830000021
和
Figure FDA0003536996830000022
得到相对于生长后构型的增量应力
Figure FDA0003536996830000028
和增量电场强度
Figure FDA0003536996830000023
其中，
Figure FDA0003536996830000029
是对应于生长后构型中的物质点x摄动的位移；
Figure FDA00035369968300000210
为增量的拉格朗日乘子，
Figure FDA0003536996830000024
为增量电位移的相对于变形后构型的分量，
Figure FDA0003536996830000025
ΓI，KI为相对于生长后构型的增量即时模量；
S2.2：将增量应力的具体表达式代入到增量应力和电场的平衡方程
Figure FDA0003536996830000026
和边界条件
Figure FDA0003536996830000027
得到增量场的控制方程η′＝Gη，这里η是增量的位移应力以及电位仪向量，G是取决于生长和电场的Stroh矩阵；
S2.3：通过求解一个数学上存在的增量位移应力向量能够恰好满足控制方程η′＝Gη，即得到软物质生长之后失稳的临界时刻以及相应的模态，从而实现通过力电耦合效应控制生物组织的生长以及自组装。
CN201911400313.5A 2019-12-30 2019-12-30 基于非线性力电耦合生长模型的软物质生长自组装方法 Active CN111210880B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911400313.5A CN111210880B (zh) 2019-12-30 2019-12-30 基于非线性力电耦合生长模型的软物质生长自组装方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911400313.5A CN111210880B (zh) 2019-12-30 2019-12-30 基于非线性力电耦合生长模型的软物质生长自组装方法

Publications (2)

Publication Number Publication Date
CN111210880A CN111210880A (zh) 2020-05-29
CN111210880B true CN111210880B (zh) 2022-05-03

Family

ID=70787668

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911400313.5A Active CN111210880B (zh) 2019-12-30 2019-12-30 基于非线性力电耦合生长模型的软物质生长自组装方法

Country Status (1)

Country Link
CN (1) CN111210880B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104114349A (zh) * 2011-11-17 2014-10-22 哈佛学院院长等 用于制造聚合物纤维的系统、装置和方法
CN108536936A (zh) * 2018-03-27 2018-09-14 南京信息工程大学 一种多重优化的无网格软组织形变模拟方法
CN109033742A (zh) * 2018-06-21 2018-12-18 哈尔滨理工大学 一种用于模拟软组织剪切变形的超弹性模型

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7810750B2 (en) * 2006-12-13 2010-10-12 Marcio Marc Abreu Biologically fit wearable electronics apparatus and methods

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104114349A (zh) * 2011-11-17 2014-10-22 哈佛学院院长等 用于制造聚合物纤维的系统、装置和方法
CN108536936A (zh) * 2018-03-27 2018-09-14 南京信息工程大学 一种多重优化的无网格软组织形变模拟方法
CN109033742A (zh) * 2018-06-21 2018-12-18 哈尔滨理工大学 一种用于模拟软组织剪切变形的超弹性模型

Also Published As

Publication number Publication date
CN111210880A (zh) 2020-05-29

Similar Documents

Publication Publication Date Title
Previdi et al. Data-driven control design for neuroprotheses: a virtual reference feedback tuning (VRFT) approach
Villoslada et al. Position control of a shape memory alloy actuator using a four-term bilinear PID controller
Ding et al. Design of soft multi-material pneumatic actuators based on principal strain field
CN109557816B (zh) 一种压电陶瓷执行器迟滞特性的抑制方法、系统及介质
Ramasubramanian et al. Computational modeling of morphogenesis regulated by mechanical feedback
RU2002107315A (ru) Система и способ изготовления
McDaid et al. Adaptive tuning of a 2DOF controller for robust cell manipulation using IPMC actuators
CN111210880B (zh) 基于非线性力电耦合生长模型的软物质生长自组装方法
Xie et al. Improved sliding mode control with time delay estimation for motion tracking of cell puncture mechanism
Zaki Time integration of a model for martensite detwinning and reorientation under nonproportional loading using Lagrange multipliers
CN109271746B (zh) 气动肌肉高阶非线性动态模型及其构建方法与系统
WO2018219193A1 (zh) 一种基于有源和无源器件的小型力触觉再现方法
Arunachalam et al. Intelligent fractional-order PID (FOPID) heart rate controller for cardiac pacemaker
Gong et al. A bionic soft tongue driven by shape memory alloy and pneumatics
Luo et al. Hysteresis modeling and compensation of a pneumatic end-effector based on Gaussian process regression
Suzuki et al. Dynamic tracking control of an SMA wire actuator based on model matching
Cheng et al. Abnormal deformation and negative pressure of a hard magnetic disc under the action of a magnet
Amiri et al. Adaptive Particle Swarm Optimization of PID Gain Tuning for Lower-Limb Hman Exoskeleton in Virtual Environment
Kinoshita et al. Design of a database-driven kansei feedback control system using a hydraulic excavators simulator
Ru et al. Design and control of a soft bending pneumatic actuator based on visual feedback
CN112589798B (zh) 基于介电弹性体执行器的软体机器人状态反馈控制方法
Melchior et al. Analysis of muscle length effect on an S type motor-unit fractional multi-model
Meattini et al. A control architecture for grasp strength regulation in myocontrolled robotic hands using vibrotactile feedback: Preliminary results
Leal et al. Experimental and computational assessment of a shape memory alloy based morphing wing incorporating linear and non-linear control
Hu et al. Model-based contact detection and accommodation for soft bending actuators: An integrated direct/indirect adaptive robust approach

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant