虽然投nature 被拒了，但是这里也记录一下时间线，按照nature 投稿系统的排序结构，这里时间上由近到远。

Decision sent to author (Rejected) 20th September 21

Manuscript under consideration 09th Sept 21

All Reviewers Assigned 25th Aug 21

这之后editor会给corresponding author发邮件问confirmation form，顺便通知准备送审。

Editor Decision Started 26th July 21

Manuscript under consideration 23rd July 21

Editor assigned 23rd July 21

Manuscript received 22nd July 21

Manuscript under submission 22nd July 21

visual studio code 与matlab的配置方法

In order to have a more decent program editor for windows, I try to use vs code to write the m file for Matlab, however, I cannot run it in vs code. Then I start to search for a way to run Matlab script in vs code. Because of the updates of the vs code and related plugins, it takes me about one hour to figure out how to do this. I have read a Chinese tutorial, which is helpful but not enough. Here, I try to list the key steps for the configuration in vs code.

First, install plugins: Code runner (publisher:”Jun Han”), Matlab (publisher:”Xavier Hahn”), Matlab-formatter (publisher:”AffenWiesel”). Just choose the most popular one. Here, matlab-formatter is not necessary.

Second, modify the settings.json of vs code. the file of the setting.json is located in “C:\Users\Rocket\AppData\Roaming\Code\User\settings.json” if you install the vs code in the default directory. The modification parts are highlighted below. You may install Matlab in other folders that you need to change the related path as well.

{
    "python.pythonPath": "C:\\Program Files\\WindowsApps\\PythonSoftwareFoundation.Python.3.9_3.9.1776.0_x64__qbz5n2kfra8p0\\python3.9.exe",
    "editor.renderControlCharacters": true,
    "workbench.editorAssociations": {
        "*.ipynb": "jupyter-notebook",
        "*.m": "default"
    },
    "python.languageServer": "Pylance",
    "auto-build.defaultEnv.name": "STM32F103RET6_creality",
    "security.workspace.trust.untrustedFiles": "open",
    "python.defaultInterpreterPath": "C:\\Program Files\\Python39\\python.exe",
    "diffEditor.renderSideBySide": false,
    "matlab.matlabpath": "C:\\Program Files\\MATLAB\\R2021a\\bin\\matlab.exe",
    "matlab.mlintpath": "C:\\Program Files\\MATLAB\\R2021a\\bin\\win64\\mlint.exe",
    "matlab.linterEncoding": "utf8",
    "code-runner.runInTerminal": true,
    "code-runner.executorMapByFileExtension": {
        ".m" :"cd $dir && matlab -nosplash -nodesktop -r $fileNameWithoutExt",
        ".vb": "cd $dir && vbc /nologo $fileName && $dir$fileNameWithoutExt",
        ".vbs": "cscript //Nologo",
        ".scala": "scala",
        ".jl": "julia",
        ".cr": "crystal",
        ".ml": "ocaml",
        ".exs": "elixir",
        ".hx": "haxe --cwd $dirWithoutTrailingSlash --run $fileNameWithoutExt",
        ".rkt": "racket",
        ".scm": "csi -script",
        ".ahk": "autohotkey",
        ".au3": "autoit3",
        ".kt": "cd $dir && kotlinc $fileName -include-runtime -d $fileNameWithoutExt.jar && java -jar $fileNameWithoutExt.jar",
        ".kts": "kotlinc -script",
        ".dart": "dart",
        ".pas": "cd $dir && fpc $fileName && $dir$fileNameWithoutExt",
        ".pp": "cd $dir && fpc $fileName && $dir$fileNameWithoutExt",
        ".d": "cd $dir && dmd $fileName && $dir$fileNameWithoutExt",
        ".hs": "runhaskell",
        ".nim": "nim compile --verbosity:0 --hints:off --run",
        ".csproj": "dotnet run --project",
        ".fsproj": "dotnet run --project",
        ".lisp": "sbcl --script",
        ".kit": "kitc --run",
        ".v": "v run",
        ".vsh": "v run",
        ".sass": "sass --style expanded",
        ".cu": "cd $dir && nvcc $fileName -o $fileNameWithoutExt && $dir$fileNameWithoutExt"
    }
}

The third step, save the JSON files and try to run the Matlab file in vs code.

To use matlab-formatter, you need to right-click the open m file and choose the command ” format document”.

the running window is as follows:

The drawback of this method is that the running result will show up in a new window rather than in the terminal of vs code. And also be careful of the name of the running script, the program (Matlab ) may run the script with the same name but in a different folder.

after using this several times, it is not convenient. it is better to editor the script with vs code but run the script in Matlab GUI, then no need to close the pop-up windows. The MathWorks company should do something, even the live script is developed, but not convenient to use.

Recently buy the Ender5 pro, a very cheap version, and only 299 euros from 3Djake.at. The frame of the 3D printer is very solid, no complaints.

Do a simple modification of the 3D printer and update the printer with a BLTouch sensor and use a Raspberry Pi (with Octoprinter) to control it.

Because the motherboard is V4.2.2 which is special and it costs some time to figure out how to connect the BLtouch with the 3D printer. We change the firmware with Marlin. Before this mode, I also install the BLtouch with CR-10S which is simpler. And be careful of the compatibility between the Marlin Version and the build environment. For the Ender 5 pro, you need to use PlatformIO. I suggest using Microsoft visual studio code with the Marlin plugin (quite easy to modify).

For the printing head, we have tried both 0.4 mm and 0.2 mm nozzles. Quality is fine. 0.1 mm nozzle is not support.

2021-8-24 12:24:19

A “like” button is added to all the pages of the website. This can help me to value my articles if visitors press a lot.

Science 投稿记录 2021年7月

submission 09-Jul-2021

Under Evaluation 09-Jul-2021

To Advisor 09-Jul-2021

Under Evaluation 11-Jul-2021

Under Evaluation 19-Jul-2021

Rejected-Transfer Offered 21-Jul-2021

• Science
• Received 2019-11-20
• Under Evaluation  2019-11-20
• To Advisor 2019-11-22
• Rejected-Transfer Offered-12-Dec-2019

• Received – Transferred-12-Dec-2019
• Under Evaluation-12-Dec-2019
• Under Evaluation 16-Dec-2019
• Rejected-Transferred 18-Dec-2019

没想到在外国，护照竟然过期了，还得从LINZ去WIEN 维也纳更换护照。

首先需要网上申请(http://www.chinaembassy.at/chn/lsyw/c/t1455618.htm)，然后准备材料。网上申请后，打印申请表，以及附带的几张表格。

换护照的地点是Neulinggasse 29, 1030 Wien， 坐电梯上三楼。因为网上预约的只是一个时间段，早上或者下午。先取号，再办理业务。不得不说，中国领事馆内的时候有点吵，开放空间，三个窗口，自带喇叭，略像小型菜市场。不过工作人员都挺有耐心的。

以下内容采集自中国驻奥地利大使馆的官方网站。

一、普通护照换发
　　（一）适用情形
　　申请人具有中国国籍、持有中华人民共和国普通护照，且有下列情形之一者，可申请换发：
　　1、护照有效期不足一年的；
　　2、护照签证页即将用完的；
　　3、护照上有携行儿童，需单独持用护照的；
　　4、容貌变化较大的；
　　5、合法定居在奥地利但护照过期；
　　6、未成年人在当地使用新的名字。
　　（二）所需材料
　　1、如实、完整填写《中华人民共和国护照/旅行证/回国证明申请表》1份（电子版需要打印出来）；电子版申请表自带《国籍状况声明书》（须当场填写并签名）；
　　2、近期（半年内）正面免冠彩色半身证件照片（光面相纸）2张（工作人员只会要求一张，因为电子版《申请表》自带一张照片）。照片尺寸为48mmx33mm，头部宽15mm-22mm、长28mm-33mm，背景为白色；
　　3、护照原件及复印件（包括带照片资料页、曾办理过换发、加注或延期页）；
　　4、有效期不少于一个月的奥地利居留证件原件和复印件；
　　5、持有有效居留证件但护照过期者还需提交护照过期的书面说明；
　　6、在奥地利出生、无国内户籍的中国籍未成年人因更名申请换照者，还需提供奥地利相关的市政府（Magistrat)和户籍部门（Standesamt)出具的使用新名的相关证明文件（Meldesbestaetigung und Geburtsurkunde)；申请人父母出具的协商同意为子女更名的书面证明。

7. 除此之外，还需要： 回邮服务申请表，贴有邮票的信封。具体参见末尾的回邮服务。
　　（三）申请方式
　　1、申请人应亲自到中国驻奥地利使馆领事部办理，不能邮寄申请或委托他人代办。
　　2、未满16周岁未成年人换发护照应由父母双方（或法定监护人）陪同，并出具父母双方（或法定监护人）同意为儿童办理护照的书面声明。如父母有一方不能陪同，应出具具有法律效力的委托书，同意并委托另一方为儿童办理换照。
　　（四）办证时间
　　由于护照由制证中心统一制作，从申请人到馆开始申请计算15个工作日，无加急业务。
　　（五）取证及收费
　　1、申请人在取证单标明的日期，或按电话通知，由本人或委托他人凭取证单原件交费并领取证件。
　　2、普通护照换发收费16欧元，中国驻奥地利使馆只接受银行储蓄卡（POS机刷卡）和信用卡，不接受现金、支票。
　　3、护照办妥一年内未取证者，证件将被作废销毁，并将可能影响申请人下次证件申请。
　　（六）注意事项
　　换发新照后，原照将被注销并退还申请人。

回邮服务：

驻奥地利大使馆关于提供护照、旅行证回邮服务的通知2018/10/02　　为进一步便利在奥中国公民申领护照、旅行证，提高领事证件服务水平，中国驻奥地利大使馆拟于2018年10月1日起（以递交申请日期为准）提供护照旅行证回邮服务。现将回邮服务有关事项通知如下：
　　一、适用范围：
　　1、回邮服务对象仅限居住在奥地利的申请人；
　　2、凡在我馆申请护照颁发、换发、补发、姓名加注、旅行证业务，递交申请时资料齐全、无需补充签名及注销旧照等情形的，均可选择回邮服务。
　　3、我馆暂不提供签证回邮服务。
　　二、具体办法：
　　请注意：我馆暂不收取回邮服务费，但申请人需自行承担邮寄资费。
　　1、申请人如选择回邮服务，需自行购买信封，用于回寄护照、旅行证。为确保证件安全，建议选择如下信封：
　　Typ:Luftpolsterumschlag
　　Maße: 120x210mm
　　Standard-Umschlagformat:B/00；（信封大一点没关系，要求是能投入邮筒。A4纸大小 都是可以的）
　　2、贴好奥地利境内挂号（Einschreiben）所需邮资的邮票；（目前是5.05欧元邮票，可以前往邮局购买，据说可以在奥地利境内寄2kg的文件）
　　3、在回邮信封上正确、清楚填写收件人姓名、地址（含邮编）。如收件人信息与申请人的个人信息不同，请在信封背面注明申请人的姓名；
　　收件人姓名
　　街道名，门牌号
　　邮编，城市名
　　（以上信息请用德语书写，字迹务必清晰可辨认）
　　4、受理申请后即凭取证单到收费窗口交费（护照换发的费用，回邮服务不收费）并将回邮信
　　封交工作人员。
　　三、特别声明：
　　1、护照办理时间通常为15个工作日（无加急），旅行证办理时间为4个工作日（可加急），使馆在取证日期次日起5个工作日内将护照、旅行证寄出。
　　2、使馆无法掌握回邮具体进度，不接受邮件状态查询。申请人可通过中国公民办证咨询电话或电子邮箱向使馆索要邮寄挂号单据，并自行向邮局进行查询。
　　3、申请人选择回邮服务时需填写“回邮服务申请表”。因邮寄而产生的证件丢失、损毁、延误等情况，使馆不承担相关责任。
　　使馆再次感谢广大中国同胞对我馆领事工作的支持和理解！今后我馆将继续创造条件，争取推出更多便利中国公民的办证措施，提升服务质量和水平。
　　中国驻奥地利大使馆
　　2018年9月24日

配置 ABAQUS的用户子程序

My platform is listed as below

Windows 10 Pro 21H1
Device name DESKTOP-KO7EL1P
Processor AMD Ryzen 9 3950X 16-Core Processor 3.50 GHz
Installed RAM 64.0 GB
System type 64-bit operating system, x64-based processor

Software platform:
Visual studio community 2015
Intel parallel studio XE 2016 Update 1 Cluster Edition
ABAQUS 2018

The steps to install the software, 1. visual studo, 2 intel, 3 ABAQUS
when installing intel, please make sure that installing all the components, especially intel 64. Generally, the intel stdio will connect the visual studio during installment (you can even specify the version of studio to link during the installment)

After installing all the software, change the file (abq2018.bat) in the fold(: C:\SIMULIA\Commands as follows:

@call "C:\Program Files (x86)\IntelSWTools\compilers_and_libraries_2016.1.146\windows\bin\ifortvars.bat" intel64 vs2015
@echo off
setlocal
set ABA_COMMAND=%~nx0
set ABA_COMMAND_FULL=%~f0
"C:\SIMULIA\CAE\2018\win_b64\code\bin\ABQLauncher.exe" %*
endlocal

The subroutine in ABAQUS should work now.

介电高弹体数值模拟，ABAQUS，有限元模拟

Finite element for the dielectric elastomer. 

When I studied the dielectric elastic, there are not so many ways to simulate the dielectric elastomer. The only way to do it is using a UMAT[1] which is written by Xuanhe Zhao [1]. However, it is not easy to reprogram it, especially since the hyperelastic model is complex, for example, the Gent model. Of course, there is a lot of ways to solve this problem. But it is so difficult for entry-level research who is not familiar with the numerical program. In 2014, I started to know the problem: wrinkle of dielectric elastomer. I have tried with Xuanhe’s UMAT, but it does not work well. I also reprogram it with the Mooney-Rivlin model, but still not work well with my wrinkle problem. I think it is because of the large strain of the dielectric elastomer. And As in the UMAT, I need to build the Jacobian matrix (DDSDDE) which is very difficult.

Then Wei hong introduces a method using a new subroutine called UANISOHYPER_INV, which is very easy to program and can be easily programmed with different hyperelastic models, also using the build-in function of ABAQUS. All the theoretical background is provided in our paper [3]. In the paper, we did not show the details, here, I would describe a little bit how to use this method.

1. Create a material, but don’t add anything

2. Model-\textgreater Edit Keywords
Find the material you just defined. “Add After” *Material, name=…
Add the following 2 lines:

*ANISOTROPIC HYPERELASTIC, USER, FORMULATION=INVARIANT, TYPE=INCOMPRESSIBLE, LOCAL DIRECTIONS=1, PROPERTIES=3

20.62e3,3.983e-11,220

The second line defines two material parameters, the modulus, and the nominal electric field (voltage/original thickness). This example is for gent model that shear modulus=20.62e3, permittivity= 3.983e-11, stretch limit (J), 220

3. Create a section, a part … and all regular stuff

4. Model-\textgreater Edit Keywords Find the part of EAP, Add the following 3 lines before *Element, type=…

*ORIENTATION, NAME=FIELD, LOCAL DIRECTIONS=1

1, 0, 0, 1, 1, 0

3,0

0, 0, 1

This defines a local coordinate system: x in the 100 direction, y in the 010 direction, and 0 rotation (second line)
The last line means the electric field is applied in the z-direction of the local system.
for cylindric system we should define the z-axis

*ORIENTATION, NAME=FIELD, system=cylindrical,LOCAL DIRECTIONS=1

0, 0, 0, 0, 1, 0

3,0

0, 0, 1

5. Find the section you created, add the following keyword to the same line of *Solid Section, …
orientation=FIELD
This would assign the direction defined in step 4 to this section

6. Create a job and add the Fortran file as the user subroutine

Please cite our paper [3], if you would like to use this method.

Reference.
1. UMAT is a subroutine in ABAQUS
2. Zhao, X., and Suo, Z., 2008, “Method to analyze programmable deformation of dielectric elastomer layers,” Applied Physics Letters, 93(25), p. 251902.
3. Mao, G., Hong, W., Kaltenbrunner, M., and Qu, S. (May 19, 2021). “A numerical approach based on finite element method for the wrinkling analysis of dielectric elastomer membranes.” ASME. J. Appl. Mech. doi: https://doi.org/10.1115/1.4051212

code: the first zip file is the subroutine (UANISOHYPER_INV) for the Mooney-Rivlin model, also an input file is added as an example. The second zip file is the subroutine (UANISOHYPER_INV ) for the Gent model and the neo_Hookean model.

Files: